The anti-cancer immune response against mutated peptides of potential immunological relevance (neoantigens) is primarily attributed to MHC-I-restricted cytotoxic CD8+ T cell responses. ...MHC-II-restricted CD4+ T cells also drive anti-tumor responses, but their relation to neoantigen selection and tumor evolution has not been systematically studied. Modeling the potential of an individual’s MHC-II genotype to present 1,018 driver mutations in 5,942 tumors, we demonstrate that the MHC-II genotype constrains the mutational landscape during tumorigenesis in a manner complementary to MHC-I. Mutations poorly bound to MHC-II are positively selected during tumorigenesis, even more than mutations poorly bound to MHC-I. This emphasizes the importance of CD4+ T cells in anti-tumor immunity. In addition, we observed less inter-patient variation in mutation presentation for MHC-II than for MHC-I. These differences were reflected by age at diagnosis, which was correlated with presentation by MHC-I only. Collectively, our results emphasize the central role of MHC-II presentation in tumor evolution.
Display omitted
•Patient- and residue-specific PHBR-II score estimates mutation presentation by MHC-II•Tumors are less likely to harbor driver mutations that bind well to MHC-II•Frequent driver mutations are universally poorly presented by MHC-II•MHC-II shows less inter-patient variability but stronger selective effects than MHC-I
Inter-individual variation in MHC-II genotype influences the evolution of patient-specific tumor mutational spectra, emphasizing the key role of CD4+ T cells in anti-tumor immunity.
Individual MHC genotype constrains the mutational landscape during tumorigenesis. Immune checkpoint inhibition reactivates immunity against tumors that escaped immune surveillance in approximately ...30% of cases. Recent studies demonstrated poorer response rates in female and younger patients. Although immune responses differ with sex and age, the role of MHC-based immune selection in this context is unknown. We find that tumors in younger and female individuals accumulate more poorly presented driver mutations than those in older and male patients, despite no differences in MHC genotype. Younger patients show the strongest effects of MHC-based driver mutation selection, with younger females showing compounded effects and nearly twice as much MHC-II based selection. This study presents evidence that strength of immune selection during tumor development varies with sex and age, and may influence the availability of mutant peptides capable of driving effective response to immune checkpoint inhibitor therapy.
The major histocompatibility complex class I (MHC-I) molecule is a protein complex that displays intracellular peptides to T cells, allowing the immune system to recognize and destroy infected or ...cancerous cells. MHC-I is composed of a highly polymorphic HLA-encoded alpha chain that binds the peptide and a Beta-2-microglobulin (B2M) protein that acts as a stabilizing scaffold. HLA mutations have been implicated as a mechanism of immune evasion during tumorigenesis, and B2M is considered a tumor suppressor gene. However, the implications of somatic HLA and B2M mutations have not been fully explored in the context of antigen presentation via the MHC-I molecule during tumor development. To understand the effect that B2M and HLA MHC-I molecule mutations have on mutagenesis, we analyzed the accumulation of mutations in patients from The Cancer Genome Atlas according to their MHC-I molecule mutation status.
Somatic B2M and HLA mutations in microsatellite stable tumors were associated with higher overall mutation burden and a larger fraction of HLA-binding neoantigens when compared to B2M and HLA wild type tumors. B2M and HLA mutations were highly enriched in patients with microsatellite instability. B2M mutations tended to occur relatively early during patients' respective tumor development, whereas HLA mutations were either early or late events. In addition, B2M and HLA mutated patients had higher levels of immune infiltration by natural killer and CD8+ T cells and higher levels of cytotoxicity.
Our findings add to a growing body of evidence that somatic B2M and HLA mutations are a mechanism of immune evasion by demonstrating that such mutations are associated with a higher load of neoantigens that should be presented via MHC-I.
Immune checkpoint blockade (ICB) with antibodies inhibiting cytotoxic T lymphocyte-associated protein-4 (CTLA-4) and programmed cell death protein-1 (PD-1) (or its ligand (PD-L1)) can stimulate ...immune responses against cancer and have revolutionized the treatment of tumors. The influence of host germline genetics and its interaction with tumor neoantigens remains poorly defined. We sought to determine the interaction between tumor mutational burden (TMB) and the ability of a patient's major histocompatibility complex class I (MHC-I) to efficiently present mutated driver neoantigens in predicting response ICB.
Comprehensive genomic profiling was performed on 83 patients with diverse cancers treated with ICB to determine TMB and human leukocyte antigen-I (HLA-I) genotype. The ability of a patient's MHC-I to efficiently present mutated driver neoantigens (defined by the Patient Harmonic-mean Best Rank (PHBR) score (with lower PHBR indicating more efficient presentation)) was calculated for each patient.
The median progression-free survival (PFS) for PHBR score < 0.5 vs. ≥ 0.5 was 5.1 vs. 4.4 months (P = 0.04). Using a TMB cutoff of 10 mutations/mb, the stable disease > 6 months/partial response/complete response rate, median PFS, and median overall survival (OS) of TMB high/PHBR high vs. TMB high/PHBR low were 43% vs. 78% (P = 0.049), 5.8 vs. 26.8 months (P = 0.03), and 17.2 months vs. not reached (P = 0.23), respectively. These findings were confirmed in an independent validation cohort of 32 patients.
Poor presentation of driver mutation neoantigens by MHC-I may explain why some tumors (even with a high TMB) do not respond to ICB.
Human leukocyte antigen loss of heterozygosity (HLA LOH) allows cancer cells to escape immune recognition by deleting HLA alleles, causing the suppressed presentation of tumor neoantigens. Despite ...its importance in immunotherapy response, few methods exist to detect HLA LOH, and their accuracy is not well understood. Here, we develop DASH (Deletion of Allele-Specific HLAs), a machine learning-based algorithm to detect HLA LOH from paired tumor-normal sequencing data. With cell line mixtures, we demonstrate increased sensitivity compared to previously published tools. Moreover, our patient-specific digital PCR validation approach provides a sensitive, robust orthogonal approach that could be used for clinical validation. Using DASH on 610 patients across 15 tumor types, we find that 18% of patients have HLA LOH. Moreover, we show inflated HLA LOH rates compared to genome-wide LOH and correlations between CD274 (encodes PD-L1) expression and microsatellite instability status, suggesting the HLA LOH is a key immune resistance strategy.
This article has been withdrawn by the authors. A publication of the manuscript with the correct figures and tables has been approved and the authors state the conclusions of the manuscript remain ...unaffected. Specifically, errors are in Figure 6A, Supplementary Figure 10B, Supplementary Figure 10C, and Supplementary Table 5. The details of the errors are as follows: the HLA types for one sample were incorrectly assigned because of a tumor/normal mislabeling from the biobank vendor. Due to the differing HLA types between the tumor and normal sample, the sequence analysis established that the HLA alleles for this patient had been deleted (HLA LOH). The authors conclude that this was an artifact caused by the normal sample mislabeling. The corrected version can be accessed (Pyke, R.M., Mellacheruvu, D., Dea, S., Abbott, C.W., Zhang, S.V., Philips, N.A., Harris, J., Bartha, G., Desai, S., McClory, R., West, J., Snyder, M,P., Chen, R., Boyle, S.M. (2022) Precision Neoantigen Discovery Using Large-Scale Immunopeptidomics and Composite Modeling of MHC Peptide Presentation. Mol. Cell. Proteomics 22, 100506
While immune checkpoint blockade (ICB) has become a pillar of cancer treatment, biomarkers that consistently predict patient response remain elusive due to the complex mechanisms driving immune ...response to tumors. We hypothesized that a multi-dimensional approach modeling both tumor and immune-related molecular mechanisms would better predict ICB response than simpler mutation-focused biomarkers, such as tumor mutational burden (TMB).
Tumors from a cohort of patients with late-stage melanoma (
= 51) were profiled using an immune-enhanced exome and transcriptome platform. We demonstrate increasing predictive power with deeper modeling of neoantigens and immune-related resistance mechanisms to ICB.
Our neoantigen burden score, which integrates both exome and transcriptome features, more significantly stratified responders and nonresponders (
= 0.016) than TMB alone (
= 0.049). Extension of this model to include immune-related resistance mechanisms affecting the antigen presentation machinery, such as HLA allele-specific LOH, resulted in a composite neoantigen presentation score (NEOPS) that demonstrated further increased association with therapy response (
= 0.002).
NEOPS proved the statistically strongest biomarker compared with all single-gene biomarkers, expression signatures, and TMB biomarkers evaluated in this cohort. Subsequent confirmation of these findings in an independent cohort of patients (
= 110) suggests that NEOPS is a robust, novel biomarker of ICB response in melanoma.
Major histocompatibility complex (MHC)–bound peptides that originate from tumor-specific genetic alterations, known as neoantigens, are an important class of anticancer therapeutic targets. ...Accurately predicting peptide presentation by MHC complexes is a key aspect of discovering therapeutically relevant neoantigens. Technological improvements in mass spectrometry–based immunopeptidomics and advanced modeling techniques have vastly improved MHC presentation prediction over the past 2 decades. However, improvement in the accuracy of prediction algorithms is needed for clinical applications like the development of personalized cancer vaccines, the discovery of biomarkers for response to immunotherapies, and the quantification of autoimmune risk in gene therapies. Toward this end, we generated allele-specific immunopeptidomics data using 25 monoallelic cell lines and created Systematic Human Leukocyte Antigen (HLA) Epitope Ranking Pan Algorithm (SHERPA), a pan-allelic MHC-peptide algorithm for predicting MHC-peptide binding and presentation. In contrast to previously published large-scale monoallelic data, we used an HLA-null K562 parental cell line and a stable transfection of HLA allele to better emulate native presentation. Our dataset includes five previously unprofiled alleles that expand MHC diversity in the training data and extend allelic coverage in underprofiled populations. To improve generalizability, SHERPA systematically integrates 128 monoallelic and 384 multiallelic samples with publicly available immunoproteomics data and binding assay data. Using this dataset, we developed two features that empirically estimate the propensities of genes and specific regions within gene bodies to engender immunopeptides to represent antigen processing. Using a composite model constructed with gradient boosting decision trees, multiallelic deconvolution, and 2.15 million peptides encompassing 167 alleles, we achieved a 1.44-fold improvement of positive predictive value compared with existing tools when evaluated on independent monoallelic datasets and a 1.17-fold improvement when evaluating on tumor samples. With a high degree of accuracy, SHERPA has the potential to enable precision neoantigen discovery for future clinical applications.
Display omitted
•Generated 25 stably transfected monoallelic cell lines and applied immunopeptidomics.•Harmonized 512 public immunopeptidomic samples through systematic reprocessing.•Developed pan-allele MHC binding algorithm (SHERPA) utilizing 167 human HLA alleles.•Employed empirically derived antigen-processing features to predict MHC presentation.•SHERPA demonstrates up to 1.44-fold increased precision over competing algorithms.
Accurately identifying neoantigens is critical for many clinical applications. We generated immunopeptidomics data from 25 stably transfected monoallelic cell lines. Then, we systematically reprocessed a large corpus of public data to improve MHC binding pocket diversity and to empirically learn the rules of antigen presentation. In applying these datasets, we trained SHERPA, an MHC binding and presentation prediction algorithm. SHERPA improves performance compared with existing tools by 1.44-fold in held-out monoallelic data and 1.11-fold for known immunogenic epitopes.