Presence of tumor-infiltrating lymphocytes (TIL) predicts survival in many cancer types. In HPV-driven cancers, cervical and oropharyngeal squamous cell carcinomas (CSCC and OPSCC, respectively), numbers of infiltrating T cells, particularly CD8+ T cells, and presentation of HPV E6/E7 epitopes are associated with improved prognosis. Endoplasmic reticulum aminopeptidase 1 (ERAP1) regulates the presented peptide repertoire, trimming peptide precursors prior to MHC I loading. ERAP1 is polymorphic, and allotypic variation of ERAP1 enzyme activity has an impact on the presented peptide repertoire. Individual SNPs are associated with incidence and outcome in a number of diseases, including CSCC. Here, we highlight the requirement for ERAP1 in the generation of HPV E6/E7 epitopes and show that the functional activity of ERAP1 allotype combinations identified in OPSCC correlate with tumor-infiltrating CD8+ T-cell (CD8)/TIL (CD8/TIL) status of the tumor. Functional analyses revealed that ERAP1 allotype combinations associated with CD8/TILlow tumors have a reduced capacity to generate both a model antigen SIINFEHL and the HPV-16 E782-90 epitope LLMGTLGIV from N-terminally extended precursor peptides. In contrast, ERAP1 allotypes from CD8/TILhigh tumors generated the epitopes efficiently. These data reveal that ERAP1 function correlates with CD8/TIL numbers and, by implication, prognosis, suggesting that the presentation of HPV-16 epitopes at the cell surface, resulting in an anti-HPV T-cell response, may depend on the ERAP1 allotype combinations expressed within an individual.

Lymphocyte infiltration of tumors, in particular by CD8+ T cells, is associated with better clinical outcome in many cancers (1–6). Although this tumor infiltration is driven by a number of factors such as expression of chemokines and their receptors, the presentation of immunogenic peptide epitopes is paramount. The identity of immunogenic peptide epitopes and hence, specificity of tumor-infiltrating CD8+ T cells, has not been determined in most instances, but it is likely that T-cell specificity and efficacy are linked, as they are in other diseases. However, in human papillomavirus (HPV)-driven cancers such as cervical squamous cell carcinoma (CSCC) and oropharyngeal squamous cell carcinoma (OPSCC), where HPV infection accounts for >95% and 40% to 80% of cases in the western world respectively, peripheral and tumor-infiltrating CD8+ T cells have been shown to recognize epitopes derived from the E6/E7 oncoproteins expressed by the virus (5, 7, 8). Furthermore, more tumor-infiltrating CD8+ T cells (tumor-infiltrating lymphocytes, TIL) in OPSCC are associated with improved survival (5, 9).

The presentation of peptide epitopes to CD8+ T cells is dependent on the antigen-processing pathway. Several components of this pathway are polymorphic, including HLA, endoplasmic reticulum aminopeptidase 1 (ERAP1), the transporter associated with antigen processing (TAP), and to a lesser extent Tapasin (10). ERAP1 trims N-terminally extended peptide epitopes to the correct length for optimal binding to HLA and therefore serves as an editor of the peptide repertoire presented at the cell surface. We have previously shown that ERAP1 allotypes arising from its polymorphism can differ in their peptide trimming function, defined by both overall activity and specificity for the N-terminal amino acid of peptide epitope precursors (11). In CSCC, genome-wide association studies (GWAS) showed nonsynonymous exonic SNPs of ERAP1 were associated with disease incidence (12), with 2 of these SNPs associated with disease-free and overall survival (13). Given that CSCC is predominantly driven by HPV infection and, similar to OPSCC, shows improved survival with more TIL (14), the ability to control HPV+ tumors is likely mediated by CD8+ T cells specific for HPV-derived antigens (4, 5). Thus, the ability to generate and present peptide epitopes derived from viral proteins, such as E6 and E7, likely drives antitumor responses.

In this study, having observed similarities in SNP variants between CSCC and a cohort of HPV+ patients with OPSCC, we tested the hypothesis that allotypic variation in ERAP1 might underpin the extent of CD8+ T-cell tumor infiltration via its impact on processing of tumor-specific (HPV-derived) epitopes. In individual patients, we found that ERAP1 function correlates with tumor infiltrating CD8+ T cell (CD8)/TIL (CD8/TIL) status in resected tumors: ERAP1 in CD8/TILlow tumors had reduced capacity to generate the HPV E7-derived epitope, LLMGTLGIV, suggesting that the lack of CD8+ T-cell infiltration may be linked to the ability to generate and present virally derived peptide epitopes. These data suggest that ERAP1 function may identify “high-risk” patients with OPSCC with poor prognosis.

Patients, HPV classification, and CD8/TIL scoring

Patients previously classified for HPV+ tumor status using HPV ISH and p16 immunochemistry were used in this study (REC reference 09/H0501/90). Tumor histology was previously characterized for each patient (9) and reviewed for this study by pathologist GJT according to the 1998 UK Royal College of Pathologist Guidelines. Lymphocyte infiltrate (TIL) was determined and scored under low-power magnification (×2.5 objective) as high (diffuse; present in >80% of tumor/stroma), moderate (patchy; present in 20%–80% of tumor/stroma), or low (weak/absent; present in <20% of tumor/stroma; ref. 9). In addition, CD8+ T-cell infiltration was scored by IHC and enumeration of CD8+ T cells was expressed as an average of 10 high-power fields (9). Images were collected using a CKX41 inverted microscope with DP-22 color camera running under CellSens (Version2) software and ×10 or ×40 achromat objective lenses (all Olympus; Supplementary Fig. S1A). Grouping was defined by the tertile values: high (>39 CD8+ T cells/field), moderate (>13 and <39 CD8+ T cells/field), and low (<13 CD8+ T cells/field; Supplementary Fig. S1B).

ERAP1 isolation and identification

ERAP1 was amplified from whole blood samples and cloned into the pcDNA3.1 expression vector (Life Technologies) as described previously (11). ERAP1 clones from each patient were sequenced (Source Bioscience Ltd.) and compared with the wild-type ERAP1 sequence (UniProt – Q9NZ08) to identify the presence of polymorphic variants using DNADynamo alignment software (Blue Tractor Software). To further validate the presence of ERAP1 polymorphisms, gDNA was extracted from whole blood using DNeasy Blood and Tissue Kit (Qiagen) and targeted ERAP1 exons amplified and sequenced using the following primers: M349V 5-GATGATTGTTTGGGAGAATG-3′ and 5′-GCCTTATATGTCCACGCT-3′; K528R 5′-CTTTCAGTCTGGTGCTATGG-3′ and 5′-GTGTGATGGCTGGGGACATC-3′; D575N 5′-GCTCTACTCAAGGAGTCCAAG-3′ and 5′-CGATTTTCAATCTGATCCC-3′; and R725Q/Q730E 5′-CCTGACTTTAAGTAGATGG-3′ and 5′-GGTTATTTAGTAAGGACTGAC-3′.

DNA constructs

The pcDNA3.1 minigenes (ES)-SHL8, (ES)-X-SHL8, H2-Kb, ERAP1*002 as well as the nonfunctional active site ERAP1 mutant, E320A, have been described previously (11). The HPV-16 E7 epitope-specific construct (ES)-LV9 was generated by the insertion of oligonucleotides into the EcoRI/XbaI sites of pcDNA3.1, encoding the ER translocation signal sequence (ES) followed by the amino acid sequence LLMGTLGIV (LV9). The (ES)-X-LV9 constructs were generated by site directed mutagenesis to incorporate a single amino acid into the (ES)-LV9 construct, using the following primers: 5′-GCGGCAGTCTGCAGCGCGNNSCTGTTAATGGGCACACTAG-3′ and 5′- CTAGTGTGCCCATTAACAGSNNCGCGCTGCAGACTGCCGC-3′, where N = any nucleotide and S = C or G. Constructs were sequence verified, and the most frequent codon for each amino chosen for use where possible.

Dextramer detection of HPV-specific T cells

Patient tumor samples were mechanically disaggregated into single-cell suspensions before being analyzed for the presence of HPV16-specific T cells by flow cytometry. Dextramers specific for HPV16 E618-26 (KLPQLCTEL; KL9), E629-38 (TIHDIILEV; TV10), and E782-90 (LLMGTLGIV; LV9) were obtained from Immudex and directly conjugated with PE, APC, and FITC fluorophores, respectively. Tumor cells were preincubated with 50 nmol/L protein kinase inhibitor dasatanib for 30 minutes at 37°C, before the addition of 10 μL dextramer and incubated for a further 30 minutes at 4°C. After two washes, surface antibodies for CD8 (anti-CD8-Pacific Blue, clone RPA-T8; BD Biosciences) and CD3 (anti-CD3-AmCyan, clone SK7; BD Biosciences) were incubated for 20 minutes at 4°C before washing. Cells were acquired on BD FACS Canto II (BD Biosciences) and analyzed using Flow Jo Software (TreeStar Inc.).

Generation of BE7A2Z hybridoma

The generation of LacZ-inducible T-cell hybridomas has been described previously (15). Briefly, the BE7A2Z hybridoma was generated by immunizing HLA-A*0201 transgenic HHD mice with LLMGTLGIV-pulsed BMDC and splenocytes fused with the BWZ.36/CD8α cell line as described previously (15). To determine peptide specificity of the hybridoma, splenocytes were pulsed with 10 μmol/L LV9 or 10 μmol/L irrelevant peptide and hybridomas assessed for IFNγ production. To determine sensitivity of subsequent LV9-specific clones, 293T cells expressing endogenous HLA-A2 were pulsed with 10 nmol/L LV9 peptide and cocultured with hybridoma clones for 24 hours. Each hybridoma clone sensitivity to LV9 was assessed by measurement of intracellular LacZ activity by chlorophenol red-β-d-galactopyranoside (CPRG; Roche).

Cell lines, transfection, and T-cell activation assay

Culture conditions for ERAP1 knockout HEK293T (293T E1KO) and T-cell hybridoma cells have been described before (16). 293T E1KO cells were cultured from freshly thawed aliquots from passage 2 (generated from original 293T cell line in 2013) and used over 2 weeks for the indicated experiments. Cells were authenticated by STR analysis and were confirmed Mycoplasma free. For assessment of ERAP1 trimming activity, 293T E1KO cells were transfected with DNA using FuGENE 6 Transfection Reagent (Promega). For assessment of single amino acid trimming specificity, 293T E1KO cells were transfected with 0.1 μg total DNA, consisting of 0.05 μg E320A nonfunctional ERAP1 or 0.025 μg of each patient ERAP1 allotype, together with 0.025 μg H2-Kb (SHL8) or HLA-A*0201 (LV9) and 0.025 μg X-SHL8 or X-LV9. For assessment of LV9 extended precursor trimming (ED/D-LV9), 293T E1KO cells were transfected with 1 μg total DNA, consisting of 0.5 μg E320A or 0.25 μg each patient ERAP1 allotype, together with 0.25 μg ED/D-LV9 or LV9 and 0.25 μg HLA-A*0201. All transfected cells were incubated for 24 hours at 37°C. Presentation of final peptide SHL8 or LV9 and activation of the respective LacZ-inducible T-cell hybridoma, B3Z or BE7A2Z, was assessed by measurement of intracellular LacZ activity by CPRG.

Statistical analysis

Two or one-way ANOVA with Dunnett posttest was performed for analysis of differences between multiple groups and controls (GraphPad Prism).

ERAP1 allotype sequences in patients with OPSCC

ERAP1 polymorphism is associated with disease risk and disease-free survival in HPV-induced CSCC (12, 13). For OPSCC, 40% to 80% of cases in the western world are also associated with HPV infection and within this group of patients, the amount of tumor-infiltrating lymphocytes (TIL) observed in resected tumors correlates with disease survival (9, 17). We therefore investigated whether patients with OPSCC shared the same ERAP1 polymorphisms observed in CSCC and whether the allotypes expressed in patients correlated with numbers of tumor-infiltrating CD8+ T cell (CD8)/TIL. We sequenced both ERAP1 copies expressed in a cohort of 26 patients with OPSCC whose tumors had been scored with high, moderate, or low numbers of CD8/TIL. Two of the CSCC-associated ERAP1 SNPs (P127 and E730) were detected in the majority of patients, regardless of CD8/TIL status, with E730 being more frequently observed (E730 = 9 of 10 CD8/TILhigh, 6 of 7 CD8/TILmod and 7 of 8 CD8/TILlow patients; P127 = 8 of 10 CD8/TILhigh, 2 of 7 CD8/TILmod, and 4 of 8 CD8/TILlow patients; Table 1). The number of patients that contained both P127 and E730 varied between CD8/TIL groups with CD8/TILhigh having the greatest proportion (8/10 vs. 2/7 CD8/TILmod and 4/8 CD8/TILlow; Table 1). Examination of the 2 individual ERAP1 allotypes from each patient that contained either of the P127 or E730 variants in a single molecule (CD8/TILhigh = 13/20 total allotypes; CD8/TILmod = 11/14; CD8/TILlow = 12/16), revealed that a proportion did not contain both of the variants (P127 and E730) in the same molecule, with CD8/TILmod and CD8/TILlow having the lowest frequency (8/13 TILhigh, 3/11 TILmod, 4/12 TILlow; Table 1).

Table 1.

ERAP1 allotype identity in patients with OPSCC

Amino acid at indicated position
56127346349528575725730
Patient numberCD8/TIL statusERAP1 allotypeERGMKDRQ
High *001 P V R N Q E 
  *002 
High *011 R E 
  *016 V Q E 
High *001 P V R N Q E 
  *018 E 
High *001 P V R N Q E 
  *002 
High *001 P V R N Q E 
  *002 
High *002 
  *002 
High *002 
  *014 K P R E 
High *001 P V R N Q E 
  *011 R E 
High *011 R E 
  *015 P R E 
10 High *001 P V R N Q E 
  *002 
11 Moderate *015 P R E 
  *015 P R E 
12 Moderate *001 P V R N Q E 
  *011 R E 
13 Moderate *002 
  *017 N 
14 Moderate *011 R E 
  *018 E 
15 Moderate *002 
  *011 R E 
16 Moderate *011 R E 
  *011 R E 
17 Moderate *016 V Q E 
  *018 E 
18 Low *007 R Q 
  *011 R E 
19 Low *001 P V R N Q E 
  *019 D R E 
20 Low *001 P V R N Q E 
  *001 P V R N Q E 
21 Low *011 R E 
  *011 R E 
22 Low *013 P 
  *019 D R E 
23 Low *002 
  *002 
24 Low *002 
  *011 R E 
25 Low *011 R E 
  *011 P R E 
Amino acid at indicated position
56127346349528575725730
Patient numberCD8/TIL statusERAP1 allotypeERGMKDRQ
High *001 P V R N Q E 
  *002 
High *011 R E 
  *016 V Q E 
High *001 P V R N Q E 
  *018 E 
High *001 P V R N Q E 
  *002 
High *001 P V R N Q E 
  *002 
High *002 
  *002 
High *002 
  *014 K P R E 
High *001 P V R N Q E 
  *011 R E 
High *011 R E 
  *015 P R E 
10 High *001 P V R N Q E 
  *002 
11 Moderate *015 P R E 
  *015 P R E 
12 Moderate *001 P V R N Q E 
  *011 R E 
13 Moderate *002 
  *017 N 
14 Moderate *011 R E 
  *018 E 
15 Moderate *002 
  *011 R E 
16 Moderate *011 R E 
  *011 R E 
17 Moderate *016 V Q E 
  *018 E 
18 Low *007 R Q 
  *011 R E 
19 Low *001 P V R N Q E 
  *019 D R E 
20 Low *001 P V R N Q E 
  *001 P V R N Q E 
21 Low *011 R E 
  *011 R E 
22 Low *013 P 
  *019 D R E 
23 Low *002 
  *002 
24 Low *002 
  *011 R E 
25 Low *011 R E 
  *011 P R E 

NOTE: Bold type indicates variants at the indicated amino acid position.

Overall, most of the allotypes observed have been previously described (*001; Hap10, *002; Hap2, *007, *011; Hap4, *013; Hap1, *014; Hap7, *015; Hap6, *018; Hap3 and *019; Hap5; Table 1; refs. 16, 18). Two allotypes that have not previously been identified (*016, *017). Comparison of allotypes from different CD8/TIL groups showed that 5 allotypes were only found in one of the 3 groups; however, almost all of these represented a single incidence (Supplementary Table S1). Three allotypes were observed in all CD8/TIL groups and three in two groups (Supplementary Table S1), showing that, for this cohort, analysis of individual ERAP1 allotypes did not correlate with CD8/TIL status.

CD8/TILlow allotype combinations have poor trimming function

We have previously shown, for an autoinflammatory disease, that the strongest association between ERAP1 and disease is achieved by considering ERAP1 allotype pairing, because this best represents the combined enzyme function of the codominantly expressed allotypes (16). We therefore investigated the peptide-editing function of ERAP1 allotype combinations in patients with OPSCC, utilizing the model system previously used to assess trimming activity (11). This system classified ERAP1 allotypes with respect to their differential ability to generate a model epitope, SIINFEHL (SHL8) a modified form of SIINFEKL that stimulates T cells to the same degree as the original peptide (Supplementary Fig. S2), from precursor peptides with single amino acid N-terminal extensions in living cells. The SHL8/Kb complexes are detected at the cell surface by the specific T-cell hybridoma, B3Z. Briefly, ERAP1 allotypes fall into three general classes of activity: those that process the epitope efficiently, those that are hypofunctional, and those that are hyperfunctional and overtrim the epitope, thus destroying it. Further differentiation of function is seen in N-terminal specificity. This assay is preferable to in vitro assays of function using small chromogenic substrates because it evaluates function in vivo, in the context of a functioning antigen-processing pathway, and provides a readout that relates directly to the impact on antigen presentation.

We analyzed the enzyme activity of ERAP1 allotype combinations that were observed exclusively in either high, moderate, or low infiltrated tumors using this assay (Fig. 1A–C). Each allotype tested gave rise to a characteristic activity signature according to their relative ability to cleave X from the precursor substrate (where X is any one of the 20 naturally occurring amino acids). We also measured the activity of allotypes observed in both high and moderate infiltrated tumors (Fig. 1D), high and low infiltrated tumors (Fig. 1E), and in both moderate and low infiltrated tumors (Fig. 1F). The activity toward each substrate was expressed as a percentage of the maximal activity observed with the SHL8 minigene (which does not require processing by ERAP1). By calculating the total activity across all substrates for each allotype combination as AUC, we found that allotype combinations that were associated exclusively with low infiltrated tumors had a significantly reduced activity compared with those associated with high or moderate infiltrated tumors (50% and 55%, respectively; Fig. 2A and B).

Figure 1.

N-terminal amino acid trimming specificity by OPSCC ERAP1 allotype combinations. E1KO 293T cells were transfected with ERAP1 allotype combinations from patients with OPSCC together with H2-Kb and X-SHL8 minigenes representing 20 amino acids and assessed for generation of SHL8 by B3Z activation. OPSCC ERAP1 allotype identity from patients with TIL status are shown in (A) CD8/TILhigh, (B) CD8/TILmoderate, (C) CD8/TILLow, (D) CD8/TILhigh and CD8/TILmoderate, (E) CD8/TILhigh and CD8/TILlow, and (F) CD8/TILmoderate and CD8/TILlow. The relative presentation of trimmed X-SHL8 was compared with that of the maximal response using SHL8, which does not require ERAP1 activity. Data pooled from 4 independent experimental repeats ± SEM.

Figure 1.

N-terminal amino acid trimming specificity by OPSCC ERAP1 allotype combinations. E1KO 293T cells were transfected with ERAP1 allotype combinations from patients with OPSCC together with H2-Kb and X-SHL8 minigenes representing 20 amino acids and assessed for generation of SHL8 by B3Z activation. OPSCC ERAP1 allotype identity from patients with TIL status are shown in (A) CD8/TILhigh, (B) CD8/TILmoderate, (C) CD8/TILLow, (D) CD8/TILhigh and CD8/TILmoderate, (E) CD8/TILhigh and CD8/TILlow, and (F) CD8/TILmoderate and CD8/TILlow. The relative presentation of trimmed X-SHL8 was compared with that of the maximal response using SHL8, which does not require ERAP1 activity. Data pooled from 4 independent experimental repeats ± SEM.

Close modal
Figure 2.

ERAP1 allotype combinations from CD8/TILlow tumors have reduced trimming activity. A, The total trimming activity for all 20 X-SHL8 substrates for each allotype combination were combined and grouped into CD8/TIL status. Overall trimming for each CD8/TIL group is represented as AUC. B, Table showing both the individual AUC for allotype pairs in each CD8/TIL group and the mean for each group. C, Specificity of ERAP1 trimming activity from those found exclusively in CD8/TILhigh, CD8/TILmoderate, or CD8/TILlow tumors towards X-SHL8 amino acids grouped based on physicochemical properties; basic, acidic, polar uncharged, aliphatic, aromatic, and special. Data from A–C, pooled from 4 independent experimental repeats ± SEM (***, P ≤ 0.001; **, P ≤ 0.01; *, P ≤ 0.05).

Figure 2.

ERAP1 allotype combinations from CD8/TILlow tumors have reduced trimming activity. A, The total trimming activity for all 20 X-SHL8 substrates for each allotype combination were combined and grouped into CD8/TIL status. Overall trimming for each CD8/TIL group is represented as AUC. B, Table showing both the individual AUC for allotype pairs in each CD8/TIL group and the mean for each group. C, Specificity of ERAP1 trimming activity from those found exclusively in CD8/TILhigh, CD8/TILmoderate, or CD8/TILlow tumors towards X-SHL8 amino acids grouped based on physicochemical properties; basic, acidic, polar uncharged, aliphatic, aromatic, and special. Data from A–C, pooled from 4 independent experimental repeats ± SEM (***, P ≤ 0.001; **, P ≤ 0.01; *, P ≤ 0.05).

Close modal

We next analyzed the specificity of each ERAP1 allotype combination by considering their activity against N-terminal extensions classified according to their physicochemical properties (aliphatic hydrophobic, aromatic hydrophobic, basic, acidic, uncharged polar). The results in Fig. 2C show that reduced activity of allotype combinations seen in CD8/TILlow tumors is observed across all classes of N-terminal amino acid, closer inspection suggests that precursors with charged, polar uncharged, or aliphatic hydrophobic N-1 amino acids are most effected (charged = 3/5 residues, polar uncharged = 3/4 residues, aliphatic = 2/4 residues were significantly reduced in CD8/TILlow vs. CD8/TILhigh).

Anti-HPV CD8+ T-cell responses in OPSCC tumors

These observations are consistent with the hypothesis that the generation of epitopes requiring N-terminal trimming recruits CD8+ T cells into the tumor mass and therefore that ERAP1 function is a determinant in this process. The data also suggest that the impact of ERAP1 polymorphism on CD8+ T-cell recruitment might be greatest for T cells recognizing epitopes generated from precursors with charged/polar uncharged/aliphatic hydrophobic amino acids. Better prognosis for HPV+ patients with OPSCC and CSCC is associated with more TIL (5, 9, 14) and there is evidence that CD8+ T cells recognizing epitopes derived from the HPV E6/E7 proteins are among TILs in these cancers (5, 8, 19). We therefore scanned the immediate upstream sequence of predicted proteasomal cleavage products of E6 and E7 that contained a putative HLA-A*0201–restricted epitope (Table 2; refs. 20–22). Of the 4 E6 and 4 E7-derived peptides that have been confirmed as epitopes in previous studies (7, 23), 3 consisted of a core sequence with a predicted IC50 for HLA-A*0201 of less than 100 nmol/L. In addition, 5 of the 8 peptides had a predicted 1 or 2 amino acid extension following proteasomal cleavage. However, only 1 of the epitopes had both an IC50 of less than 100 nmol/L and contained an N-terminal extension, HPV E7 82-90, LLMGTLGIV (LV9; Table 2). The N-terminal extension of LV9 is comprised of acidic amino acids (Glu, Asp; ED-LV9) both of which were trimmed less efficiently by CD8/TILlow allotype pairs when analyzed with the X-SHL8 substrate (Fig. 2C). In addition, the orthologue of this sequence from HPV-type 11 has been confirmed as an epitope in infected humans and so we selected LV9 for further analysis (24).

Table 2.

HLA-A*0201–restricted HPV-16 E6/E7 peptide epitopes

ProteinPosition (amino acid)Epitope sequenceMHC IC50 (predicted nmol/L) NetMHCPredicted proteasome cleaved peptide (N-terminal extension) NetCHOP, MAPPP
E6 7–15 AMFQDPQER — RTAMFQDPQER 
E6 18–26 KLPQLCTEL 134.4 PRKLPQLCTEL 
E6 29–38 TIHDIILECV 200.6 QTTIHDIILECV/TIHDIILECV 
E6 42–50 FAFRDLCIV 100.4 DFAFRDLCIV 
E7 7–15 TLHEYMLDL 94.1 TLHEYMLDL 
E7 11–19 (20) YMLDLQPET(T) 10.58 (33.3) YMLDLQPET(T) 
E7 82–90 LLMGTLGIV 37.11 EDLLMGTLGIV/DLLMGTLGIV 
E7 86–93 TLGIVCPI 409.85 TLGIVCPI 
ProteinPosition (amino acid)Epitope sequenceMHC IC50 (predicted nmol/L) NetMHCPredicted proteasome cleaved peptide (N-terminal extension) NetCHOP, MAPPP
E6 7–15 AMFQDPQER — RTAMFQDPQER 
E6 18–26 KLPQLCTEL 134.4 PRKLPQLCTEL 
E6 29–38 TIHDIILECV 200.6 QTTIHDIILECV/TIHDIILECV 
E6 42–50 FAFRDLCIV 100.4 DFAFRDLCIV 
E7 7–15 TLHEYMLDL 94.1 TLHEYMLDL 
E7 11–19 (20) YMLDLQPET(T) 10.58 (33.3) YMLDLQPET(T) 
E7 82–90 LLMGTLGIV 37.11 EDLLMGTLGIV/DLLMGTLGIV 
E7 86–93 TLGIVCPI 409.85 TLGIVCPI 

The precise immunodominance hierarchy of HLA-A*0201–restricted anti-HPV CD8+ T-cell responses has not been well characterized. To further confirm the role of LV9-CD8+ T cells in antitumor immunity, we assessed the presence of anti-HPV responses in 5 HPV+ OPSCC tumors. Investigation of 3 E6/E7 peptide-specific responses, LV9, TV10, and KL9, revealed that responses to all peptides were observed (Fig. 3A and B), with those directed to LV9 and TV10 the most prevalent. Indeed, LV9- and TV10-specific responses were dominant in 2 tumors (#1 and #4 vs. #3 and #5, respectively), and codominant in 1 tumor (#2) (Fig. 3B). These responses confirmed LV9 as a good candidate peptide to investigate the role of ERAP1 in generation of HPV-derived antigenic peptides.

Figure 3.

Presence of anti-HPV E6/E7 epitope CD8+ T-cell responses in HPV+ OPSCC tumors. Tumors from HPV+ patients with OPSCC were assessed for the presence of CD8+ T-cell responses to the HLA-A*0201-restricted HPV E6/E7 epitopes, LLMGTLGIV (LV9), TIHDIILECV (TV10), and KLPQLCTEL (KL9) using dextramers. A, A representative dot plot showing the presence of all 3 reactivities in a tumor. B, Quantification of responses to all three epitopes (LV9, TV10, and KL9) from 5 HPV+ OPSCC tumors.

Figure 3.

Presence of anti-HPV E6/E7 epitope CD8+ T-cell responses in HPV+ OPSCC tumors. Tumors from HPV+ patients with OPSCC were assessed for the presence of CD8+ T-cell responses to the HLA-A*0201-restricted HPV E6/E7 epitopes, LLMGTLGIV (LV9), TIHDIILECV (TV10), and KLPQLCTEL (KL9) using dextramers. A, A representative dot plot showing the presence of all 3 reactivities in a tumor. B, Quantification of responses to all three epitopes (LV9, TV10, and KL9) from 5 HPV+ OPSCC tumors.

Close modal

Cell surface presentation of LV9 epitope is ERAP1 dependent

To investigate the effect of ERAP1 on the generation of LV9, we created a T-cell hybridoma recognizing LV9 restricted by HLA-A*0201 by immunizing HLA-A*0201 transgenic (HHD) mice with peptide-pulsed BMDC and fusing splenocytes harvested at day 7 after the final immunization with BWZ.36 CD8α+ cells (Supplementary Fig. S3A; ref. 15). Several hybridomas were generated and screened on LV9-pulsed splenocytes (Supplementary Fig. S3B). Four hybridomas were found to be peptide-specific and were stimulated by peptide-pulsed HEK293T cells at the nanomolar quantities, with hybridoma #4 (called BE7A2Z) showing the greatest sensitivity and used in subsequent experiments (Supplementary Fig. S3C). To assess the in vivo processing of the LV9 peptide, we transfected ERAP1 knockout HEK293T (E1KO 293T) cells with a minigene encoding the final (LV9) peptide. This construct sensitized both ERAP1-sufficient and E1KO 293T cells, whereas constructs encoding N-terminally extended LV9 (ED-LV9 and D-LV9) were generated only by ERAP1-sufficient cells, showing that in vivo, the processing of LV9 from these 2 precursor peptides is ERAP1 dependent (Supplementary Fig. S3D).

CD8/TILlow allotype combinations have reduced functionality and different peptide specificity

We next measured the activity of ERAP1 allotype pairs that we identified exclusively in CD8/TIL high, moderate, or low groups using X-LV9 minigene encoded substrates, as described above for X-SHL8. Examination of the allotype pairs showed a similar distribution of activities as for X-SHL8, with allotype pairs from CD8/TILhigh and CD8/TILmod patients having a significantly higher overall activity than those from CD8/TILlow when the AUC was calculated (CD8/TILhigh = 45% and CD8/TILmod = 35% greater; Fig. 4A and B). Original data for the 9 allotype combinations from CD8/TILhigh, CD8/TILmod, and CD8/TILlow groups are shown in Supplementary Fig. S4A–S4C, respectively. Ala-, Cys-, Met-, Thr-, and Tyr-LV9 trimming by ERAP1 had a high background for E320A ERAP1, making it difficult to confirm the amount of trimming of each allotype pair, consequently these data are not included. When we analyzed the specificity of N-terminus preference, we found significantly reduced responses for CD8/TILlow compared with the other two CD8/TIL groups among the charged, polar uncharged, and special N-terminal amino acid extensions (Fig. 4C). The trimming of Lys, Asp, and Glu was more efficient in the context of the LV9 peptide backbone compared with SHL8 for all CD8/TIL groups (Lys = 15%–50% for SHL8, 60%–90% for LV9, Asp = 15%–30% for SHL8, 50%–70% for LV9 and Glu = 30%–50% for SHL8, 85%–90% for LV9; Figs. 2 and 4). Conversely, the aliphatic hydrophobic amino acids Ile and Leu were trimmed less efficiently in the context of LV9 than that observed for SHL8 (Ile = ∼15%–40% SHL8, 5%–10% LV9, and Leu = 30%–60% SHL8, 10%–20% LV9; Figs. 2 and 4). These results indicate that although there are differences in amino acid specificity between substrates X-SHL8 and X-LV9, the total activity of CD8/TILlow ERAP1 allotype pairs is reduced compared with that from CD8/TILhigh patients.

Figure 4.

Amino acid specificity of ERAP1 allotype combinations generating the HPV E782-90 epitope. E1KO 293T cells were transfected with ERAP1 allotype combinations found in CD8/TILhigh, CD8/TILmoderate, and CD8/TILlow tumors, together with X-LV9 minigenes representing 15 amino acids and assessed for generation of LV9 by BE7A2Z activation. The relative presentation of trimmed X-LV9 was compared with that of the maximal response using LV9, which does not require ERAP1 activity. A, Heatmap representation of % maximal LV9 response towards each amino acid from each allotype combination. Blue represents poor trimming and red represents efficient trimming activity. B, Total trimming activity of all 15 X-LV9 substrates from each ERAP1 allotype combination grouped into CD8/TIL status and represented as AUC. C, Specificity of ERAP1 trimming activity from those found exclusively in CD8/TILhigh, CD8/TILmoderate, or CD8/TILlow tumors towards X-LV9 amino acids grouped on the basis of physicochemical properties; basic, acidic, polar uncharged, aliphatic, aromatic, and special. Data pooled from 5 independent experimental repeats ± SEM (****, P ≤ 0.0001; **, P ≤ 0.01; *, P ≤ 0.05).

Figure 4.

Amino acid specificity of ERAP1 allotype combinations generating the HPV E782-90 epitope. E1KO 293T cells were transfected with ERAP1 allotype combinations found in CD8/TILhigh, CD8/TILmoderate, and CD8/TILlow tumors, together with X-LV9 minigenes representing 15 amino acids and assessed for generation of LV9 by BE7A2Z activation. The relative presentation of trimmed X-LV9 was compared with that of the maximal response using LV9, which does not require ERAP1 activity. A, Heatmap representation of % maximal LV9 response towards each amino acid from each allotype combination. Blue represents poor trimming and red represents efficient trimming activity. B, Total trimming activity of all 15 X-LV9 substrates from each ERAP1 allotype combination grouped into CD8/TIL status and represented as AUC. C, Specificity of ERAP1 trimming activity from those found exclusively in CD8/TILhigh, CD8/TILmoderate, or CD8/TILlow tumors towards X-LV9 amino acids grouped on the basis of physicochemical properties; basic, acidic, polar uncharged, aliphatic, aromatic, and special. Data pooled from 5 independent experimental repeats ± SEM (****, P ≤ 0.0001; **, P ≤ 0.01; *, P ≤ 0.05).

Close modal

Generation of LV9 peptide epitope is impaired in CD8/TILlow ERAP1 allotype pairs

The best prognosis for HPV+ OPSCC is associated with high infiltration of CD8+ T cells into the tumor (9), many of which are HPV specific (Fig. 3; refs. 5, 19, 25). With the reduced trimming capacity observed for CD8/TILlow ERAP1 pairs, the ability to generate HPV-specific epitopes may be impaired preventing the activation and subsequent migration and infiltration of CD8+ T cells into the tumor. To investigate the impact of ERAP1 on the generation of HPV-derived epitopes, we tested the ability of ERAP1 allotype combinations to generate the LV9 epitope from the minigene-encoded putative natural precursor (ED-LV9) across a range of effector:target ratios. Example results obtained from the allotype pair *001 + *002 found in CD8/TILhigh patients (#1, 4, 5, 10), *002 + *017 found in a CD8/TILmod patient (#13), and *007 + *011 from a CD8/TILlow patient (#18) are shown in Fig. 5A. As a negative control, cells were transfected with the test precursor sequence together with an inactive mutant of *002 ERAP1 allotype (E320A). These experiments show a lower efficiency of CD8/TILlow-associated allotype pairs in generating LV9 from both the single (D-) and double (ED-) N-terminal extended precursors compared with CD8/TILhigh, with ED-LV9 having the greatest reduction. When activity data from all allotype pairs observed were grouped according to the CD8/TIL status (Fig. 5B–E), it was evident that allotypes from CD8/TILlow patients were less able to generate LV9 from either precursor, D-LV9 or ED-LV9, than allotypes from CD8/TILhigh and from ED-LV9 compared with CD8/TILmod. CD8/TILhigh pairs showed a 33% and 66% increase in the generation of LV9 from D-LV9 and ED-LV9 precursors, respectively (Fig. 5B–E).

Figure 5.

Generation of HPV E782-90 from natural precursors is impaired in CD8/TILlow ERAP1 allotype combinations. E1KO 293T cells were transfected with ERAP1 allotype combinations from CD8/TILhigh, CD8/TILmoderate, and CD8/TILlow or the nonfunctional active site mutant, E320A, together with minigenes encoding final LV9 or extended precursors D-LV9 or ED-LV9 and assessed for trimming activity by the activation of BE7A2Z. A, Representative line graphs showing the trimming activity of selected ERAP1 allotype combinations from the 3 CD8/TIL groups. B and C, The relative maximum LV9 response of ERAP1 allotype combinations from each group, where each symbol represents an individual ERAP1 allotype combination transfection. Trimming activity toward a single amino acid extension (D-LV9; B) and a double amino acid extension (ED-LV9; ****, P ≤ 0.0001; **, P ≤ 0.01; C). D and E, The relative maximum LV9 response of each individual ERAP1 allotype combination from CD8/TILhigh (black), CD8/TILmoderate (gray), and CD8/TILlow (white) toward D-LV9 (D) and ED-LV9 (E). Data pooled from three independent experimental repeats ± SEM.

Figure 5.

Generation of HPV E782-90 from natural precursors is impaired in CD8/TILlow ERAP1 allotype combinations. E1KO 293T cells were transfected with ERAP1 allotype combinations from CD8/TILhigh, CD8/TILmoderate, and CD8/TILlow or the nonfunctional active site mutant, E320A, together with minigenes encoding final LV9 or extended precursors D-LV9 or ED-LV9 and assessed for trimming activity by the activation of BE7A2Z. A, Representative line graphs showing the trimming activity of selected ERAP1 allotype combinations from the 3 CD8/TIL groups. B and C, The relative maximum LV9 response of ERAP1 allotype combinations from each group, where each symbol represents an individual ERAP1 allotype combination transfection. Trimming activity toward a single amino acid extension (D-LV9; B) and a double amino acid extension (ED-LV9; ****, P ≤ 0.0001; **, P ≤ 0.01; C). D and E, The relative maximum LV9 response of each individual ERAP1 allotype combination from CD8/TILhigh (black), CD8/TILmoderate (gray), and CD8/TILlow (white) toward D-LV9 (D) and ED-LV9 (E). Data pooled from three independent experimental repeats ± SEM.

Close modal

Examination of individual allotype pairs shows that the range of responses to both precursors observed for each CD8/TIL group overlap for both CD8/TILhigh and CD8/TILmod, but not CD8/TILlow (D-LV9 = CD8/TILhigh 70%–90%, CD8/TILmod 60%–80% and CD8/TILlow 50%–60%; ED-LV9 = CD8/TILhigh 60%–80%, CD8/TILmod 55%–65%, and TILlow 40%–45%; Fig. 4D and E). In particular, trimming of ED-LV9 showed stratification of allotype pairs based of CD8/TIL numbers (Fig. 5E). The *007 allotype has been previously reported to be a hyperactive variant (16), therefore the reduced overall trimming activity in this CD8/TILlow allotype pair is likely due to the dominant negative effect of allotype *007 (Fig. 5D and E). These results indicate that there is a correlation between CD8/TIL numbers and the ability to generate the HPV-E7 LV9 peptide epitope. Thus, the decreased cell surface expression of this HPV epitope in HPV-infected cells and tumors from CD8/TILlow patients may compromise the ability of the immune system to recognize and kill HPV-infected cells.

In this study, we investigated the correlation between ERAP1 allotype sequence and the amount of TILs in 22 patients with HPV+ OPSCC. CD8+ T-cell tumor infiltration has been linked with improved disease prognosis (5, 9). Examination of ERAP1 allotypes revealed that although many allotypes were shared between the CD8/TIL groups, some differentiation between patients with high versus low CD8/TIL numbers was seen even in this small sample. Moreover, we observed that ERAP1 allotype pairs expressed in the high CD8/TIL group had a greater trimming capacity compared with those identified from patients with low CD8/TIL numbers. This was further evident from the ability of the ERAP1 allotype pairs from different CD8/TIL groups to generate an E7 derived epitope, with low CD8/TIL ERAP1 pairs being poorer at generating the epitope and stimulating an LV9-specific T-cell hybridoma. The anti-HPV T-cell response in patients with HPV+ OPSCC, a cohort linked to a better response to therapy and prognosis (9, 26–29), is primarily directed towards antigens derived from the E6 and E7 oncoproteins (5, 7, 19). Indeed, we observed responses to 3 epitopes derived from these proteins in tumors. Thus, the reduced capacity to generate the E7 derived LV9 epitope by CD8/TIL low ERAP1 pairs suggests that ERAP1 activity may control presentation of epitopes derived from E6 and E7 that provoke antitumor T-cell responses. In addition, within each CD8/TIL group, the ERAP1 pairs identified in patients had different peptide trimming specificities. This may explain why the correlation between ERAP1 identity and CD8/TIL level is not absolute.

ERAP1 polymorphism is associated with CSCC in GWAS (12). Despite the lack of GWAS of HPV+ OPSCC, our observations of ERAP1 allotypes, both their sequence and functionality, are consistent with a common immunologic mechanism between HPV+ OPSCC and CSCC. SNPs in ERAP1 together with TAP2 and LMP7 (also part of the antigen processing pathway) are associated with a 3-fold increase in cervical carcinoma risk (12). In addition, downregulation of ERAP1 expression is an independent predictor of decreased overall survival and disease-free survival in cervical carcinoma (30). In contrast, a study revealed ERAP1 expression was upregulated in head and neck squamous cell carcinoma (HNSCC) and HPV16+ CSCC cell lines and cervical cancer (31). Although the identity of the overexpressed ERAP1 allotypes was not characterized in all cell lines in that study, to assess whether they expressed SNPs associated with CSCC risk, knockdown of ERAP1 in 1 of 2 cervical cancer cell lines in which ERAP1 alleles were identified resulted in increased killing by E781-91–specific T cells. This cell line expressed ERAP1 that is known to efficiently trim peptides, but for this epitope overtrims instead, indicating that ERAP1 can affect the presentation of HPV-derived epitopes. Therefore, the differences in function we have observed between ERAP1 allotypes could be relevant for controlling HPV-associated cancer. We did not examine expression of the different ERAP1 allotypes in tumors in this study, but altering expression of specific ERAP1 allotypes that are efficient trimmers, such as *002, when coexpressed with a poor trimmer, *001, within the tumor may alter the ability to generate of tumor epitopes and diminish CD8+ T-cell activation. How the identity and expression of different ERAP1 allotype pairs affect generation of HPV-derived E6/E7 epitopes and their correlation with the antigen-specificity of HPV-specific CD8+ T cells in patient cohorts, and whether these correlate with positive outcome and response to checkpoint blockade therapy, remains unclear. ERAP1 allotypes expressed in germline DNA may provide a biomarker for prognosis that would guide the patient's care pathway at the time of diagnosis.

The mechanism by which ERAP1 trims peptides and the factors that determine the efficacy of trimming such as the backbone sequence and the amino acid preceding the peptide (N-1) is poorly understood. Assessment of the trimming efficacy of ERAP1 containing single SNP variants showed that the catalytic rate of the variant enzymes was affected to a similar degree for a number of peptide substrates with different backbones (32). The identity of the N-1 amino acid, however, influenced the ability of ERAP1 to generate the final peptide epitope (11, 33) suggesting that the N-1 amino acid affects trimming to a greater extent than the backbone. Consistent with this, we found that for 2 different peptide backbones, SHL8 and LV9 (which bind to different MHC molecules), the general trimming activity of the different ERAP1 pairs was similar. However, we found that the precise amino acid trimming specificity profiles of ERAP1 allotypes/pairs were different between X-SHL8 and X-LV9. For example, the “top” 3 specificities for the single allotype *015 (CD8/TILmod) were E, H, and Q for X-SHL8 and E, H, and K for X-LV9. For single allotype *001 (CD8/TILlow), we observed a greater difference where the “top” 3 specificities were C, A, and E for X-SHL8, but E, K, and N for X-LV9. These differences indicate that compared with overall trimming phenotype, amino acid specificity profiles may depend on the peptide backbone. A greater understanding of how these specificity profiles relate to the core peptide sequence could enable the development of ERAP1 processing algorithms for individual and combined allotypes. This could be integrated with algorithms that predict proteasome cleavage (20–22), peptide binding affinity (34), and the effects of peptide abundance and intracellular competition (35) to improve the ability to predict whether candidate peptides are likely to be presented at the cell surface in sufficient abundance to prime an effective CTL response.

ERAP2, a homolog of ERAP1 that is less polymorphic, may also play a role in the generation of HPV-derived epitopes. ERAP2 has altered substrate handling and amino acid trimming preferences compared with ERAP1, showing a preference for Arg (R) and less enzymatic activity (36). Despite its reduced activity, ERAP2 likely also contributes to the generation of HPV-derived epitopes.

Currently, three commercially available prophylactic HPV vaccines target high-risk HPV types (37). These vary in efficacy for prevention of both CIN2+ and CIN3+ (37). Neither of these vaccines are effective at eliminating preexisting HPV infections (38, 39) and therefore infected individuals remain at risk. The clearance of established HPV infections is associated with CD4+ T-helper 1 and CD8+ T cells (40), with the E6 and E7 oncoproteins targets for the response. No therapeutic vaccines are yet available despite several promising clinical trials (41–43). One contributing factor to this failure may be related to ERAP1, as the expression of particular ERAP1 allotype pairs may prevent or reduce the capacity to generate immunogenic E6/E7 peptide epitopes. Thus, affecting the ability of antigen-presenting cells to prime HPV-specific CD8+ T-cell responses in tumor draining lymph node and, for T cells that traffic to the tumor, the ability of tumor cells to activate them. Future vaccine design may need to incorporate the effect different ERAP1 allotypes have on the generation of candidate HPV-derived epitopes. In addition, the use of ERAP1 inhibitors (44, 45) alone or in combination with vaccines may prevent destruction of immunogenic epitopes through overtrimming by ERAP1 allotypes (46, 47), enabling activation of anti-HPV CD8+ T-cell responses and tumor protection.

In conclusion, this study demonstrates that the function of expressed ERAP1 allotypes correlates with the amount of infiltrating lymphocytes in tumors of patients with OPSCC and thus, prognosis. CD8/TILlow ERAP1 allotype pairs were poor at generating an HPV E7 epitope, indicating that the lack of T-cell infiltration and subsequent poorer antitumor response may be related to the ability of HPV-infected/tumor cells to generate HPV-specific peptide epitopes.

T. Elliott is a consultant/advisory board member of Avvinity Therapeutics. No potential conflicts of interest were disclosed by the other authors.

Conception and design: E. Reeves, E.V. King, T. Elliott, E. James

Development of methodology: E. Reeves, O. Wood, G.J. Thomas, T. Elliott, E. James

Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): E. Reeves, O. Wood, C.H. Ottensmeier, E.V. King, G.J. Thomas

Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): E. Reeves, T. Elliott, E. James

Writing, review, and/or revision of the manuscript: E. Reeves, C.H. Ottensmeier, E.V. King, G.J. Thomas, T. Elliott, E. James

Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): E. James

Study supervision: T. Elliott, E. James

The authors would like to thank Pandurangan Vijayanand for his support in collating samples and Nasia Kontouli for technical support. This work was supported by Cancer Research UK Programme Grant A16997 awarded to E. James and T. Elliott.

The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

1.
Galon
J
,
Costes
A
,
Sanchez-Cabo
F
,
Kirilovsky
A
,
Mlecnik
B
,
Lagorce-Pages
C
, et al
Type, density, and location of immune cells within human colorectal tumors predict clinical outcome
.
Science
2006
;
313
:
1960
4
.
2.
Leffers
N
,
Gooden
MJ
,
de Jong
RA
,
Hoogeboom
BN
,
ten Hoor
KA
,
Hollema
H
, et al
Prognostic significance of tumor-infiltrating T-lymphocytes in primary and metastatic lesions of advanced stage ovarian cancer
.
Cancer Immunol Immunother
2009
;
58
:
449
59
.
3.
Noble
F
,
Mellows
T
,
McCormick Matthews
LH
,
Bateman
AC
,
Harris
S
,
Underwood
TJ
, et al
Tumor infiltrating lymphocytes correlate with improved survival in patients with oesophageal adenocarcinoma
.
Cancer Immunol Immunother
2016
;
65
:
651
62
.
4.
Piersma
SJ
,
Jordanova
ES
,
van Poelgeest
MI
,
Kwappenberg
KM
,
van der Hulst
JM
,
Drijfhout
JW
, et al
High number of intraepithelial CD8+ tumor-infiltrating lymphocytes is associated with the absence of lymph node metastases in patients with large early-stage cervical cancer
.
Cancer Res
2007
;
67
:
354
61
.
5.
Welters
MJP
,
Ma
W
,
Santegoets
SJAM
,
Goedemans
R
,
Ehsan
I
,
Jordanova
ES
, et al
Intratumoral HPV16-specific T cells constitute a type I-oriented tumor microenvironment to improve survival in HPV16-driven oropharyngeal cancer
.
Clin Cancer Res
2018
;
24
:
634
47
.
6.
Brambilla
E
,
Le Teuff
G.
,
Marguet
S
,
Lantuejoul
S
,
Dunant
A
,
Graziano
S
, et al
Prognostic effect of tumor lymphocytic infiltration in resectable non-small-cell lung cancer
.
J Clin Oncol
2016
;
34
:
1223
30
.
7.
Albers
A
,
Abe
K
,
Hunt
J
,
Wang
J
,
Lopez-Albaitero
A
,
Schaefer
C
, et al
Antitumor activity of human papillomavirus type 16 E7-specific T cells against virally infected squamous cell carcinoma of the head and neck
.
Cancer Res
2005
;
65
:
11146
55
.
8.
de Vos van Steenwijk
PJ
,
Heusinkveld
M
,
Ramwadhdoebe
TH
,
Lowik
MJ
,
van der Hulst
JM
,
Goedemans
R
, et al
An unexpectedly large polyclonal repertoire of HPV-specific T cells is poised for action in patients with cervical cancer
.
Cancer Res
2010
;
70
:
2707
17
.
9.
Ward
MJ
,
Thirdborough
SM
,
Mellows
T
,
Riley
C
,
Harris
S
,
Suchak
K
, et al
Tumor-infiltrating lymphocytes predict for outcome in HPV-positive oropharyngeal cancer
.
Br J Cancer
2014
;
110
:
489
500
.
10.
Blum
JS
,
Wearsch
PA
,
Cresswell
P
. 
Pathways of antigen processing
.
Annu Rev Immunol
2013
;
31
:
443
73
.
11.
Reeves
E
,
Edwards
CJ
,
Elliott
T
,
James
E
. 
Naturally occurring ERAP1 haplotypes encode functionally distinct alleles with fine substrate specificity
.
J Immunol
2013
;
191
:
35
43
.
12.
Mehta
AM
,
Jordanova
ES
,
van Wezel
T
,
Uh
HW
,
Corver
WE
,
Kwappenberg
KM
, et al
Genetic variation of antigen processing machinery components and association with cervical carcinoma
.
Genes Chromosomes Cancer
2007
;
46
:
577
86
.
13.
Mehta
AM
,
Jordanova
ES
,
Corver
WE
,
van Wezel
T
,
Uh
HW
,
Kenter
GG
, et al
Single nucleotide polymorphisms in antigen processing machinery component ERAP1 significantly associate with clinical outcome in cervical carcinoma
.
Genes Chromosomes Cancer
2009
;
48
:
410
8
.
14.
Komdeur
FL
,
Prins
TM
,
van de Wall
S
,
Plat
A
,
Wisman
GBA
,
Hollema
H
, et al
CD103+ tumor-infiltrating lymphocytes are tumor-reactive intraepithelial CD8+ T cells associated with prognostic benefit and therapy response in cervical cancer
.
Oncoimmunology
2017
;
6
:
e1338230
.
15.
Sanderson
S
,
Shastri
N
. 
LacZ inducible, antigen/MHC-specific T cell hybrids
.
Int Immunol
1994
;
6
:
369
76
.
16.
Reeves
E
,
Colebatch-Bourn
A
,
Elliott
T
,
Edwards
CJ
,
James
E
. 
Functionally distinct ERAP1 allotype combinations distinguish individuals with Ankylosing Spondylitis
.
Proc Natl Acad Sci U S A
2014
;
111
:
17594
9
.
17.
Wansom
D
,
Light
E
,
Worden
F
,
Prince
M
,
Urba
S
,
Chepeha
DB
, et al
Correlation of cellular immunity with human papillomavirus 16 status and outcome in patients with advanced oropharyngeal cancer
.
Arch Otolaryngol Head Neck Surg
2010
;
136
:
1267
73
.
18.
Ombrello
MJ
,
Kastner
DL
,
Remmers
EF
. 
Endoplasmic reticulum-associated amino-peptidase 1 and rheumatic disease: genetics
.
Curr Opin Rheumatol
2015
;
27
:
349
56
.
19.
Heusinkveld
M
,
Goedemans
R
,
Briet
RJ
,
Gelderblom
H
,
Nortier
JW
,
Gorter
A
, et al
Systemic and local human papillomavirus 16-specific T-cell immunity in patients with head and neck cancer
.
Int J Cancer
2012
;
131
:
E74
85
.
20.
Hakenberg
J
,
Nussbaum
AK
,
Schild
H
,
Rammensee
HG
,
Kuttler
C
,
Holzhutter
HG
, et al
MAPPP: MHC class I antigenic peptide processing prediction
.
Appl Bioinformatics
2003
;
2
:
155
8
.
21.
Kuttler
C
,
Nussbaum
AK
,
Dick
TP
,
Rammensee
HG
,
Schild
H
,
Hadeler
KP
. 
An algorithm for the prediction of proteasomal cleavages
.
J Mol Biol
2000
;
298
:
417
29
.
22.
Nielsen
M
,
Lundegaard
C
,
Lund
O
,
Kesmir
C
. 
The role of the proteasome in generating cytotoxic T-cell epitopes: insights obtained from improved predictions of proteasomal cleavage
.
Immunogenetics
2005
;
57
:
33
41
.
23.
Ressing
ME
,
Sette
A
,
Brandt
RM
,
Ruppert
J
,
Wentworth
PA
,
Hartman
M
, et al
Human CTL epitopes encoded by human papillomavirus type 16 E6 and E7 identified through in vivo and in vitro immunogenicity studies of HLA-A*0201-binding peptides
.
J Immunol
1995
;
154
:
5934
43
.
24.
Xu
Y
,
Zhu
KJ
,
Chen
XZ
,
Zhao
KJ
,
Lu
ZM
,
Cheng
H
. 
Mapping of cytotoxic T lymphocytes epitopes in E7 antigen of human papillomavirus type 11
.
Arch Dermatol Res
2008
;
300
:
235
42
.
25.
Krishna
S
,
Ulrich
P
,
Wilson
E
,
Parikh
F
,
Narang
P
,
Yang
S
, et al
Human papilloma virus specific immunogenicity and dysfunction of CD8(+) T cells in head and neck cancer
.
Cancer Res
2018
;
78
:
6159
70
.
26.
Badoual
C
,
Hans
S
,
Merillon
N
,
Van Ryswick
C
,
Ravel
P
,
Benhamouda
N
, et al
PD-1-expressing tumor-infiltrating T cells are a favorable prognostic biomarker in HPV-associated head and neck cancer
.
Cancer Res
2013
;
73
:
128
38
.
27.
Oguejiofor
K
,
Galletta-Williams
H
,
Dovedi
SJ
,
Roberts
DL
,
Stern
PL
,
West
CM.
Distinct patterns of infiltrating CD8+ T cells in HPV+ and CD68 macrophages in HPV- oropharyngeal squamous cell carcinomas are associated with better clinical outcome but PD-L1 expression is not prognostic
.
Oncotarget
2017
;
8
:
14416
27
.
28.
Ang
KK
,
Harris
J
,
Wheeler
R
,
Weber
R
,
Rosenthal
DI
,
Nguyen-Tan
PF
, et al
Human papillomavirus and survival of patients with oropharyngeal cancer
.
N Engl J Med
2010
;
363
:
24
35
.
29.
Fakhry
C
,
Westra
WH
,
Li
S
,
Cmelak
A
,
Ridge
JA
,
Pinto
H
, et al
Improved survival of patients with human papillomavirus-positive head and neck squamous cell carcinoma in a prospective clinical trial
.
J Natl Cancer Inst
2008
;
100
:
261
9
.
30.
Mehta
AM
,
Jordanova
ES
,
Kenter
GG
,
Ferrone
S
,
Fleuren
GJ
. 
Association of antigen processing machinery and HLA class I defects with clinicopathological outcome in cervical carcinoma
.
Cancer Immunol Immunother
2008
;
57
:
197
206
.
31.
Steinbach
A
,
Winter
J
,
Reuschenbach
M
,
Blatnik
R
,
Klevenz
A
,
Bertrand
M
, et al
ERAP1 overexpression in HPV-induced malignancies: a possible novel immune evasion mechanism
.
Oncoimmunology
2017
;
6
:
e1336594
.
32.
Evnouchidou
I
,
Kamal
RP
,
Seregin
SS
,
Goto
Y
,
Tsujimoto
M
,
Hattori
A
, et al
Cutting edge: coding single nucleotide polymorphisms of endoplasmic reticulum aminopeptidase 1 can affect antigenic peptide generation in vitro by influencing basic enzymatic properties of the enzyme
.
J Immunol
2011
;
186
:
1909
13
.
33.
Hearn
A
,
York
IA
,
Rock
KL
. 
The specificity of trimming of MHC class I-presented peptides in the endoplasmic reticulum
.
J Immunol
2009
;
183
:
5526
36
.
34.
Andreatta
M
,
Nielsen
M
. 
Gapped sequence alignment using artificial neural networks: application to the MHC class I system
.
Bioinformatics
2016
;
32
:
511
7
.
35.
Boulanger
DSM
,
Eccleston
RC
,
Phillips
A
,
Coveney
PV
,
Elliott
T
,
Dalchau
N
. 
A mechanistic model for predicting cell surface presentation of competing peptides by MHC class I molecules
.
Front Immunol
2018
;
9
:1538.
36.
Evnouchidou
I
,
Birtley
J
,
Seregin
S
,
Papakyriakou
A
,
Zervoudi
E
,
Samiotaki
M
, et al
A common single nucleotide polymorphism in endoplasmic reticulum aminopeptidase 2 induces a specificity switch that leads to altered antigen processing
.
J Immunol
2012
;
189
:
2383
92
.
37.
Harper
DM
,
DeMars
LR
. 
HPV vaccines - A review of the first decade
.
Gynecol Oncol
2017
;
146
:
196
204
.
38.
Olsson
SE
,
Kjaer
SK
,
Sigurdsson
K
,
Iversen
OE
,
Hernandez-Avila
M
,
Wheeler
CM
, et al
Evaluation of quadrivalent HPV 6/11/16/18 vaccine efficacy against cervical and anogenital disease in subjects with serological evidence of prior vaccine type HPV infection
.
Hum Vaccin
2009
;
5
:
696
704
.
39.
Paavonen
J
,
Naud
P
,
Salmeron
J
,
Wheeler
CM
,
Chow
SN
,
Apter
D
, et al
Efficacy of human papillomavirus (HPV)-16/18 AS04-adjuvanted vaccine against cervical infection and precancer caused by oncogenic HPV types (PATRICIA): final analysis of a double-blind, randomised study in young women
.
Lancet
2009
;
374
:
301
14
.
40.
van der Burg
SH
,
Melief
CJ
. 
Therapeutic vaccination against human papilloma virus induced malignancies
.
Curr Opin Immunol
2011
;
23
:
252
7
.
41.
Kim
HJ
,
Kim
HJ.
Current status and future prospects for human papillomavirus vaccines
.
Arch Pharm Res
2017
;
40
:
1050
63
.
42.
Vici
P
,
Pizzuti
L
,
Mariani
L
,
Zampa
G
,
Santini
D
,
Di Lauro
L
, et al
Targeting immune response with therapeutic vaccines in premalignant lesions and cervical cancer: hope or reality from clinical studies
.
Expert Rev Vaccines
2016
;
15
:
1327
36
.
43.
Yang
A
,
Farmer
E
,
Wu
TC
,
Hung
CF
. 
Perspectives for therapeutic HPV vaccine development
.
J Biomed Sci
2016
;
23
:
75
.
44.
Weglarz-Tomczak
E
,
Vassiliou
S
,
Mucha
A
. 
Discovery of potent and selective inhibitors of human aminopeptidases ERAP1 and ERAP2 by screening libraries of phosphorus-containing amino acid and dipeptide analogues
.
Bioorg Med Chem Lett
2016
;
26
:
4122
6
.
45.
Zervoudi
E
,
Saridakis
E
,
Birtley
JR
,
Seregin
SS
,
Reeves
E
,
Kokkala
P
, et al
Rationally designed inhibitor targeting antigen-trimming aminopeptidases enhances antigen presentation and cytotoxic T-cell responses
.
Proc Natl Acad Sci U S A
2013
;
110
:
19890
5
.
46.
James
E
,
Bailey
I
,
Sugiyarto
G
,
Elliott
T
. 
Induction of protective antitumor immunity through attenuation of ERAAP function
.
J Immunol
2013
;
190
:
5839
46
.
47.
Keller
M
,
Ebstein
F
,
Burger
E
,
Textoris-Taube
K
,
Gorny
X
,
Urban
S
, et al
The proteasome immunosubunits, PA28 and ER-aminopeptidase 1 protect melanoma cells from efficient MART-126-35 -specific T-cell recognition
.
Eur J Immunol
2015
;
45
:
3257
68
.