The Choice of Search Engine Affects Sequencing Depth and HLA Class I Allele-Specific Peptide Repertoires

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.

Immunopeptidomics-based identification of naturally presented non-canonical circRNA-derived peptides

Circular RNAs (circRNAs) are covalently closed non-coding RNAs lacking the 5' cap and the poly-A tail. Nevertheless, it has been demonstrated that certain circRNAs can undergo active translation. Therefore, aberrantly expressed circRNAs in human cancers could be an unexplored source of tumor-specific antigens, potentially mediating anti-tumor T cell responses. This study presents an immunopeptidomics workflow with a specific focus on generating a circRNA-specific protein fasta reference. The main goal of this workflow is to streamline the process of identifying and validating human leukocyte antigen (HLA) bound peptides potentially originating from circRNAs. We increase the analytical stringency of our workflow by retaining peptides identified independently by two mass spectrometry search engines and/or by applying a group-specific FDR for canonical-derived and circRNA-derived peptides. A subset of circRNA-derived peptides specifically encoded by the region spanning the back-splice junction (BSJ) are validated with targeted MS, and with direct Sanger sequencing of the respective source transcripts. Our workflow identifies 54 unique BSJ-spanning circRNA-derived peptides in the immunopeptidome of melanoma and lung cancer samples. Our approach enlarges the catalog of source proteins that can be explored for immunotherapy.

UniSpec: Deep Learning for Predicting the Full Range of Peptide Fragment Ion Series to Enhance the Proteomics Data Analysis Workflow

We present UniSpec, an attention-driven deep neural network designed to predict comprehensive collision-induced fragmentation spectra, thereby improving peptide identification in shotgun proteomics. Utilizing a training data set of 1.8 million unique high-quality tandem mass spectra (MS2) from 0.8 million unique peptide ions, UniSpec learned with a peptide fragmentation dictionary encompassing 7919 fragment peaks. Among these, 5712 are neutral loss peaks, with 2310 corresponding to modification-specific neutral losses. Remarkably, UniSpec can predict 73%-77% of fragment intensities based on our NIST reference library spectra, a significant leap from the 35%-45% coverage of only b and y ions. Comparative studies with Prosit elucidate that while both models are strong at predicting their respective fragment ion series, UniSpec particularly shines in generating more complex MS2 spectra with diverse ion annotations. The integration of UniSpec's predictions into shotgun proteomics data analysis boosts the identification rate of tryptic peptides by 48% at a 1% false discovery rate (FDR) and 60% at a more confident 0.1% FDR. Using UniSpec's predicted in-silico spectral library, the search results closely matched those from search engines and experimental spectral libraries used in peptide identification, highlighting its potential as a stand-alone identification tool. The source code and Python scripts are available on GitHub (https://github.com/usnistgov/UniSpec) and Zenodo (https://zenodo.org/records/10452792), and all data sets and analysis results generated in this work were deposited in Zenodo (https://zenodo.org/records/10052268).

Protein degradation by human 20S proteasomes elucidates the interplay between peptide hydrolysis and splicing

If and how proteasomes catalyze not only peptide hydrolysis but also peptide splicing is an open question that has divided the scientific community. The debate has so far been based on immunopeptidomics, in vitro digestions of synthetic polypeptides as well as ex vivo and in vivo experiments, which could only indirectly describe proteasome-catalyzed peptide splicing of full-length proteins. Here we develop a workflow-and cognate software - to analyze proteasome-generated non-spliced and spliced peptides produced from entire proteins and apply it to in vitro digestions of 15 proteins, including well-known intrinsically disordered proteins such as human tau and α-Synuclein. The results confirm that 20S proteasomes produce a sizeable variety of cis-spliced peptides, whereas trans-spliced peptides are a minority. Both peptide hydrolysis and splicing produce peptides with well-defined characteristics, which hint toward an intricate regulation of both catalytic activities. At protein level, both non-spliced and spliced peptides are not randomly localized within protein sequences, but rather concentrated in hotspots of peptide products, in part driven by protein sequence motifs and proteasomal preferences. At sequence level, the different peptide sequence preference of peptide hydrolysis and peptide splicing suggests a competition between the two catalytic activities of 20S proteasomes during protein degradation.

MARS an improved de novo peptide candidate selection method for non-canonical antigen target discovery in cancer

Understanding the nature and extent of non-canonical human leukocyte antigen (HLA) presentation in tumour cells is a priority for target antigen discovery for the development of next generation immunotherapies in cancer. We here employ a de novo mass spectrometric sequencing approach with a refined, MHC-centric analysis strategy to detect non-canonical MHC-associated peptides specific to cancer without any prior knowledge of the target sequence from genomic or RNA sequencing data. Our strategy integrates MHC binding rank, Average local confidence scores, and peptide Retention time prediction for improved de novo candidate Selection; culminating in the machine learning model MARS. We benchmark our model on a large synthetic peptide library dataset and reanalysis of a published dataset of high-quality non-canonical MHC-associated peptide identifications in human cancer. We achieve almost 2-fold improvement for high quality spectral assignments in comparison to de novo sequencing alone with an estimated accuracy of above 85.7% when integrated with a stepwise peptide sequence mapping strategy. Finally, we utilize MARS to detect and validate lncRNA-derived peptides in human cervical tumour resections, demonstrating its suitability to discover novel, immunogenic, non-canonical peptide sequences in primary tumour tissue.

Deep Learning-Assisted Analysis of Immunopeptidomics Data

Liquid chromatography-coupled mass spectrometry (LC-MS/MS) is the primary method to obtain direct evidence for the presentation of disease- or patient-specific human leukocyte antigen (HLA). However, compared to the analysis of tryptic peptides in proteomics, the analysis of HLA peptides still poses computational and statistical challenges. Recently, fragment ion intensity-based matching scores assessing the similarity between predicted and observed spectra were shown to substantially increase the number of confidently identified peptides, particularly in use cases where non-tryptic peptides are analyzed. In this chapter, we describe in detail three procedures on how to benefit from state-of-the-art deep learning models to analyze and validate single spectra, single measurements, and multiple measurements in mass spectrometry-based immunopeptidomics. For this, we explain how to use the Universal Spectrum Explorer (USE), online Oktoberfest, and offline Oktoberfest. For intensity-based scoring, Oktoberfest uses fragment ion intensity and retention time predictions from the deep learning framework Prosit, a deep neural network trained on a very large number of synthetic peptides and tandem mass spectra generated within the ProteomeTools project. The examples shown highlight how deep learning-assisted analysis can increase the number of identified HLA peptides, facilitate the discovery of confidently identified neo-epitopes, or provide assistance in the assessment of the presence of cryptic peptides, such as spliced peptides.

Current perspectives on mass spectrometry-based immunopeptidomics: the computational angle to tumor antigen discovery

Identification of tumor antigens presented by the human leucocyte antigen (HLA) molecules is essential for the design of effective and safe cancer immunotherapies that rely on T cell recognition and killing of tumor cells. Mass spectrometry (MS)-based immunopeptidomics enables high-throughput, direct identification of HLA-bound peptides from a variety of cell lines, tumor tissues, and healthy tissues. It involves immunoaffinity purification of HLA complexes followed by MS profiling of the extracted peptides using data-dependent acquisition, data-independent acquisition, or targeted approaches. By incorporating DNA, RNA, and ribosome sequencing data into immunopeptidomics data analysis, the proteogenomic approach provides a powerful means for identifying tumor antigens encoded within the canonical open reading frames of annotated coding genes and non-canonical tumor antigens derived from presumably non-coding regions of our genome. We discuss emerging computational challenges in immunopeptidomics data analysis and tumor antigen identification, highlighting key considerations in the proteogenomics-based approach, including accurate DNA, RNA and ribosomal sequencing data analysis, careful incorporation of predicted novel protein sequences into reference protein database, special quality control in MS data analysis due to the expanded and heterogeneous search space, cancer-specificity determination, and immunogenicity prediction. The advancements in technology and computation is continually enabling us to identify tumor antigens with higher sensitivity and accuracy, paving the way toward the development of more effective cancer immunotherapies.

MSBooster: improving peptide identification rates using deep learning-based features

Peptide identification in liquid chromatography-tandem mass spectrometry (LC-MS/MS) experiments relies on computational algorithms for matching acquired MS/MS spectra against sequences of candidate peptides using database search tools, such as MSFragger. Here, we present a new tool, MSBooster, for rescoring peptide-to-spectrum matches using additional features incorporating deep learning-based predictions of peptide properties, such as LC retention time, ion mobility, and MS/MS spectra. We demonstrate the utility of MSBooster, in tandem with MSFragger and Percolator, in several different workflows, including nonspecific searches (immunopeptidomics), direct identification of peptides from data independent acquisition data, single-cell proteomics, and data generated on an ion mobility separation-enabled timsTOF MS platform. MSBooster is fast, robust, and fully integrated into the widely used FragPipe computational platform.

Novel canonical and non-canonical viral antigens extend current targets for immunotherapy of HPV-driven cervical cancer

Current immunotherapeutic approaches for human papillomavirus (HPV)-driven cervical cancer target the viral oncogenes E6 and E7. We report viral canonical and alternative reading frame (ARF)-derived sequences presented on cervical tumor cells, including antigens encoded by the conserved viral gene E1. We confirm immunogenicity of the identified viral peptides in HPV-positive women, and women with cervical intraepithelial neoplasia. We observe consistent transcription of the E1, E6, and E7 genes in 10 primary cervical tumor resections from the four most common high-risk HPV subtypes (HPV16, 18, 31, and 45), suggesting the suitability of E1 as therapeutic target. We finally confirm HLA presentation of canonical peptides derived from E6 and E7, and ARF-derived viral peptides from a reverse-strand transcript spanning the HPV E1 and E2 genes in primary human cervical tumor tissue. Our results extend currently known viral immunotherapeutic targets in cervical cancer and highlight E1 as an important cervical cancer antigen.

Antigen discovery for the development of cancer immunotherapy

Central to successful cancer immunotherapy is effective T cell antitumor immunity. Multiple targeted immunotherapies engineered to invigorate T cell-driven antitumor immunity rely on identifying the repertoire of T cell antigens expressed on the tumor cell surface. Mass spectrometry-based survey of such antigens ("immunopeptidomics") combined with other omics platforms and computational algorithms has been instrumental in identifying and quantifying tumor-derived T cell antigens. In this review, we discuss the types of tumor antigens that have emerged for targeted cancer immunotherapy and the immunopeptidomics methods that are central in MHC peptide identification and quantification. We provide an overview of the strength and limitations of mass spectrometry-driven approaches and how they have been integrated with other technologies to discover targetable T cell antigens for cancer immunotherapy. We highlight some of the emerging cancer immunotherapies that successfully capitalized on immunopeptidomics, their challenges, and mass spectrometry-based strategies that can support their development.

Ionizing Radiation Drives Key Regulators of Antigen Presentation and a Global Expansion of the Immunopeptidome

Little is known about the pathways regulating MHC antigen presentation and the identity of treatment-specific T cell antigens induced by ionizing radiation. For this reason, we investigated the radiation-specific changes in the colorectal tumor cell proteome. We found an increase in DDX58 and ZBP1 protein expression, two nucleic acid sensing molecules likely involved in induction of the dominant interferon response signature observed after genotoxic insult. We further observed treatment-induced changes in key regulators and effector proteins of the antigen processing and presentation machinery. Differential regulation of MHC allele expression was further driving the presentation of a significantly broader MHC-associated peptidome postirradiation, defining a radiation-specific peptide repertoire. Interestingly, treatment-induced peptides originated predominantly from proteins involved in catecholamine synthesis and metabolic pathways. A nuanced relationship between protein expression and antigen presentation was observed where radiation-induced changes in proteins do not correlate with increased presentation of associated peptides. Finally, we detected an increase in the presentation of a tumor-specific neoantigen derived from Mtch1. This study provides new insights into how radiation enhances antigen processing and presentation that could be suitable for the development of combinatorial therapies. Data are available via ProteomeXchange with identifier PXD032003.

New insights and approaches for analyses of immunopeptidomes

Human leucocyte antigen (HLA) molecules play a key role in health and disease by presenting antigen to T-lymphocytes for immunosurveillance. Immunopeptidomics involves the study of the collection of peptides presented within the antigen-binding groove of HLA molecules. Identifying their nature and diversity is crucial to understanding immunosurveillance especially during infection or for the recognition and potential eradication of tumours. This review discusses recent advances in the isolation, identification, and quantitation of these peptide antigens. New informatics approaches and databases have shed light on the extent of peptide antigens derived from unconventional sources including peptides derived from transcripts associated with frame shifts, long noncoding RNA, incorrectly annotated untranslated regions, post-translational modifications, and proteasomal splicing. Several challenges remain in successful analysis of immunopeptides, yet recent developments point to unexplored biology waiting to be unravelled.

T Cell Epitope Discovery in the Context of Distinct and Unique Indigenous HLA Profiles

CD8⁺ T cells are a pivotal part of the immune response to viruses, playing a key role in disease outcome and providing long-lasting immunity to conserved pathogen epitopes. Understanding CD8⁺ T cell immunity in humans is complex due to CD8⁺ T cell restriction by highly polymorphic Human Leukocyte Antigen (HLA) proteins, requiring T cell epitopes to be defined for different HLA allotypes across different ethnicities. Here we evaluate strategies that have been developed to facilitate epitope identification and study immunogenic T cell responses. We describe an immunopeptidomics approach to sequence HLA-bound peptides presented on virus-infected cells by liquid chromatography with tandem mass spectrometry (LC-MS/MS). Using antigen presenting cell lines that stably express the HLA alleles characteristic of Indigenous Australians, this approach has been successfully used to comprehensively identify influenza-specific CD8⁺ T cell epitopes restricted by HLA allotypes predominant in Indigenous Australians, including HLA-A*24:02 and HLA-A*11:01. This is an essential step in ensuring high vaccine coverage and efficacy in Indigenous populations globally, known to be at high risk from influenza disease and other respiratory infections.

Database search engines and target database features impinge upon the identification of post-translationally cis-spliced peptides in HLA class I immunopeptidomes

Unconventional epitopes presented by HLA class I complexes are emerging targets for T cell targeted immunotherapies. Their identification by mass spectrometry (MS) required development of novel methods to cope with the large number of theoretical candidates. Methods to identify post-translationally spliced peptides led to a broad range of outcomes. We here investigated the impact of three common database search engines - that is, Mascot, Mascot+Percolator, and PEAKS DB - as final identification step, as well as the features of target database on the ability to correctly identify non-spliced and cis-spliced peptides. We used ground truth datasets measured by MS to benchmark methods' performance and extended the analysis to HLA class I immunopeptidomes. PEAKS DB showed better precision and recall of cis-spliced peptides and larger number of identified peptides in HLA class I immunopeptidomes than the other search engine strategies. The better performance of PEAKS DB appears to result from better discrimination between target and decoy hits and hence a more robust FDR estimation, and seems independent to peptide and spectrum features here investigated.

The interdependence of machine learning and LC-MS approaches for an unbiased understanding of the cellular immunopeptidome

INTRODUCTION

The comprehensive collection of peptides presented by Major Histocompatibility Complex (MHC) molecules on the cell surface is collectively known as the immunopeptidome. The analysis and interpretation of such data sets holds great promise for furthering our understanding of basic immunology and adaptive immune activation and regulation, and for direct rational discovery of T cell antigens and the design of T-cell based therapeutics and vaccines. These applications are however challenged by the complex nature of immunopeptidome data.

AREAS COVERED

Here, we describe the benefits and shortcomings of applying liquid chromatography-tandem mass spectrometry (MS) to obtain large scale immunopeptidome data sets and illustrate how the accurate analysis and optimal interpretation of such data is reliant on the availability of refined and highly optimized machine learning approaches.

EXPERT OPINION

Further we demonstrate how the accuracy of immunoinformatics prediction methods within the field of MHC antigen presentation has benefited greatly from the availability of MS-immunopeptidomics data, and exemplify how optimal antigen discovery is best performed in a synergistic combination of MS experiments and such in silico models trained on large scale immunopeptidomics data.

Immunopeptidogenomics: Harnessing RNA-Seq to Illuminate the Dark Immunopeptidome

Human leukocyte antigen (HLA) molecules are cell-surface glycoproteins that present peptide antigens on the cell surface for surveillance by T lymphocytes, which contemporaneously seek signs of disease. Mass spectrometric analysis allows us to identify large numbers of these peptides (the immunopeptidome) following affinity purification of solubilized HLA-peptide complexes. However, in recent years, there has been a growing awareness of the "dark side" of the immunopeptidome: unconventional peptide epitopes, including neoepitopes, which elude detection by conventional search methods because their sequences are not present in reference protein databases (DBs). Here, we establish a bioinformatics workflow to aid identification of peptides generated by noncanonical translation of mRNA or by genome variants. The workflow incorporates both standard transcriptomics software and novel computer programs to produce cell line-specific protein DBs based on three-frame translation of the transcriptome. The final protein DB also includes sequences resulting from variants determined by variant calling on the same RNA-Seq data. We then searched our experimental data against both transcriptome-based and standard DBs using PEAKS Studio (Bioinformatics Solutions, Inc). Finally, further novel software helps to compare the various result sets arising for each sample, pinpoint putative genomic origins for unconventional sequences, and highlight potential neoepitopes. We applied the workflow to study the immunopeptidome of the acute myeloid leukemia cell line THP-1, using RNA-Seq and immunopeptidome data. We confidently identified over 14,000 peptides from three replicates of purified HLA peptides derived from THP-1 cells using the conventional UniProt human proteome. Using the transcriptome-based DB generated using our workflow, we recapitulated >85% of these and also identified 1029 unconventional peptides not explained by UniProt, including 16 sequences caused by nonsynonymous variants. Our workflow, which we term "immunopeptidogenomics," can provide DBs, which include pertinent unconventional sequences and allow neoepitope discovery, without becoming too large to search. Immunopeptidogenomics is a step toward unbiased search approaches that are needed to illuminate the dark side of the immunopeptidome.