Deep neural networks predict class I major histocompatibility complex epitope presentation and transfer learn neoepitope immunogenicity.

Identifying neoepitopes that elicit an adaptive immune response is a major bottleneck to developing personalized cancer vaccines. Experimental validation of candidate neoepitopes is extremely resource intensive and the vast majority of candidates are non-immunogenic, creating a needle-in-a-haystack problem. Here we address this challenge, presenting computational methods for predicting class I major histocompatibility complex (MHC-I) epitopes and identifying immunogenic neoepitopes with improved precision. The BigMHC method comprises an ensemble of seven pan-allelic deep neural networks trained on peptide-MHC eluted ligand data from mass spectrometry assays and transfer learned on data from assays of antigen-specific immune response. Compared with four state-of-the-art classifiers, BigMHC significantly improves the prediction of epitope presentation on a test set of 45,409 MHC ligands among 900,592 random negatives (area under the receiver operating characteristic = 0.9733; area under the precision-recall curve = 0.8779). After transfer learning on immunogenicity data, BigMHC yields significantly higher precision than seven state-of-the-art models in identifying immunogenic neoepitopes, making BigMHC effective in clinical settings.

Chimeric RNAs reveal putative neoantigen peptides for developing tumor vaccines for breast cancer.

Introduction: We present here a strategy to identify immunogenic neoantigen candidates from unique amino acid sequences at the junctions of fusion proteins which can serve as targets in the development of tumor vaccines for the treatment of breastcancer. Method: We mined the sequence reads of breast tumor tissue that are usually discarded as discordant paired-end reads and discovered cancer specific fusion transcripts using tissue from cancer free controls as reference. Binding affinity predictions of novel peptide sequences crossing the fusion junction were analyzed by the MHC Class I binding predictor, MHCnuggets. CD8+ T cell responses against the 15 peptides were assessed through in vitro Enzyme Linked Immunospot (ELISpot). Results: We uncovered 20 novel fusion transcripts from 75 breast tumors of 3 subtypes: TNBC, HER2+, and HR+. Of these, the NSFP1-LRRC37A2 fusion transcript was selected for further study. The 3833 bp chimeric RNA predicted by the consensus fusion junction sequence is consistent with a read-through transcription of the 5'-gene NSFP1-Pseudo gene NSFP1 (NSFtruncation at exon 12/13) followed by trans-splicing to connect withLRRC37A2 located immediately 3' through exon 1/2. A total of 15 different 8-mer neoantigen peptides discovered from the NSFP1 and LRRC37A2 truncations were predicted to bind to a total of 35 unique MHC class I alleles with a binding affinity of IC50<500nM.); 1 of which elicited a robust immune response. Conclusion: Our data provides a framework to identify immunogenic neoantigen candidates from fusion transcripts and suggests a potential vaccine strategy to target the immunogenic neopeptides in patients with tumors carrying the NSFP1-LRRC37A2 fusion.

Persistent mutation burden drives sustained anti-tumor immune responses.

Tumor mutation burden is an imperfect proxy of tumor foreignness and has therefore failed to consistently demonstrate clinical utility in predicting responses in the context of immunotherapy. We evaluated mutations in regions of the genome that are unlikely to undergo loss in a pan-cancer analysis across 31 tumor types (n = 9,242) and eight immunotherapy-treated cohorts of patients with non-small-cell lung cancer, melanoma, mesothelioma, and head and neck cancer (n = 524). We discovered that mutations in single-copy regions and those present in multiple copies per cell constitute a persistent tumor mutation burden (pTMB) which is linked with therapeutic response to immune checkpoint blockade. Persistent mutations were retained in the context of tumor evolution under selective pressure of immunotherapy and tumors with a high pTMB content were characterized by a more inflamed tumor microenvironment. pTMB imposes an evolutionary bottleneck that cancer cells cannot overcome and may thus drive sustained immunologic tumor control in the context of immunotherapy.

Sex-specific differences in immunogenomic features of response to immune checkpoint blockade.

Introduction: The magnitude of response to immune checkpoint inhibitor (ICI) therapy may be sex-dependent, as females have lower response rates and decreased survival after ICI monotherapy. The mechanisms underlying this sex dimorphism in ICI response are unknown, and may be related to sex-driven differences in the immunogenomic landscape of tumors that shape anti-tumor immune responses in the context of therapy. Methods: To investigate the association of immunogenic mutations with HLA haplotypes, we leveraged whole exome sequence data and HLA genotypes from 482 non-small cell lung cancer (NSCLC) tumors from The Cancer Genome Atlas (TCGA). To explore sex-specific genomic features linked with ICI response, we analyzed whole exome sequence data from patients with NSCLC treated with ICI. Tumor mutational burden (TMB), HLA class I and II restricted immunogenic missense mutation (IMM) load, and mutational smoking signature were defined for each tumor. IMM load was combined with HLA class I and II haplotypes and correlated with therapeutic response and survival following ICI treatment. We examined rates of durable clinical benefit (DCB) for at least six months from ICI treatment initiation. Findings were validated utilizing whole exome sequence data from an independent cohort of ICI treated NSCLC. Results: Analysis of whole exome sequence data from NSCLC tumors of females and males revealed that germline HLA class II diversity (≥9 unique HLA alleles) was associated with higher tumor class II IMM load in females (p=0.01) and not in males (p=0.64). Similarly, in tumors of female patients, somatic HLA class II loss of heterozygosity was associated with increased IMM load (p=0.01) while this association was not observed in tumors in males (p=0.20). In females, TMB (p=0.005), class I IMM load (p=0.005), class II IMM load (p=0.004), and mutational smoking signature (p<0.001) were significantly higher in tumors responding to ICI as compared to non-responding tumors. In contrast, among males, there was no significant association between DCB and any of these features. When IMM was considered in the context of HLA zygosity, high MHC-II restricted IMM load and high HLA class II diversity was significantly associated with overall survival in males (p=0.017). Conclusions: Inherent sex-driven differences in immune surveillance affect the immunogenomic determinants of response to ICI and likely mediate the dimorphic outcomes with ICI therapy. Deeper understanding of the selective pressures and mechanisms of immune escape in tumors in males and females can inform patient selection strategies and can be utilized to further hone immunotherapy approaches in cancer.

Durvalumab with platinum-pemetrexed for unresectable pleural mesothelioma: survival, genomic and immunologic analyses from the phase 2 PrE0505 trial.

Mesothelioma is a rare and fatal cancer with limited therapeutic options until the recent approval of combination immune checkpoint blockade. Here we report the results of the phase 2 PrE0505 trial ( NCT02899195 ) of the anti-PD-L1 antibody durvalumab plus platinum-pemetrexed chemotherapy for 55 patients with previously untreated, unresectable pleural mesothelioma. The primary endpoint was overall survival compared to historical control with cisplatin and pemetrexed chemotherapy; secondary and exploratory endpoints included safety, progression-free survival and biomarkers of response. The combination of durvalumab with chemotherapy met the pre-specified primary endpoint, reaching a median survival of 20.4 months versus 12.1 months with historical control. Treatment-emergent adverse events were consistent with known side effects of chemotherapy, and all adverse events due to immunotherapy were grade 2 or lower. Integrated genomic and immune cell repertoire analyses revealed that a higher immunogenic mutation burden coupled with a more diverse T cell repertoire was linked to favorable clinical outcome. Structural genome-wide analyses showed a higher degree of genomic instability in responding tumors of epithelioid histology. Patients with germline alterations in cancer predisposing genes, especially those involved in DNA repair, were more likely to achieve long-term survival. Our findings indicate that concurrent durvalumab with platinum-based chemotherapy has promising clinical activity and that responses are driven by the complex genomic background of malignant pleural mesothelioma.

High-Throughput Prediction of MHC Class I and II Neoantigens with MHCnuggets.

Computational prediction of binding between neoantigen peptides and major histocompatibility complex (MHC) proteins can be used to predict patient response to cancer immunotherapy. Current neoantigen predictors focus on in silico estimation of MHC binding affinity and are limited by low predictive value for actual peptide presentation, inadequate support for rare MHC alleles, and poor scalability to high-throughput data sets. To address these limitations, we developed MHCnuggets, a deep neural network method that predicts peptide-MHC binding. MHCnuggets can predict binding for common or rare alleles of MHC class I or II with a single neural network architecture. Using a long short-term memory network (LSTM), MHCnuggets accepts peptides of variable length and is faster than other methods. When compared with methods that integrate binding affinity and MHC-bound peptide (HLAp) data from mass spectrometry, MHCnuggets yields a 4-fold increase in positive predictive value on independent HLAp data. We applied MHCnuggets to 26 cancer types in The Cancer Genome Atlas, processing 26.3 million allele-peptide comparisons in under 2.3 hours, yielding 101,326 unique predicted immunogenic missense mutations (IMM). Predicted IMM hotspots occurred in 38 genes, including 24 driver genes. Predicted IMM load was significantly associated with increased immune cell infiltration (P < 2 × 10-16), including CD8+ T cells. Only 0.16% of predicted IMMs were observed in more than 2 patients, with 61.7% of these derived from driver mutations. Thus, we describe a method for neoantigen prediction and its performance characteristics and demonstrate its utility in data sets representing multiple human cancers.

Integrative Tumor and Immune Cell Multi-omic Analyses Predict Response to Immune Checkpoint Blockade in Melanoma.

In this study, we incorporate analyses of genome-wide sequence and structural alterations with pre- and on-therapy transcriptomic and T cell repertoire features in immunotherapy-naive melanoma patients treated with immune checkpoint blockade. Although tumor mutation burden is associated with improved treatment response, the mutation frequency in expressed genes is superior in predicting outcome. Increased T cell density in baseline tumors and dynamic changes in regression or expansion of the T cell repertoire during therapy distinguish responders from non-responders. Transcriptome analyses reveal an increased abundance of B cell subsets in tumors from responders and patterns of molecular response related to expressed mutation elimination or retention that reflect clinical outcome. High-dimensional genomic, transcriptomic, and immune repertoire data were integrated into a multi-modal predictor of response. These findings identify genomic and transcriptomic characteristics of tumors and immune cells that predict response to immune checkpoint blockade and highlight the importance of pre-existing T and B cell immunity in therapeutic outcomes.