Use of “one-pot, mix-and-read” peptide-MHC class I tetramers and predictive algorithms to improve detection of cytotoxic T lymphocyte responses in cattle
© Svitek et al.; licensee BioMed Central Ltd. 2014
Received: 15 October 2013
Accepted: 7 April 2014
Published: 28 April 2014
Peptide-major histocompatibility complex (p-MHC) class I tetramer complexes have facilitated the early detection and functional characterisation of epitope specific CD8+ cytotoxic T lymphocytes (CTL). Here, we report on the generation of seven recombinant bovine leukocyte antigens (BoLA) and recombinant bovine β2-microglobulin from which p-MHC class I tetramers can be derived in ~48 h. We validated a set of p-MHC class I tetramers against a panel of CTL lines specific to seven epitopes on five different antigens of Theileria parva, a protozoan pathogen causing the lethal bovine disease East Coast fever. One of the p-MHC class I tetramers was tested in ex vivo assays and we detected T. parva specific CTL in peripheral blood of cattle at day 15-17 post-immunization with a live parasite vaccine. The algorithm NetMHCpan predicted alternative epitope sequences for some of the T. parva CTL epitopes. Using an ELISA assay to measure peptide-BoLA monomer formation and p-MHC class I tetramers of new specificity, we demonstrate that a predicted alternative epitope Tp229-37 rather than the previously reported Tp227-37 epitope is the correct Tp2 epitope presented by BoLA-6*04101. We also verified the prediction by NetMHCpan that the Tp587-95 epitope reported as BoLA-T5 restricted can also be presented by BoLA-1*02301, a molecule similar in sequence to BoLA-T5. In addition, Tp587-95 specific bovine CTL were simultaneously stained by Tp5-BoLA-1*02301 and Tp5-BoLA-T5 tetramers suggesting that one T cell receptor can bind to two different BoLA MHC class I molecules presenting the Tp587-95 epitope and that these BoLA molecules fall into a single functional supertype.
CD8+ cytotoxic T lymphocytes (CTL) play a significant role in controlling intracellular pathogens such as viruses and parasites through a perforin/granule mediated cytolytic mechanism, Fas-induced apoptosis or by release of IFN-γ . The specificity of CTL is mediated by their T-cell receptor (TCR) which recognize foreign peptide epitopes associated with major histocompatibility complex (MHC) class I molecules at the surface of infected cells (reviewed in ). The antigen specificity of CTL is commonly measured using synthetic peptide epitopes in IFN-γ ELISpot or cytolytic assays. However, the demonstration that soluble peptide-MHC (p-MHC) class I complexes will bind to TCR has led to generation of multimeric forms of p-MHC class I tetramers as diagnostic reagents that can be used to detect and purify epitope specific CTL. The importance of these reagents in research is highlighted by the creation of a NIH Tetramer Core Facility, which provides to the scientific community pre-made and made-to-order p-MHC class I tetramers for a number of human, mouse, macaque and chimpanzee MHC alleles, including class I, class II and non-classical MHC molecules. Many p-MHC class I tetramers are now also commercially available. However, such tools are in general not available for livestock species.
Peptide-MHC class I complexes are hetero-trimeric in nature consisting of the MHC class I molecule, often referred to as the heavy chain (HC), a peptide of eight to eleven amino acids in length and an invariant molecule called β2-microglobulin (β2m) (reviewed in [3, 4]). Generation of the p-MHC class I complex in vivo involves many steps and is tightly regulated by chaperones . Peptides are generated by the proteasome, a multi-catalytic complex in the cytoplasm that cleaves proteins into polypeptides , which are further processed by ER-luminal peptidases  during peptide loading onto MHC class I molecules by the peptide-loading-complex  prior to export to the cell surface. However, it is possible to generate p-MHC class I complexes at high efficiency in vitro [9–13] by co-incubation of recombinant MHC class I HC, synthetic peptide and recombinant β2m. The affinity of peptides that bind to the peptide-binding groove in the α1 − α2 domain of the HC is primarily dictated by the identity of amino acids at “anchor positions” and sequence motifs describe peptides that fit into a particular peptide-binding groove . Sequence polymorphisms in this region of the HC result in different peptide binding specificities and motifs, which mediate the functional differences between different MHC class I molecules . These observations have led to the generation of several algorithms ([16–18], for a detailed review see ) that attempt to predict the presence of CTL epitopes in proteins.
Different methods have been developed for generating multimeric peptide MHC complexes [13, 20, 21]. Here we describe the use of the system of Leisner et al.  for production of bovine leukocyte antigen (BoLA) p-MHC class I tetramers to facilitate the study of immune bovine CTL as part of our subunit vaccine development program against East Coast fever (ECF), a lethal disease of cattle in sub-Saharan Africa. The disease kills about one million animals each year with severe negative socio-economic impact on the smallholder farmers and pastoralists whose livelihoods depend on their animals [22, 23]. Animals that recover from infection are solidly immune to re-infection and this observation led to development of a live parasite based infection and treatment method (ITM) of vaccination against ECF . The immunity induced by ITM is primarily dependent on parasite-specific CTL and is parasite strain specific. Vaccination with mixtures of parasite broadens the spectrum of immunity [24, 25] and such a vaccine, called the Muguga cocktail, which consists of three different parasite isolates, is available on a commercial basis in some countries . However, there are several problems with the live vaccine as it requires a liquid nitrogen cold chain for delivery, oxy-tetracycline treatment and is expensive ($8 ~ 12/dose). Hence, the availability of a subunit vaccine with a better product profile than the current vaccine is likely to benefit programs targeting the control of ECF.
A feature of the protective CTL response generated by ITM is that it is focused on a limited number of epitopes, the identity of which is influenced by the BoLA genotype of individual animals and the strains of the parasite used for immunisation . The antigen and epitope specificity of CTL induced by the Muguga isolate of T. parva and restricted by ten different BoLA class I molecules in cattle has been reported [28, 29], and unpublished data]. Identification of these T. parva CTL epitopes was carried out in a two stage process, first through screening of parasite cDNA to define relevant antigens and then through screening of synthetic overlapping peptides to define epitopes within the predicted antigen sequence . Although 12 CTL epitope sequences were defined they were primarily evaluated in IFN-γ ELISpot assays and most were not subjected to systematic analyses to determine minimal epitope sequences. We previously evaluated NetMHCpan, a pan-predictor of CTL epitopes to assess the T. parva epitopes. In general we found that the predictive power of the algorithm for the known CTL epitopes was high, although the method in some cases predicted the presence of alternative epitope sequences of stronger affinity .
Here we report on generation of p-MHC class I tetramers based on a set of seven different recombinant BoLA MHC class I molecules for detection of CTL responses to epitopes located in five different T. parva antigens and the use of one p-MHC class I tetramer in ex vivo assays. A critical feature of the Leisner system is that new p-MHC class I tetramers can be made on demand with a turnaround time of ~48 h, which as described herein accelerates the ability to test candidate epitopes and define new CTL specificities. Based on predictions from NetMHCpan and MHCcluster we provide experimental evidence that redefines a CTL epitope on antigen Tp2 and that an epitope on antigen Tp5 can be presented by two closely related but distinct BoLA MHC class I molecules.
Material and methods
Expression and purification of BoLA class I HC and β2m molecules
The expression and purification of biotinylated BoLA class I HC molecules was performed essentially as described in . Briefly, truncated versions of the BoLA class I HC without their trans-membrane domain and containing a biotinylation tag (BSP) were cloned into pET28a+. The proteins were expressed in E. coli BL21 (DE3) co-transformed with a pASYC vector encoding the BirA enzyme, which increases the efficiency of biotinylation of the BSP tag. Expression of HC and BirA was induced by adding IPTG (1 mM) during the last 3 h of culture. Biotin (Sigma, 100 mg/L) was added to the culture prior to inducing protein expression. Large-scale production and purification were performed as previously described . Briefly, this involves culture of E. coli BL21 (DE3) in ECPM 1 media  using 2-L Labfors fermenters (Infors AG). Cell lysis was performed by high-pressure cell disruption (Z PLUS, Constant Systems, Daventry, UK). Inclusion bodies were collected by centrifugation at 17 000 g for 10 min at 4 °C. Washed inclusion bodies were solubilized in 20 mM Tris-HCl (pH 8.0) and 8 M urea. HC molecules were purified in 8 M Urea first on Q-Sepharose (Amersham Biosciences) then on phenyl-Sepharose column (Amersham Biosciences) and high molecular complex aggregates were removed by size exclusion chromatography on Sephacryl columns. Bovine β2m molecule was purified using the same method as for human and porcine β2m [31, 33, 34]. Briefly, recombinant bovine Met-β2m was dissolved in 8 M urea and purified on Q-Sepharose; refolded overnight by drop-wise dilution into 25 mM Tris, 300 mM Urea, pH 8; concentrated by tangential flow ultrafiltration (Vivaflow); and buffer-exchanged into PBS by gel filtration (Superdex 200, Amersham Biosciences). Finally, the product was concentrated and frozen at −20 °C until use.
Production of peptide-MHC class I tetramers
Peptide-MHC class I tetramers were generated as described previously . Briefly, purified BoLA HC at a final concentration of 0.25 μM was incubated with the following: β2m [1.5 μM], peptide [4 μM], EDTA (1 mM), pepstatin A (1 μM), 1.10 phenantroline (100 μM), TLCK (50 μg/mL), TPCK (100 μM), PMSF (100 μM), Lutrol (0.03%), 50 mM of Tris-maleate (pH6.6), in a final volume of 500 μL. Incubation was performed at 18 °C for 48 h and then high molecular weight complexes were removed by centrifugation with an Amicon® Ultra Centrifugal Filter Units (0.5 mL/100 K) [cat no. UFC510024] at 14 000 g for 10 min. The flow through was recovered and the p-MHC class I complexes were tetramerized at room temperature by incubating them with 1.5 μg of phycoerythrin (PE)-streptavidin (BD Pharmingen, cat no. 554061) or allophycocyanin (APC)-streptavidin (BD Pharmingen, cat no. 554067) starting with 0.45 μg and then adding 0.21 μg every ten min over a 50 min period at room temperature and then stored at 4 °C.
List of bovine immunological reagents
Restricting BoLA class I molecule
T. parva antigen
BB007 & BV115
BW02 (Clones 5 & 7)
ELISA assay to measure peptide-MHC class I monomer formation
Folding of p-MHC class I complexes was assessed with an ELISA assay as described in . The BoLA-6*04101 molecules [0.25 μM] and bovine β2m molecule [1.5 μM] was incubated in a 96-well plate with Tp227-37 and Tp229-37 peptides at various peptide concentrations ranging from 0 to 40 μM. After an incubation of 48 h at 18 °C, the complexes were transferred to another 96-well plate pre-coated with streptavidin (Nunc, cat# 436014), which captures biotinylated HC, and incubated for three hours at 4 °C. After washing, the W6/32 monoclonal antibody (Santa Cruz Biotechnology, Texas, USA,, cat # SC-32235), which binds to a monomorphic epitope on β2m but only when it is incorporated in p-MHC class I complexes, was added for 1 h at 4 °C. After washing, anti-mouse IgG coupled to peroxidase (Sigma Aldrich, cat # A9917-1ML) was added to the plate followed by incubation for an hour at room temperature and further rounds of washes. Colorimetric change was performed by adding the TMB Plus2 “Ready to Use” Substrate (Kem En Tec, cat # 4395H) for 10 min at room temperature. The reaction was stopped by adding H2SO4 (0.3 M) and optical density was measured using an ELISA plate reader at 450 nm. The concentration of MHC class I and β2m molecules used in such assays is normally in the nanomolar range . However, the bovine β2m molecule used in this assay contains an extra methionine at the N-terminus that diminishes binding of the W6/32 antibody to the mature p-MHC class I complex, hence higher concentrations of recombinant molecules were used in the assay. We are in the process of generating a β2m without the extra start methionine, which should allow use of recombinant molecules in the nanomolar range.
Bovine CTL and parasite-infected cell lines
CTL and autologous T. parva Muguga infected lymphocyte lines (TpM) were established as previously described  and maintained in RPMI containing 10% fetal bovine serum [FBS] (decomplemented; Sigma), 1% penicillin/streptomycin, 1% L-Glutamine, 1% gentamycin, 0.1% beta-mercaptoethanol. CTL lines (Table 1) were propagated by stimulating them with irradiated autologous TpM at a ratio of 1:5 to 1:10 (CTL:TpM) and supplementing the media with 2 ng/mL of recombinant human interleukin (IL)-2 (Sigma) or 15% T-cell growth factor (TCGF). Cells were incubated at 37 °C in a 5% CO2 humidified atmosphere. Clones 5 and 7 of the BW02 CTL line were obtained by limiting dilution.
IFNγ ELISpot assays were conducted as described in . Briefly, a monoclonal anti-bovine IFNγ antibody (Serotec, Oxford, UK, cat.no. MCA1783) was incubated overnight at 4 °C on ELISpot plates (Millipore, Billerica, MA, USA) and then blocked with RPMI containing 10% heat-inactivated FBS for 2 h at 37 °C. Peptides were added at 1 μM concentration and CTL lines at 2.5 × 104 cells per well. The plates were incubated at 37 °C for 20 h. Release of IFNγ was monitored with primary rabbit polyclonal anti-bovine IFNγ antibody (Sigma–Aldrich, St. Louis, MO, USA) and secondary AP-conjugated monoclonal anti-rabbit antibody (Sigma–Aldrich, St. Louis, MO, USA, cat # A2556). Development of plates was done by addition of the substrate solution Sigma Fast (BCIP/NBT, Sigma–Aldrich, St. Louis, MO, USA).
Flow cytometry data were acquired on a BD FACSCanto™ II instrument (BD Biosciences). Compensation controls for PE, PerCP, FITC and APC were included for automatic compensation by the FACS DIVA software. 10 μL BoLA class I tetramer was used for staining cells (between 2 × 105 and 5 × 105 cells per sample). CD8 staining was done by incubating cells with ILRI monoclonal antibody ILA51 (IgG1) and using rat anti-mouse IgG1-PerCP (BD; cat # 340272) as the secondary antibody. Staining for perforin was done using anti-human-perforin-FITC (BD Pharmingen, San diego, CA, USA, cat# 556577) and staining for Fas-ligand was done using primary rabbit anti-Fas ligand (Santa Cruz biotechnology, cat no. SC-957) and secondary FITC-coupled anti-rabbit IgG antibody (Sigma, MO, USA, cat# F-0382). All staining were conducted in PBS-0.5% BSA except for intracellular staining with anti-perforin where cells were incubated and stained in PBS-0.1%Saponin-10% FBS. Negative control for intracellular staining was done with anti-dansyl-FITC antibody (BD Pharmingen, San diego, CA, USA, cat# 556577). The data were analyzed using the FlowJo software version 7.6.5.
Three cattle of BoLA-A18 haplotype (typed with PCR primers [forward] 5’-CCGGGATCCGAGGACT-3’ and [reverse] 5’-CTCCATCTTGCGTTTGGA-3’ specific for BoLA-6*01301) were vaccinated by ITM  using the T. parva strain Muguga 3308. Whole blood was collected from these cattle at the day of vaccination (day 0) and every two days till appearance of tetramer positive cells. PBMC were purified on Ficoll paque density gradients at 2500 rpm for 25 min and 5 × 105 PBMC were stained with the Tp1214-224-BoLA-6*01301 tetramer as described in the flow cytometry procedure. This animal experiment was approved by the ILRI IACUC committee [# 2012.10].
NetMHCpan and MHCcluster
We used NetMHCpan, a pan-specific predictor trained on known peptide epitope and MHC class I sequence data and the affinity of p-MHC class I interactions to predict T. parva epitopes. The binding affinity of selected epitopes was predicted using the latest version of NetMHCpan (version 2.8, ). Predictions of strong peptide binders was done by evaluating all possible eight to eleven-mers from antigens Tp2 and Tp5 with BoLA-6*04101, and BoLA-T5 and BoLA-1*02301 respectively. We also used MHCcluster to compare predicted peptide-binding specificities between the different BoLA class I molecules involved in this study.
Comparison of sequence similarities between BoLA-T5 and BoLA-1*02301
BoLA-1*02301 and BoLA-T5 amino acid sequences were downloaded from the Immuno-Polymorphism Database [40–42] and NCBI respectively, and a MUSCLE alignment was done using clcSequence Viewer 6.0. Visualization of the position of the amino acid differences between these molecules was done on the crystal structure of BoLA-2*01801  using RasWin.
Generation of peptide-MHC class I tetramers and assay on defined CTL lines
The “on-demand” Leisner system for tetramer generation involves p-MHC class I monomer formation in the presence of excess peptide (concentration of 4 μM), followed by tetramerization using either PE or APC labelled streptavidin (see Materials and methods). This method was used to generate a set of seven p-MHC class I tetramers from recombinant BoLA class I HC and bovine β2m and synthetic peptides as listed in Table 1. This system does not routinely measure monomer formation but this became necessary for some experiments (see below).
A role for NetMHCpan in T. parva CTL epitope identification
We have previously reported that the NetMHCpan algorithm  predicted an alternative epitope sequence for six of the 12 epitopes . Since that time NetMHCpan has been improved and version 2.8 predicts alternative epitopes for only 4 of the 12 epitopes: Tp7207-214 instead of Tp7206-214, Tp241-48 instead of Tp240-48, Tp2138-146 instead of Tp2138-147 and Tp229-37 instead of Tp227-37. As a measure of predictive accuracy, we use the false positive fraction score. For a given CTL epitope, this score is defined as the number of 8-11mer peptides with a predicted binding affinity stronger than the epitope divided with the total number of 8-11mer peptides found within the source protein sequence. A false positive fraction score closer to 0 is indicative of a high predictive performance.
NetMHCpan predicted that Tp229-37 (EELKKLGML) has a false positive fraction score of zero for BoLA-6*04101 while that for the peptide Tp227-37 (SHEELKKLGML) we have been working with (Table 1) was 0.039 (this value has changed from previous predictions with older version of Net MHCpan but is still indicative of a low binding affinity ). To experimentally verify this prediction we developed a semi-quantitative ELISA assay as described in  to measure p-MHC class I monomer formation (see Materials and Methods). The W6/32 monoclonal antibody used in the ELISA assay binds to β2m from different species but only when β2m is in a complex with HC and peptide. If a stable p-MHC class I complex is not formed, for example in the absence of peptides, β2m remains free in solution and will not bind to the streptavidin-bound HC on the plate. Figure 3A clearly demonstrates that Tp227-37 does not bind to BoLA-6*04101, whereas Tp229-37 demonstrates a good binding curve.
We then generated p-MHC class I tetramers made with each one of the two Tp2 sequences to stain BW02 CTL line clones that in ELISpot assay respond to both peptide sequences (Figure 3B). Figure 3C shows a clear difference in staining profile between the p-MHC class I tetramer generated with Tp229-37 and the p-MHC class I tetramer generated with Tp227-37, with the latter showing a profile similar to that seen for non-specific staining of cells. Thus, our results strongly support the NetMHCpan prediction that Tp229-37 and not Tp227-37 is the correct epitope sequence for BoLA-6*04101.
Tp1-tetramer positive CTL are also positive for CD8, FasL, and perforin
Ex vivo analysis of CTL following ITM vaccination of cattle
Predicted peptide binding specificity as a means of clustering of the seven BoLA MHC class I molecules
Peptide-MHC class I tetramer complexes have become essential tools in cellular immunology to decipher and characterize epitope-specific CD8+ T cell immune responses [12, 13, 21]. They are especially useful for detection of CTL that occur at low frequency and are finding many applications in ex vivo assays for direct assessment of T cell responses. Such reagents have become indispensable tools in murine and human immunology but are generally lacking for many important livestock species. In an effort to provide reagents to fill this gap, we recently described initial efforts to develop p-MHC class I tetramers for use in swine and bovine research [30, 34] using a simple in-house system that permits generation of p-MHC class I tetramers with a turnaround time of 48 h . Here, we report on p-MHC class I tetramers made from seven BoLA MHC class I genes and the power of integrating the on demand tetramer system with CTL epitope prediction data from NetMHCpan, a pan specific algorithm, and MHCcluster, a method for clustering of MHC molecules based on function, i.e., their predicted peptide binding specificities, rather than sequence.
We used the novel immunological reagents we generated to facilitate our research on ECF, a lethal disease of cattle. Parasite-specific CTL induced by the ITM vaccine play a major role in mediating immunity to ECF [25, 28, 29]. A number of T. parva antigens that are the targets of CTL induced by ITM have been identified using high-throughput IFN-γ ELISpot immunoscreens of antigen presenting cells transfected with either candidate T. parva genes or pools of parasite cDNA libraries . Positive genes were then tested in cytolytic assays and CTL epitopes mapped by use of synthetic overlapping peptides. Based on such data a panel of 12 experimentally verified epitope sequences derived from seven T. parva antigens and restricted by ten BoLA class I molecules have been reported  and unpublished data]. However, these peptide sequences were not subjected to systematic analyses to determine minimal epitope sequence data and due to a lack of appropriate reagents the epitopes were not verified in independent assays, e.g., in peptide-BoLA binding assays. The prediction of alternative sequences for some of the CTL epitopes by NetMHCpan has raised questions regarding the performance of the algorithm as well as the veracity of CTL epitope sequences. We have addressed some of these issues.
The method of Leisner et al. works well for production of seven recombinant BoLA MHC class I molecules, adding to the list of mammalian MHC molecules generated by this system. We have made and validated the use of p-MHC class I tetramers that allow us to rapidly identify T. parva specific CTL responses directed to one of seven epitopes on five different antigens (Figures 1, 2, and 3). We used these reagents to demonstrate that all Tp1-specific CTL co-express CD8, FasL or perforin (Figure 4). Such assays facilitate rapid typing of the predicted functionality of T. parva specific CTL. This is important for subunit vaccine development studies since in the past we have demonstrated priming of peptide specific IFN-γ ELISpot CD8+ T cell responses in the absence of cytolytic activity , and only the latter activity correlates with immunity to ECF. The ability to detect cytolytic activity in peripheral blood post-vaccination is the range of 14 to 21 days [48–50]. However, several rounds of CTL in vitro stimulation are required for detection of cytolytic activity. Using the Tp1-specific tetramer we were able to detect CTL in BoLA-A18 cattle as early as day 13 post-vaccination by the live parasite based ITM vaccine (data not shown). The CTL response peaked at about day 15 (Figure 5). Thus, the Tp1-specific tetramers demonstrate high sensitivity and great potential as diagnostic tools for detecting ECF immunity in cattle with the ability to stain cells isolated ex vivo from infected animals without further stimulation and expansion in vitro. While the remaining p-MHC class I tetramers remain to be validated in ex vivo assays, we expect that our collection of p-MHC class I tetramers will permit studies on the kinetics and functionality of antigen specific CTL responses to T. parva during infection and facilitate immunogenicity studies during subunit vaccine trials. Moreover, as new CTL epitopes restricted by these BoLA molecules are identified, in future studies we will be able to rapidly expand the repertoire of p-MHC class I tetramers available for ECF research.
Several unexpected results were predicted through the use of NetMHCpan and MHCcluster. The first is that NetMHCpan predicted that Tp229-37 rather than Tp227-37 is likely to be the correct peptide sequence that binds to BoLA-6*04101 (Figure 3). The longer Tp2 peptide sequence works in both ELISpot and cytolytic assays (our data and ) but this peptide does not bind to BoLA-6*04101 (Figure 3). We presume that the peptide is processed to the functional epitope during the ELISpot assays, which does not occur during the peptide MHC binding assays or formation of p-MHC class I tetramers. NetMHCpan also predicted three alternative epitopes from those we have previously described: Tp7207-214 instead of Tp7206-214, Tp241-48 instead of Tp240-48, and Tp2138-146 instead of Tp2138-147. Due to a lack of CTL lines of appropriate specificity, we have not been able to test these predictions. However, Tp7207-214 has been previously shown to be negative in ELISpot assays . Preliminary data from our group (not shown) has determined a higher binding affinity of Tp7207-214 for BoLA-T7 (Kd value of 142 nM) compared to the binding affinity of Tp7206-214 (Kd value of 252 nM). This discrepancy and the alternative epitopes predicted for Tp2 remains to be resolved.
One of the outputs of the MHCcluster server is a tree-based (Figure 6A) and heat-map (Figure 6B) visualization of the functional relationship between MHC molecules based on their predicted peptide binding specificity rather than sequence alignment of the MHC molecules themselves. For example, it is known that HLA-A*0201 and HLA A*0301 share 93% sequence identity but exhibit completely different peptide binding specificities . Thus, clustering on function rather than sequence of MHC molecules is a more important parameter in CTL biology and immunology. NetMHCpan predicted that the epitope Tp587-95 should bind to BoLA-1*02301 (with a false positive fraction of 0.017) in addition to BoLA-T5, the original molecule that was defined to present the Tp587-95 epitope by parasite infected cells (Table 1). Using both peptide-BoLA binding assays and p-MHC class I tetramers we demonstrate that the Tp587-95 epitope does indeed bind both BoLA molecules and that a single TCR can recognize the Tp587-95 epitope presented by either BoLA molecule (Figure 7). This suggests that as for HLA molecules, it will be possible to categorize BoLA molecules into supertypes based on function. More importantly the functional overlap of BoLA molecules suggest that it should be possible to select a repertoire of T. parva antigens to cover a broader population of cattle than indicted by the limited specificity of CTL responses generated by the live ITM vaccine.
The finding that the Tp587-95 epitope binds to two BoLA molecules was unexpected. What is more surprising was the observation that three of 15 amino acid sequence polymorphisms between the BoLA molecules map to the peptide binding groove (Figure 8B). More specifically they map to the pseudo-sequence of 34 amino acids in HLA molecules, which form the basis of the pan prediction capacity of NetMHCpan. These 34 polymorphic residues have been defined from crystal structure studies as being within 4.0 Angstrom (Å) of the peptide in any of a representative set of HLA-A and HLA-B structures binding a nonameric peptide  and hence are predicted to play a critical role in peptide binding. Leucine in BoLA-T5 is replaced by isoleucine in BoLA-1*02301 at position 94; leucine by alanine at position 96 and arginine by serine at position 113 in the mature protein (highlighted with red asterisk in Figure 8A). Of these substitutions two are conservative in nature. Our data indicates that these changes and the non-conservative amino acid substitution do not have a dramatic influence on the specificity of peptide binding. The NetMHCpan method has, in several benchmark studies, proven to accurately predict the peptide binding motif for molecules characterised by limited or no experimental peptide binding data [16, 53, 54]. In addition, NetMHCpan can differentiate between critical and less critical amino acid substitutions in the MHC pseudo-sequence in situations where the amino acid variations have been observed and used in training the method . As the training data for NetMHCpan increases it should become possible to predict such subtleties in peptide binding to MHC molecules.
An essential feature of the system we have described is that p-MHC class I tetramers are generated over a 48 h time frame using an easy protocol starting from the raw ingredients of recombinant MHC class I HC, synthetic peptide and recombinant β2m. We have also used the BoLA molecules to develop a simple ELISA based method to determine peptide-MHC complex formation. The power of combining such experimental verifications with predictions from NetMHCpan is highlighted by our results and have important implications for research on infectious diseases affecting not only cattle but other animal species including humans. We have continued to improve NetMHCpan using peptide-BoLA binding data (manuscript in preparation) and expect similar improvements to occur through incorporation of data from other MHC projects. These new bovine MHC class I tetramers are sensitive tools that can be used to study the immune responses in cattle e.g. the CD8+ T cells that play a role in Theileria parva immune responses or in other parasitic or viral diseases where CTL contribute to immunity [43, 55]. More BoLA class I alleles could be included in the pipeline as more information about MHC class I expression in cattle of eastern, central and southern Africa is generated.
Bovine leukocyte antigen
Basic research to enable agricultural development
Bovine serum albumin
Cytotoxic T lymphocytes
- E. coli:
East Coast fever
Enzyme linked immunosorbent assay
Fetal bovine serum
Infection and treatment method
Major histocompatibility complex
National Science foundation
Peripheral blood mononuclear cells
Phosphate buffered saline
Peridinin chlorophyll protein
- T. parva:
T-cell growth factor
Tosyl phenylalanyl chloromethyl ketone
T. parva Muguga transformed lymphocytes
Particular thanks go to John Wasilwa for technical assistance with the flow cytometer. Thanks also go to FlowJo for the free licence to run the flow cytometry analysis software given to scientists in Africa. This research was funded through an award (No. 0965346) from the BREAD program of the National Science Foundation (USA); through the CGIAR Research Program on Livestock and Fish; and through the Norman Borlaug Commemorative Research Initiative, an initiative between the Feed the Future program of USAID and USDA-ARS (58-5348-2-117F).
- Chavez-Galan L, Arenas-Del Angel MC, Zenteno E, Chavez R, Lascurain R: Cell death mechanisms induced by cytotoxic lymphocytes. Cell Mol Immunol. 2009, 6: 15-25. 10.1038/cmi.2009.3.PubMed CentralView ArticlePubMedGoogle Scholar
- Rudolph MG, Stanfield RL, Wilson IA: How TCRs bind MHCs, peptides, and coreceptors. Annu Rev Immunol. 2006, 24: 419-466. 10.1146/annurev.immunol.23.021704.115658.View ArticlePubMedGoogle Scholar
- Apostolopoulos V, Lazoura E, Yu M: MHC and MHC-like molecules: structural perspectives on the design of molecular vaccines. Adv Exp Med Biol. 2008, 640: 252-267. 10.1007/978-0-387-09789-3_19.View ArticlePubMedGoogle Scholar
- Pamer E, Cresswell P: Mechanisms of MHC class I–restricted antigen processing. Annu Rev Immunol. 1998, 16: 323-358. 10.1146/annurev.immunol.16.1.323.View ArticlePubMedGoogle Scholar
- Donaldson JG, Williams DB: Intracellular assembly and trafficking of MHC class I molecules. Traffic. 2009, 10: 1745-1752. 10.1111/j.1600-0854.2009.00979.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Kisselev AF, Akopian TN, Woo KM, Goldberg AL: The sizes of peptides generated from protein by mammalian 26 and 20 S proteasomes. Implications for understanding the degradative mechanism and antigen presentation. J Biol Chem. 1999, 274: 3363-3371. 10.1074/jbc.274.6.3363.View ArticlePubMedGoogle Scholar
- Serwold T, Gonzalez F, Kim J, Jacob R, Shastri N: ERAAP customizes peptides for MHC class I molecules in the endoplasmic reticulum. Nature. 2002, 419: 480-483. 10.1038/nature01074.View ArticlePubMedGoogle Scholar
- Ortmann B, Copeman J, Lehner PJ, Sadasivan B, Herberg JA, Grandea AG, Riddell SR, Tampe R, Spies T, Trowsdale J, Cresswell P: A critical role for tapasin in the assembly and function of multimeric MHC class I-TAP complexes. Science. 1997, 277: 1306-1309. 10.1126/science.277.5330.1306.View ArticlePubMedGoogle Scholar
- Whelan JA, Dunbar PR, Price DA, Purbhoo MA, Lechner F, Ogg GS, Griffiths G, Phillips RE, Cerundolo V, Sewell AK: Specificity of CTL interactions with peptide-MHC class I tetrameric complexes is temperature dependent. J Immunol. 1999, 163: 4342-4348.PubMedGoogle Scholar
- Wooldridge L, Lissina A, Cole DK, van den Berg HA, Price DA, Sewell AK: Tricks with tetramers: how to get the most from multimeric peptide-MHC. Immunology. 2009, 126: 147-164. 10.1111/j.1365-2567.2008.02848.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Leisner C, Loeth N, Lamberth K, Justesen S, Sylvester-Hvid C, Schmidt EG, Claesson M, Buus S, Stryhn A: One-pot, mix-and-read peptide-MHC tetramers. PLoS One. 2008, 3: e1678-10.1371/journal.pone.0001678.PubMed CentralView ArticlePubMedGoogle Scholar
- Sims S, Willberg C, Klenerman P: MHC-peptide tetramers for the analysis of antigen-specific T cells. Expert Rev Vaccine. 2010, 9: 765-774. 10.1586/erv.10.66.View ArticleGoogle Scholar
- Xu XN, Screaton GR: MHC/peptide tetramer-based studies of T cell function. J Immunol Methods. 2002, 268: 21-28. 10.1016/S0022-1759(02)00196-5.View ArticlePubMedGoogle Scholar
- Deres K, Beck W, Faath S, Jung G, Rammensee HG: MHC/peptide binding studies indicate hierarchy of anchor residues. Cell Immunol. 1993, 151: 158-167. 10.1006/cimm.1993.1228.View ArticlePubMedGoogle Scholar
- Neefjes J, Jongsma ML, Paul P, Bakke O: Towards a systems understanding of MHC class I and MHC class II antigen presentation. Nat Rev Immunol. 2011, 11: 823-836.PubMedGoogle Scholar
- Karosiene E, Lundegaard C, Lund O, Nielsen M: NetMHCcons: a consensus method for the major histocompatibility complex class I predictions. Immunogenetics. 2012, 64: 177-186. 10.1007/s00251-011-0579-8.View ArticlePubMedGoogle Scholar
- Nielsen M, Lundegaard C, Blicher T, Lamberth K, Harndahl M, Justesen S, Roder G, Peters B, Sette A, Lund O, Buus S: NetMHCpan, a method for quantitative predictions of peptide binding to any HLA-A and -B locus protein of known sequence. PLoS One. 2007, 2: e796-10.1371/journal.pone.0000796.PubMed CentralView ArticlePubMedGoogle Scholar
- Nielsen M, Lundegaard C, Worning P, Lauemoller SL, Lamberth K, Buus S, Brunak S, Lund O: Reliable prediction of T-cell epitopes using neural networks with novel sequence representations. Protein Sci. 2003, 12: 1007-1017. 10.1110/ps.0239403.PubMed CentralView ArticlePubMedGoogle Scholar
- Lundegaard C, Lund O, Buus S, Nielsen M: Major histocompatibility complex class I binding predictions as a tool in epitope discovery. Immunology. 2010, 130: 309-318. 10.1111/j.1365-2567.2010.03300.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Kurtulus S, Hildeman D: Assessment of CD4(+) and CD8 (+) T cell responses using MHC class I and II tetramers. Methods Mol Biol. 2013, 979: 71-79. 10.1007/978-1-62703-290-2_8.PubMed CentralView ArticlePubMedGoogle Scholar
- He XH, Xu LH, Liu Y: Procedure for preparing peptide-major histocompatibility complex tetramers for direct quantification of antigen-specific cytotoxic T lymphocytes. World J Gastroenterol. 2005, 11: 4180-4187.PubMedGoogle Scholar
- Bishop R, Musoke A, Morzaria S, Gardner M, Nene V: Theileria: intracellular protozoan parasites of wild and domestic ruminants transmitted by ixodid ticks. Parasitology. 2004, 129 (Suppl): S271-S283.View ArticlePubMedGoogle Scholar
- Gachohi J, Skilton R, Hansen F, Ngumi P, Kitala P: Epidemiology of east coast fever (Theileria parva infection) in Kenya: past, present and the future. Parasit Vectors. 2012, 5: 194-10.1186/1756-3305-5-194.PubMed CentralView ArticlePubMedGoogle Scholar
- Di Giulio G, Lynen G, Morzaria S, Oura C, Bishop R: Live immunization against East Coast fever–current status. Trends Parasitol. 2009, 25: 85-92. 10.1016/j.pt.2008.11.007.View ArticlePubMedGoogle Scholar
- Morrison WI: Progress towards understanding the immunobiology of Theileria parasites. Parasitology. 2009, 136: 1415-1426. 10.1017/S0031182009990916.View ArticlePubMedGoogle Scholar
- Global Alliance for Livestock Veterinary Medicines (GALVmed).http://www.galvmed.org/2012/03/east-cost-fever-vaccine-launch-arusha-tanzania/,
- MacHugh ND, Connelley T, Graham SP, Pelle R, Formisano P, Taracha EL, Ellis SA, McKeever DJ, Burrells A, Morrison WI: 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. 10.1002/eji.200939227.PubMed CentralView ArticlePubMedGoogle Scholar
- Graham SP, Pelle R, Honda Y, Mwangi DM, Tonukari NJ, Yamage M, Glew EJ, de Villiers EP, Shah T, Bishop R, Abuya E, Awino E, Gachanja J, Luyai AE, Mbwika F, Muthiani AM, Ndegwa DM, Njahira M, Nyanjui JK, Onono FO, Osaso J, Saya RM, Wildmann C, Fraser CM, Maudlin I, Gardner MJ, Morzaria SP, Loosmore S, Gilbert SC, Audonnet JC: Theileria parva candidate vaccine antigens recognized by immune bovine cytotoxic T lymphocytes. Proc Natl Acad Sci USA. 2006, 103: 3286-3291. 10.1073/pnas.0511273103.PubMed CentralView ArticlePubMedGoogle Scholar
- Graham SP, Pelle R, Yamage M, Mwangi DM, Honda Y, Mwakubambanya RS, de Villiers EP, Abuya E, Awino E, Gachanja J, Mbwika F, Muthiani AM, Muriuki C, Nyanjui JK, Onono FO, Osaso J, Riitho V, Saya RM, Ellis SA, McKeever DJ, MacHugh ND, Gilbert SC, Audonnet JC, Morrison WI, van der Bruggen P, Taracha EL: Characterization of the fine specificity of bovine CD8 T-cell responses to defined antigens from the protozoan parasite Theileria parva. Infect Immun. 2008, 76: 685-694. 10.1128/IAI.01244-07.PubMed CentralView ArticlePubMedGoogle Scholar
- Nene V, Svitek N, Toye P, Golde WT, Barlow J, Harndahl M, Buus S, Nielsen M: Designing bovine T cell vaccines via reverse immunology. Ticks Tick Borne Dis. 2012, 3: 188-192. 10.1016/j.ttbdis.2011.12.001.View ArticlePubMedGoogle Scholar
- Ferre H, Ruffet E, Blicher T, Sylvester-Hvid C, Nielsen LL, Hobley TJ, Thomas OR, Buus S: Purification of correctly oxidized MHC class I heavy-chain molecules under denaturing conditions: a novel strategy exploiting disulfide assisted protein folding. Protein Sci. 2003, 12: 551-559. 10.1110/ps.0233003.PubMed CentralView ArticlePubMedGoogle Scholar
- Bernard A, Payton M: Fermentation and growth of Escherichia coli for optimal protein production. Curr Protoc Protein Sci. 2001, Chapter 5: 3-PubMedGoogle Scholar
- Ostergaard Pedersen L, Nissen MH, Hansen NJ, Nielsen LL, Lauenmoller SL, Blicher T, Nansen A, Sylvester-Hvid C, Thromsen AR, Buus S: Efficient assembly of recombinant major histocompatibility complex class I molecules with preformed disulfide bonds. Eur J Immunol. 2001, 31: 2986-2996. 10.1002/1521-4141(2001010)31:10<2986::AID-IMMU2986>3.0.CO;2-R.View ArticlePubMedGoogle Scholar
- Pedersen LE, Harndahl M, Rasmussen M, Lamberth K, Golde WT, Lund O, Nielsen M, Buus S: Porcine major histocompatibility complex (MHC) class I molecules and analysis of their peptide-binding specificities. Immunogenetics. 2011, 63: 821-834. 10.1007/s00251-011-0555-3.PubMed CentralView ArticlePubMedGoogle Scholar
- Sylvester-Hvid C, Kristensen N, Blicher T, Ferre H, Lauemoller SL, Wolf XA, Lamberth K, Nissen MH, Pedersen LO, Buus S: Establishment of a quantitative ELISA capable of determining peptide - MHC class I interaction. Tissue Antigens. 2002, 59: 251-258. 10.1034/j.1399-0039.2002.590402.x.View ArticlePubMedGoogle Scholar
- Steinaa L, Saya R, Awino E, Toye P: Cytotoxic T lymphocytes from cattle immunized against Theileria parva exhibit pronounced cross-reactivity among different strain-specific epitopes of the Tp1 antigen. Vet Immunol Immunopathol. 2012, 145: 571-581. 10.1016/j.vetimm.2011.12.007.View ArticlePubMedGoogle Scholar
- Hoof I, Peters B, Sidney J, Pedersen LE, Sette A, Lund O, Buus S, Nielsen M: NetMHCpan, a method for MHC class I binding prediction beyond humans. Immunogenetics. 2009, 61: 1-13. 10.1007/s00251-008-0341-z.PubMed CentralView ArticlePubMedGoogle Scholar
- NetMHCpan 2.8. [http://www.cbs.dtu.dk/services/NetMHCpan-2.8/]
- Thomsen M, Lundegaard C, Buus S, Lund O, Nielsen M: MHCcluster, a method for functional clustering of MHC molecules. Immunogenetics. 2013, 65: 655-665. 10.1007/s00251-013-0714-9.PubMed CentralView ArticlePubMedGoogle Scholar
- Robinson J, Halliwell JA, McWilliam H, Lopez R, Parham P, Marsh SG: The IMGT/HLA database. Nucleic Acids Res. 2013, 41: D1222-D1227. 10.1093/nar/gks949.PubMed CentralView ArticlePubMedGoogle Scholar
- Robinson J, Halliwell JA, McWilliam H, Lopez R, Marsh SG: IPD–the Immuno Polymorphism Database. Nucleic Acids Res. 2013, 41: D1234-D1240. 10.1093/nar/gks1140.PubMed CentralView ArticlePubMedGoogle Scholar
- Robinson J, Waller MJ, Fail SC, Marsh SG: The IMGT/HLA and IPD databases. Hum Mutat. 2006, 27: 1192-1199. 10.1002/humu.20406.View ArticlePubMedGoogle Scholar
- Li X, Liu J, Qi J, Gao F, Li Q, Li X, Zhang N, Xia C, Gao GF: Two distinct conformations of a rinderpest virus epitope presented by bovine major histocompatibility complex class I N*01801: a host strategy to present featured peptides. J Virol. 2011, 85: 6038-6048. 10.1128/JVI.00030-11.PubMed CentralView ArticlePubMedGoogle Scholar
- Akoolo L, Pelle R, Saya R, Awino E, Nyanjui J, Taracha EL, Kanyari P, Mwangi DM, Graham SP: Evaluation of the recognition of Theileria parva vaccine candidate antigens by cytotoxic T lymphocytes from Zebu cattle. Vet Immunol Immunopathol. 2008, 121: 216-221. 10.1016/j.vetimm.2007.09.012.View ArticlePubMedGoogle Scholar
- NetMHCpan 2.4. [http://www.galvmed.org/2012/03/east-cost-fever-vaccine-launch-arusha-tanzania/]
- Hirano A, Brown WC, Trigona W, Tuo W, Estes DM: Kinetics of expression and subset distribution of the TNF superfamily members CD40 ligand and Fas ligand on T lymphocytes in cattle. Vet Immunol Immunopathol. 1998, 61: 251-263. 10.1016/S0165-2427(97)00155-4.View ArticlePubMedGoogle Scholar
- Hogg AE, Parsons K, Taylor G, Worth A, Beverley P, Howard CJ, Villarreal-Ramos B: Characterization of age-related changes in bovine CD8+ T-cells. Vet Immunol Immunopathol. 2011, 140: 47-54. 10.1016/j.vetimm.2010.11.012.View ArticlePubMedGoogle Scholar
- Emery DL, Eugui EM, Nelson RT, Tenywa T: Cell-mediated immune responses to Theileria parva (East Coast fever) during immunization and lethal infections in cattle. Immunology. 1981, 43: 323-336.PubMed CentralPubMedGoogle Scholar
- Emery DL, Kar SK: Immune responses of cattle to Theileria parva (East Coast fever): specificity of cytotoxic cells generated in vivo and in vitro. Immunology. 1983, 48: 723-731.PubMed CentralPubMedGoogle Scholar
- McKeever DJ, Taracha EL, Innes EL, MacHugh ND, Awino E, Goddeeris BM, Morrison WI: Adoptive transfer of immunity to Theileria parva in the CD8+ fraction of responding efferent lymph. Proc Natl Acad Sci USA. 1994, 91: 1959-1963. 10.1073/pnas.91.5.1959.PubMed CentralView ArticlePubMedGoogle Scholar
- MHCcluster. [http://www.cbs.dtu.dk/services/MHCcluster-2.0]
- Thomsen MC, Nielsen M: Seq2Logo: a method for construction and visualization of amino acid binding motifs and sequence profiles including sequence weighting, pseudo counts and two-sided representation of amino acid enrichment and depletion. Nucleic Acids Res. 2012, 40: W281-W287. 10.1093/nar/gks469.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang L, Udaka K, Mamitsuka H, Zhu S: Toward more accurate pan-specific MHC-peptide binding prediction: a review of current methods and tools. Brief Bioinform. 2012, 13: 350-364. 10.1093/bib/bbr060.View ArticlePubMedGoogle Scholar
- Zhang H, Lund O, Nielsen M: The PickPocket method for predicting binding specificities for receptors based on receptor pocket similarities: application to MHC-peptide binding. Bioinformatics. 2009, 25: 1293-1299. 10.1093/bioinformatics/btp137.PubMed CentralView ArticlePubMedGoogle Scholar
- Tenaya IW, Heel K, Stumbles PA, Wilcox GE: Flow cytometric analysis of lymphocyte subset kinetics in Bali cattle experimentally infected with Jembrana disease virus. Vet Immunol Immunopathol. 2012, 149: 167-176. 10.1016/j.vetimm.2012.06.013.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.