Proteogenomic Analysis Unveils the HLA Class I-Presented Immunopeptidome in Melanoma and EGFR-Mutant Lung Adenocarcinoma

ImmuneApp for HLA-I epitope prediction and immunopeptidome analysis

Advances in mass spectrometry accelerates the characterization of HLA ligandome, necessitating the development of efficient methods for immunopeptidomics analysis and (neo)antigen prediction. We develop ImmuneApp, an interpretable deep learning framework trained on extensive HLA ligand datasets, which improves the prediction of HLA-I epitopes, prioritizes neoepitopes, and enhances immunopeptidomics deconvolution. ImmuneApp extracts informative embeddings and identifies key residues for pHLA binding. We also present a more accurate model-based deconvolution approach and systematically analyzed 216 multi-allelic immunopeptidomics samples, identifying 835,551 ligands restricted to over 100 HLA-I alleles. Our investigation reveals the effectiveness of the composite model, denoted as ImmuneApp-MA, which integrates mono- and multi-allelic data to enhance predictive performance. Leveraging ImmuneApp-MA as a pre-trained model, we built ImmuneApp-Neo, an immunogenicity predictor that outperforms existing methods for prioritizing immunogenic neoepitope. ImmuneApp demonstrates its utility across various immunopeptidomics datasets, which will promote the discovery of novel neoantigens and the development of new immunotherapies.

Exploring the dynamic landscape of immunopeptidomics: Unravelling posttranslational modifications and navigating bioinformatics terrain

Immunopeptidomics is becoming an increasingly important field of study. The capability to identify immunopeptides with pivotal roles in the human immune system is essential to shift the current curative medicine towards personalized medicine. Throughout the years, the field has matured, giving insight into the current pitfalls. Nowadays, it is commonly accepted that generalizing shotgun proteomics workflows is malpractice because immunopeptidomics faces numerous challenges. While many of these difficulties have been addressed, the road towards the ideal workflow remains complicated. Although the presence of Posttranslational modifications (PTMs) in the immunopeptidome has been demonstrated, their identification remains highly challenging despite their significance for immunotherapies. The large number of unpredictable modifications in the immunopeptidome plays a pivotal role in the functionality and these challenges. This review provides a comprehensive overview of the current advancements in immunopeptidomics. We delve into the challenges associated with identifying PTMs within the immunopeptidome, aiming to address the current state of the field.

pXg: Comprehensive Identification of Noncanonical MHC-I-Associated Peptides From De Novo Peptide Sequencing Using RNA-Seq Reads

Discovering noncanonical peptides has been a common application of proteogenomics. Recent studies suggest that certain noncanonical peptides, known as noncanonical major histocompatibility complex-I (MHC-I)-associated peptides (ncMAPs), that bind to MHC-I may make good immunotherapeutic targets. De novo peptide sequencing is a great way to find ncMAPs since it can detect peptide sequences from their tandem mass spectra without using any sequence databases. However, this strategy has not been widely applied for ncMAP identification because there is not a good way to estimate its false-positive rates. In order to completely and accurately identify immunopeptides using de novo peptide sequencing, we describe a unique pipeline called proteomics X genomics. In contrast to current pipelines, it makes use of genomic data, RNA-Seq abundance and sequencing quality, in addition to proteomic features to increase the sensitivity and specificity of peptide identification. We show that the peptide-spectrum match quality and genetic traits have a clear relationship, showing that they can be utilized to evaluate peptide-spectrum matches. From 10 samples, we found 24,449 canonical MHC-I-associated peptides and 956 ncMAPs by using a target-decoy competition. Three hundred eighty-seven ncMAPs and 1611 canonical MHC-I-associated peptides were new identifications that had not yet been published. We discovered 11 ncMAPs produced from a squirrel monkey retrovirus in human cell lines in addition to the two ncMAPs originating from a complementarity determining region 3 in an antibody thanks to the unrestricted search space assumed by de novo sequencing. These entirely new identifications show that proteomics X genomics can make the most of de novo peptide sequencing's advantages and its potential use in the search for new immunotherapeutic targets.

COSMIC-based mutation database enhances identification efficiency of HLA-I immunopeptidome

BACKGROUND

Neoantigens have emerged as a promising area of focus in tumor immunotherapy, with several established strategies aiming to enhance their identification. Human leukocyte antigen class I molecules (HLA-I), which present intracellular immunopeptides to T cells, provide an ideal source for identifying neoantigens. However, solely relying on a mutation database generated through commonly used whole exome sequencing (WES) for the identification of HLA-I immunopeptides, may result in potential neoantigens being missed due to limitations in sequencing depth and sample quality.

METHOD

In this study, we constructed and evaluated an extended database for neoantigen identification, based on COSMIC mutation database. This study utilized mass spectrometry-based proteogenomic profiling to identify the HLA-I immunopeptidome enriched from HepG2 cell. HepG2 WES-based and the COSMIC-based mutation database were generated and utilized to identify HepG2-specific mutant immunopeptides.

RESULT

The results demonstrated that COSMIC-based database identified 5 immunopeptides compared to only 1 mutant peptide identified by HepG2 WES-based database, indicating its effectiveness in identifying mutant immunopeptides. Furthermore, HLA-I affinity of the mutant immunopeptides was evaluated through NetMHCpan and peptide-docking modeling to validate their binding to HLA-I molecules, demonstrating the potential of mutant peptides identified by the COSMIC-based database as neoantigens.

CONCLUSION

Utilizing the COSMIC-based mutation database is a more efficient strategy for identifying mutant peptides from HLA-I immunopeptidome without significantly increasing the false positive rate. HepG2 specific WES-based database may exclude certain mutant peptides due to WES sequencing depth or sample heterogeneity. The COSMIC-based database can effectively uncover potential neoantigens within the HLA-I immunopeptidomes.

Profound N-glycan remodelling accompanies MHC-II immunopeptide presentation

Immunopeptidomics, the study of peptide antigens presented on the cell surface by the major histocompatibility complex (MHC), offers insights into how our immune system recognises self/non-self in health and disease. We recently discovered that hyper-processed (remodelled) N-glycans are dominant features decorating viral spike immunopeptides presented via MHC-class II (MHC-II) molecules by dendritic cells pulsed with SARS-CoV-2 spike protein, but it remains unknown if endogenous immunopeptides also undergo N-glycan remodelling. Taking a multi-omics approach, we here interrogate published MHC-II immunopeptidomics datasets of cultured monocyte-like (THP-1) and breast cancer-derived (MDA-MB-231) cell lines for overlooked N-glycosylated peptide antigens, which we compare to their source proteins in the cellular glycoproteome using proteomics and N-glycomics data from matching cell lines. Hyper-processed chitobiose core and paucimannosidic N-glycans alongside under-processed oligomannosidic N-glycans were found to prevalently modify MHC-II-bound immunopeptides isolated from both THP-1 and MDA-MB-231, while complex/hybrid-type N-glycans were (near-)absent in the immunopeptidome as supported further by new N-glycomics data generated from isolated MHC-II-bound peptides derived from MDA-MB-231 cells. Contrastingly, the cellular proteomics and N-glycomics data from both cell lines revealed conventional N-glycosylation rich in complex/hybrid-type N-glycans, which, together with the identification of key lysosomal glycosidases, suggest that MHC-II peptide antigen processing is accompanied by extensive N-glycan trimming. N-glycan remodelling appeared particularly dramatic for cell surface-located glycoproteins while less remodelling was observed for lysosomal-resident glycoproteins. Collectively, our findings indicate that both under- and hyper-processed N-glycans are prevalent features of endogenous MHC-II immunopeptides, an observation that demands further investigation to enable a better molecular-level understanding of immune surveillance.

Identification of neoepitope reactive T-cell receptors guided by HLA-A*03:01 and HLA-A*11:01 immunopeptidomics

BACKGROUND

Tumor-specific mutated proteins can create immunogenic non-self, mutation-containing 'neoepitopes' that are attractive targets for adoptive T-cell therapies. To avoid the complexity of defining patient-specific, private neoepitopes, there has been major interest in targeting common shared mutations in driver genes using off-the-shelf T-cell receptors (TCRs) engineered into autologous lymphocytes. However, identifying the precise naturally processed neoepitopes to pursue is a complex and challenging process. One method to definitively demonstrate whether an epitope is presented at the cell surface is to elute peptides bound to a specific major histocompatibility complex (MHC) allele and analyze them by mass spectrometry (MS). These MS data can then be prospectively applied to isolate TCRs specific to the neoepitope.

METHODS

We created mono-allelic cell lines expressing one class I HLA allele and one common mutated oncogene in order to eliminate HLA deconvolution requirements and increase the signal of recovered peptides. MHC-bound peptides on the surface of these cell lines were immunoprecipitated, purified, and analyzed using liquid chromatography-tandem mass spectrometry, producing a list of mutation-containing minimal epitopes. To validate the immunogenicity of these neoepitopes, HLA-transgenic mice were vaccinated using the minimal peptides identified by MS in order to generate neoepitope-reactive TCRs. Specificity of these candidate TCRs was confirmed by peptide titration and recognition of transduced targets.

RESULTS

We identified precise neoepitopes derived from mutated isoforms of KRAS, EGFR, BRAF, and PIK3CA presented by HLA-A*03:01 and/or HLA-A*11:01 across multiple biological replicates. From our MS data, we were able to successfully isolate murine TCRs that specifically recognize four HLA-A*11:01 restricted neoepitopes (KRAS G13D, PIK3CA E545K, EGFR L858R and BRAF V600E) and three HLA-A*03:01 restricted neoepitopes (KRAS G12V, EGFR L858R and BRAF V600E).

CONCLUSIONS

Our data show that an MS approach can be used to demonstrate which shared oncogene-derived neoepitopes are processed and presented by common HLA alleles, and those MS data can rapidly be used to develop TCRs against these common tumor-specific antigens. Although further characterization of these neoepitope-specific murine TCRs is required, ultimately, they have the potential to be used clinically for adoptive cell therapy.

Variant biomarker discovery using mass spectrometry-based proteogenomics

Genomic diversity plays critical roles in risk of disease pathogenesis and diagnosis. While genomic variants-including single nucleotide variants, frameshift variants, and mis-splicing isoforms-are commonly detected at the DNA or RNA level, their translated variant protein or polypeptide products are ultimately the functional units of the associated disease. These products are often released in biofluids and could be leveraged for clinical diagnosis and patient stratification. Recent emergence of integrated analysis of genomics with mass spectrometry-based proteomics for biomarker discovery, also known as proteogenomics, have significantly advanced the understanding disease risk variants, precise medicine, and biomarker discovery. In this review, we discuss variant proteins in the context of cancers and neurodegenerative diseases, outline current and emerging proteogenomic approaches for biomarker discovery, and provide a comprehensive proteogenomic strategy for detection of putative biomarker candidates in human biospecimens. This strategy can be implemented for proteogenomic studies in any field of enquiry. Our review timely addresses the need of biomarkers for aging related diseases.

Workflow enabling deepscale immunopeptidome, proteome, ubiquitylome, phosphoproteome, and acetylome analyses of sample-limited tissues

Serial multi-omic analysis of proteome, phosphoproteome, and acetylome provides insights into changes in protein expression, cell signaling, cross-talk and epigenetic pathways involved in disease pathology and treatment. However, ubiquitylome and HLA peptidome data collection used to understand protein degradation and antigen presentation have not together been serialized, and instead require separate samples for parallel processing using distinct protocols. Here we present MONTE, a highly sensitive multi-omic native tissue enrichment workflow, that enables serial, deep-scale analysis of HLA-I and HLA-II immunopeptidome, ubiquitylome, proteome, phosphoproteome, and acetylome from the same tissue sample. We demonstrate that the depth of coverage and quantitative precision of each 'ome is not compromised by serialization, and the addition of HLA immunopeptidomics enables the identification of peptides derived from cancer/testis antigens and patient specific neoantigens. We evaluate the technical feasibility of the MONTE workflow using a small cohort of patient lung adenocarcinoma tumors.

Long non-coding RNA-derived peptides are immunogenic and drive a potent anti-tumour response

Protein arginine methyltransferase (PRMT) 5 is over-expressed in a variety of cancers and the master transcription regulator E2F1 is an important methylation target. We have explored the role of PRMT5 and E2F1 in regulating the non-coding genome and report here a striking effect on long non-coding (lnc) RNA gene expression. Moreover, many MHC class I protein-associated peptides were derived from small open reading frames in the lncRNA genes. Pharmacological inhibition of PRMT5 or adjusting E2F1 levels qualitatively altered the repertoire of lncRNA-derived peptide antigens displayed by tumour cells. When presented to the immune system as either ex vivo-loaded dendritic cells or expressed from a viral vector, lncRNA-derived peptides drove a potent antigen-specific CD8 T lymphocyte response, which translated into a significant delay in tumour growth. Thus, lncRNA genes encode immunogenic peptides that can be deployed as a cancer vaccine.

Correlation of tumor mutational burden with prognosis and immune infiltration in lung adenocarcinoma

Background

Tumor mutational burden (TMB) plays an important role in the evaluation of immunotherapy efficacy in lung adenocarcinoma (LUAD).

Objective

To improve the clinical management of LUAD by investigating the prognostic value of TMB and the relationship between TMB and immune infiltration.

Methods

TMB scores were calculated from the mutation data of 587 LUAD samples from The Cancer Genome Atlas (TCGA), and patients were divided into low-TMB and high-TMB groups based on the quartiles of the TMB score. Differentially expressed genes (DEGs), immune cell infiltration and survival analysis were compared between the low-TMB and high-TMB groups. We queried the expression of genes in lung cancer tissues through the GEPIA online database and performed experimental validation of the function of aberrant genes expressed in lung cancer tissues.

Results

We obtained sample information from TCGA for 587 LUAD patients, and the results of survival analysis for the high- and low- TMB groups suggested that patients in the high-TMB group had lower survival rates than those in the low-TMB group. A total of 756 DEGs were identified in the study, and gene set enrichment analysis (GSEA) showed that DEGs in the low-TMB group were enriched in immune-related pathways. Among the differentially expressed genes obtained, 15 immune-related key genes were screened with the help of ImmPort database, including 5 prognosis-related genes (CD274, PDCD1, CTLA4, LAG3, TIGIT). No difference in the expression of PDCD1, CTLA4, LAG3, TIGIT in lung cancer tissues and differential expression of CD274 in lung cancer tissues.

Conclusions

The survival rate of LUAD patients with low TMB was better than that of LUAD patients with high TMB. CD274 expression was down regulated in human LUAD cell lines H1299, PC-9, A549 and SPC-A1, which inhibited malignant progression of A549 cells.

Proteogenomic identification of Hepatitis B virus (HBV) genotype-specific HLA-I restricted peptides from HBV-positive patient liver tissues

The presentation of virus-derived peptides by HLA class I molecules on the surface of an infected cell and the recognition of these HLA-peptide complexes by, and subsequent activation of, CD8⁺ cytotoxic T cells provides an important mechanism for immune protection against viruses. Recent advances in proteogenomics have allowed researchers to discover a growing number of unique HLA-restricted viral peptides, resulting in a rapidly expanding repertoire of targets for immunotherapeutics (i.e. bispecific antibodies, engineered T-cell receptors (TCRs), chimeric antigen receptor T-cells (CAR-Ts)) to infected tissues. However, genomic variability between viral strains, such as Hepatitis-B virus (HBV), in combination with differences in patient HLA alleles, make it difficult to develop therapeutics against these targets. To address this challenge, we developed a novel proteogenomics approach for generating patient-specific databases that enable the identification of viral peptides based on the viral transcriptomes sequenced from individual patient liver samples. We also utilized DNA sequencing of patient samples to identify HLA genotypes and assist in target selection. Liver samples from 48 HBV infected patients, primarily from Asia, were examined to reconstruct patient-specific HBV genomes, identify regions within the human chromosomes targeted by HBV integrations and obtain a comprehensive view of HBV peptide epitopes using our HLA class-I (HLA-I) immunopeptidomics discovery platform. Two previously reported HLA associated HBV-derived peptides, HLA-A02 binder FLLTRILTI (S_194-202) from the large surface antigen and HLA-A11 binder STLPETTVVRR (C_141-151) from the capsid protein were validated by our discovery platform, but both were detected at very low frequencies. In addition, we identified and validated, using heavy peptide analogues, novel strain-specific HBV-HLA associated peptides, such as GSLPQEHIVQK (P_606-616) and variants. Overall, our novel approach can guide the development of bispecific antibody, TCR-T, or CAR-T based therapeutics for the treatment of HBV-related HCC and inform vaccine development.

Scoring model based on the signature of non-m6A-related neoantigen-coding lncRNAs assists in immune microenvironment analysis and TCR-neoantigen pair selection in gliomas

Background

Small peptides encoded by long non-coding RNAs (lncRNAs) have attracted attention for their various functions. Recent studies indicate that these small peptides participate in immune responses and antigen presentation. However, the significance of RNA modifications remains unclear.

Methods

Thirteen non-m6A-related neoantigen-coding lncRNAs were selected for analysis from the TransLnc database. Next, a neoantigen activation score (NAS) model was established based on the characteristics of the lncRNAs. Machine learning was employed to expand the model to two additional RNA-seq and two single-cell sequencing datasets for further validation. The DLpTCR algorithm was used to predict T cell receptor (TCR)-peptide binding probability.

Results

The non-m6A-related NAS model predicted patients’ overall survival outcomes more precisely than the m6A-related NAS model. Furthermore, the non-m6A-related NAS was positively correlated with tumor cells’ evolutionary level, immune infiltration, and antigen presentation. However, high NAS gliomas also showed more PD-L1 expression and high mutation frequencies of T-cell positive regulators. Interestingly, results of intercellular communication analysis suggest that T cell-high neoplastic cell interaction is weaker in both of the NAS groups which might arise from decreased IFNGR1 expression. Moreover, we identified unique TCR-peptide pairs present in all glioma samples based on peptides encoded by the 13 selected lncRNAs. And increased levels of neoantigen-active TCR patterns were found in high NAS gliomas.

Conclusions

Our work suggests that non-m6A-related neoantigen-coding lncRNAs play an essential role in glioma progression and that screened TCR clonotypes might provide potential avenues for chimeric antigen receptor T cell (CAR-T) therapy for gliomas.

Soluble HLA peptidome of pleural effusions is a valuable source for tumor antigens

Background

Soluble human leucocyte antigen (sHLA) molecules, released into the plasma, carry their original peptide cargo and provide insight into the protein synthesis and degradation schemes of their source cells and tissues. Other body fluids, such as pleural effusions, may also contain sHLA-peptide complexes, and can potentially serve as a source of tumor antigens since these fluids are drained from the tumor microenvironment. We explored this possibility by developing a methodology for purifying and analyzing large pleural effusion sHLA class I peptidomes of patients with malignancies or benign diseases.

Methods

Cleared pleural fluids, cell pellets present in the pleural effusions, and the primary tumor cells cultured from cancer patients’ effusions, were used for immunoaffinity purification of the HLA molecules. The recovered HLA peptides were analyzed by liquid chromatography coupled to tandem mass spectrometry (LC-MS/MS) and the resulting LC-MS/MS data were analyzed with the MaxQuant software tool. Selected tumor antigen peptides were tested for their immunogenicity potential with donor peripheral blood mononuclear cells (PBMCs) in an in vitro assay.

Results

Mass spectrometry analysis of the pleural effusions revealed 39,669 peptides attributable to 11,305 source proteins. The majority of peptides identified from the pleural effusions were defined as HLA ligands that fit the patients’ HLA consensus sequence motifs. The membranal and soluble HLA peptidomes of each individual patient correlated to each other. Additionally, soluble HLA peptidomes from the same patient, obtained at different visits to the clinic, were highly similar. Compared with benign effusions, the soluble HLA peptidomes of malignant pleural effusions were larger and included HLA peptides derived from known tumor-associated antigens, including cancer/testis antigens, lung-related proteins, and vascular endothelial growth factor pathway proteins. Selected tumor-associated antigens that were identified by the immunopeptidomics were able to successfully prime CD8⁺ T cells.

Conclusions

Pleural effusions contain sHLA-peptide complexes, and the pleural effusion HLA peptidome of patients with malignant tumors can serve as a rich source of biomarkers for tumor diagnosis and potential candidates for personalized immunotherapy.

Alterations in HLA Class I-Presented Immunopeptidome and Class I-Interactome upon Osimertinib Resistance in EGFR Mutant Lung Adenocarcinoma

Simple Summary

We sought to identify molecular mechanisms of lower efficacy of immunotherapy in epidermal growth factor receptor (EGFR) mutant lung adenocarcinoma and the differences in those mechanisms with the emergence of tyrosine kinase inhibitor (TKI)-resistance. To this end, we conducted affinity purification and quantitative mass spectrometry-based proteomic profiling of human leukocyte antigen (HLA) Class I-presented immunopeptides and Class I-interacting proteins. This large-scale dataset revealed that the Class I-presented immunopeptidome was suppressed in two third-generation EGFR TKI, osimertinib-resistant lung adenocarcinoma cell lines compared to their isogenic TKI-sensitive counterparts. The whole-cell proteomic profiling show that antigen presentation complex proteins and immunoproteasome were downregulated upon EGFR TKI resistance. Furthermore, HLA class I-interactome profiling demonstrated altered interaction with key apoptosis and autophagy pathway proteins. In summary, our comprehensive multi-proteomic characterization in antigen presentation machinery provides potentially novel evidence of poor immune response in osimertinib-resistant lung adenocarcinoma.

Abstract

Immune checkpoint inhibitor (ICI) therapy has been a paradigm shift in the treatment of cancer. ICI therapy results in durable responses and survival benefit for a large number of tumor types. Osimertinib, a third-generation epidermal growth factor receptor (EGFR) tyrosine kinase inhibitor (TKI) has shown great efficacy treating EGFR mutant lung cancers; however, all patients eventually develop resistance. ICI therapy has not benefitted EGFR mutant lung cancer. Herein, we employed stable isotope labeling by amino acids in cell culture (SILAC) quantitative mass spectrometry-based proteomics to investigate potential immune escape molecular mechanisms in osimertinib resistant EGFR mutant lung adenocarcinoma by interrogating the alterations in the human leukocyte antigen (HLA) Class I-presented immunopeptidome, Class I-interactome, and the whole cell proteome between isogenic osimertinib-sensitive and -resistant human lung adenocarcinoma cells. Our study demonstrates an overall reduction in HLA class I-presented immunopeptidome and downregulation of antigen presentation core complex (e.g., TAP1 and ERAP1/2) and immunoproteasome in osimertinib resistant lung adenocarcinoma cells. Several key components in autophagy pathway are differentially altered. S100 proteins and SLC3A2 may play critical roles in reduced antigen presentation. Our dataset also includes ~1000 novel HLA class I interaction partners and hundreds of Class I-presented immunopeptides in EGFR mutant lung adenocarcinoma. This large-scale unbiased proteomics study provides novel insights and potential mechanisms of immune evasion of EGFR mutant lung adenocarcinoma.