Abstract
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.
Introduction
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.
Materials and Methods
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).
Results
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).
ERAP1 allotype identity in patients with OPSCC
. | . | . | Amino acid at indicated position . | |||||||
---|---|---|---|---|---|---|---|---|---|---|
. | . | . | 56 . | 127 . | 346 . | 349 . | 528 . | 575 . | 725 . | 730 . |
Patient number . | CD8/TIL status . | ERAP1 allotype . | E . | R . | G . | M . | K . | D . | R . | Q . |
1 | High | *001 | E | P | G | V | R | N | Q | E |
*002 | E | R | G | M | K | D | R | Q | ||
2 | High | *011 | E | R | G | M | R | D | R | E |
*016 | E | R | G | V | K | D | Q | E | ||
3 | High | *001 | E | P | G | V | R | N | Q | E |
*018 | E | R | G | M | K | D | R | E | ||
4 | High | *001 | E | P | G | V | R | N | Q | E |
*002 | E | R | G | M | K | D | R | Q | ||
5 | High | *001 | E | P | G | V | R | N | Q | E |
*002 | E | R | G | M | K | D | R | Q | ||
6 | High | *002 | E | R | G | M | K | D | R | Q |
*002 | E | R | G | M | K | D | R | Q | ||
7 | High | *002 | E | R | G | M | K | D | R | Q |
*014 | K | P | G | M | R | D | R | E | ||
8 | High | *001 | E | P | G | V | R | N | Q | E |
*011 | E | R | G | M | R | D | R | E | ||
9 | High | *011 | E | R | G | M | R | D | R | E |
*015 | E | P | G | M | R | D | R | E | ||
10 | High | *001 | E | P | G | V | R | N | Q | E |
*002 | E | R | G | M | K | D | R | Q | ||
11 | Moderate | *015 | E | P | G | M | R | D | R | E |
*015 | E | P | G | M | R | D | R | E | ||
12 | Moderate | *001 | E | P | G | V | R | N | Q | E |
*011 | E | R | G | M | R | D | R | E | ||
13 | Moderate | *002 | E | R | G | M | K | D | R | Q |
*017 | E | R | G | M | K | N | R | Q | ||
14 | Moderate | *011 | E | R | G | M | R | D | R | E |
*018 | E | R | G | M | K | D | R | E | ||
15 | Moderate | *002 | E | R | G | M | K | D | R | Q |
*011 | E | R | G | M | R | D | R | E | ||
16 | Moderate | *011 | E | R | G | M | R | D | R | E |
*011 | E | R | G | M | R | D | R | E | ||
17 | Moderate | *016 | E | R | G | V | K | D | Q | E |
*018 | E | R | G | M | K | D | R | E | ||
18 | Low | *007 | E | R | G | M | R | D | Q | Q |
*011 | E | R | G | M | R | D | R | E | ||
19 | Low | *001 | E | P | G | V | R | N | Q | E |
*019 | E | R | D | M | R | D | R | E | ||
20 | Low | *001 | E | P | G | V | R | N | Q | E |
*001 | E | P | G | V | R | N | Q | E | ||
21 | Low | *011 | E | R | G | M | R | D | R | E |
*011 | E | R | G | M | R | D | R | E | ||
22 | Low | *013 | E | P | G | M | K | D | R | Q |
*019 | E | R | D | M | R | D | R | E | ||
23 | Low | *002 | E | R | G | M | K | D | R | Q |
*002 | E | R | G | M | K | D | R | Q | ||
24 | Low | *002 | E | R | G | M | K | D | R | Q |
*011 | E | R | G | M | R | D | R | E | ||
25 | Low | *011 | E | R | G | M | R | D | R | E |
*011 | E | P | G | M | R | D | R | E |
. | . | . | Amino acid at indicated position . | |||||||
---|---|---|---|---|---|---|---|---|---|---|
. | . | . | 56 . | 127 . | 346 . | 349 . | 528 . | 575 . | 725 . | 730 . |
Patient number . | CD8/TIL status . | ERAP1 allotype . | E . | R . | G . | M . | K . | D . | R . | Q . |
1 | High | *001 | E | P | G | V | R | N | Q | E |
*002 | E | R | G | M | K | D | R | Q | ||
2 | High | *011 | E | R | G | M | R | D | R | E |
*016 | E | R | G | V | K | D | Q | E | ||
3 | High | *001 | E | P | G | V | R | N | Q | E |
*018 | E | R | G | M | K | D | R | E | ||
4 | High | *001 | E | P | G | V | R | N | Q | E |
*002 | E | R | G | M | K | D | R | Q | ||
5 | High | *001 | E | P | G | V | R | N | Q | E |
*002 | E | R | G | M | K | D | R | Q | ||
6 | High | *002 | E | R | G | M | K | D | R | Q |
*002 | E | R | G | M | K | D | R | Q | ||
7 | High | *002 | E | R | G | M | K | D | R | Q |
*014 | K | P | G | M | R | D | R | E | ||
8 | High | *001 | E | P | G | V | R | N | Q | E |
*011 | E | R | G | M | R | D | R | E | ||
9 | High | *011 | E | R | G | M | R | D | R | E |
*015 | E | P | G | M | R | D | R | E | ||
10 | High | *001 | E | P | G | V | R | N | Q | E |
*002 | E | R | G | M | K | D | R | Q | ||
11 | Moderate | *015 | E | P | G | M | R | D | R | E |
*015 | E | P | G | M | R | D | R | E | ||
12 | Moderate | *001 | E | P | G | V | R | N | Q | E |
*011 | E | R | G | M | R | D | R | E | ||
13 | Moderate | *002 | E | R | G | M | K | D | R | Q |
*017 | E | R | G | M | K | N | R | Q | ||
14 | Moderate | *011 | E | R | G | M | R | D | R | E |
*018 | E | R | G | M | K | D | R | E | ||
15 | Moderate | *002 | E | R | G | M | K | D | R | Q |
*011 | E | R | G | M | R | D | R | E | ||
16 | Moderate | *011 | E | R | G | M | R | D | R | E |
*011 | E | R | G | M | R | D | R | E | ||
17 | Moderate | *016 | E | R | G | V | K | D | Q | E |
*018 | E | R | G | M | K | D | R | E | ||
18 | Low | *007 | E | R | G | M | R | D | Q | Q |
*011 | E | R | G | M | R | D | R | E | ||
19 | Low | *001 | E | P | G | V | R | N | Q | E |
*019 | E | R | D | M | R | D | R | E | ||
20 | Low | *001 | E | P | G | V | R | N | Q | E |
*001 | E | P | G | V | R | N | Q | E | ||
21 | Low | *011 | E | R | G | M | R | D | R | E |
*011 | E | R | G | M | R | D | R | E | ||
22 | Low | *013 | E | P | G | M | K | D | R | Q |
*019 | E | R | D | M | R | D | R | E | ||
23 | Low | *002 | E | R | G | M | K | D | R | Q |
*002 | E | R | G | M | K | D | R | Q | ||
24 | Low | *002 | E | R | G | M | K | D | R | Q |
*011 | E | R | G | M | R | D | R | E | ||
25 | Low | *011 | E | R | G | M | R | D | R | E |
*011 | E | P | G | M | R | D | 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).
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.
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.
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).
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).
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).
HLA-A*0201–restricted HPV-16 E6/E7 peptide epitopes
Protein . | Position (amino acid) . | Epitope sequence . | MHC IC50 (predicted nmol/L) NetMHC . | Predicted 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 |
Protein . | Position (amino acid) . | Epitope sequence . | MHC IC50 (predicted nmol/L) NetMHC . | Predicted 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.
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.
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.
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.
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).
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).
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).
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.
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.
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.
Discussion
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.
Disclosure of Potential Conflicts of Interest
T. Elliott is a consultant/advisory board member of Avvinity Therapeutics. No potential conflicts of interest were disclosed by the other authors.
Authors' Contributions
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
Acknowledgments
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.