ProGeo-neo: a customized proteogenomic workflow for neoantigen prediction and selection

A modular, Galaxy-based immunopeptidogenomic (iPepGen) analysis pipeline for discovery, verification, and prioritization of candidate cancer neoantigen peptides

Background

Characterizing peptide antigens, processed from tumor-specific proteoforms, and bound to the major histocompatibility complex, is critical for immuno-oncology research. Next-generation sequencing predicts candidate neoantigen peptides derived from DNA mutations and/or RNA transcripts coding proteoform sequences that differ from the reference proteome. Mass spectrometry (MS)-based immunopeptidomics identifies predicted, MHC-bound neoantigen peptides and other tumor antigens. This "immunopeptidogenomic" approach requires multi-omic software integration, challenging researchers with limited bioinformatics expertise and resources. As a solution, we developed the immunopeptidogenomic (iPepGen) pipeline in the Galaxy ecosystem. iPepGen is composed of five core workflow modules, available via publicly accessible, scalable Galaxy instances, accompanied by training resources to empower community adoption.

Findings

Using representative multi-omic data from malignant peripheral nerve sheath tumors, we demonstrate the operation of iPepGen modules with these functions: 1) Predict neoantigen candidates from sequencing data and generate customized protein sequence databases, including reference and non-reference neoantigen candidate sequences; 2) Discover neoantigen peptide candidates by sequence database searching of tandem mass spectrometry (MS/MS) immunopeptidomics data; 3) Verify discovered peptide candidates through a secondary peptide-centric evaluation method against the MS/MS dataset; 4) Visualize and classify the nature of verified neoantigen peptides encoded by the genome and/or transcriptome; 5) Prioritize neoantigens for further exploration and empirical validation.

Conclusions

We demonstrate the effectiveness of the iPepGen pipeline for candidate neoantigen discovery and characterization. With tools, workflows, and training resources available in the open Galaxy ecosystem, iPepGen should provide cancer researchers with a flexible and accessible informatics resource tailored to accelerating immuno-oncology studies.

Phenotyping Tumor Heterogeneity through Proteogenomics: Study Models and Challenges

Tumor heterogeneity refers to the diversity observed among tumor cells: both between different tumors (inter-tumor heterogeneity) and within a single tumor (intra-tumor heterogeneity). These cells can display distinct morphological and phenotypic characteristics, including variations in cellular morphology, metastatic potential and variability treatment responses among patients. Therefore, a comprehensive understanding of such heterogeneity is necessary for deciphering tumor-specific mechanisms that may be diagnostically and therapeutically valuable. Innovative and multidisciplinary approaches are needed to understand this complex feature. In this context, proteogenomics has been emerging as a significant resource for integrating omics fields such as genomics and proteomics. By combining data obtained from both Next-Generation Sequencing (NGS) technologies and mass spectrometry (MS) analyses, proteogenomics aims to provide a comprehensive view of tumor heterogeneity. This approach reveals molecular alterations and phenotypic features related to tumor subtypes, potentially identifying therapeutic biomarkers. Many achievements have been made; however, despite continuous advances in proteogenomics-based methodologies, several challenges remain: in particular the limitations in sensitivity and specificity and the lack of optimal study models. This review highlights the impact of proteogenomics on characterizing tumor phenotypes, focusing on the critical challenges and current limitations of its use in different clinical and preclinical models for tumor phenotypic characterization.

Transformers meets neoantigen detection: a systematic literature review

Cancer immunology offers a new alternative to traditional cancer treatments, such as radiotherapy and chemotherapy. One notable alternative is the development of personalized vaccines based on cancer neoantigens. Moreover, Transformers are considered a revolutionary development in artificial intelligence with a significant impact on natural language processing (NLP) tasks and have been utilized in proteomics studies in recent years. In this context, we conducted a systematic literature review to investigate how Transformers are applied in each stage of the neoantigen detection process. Additionally, we mapped current pipelines and examined the results of clinical trials involving cancer vaccines.

Neoantigen-targeted TCR-engineered T cell immunotherapy: current advances and challenges

Adoptive cell therapy using T cell receptor-engineered T cells (TCR-T) is a promising approach for cancer therapy with an expectation of no significant side effects. In the human body, mature T cells are armed with an incredible diversity of T cell receptors (TCRs) that theoretically react to the variety of random mutations generated by tumor cells. The outcomes, however, of current clinical trials using TCR-T cell therapies are not very successful especially involving solid tumors. The therapy still faces numerous challenges in the efficient screening of tumor-specific antigens and their cognate TCRs. In this review, we first introduce TCR structure-based antigen recognition and signaling, then describe recent advances in neoantigens and their specific TCR screening technologies, and finally summarize ongoing clinical trials of TCR-T therapies against neoantigens. More importantly, we also present the current challenges of TCR-T cell-based immunotherapies, e.g., the safety of viral vectors, the mismatch of T cell receptor, the impediment of suppressive tumor microenvironment. Finally, we highlight new insights and directions for personalized TCR-T therapy.

PGNneo: A Proteogenomics-Based Neoantigen Prediction Pipeline in Noncoding Regions

The development of a neoantigen-based personalized vaccine has promise in the hunt for cancer immunotherapy. The challenge in neoantigen vaccine design is the need to rapidly and accurately identify, in patients, those neoantigens with vaccine potential. Evidence shows that neoantigens can be derived from noncoding sequences, but there are few specific tools for identifying neoantigens in noncoding regions. In this work, we describe a proteogenomics-based pipeline, namely PGNneo, for use in discovering neoantigens derived from the noncoding region of the human genome with reliability. In PGNneo, four modules are included: (1) noncoding somatic variant calling and HLA typing; (2) peptide extraction and customized database construction; (3) variant peptide identification; (4) neoantigen prediction and selection. We have demonstrated the effectiveness of PGNneo and applied and validated our methodology in two real-world hepatocellular carcinoma (HCC) cohorts. TP53, WWP1, ATM, KMT2C, and NFE2L2, which are frequently mutating genes associated with HCC, were identified in two cohorts and corresponded to 107 neoantigens from non-coding regions. In addition, we applied PGNneo to a colorectal cancer (CRC) cohort, demonstrating that the tool can be extended and verified in other tumor types. In summary, PGNneo can specifically detect neoantigens generated by noncoding regions in tumors, providing additional immune targets for cancer types with a low tumor mutational burden (TMB) in coding regions. PGNneo, together with our previous tool, can identify coding and noncoding region-derived neoantigens and, thus, will contribute to a complete understanding of the tumor immune target landscape. PGNneo source code and documentation are available at Github. To facilitate the installation and use of PGNneo, we provide a Docker container and a GUI.

Systems Biology Approaches for the Improvement of Oncolytic Virus-Based Immunotherapies

Oncolytic virus (OV)-based immunotherapy is mainly dependent on establishing an efficient cell-mediated antitumor immunity. OV-mediated antitumor immunity elicits a renewed antitumor reactivity, stimulating a T-cell response against tumor-associated antigens (TAAs) and recruiting natural killer cells within the tumor microenvironment (TME). Despite the fact that OVs are unspecific cancer vaccine platforms, to further enhance antitumor immunity, it is crucial to identify the potentially immunogenic T-cell restricted TAAs, the main key orchestrators in evoking a specific and durable cytotoxic T-cell response. Today, innovative approaches derived from systems biology are exploited to improve target discovery in several types of cancer and to identify the MHC-I and II restricted peptide repertoire recognized by T-cells. Using specific computation pipelines, it is possible to select the best tumor peptide candidates that can be efficiently vectorized and delivered by numerous OV-based platforms, in order to reinforce anticancer immune responses. Beyond the identification of TAAs, system biology can also support the engineering of OVs with improved oncotropism to reduce toxicity and maintain a sufficient portion of the wild-type virus virulence. Finally, these technologies can also pave the way towards a more rational design of armed OVs where a transgene of interest can be delivered to TME to develop an intratumoral gene therapy to enhance specific immune stimuli.

NAIRscore as a biomarker for the quality of immune response to neoantigens is related with an increased overall survival in multiple myeloma

Neoantigen provides a promising breakthrough in tumor immunotherapy, although only a subset of patients responds well due to the quality of their immune response. However, few biomarkers have been reported to measure the quality of immune response to neoantigens and to predict prognosis of patients with multiple myeloma (MM). Here, we first developed a neoantigen-prediction pipeline starting from outcomes of somatic mutations and gene-expression profiles. Given the expression of some specific marker genes, the human leukocyte antigen (HLA)-I score and the cytolytic score were evaluated respectively to reflect HLA-I molecular expression and CD8+ T/natural killer (NK) cell abundance. According to the process of the immune response to neoantigens, we comprehensively took neoantigen load, cytolytic score, and HLA-I score to construct a neoantigen immune response score (NAIRscore), in which the HLA-I score presented a hazard ratio (HR) of less than 1, while the cytolytic score and neoantigen load presented a HR of greater than 1. Meanwhile, NAIRscore presented a competitive advantage to stratify MM samples. Especially, those exhibiting high NAIRscore correlated with an increased overall survival (OS), echoing the underlying molecular signatures of lower driver-gene mutations and down-regulated immune response. Notably, an online tool based on this study is provided to identify neoantigens and predict OS.

dbPepNeo2.0: A Database for Human Tumor Neoantigen Peptides From Mass Spectrometry and TCR Recognition

Neoantigens are widely reported to induce T-cell response and lead to tumor regression, indicating a promising potential to immunotherapy. Previously, we constructed an open-access database, i.e., dbPepNeo, providing a systematic resource for human tumor neoantigens to storage and query. In order to expand data volume and application scope, we updated dbPepNeo to version 2.0 (http://www.biostatistics.online/dbPepNeo2). Here, we provide about 801 high-confidence (HC) neoantigens (increased by 170%) and 842,289 low-confidence (LC) HLA immunopeptidomes (increased by 107%). Notably, 55 class II HC neoantigens and 630 neoantigen-reactive T-cell receptor-β (TCRβ) sequences were firstly included. Besides, two new analytical tools are developed, DeepCNN-Ineo and BLASTdb. DeepCNN-Ineo predicts the immunogenicity of class I neoantigens, and BLASTdb performs local alignments to look for sequence similarities in dbPepNeo2.0. Meanwhile, the web features and interface have been greatly improved and enhanced.

ProGeo-Neo v2.0: A One-Stop Software for Neoantigen Prediction and Filtering Based on the Proteogenomics Strategy

A proteogenomics-based neoantigen prediction pipeline, namely ProGeo-neo, was previously developed by our team to predict neoantigens, allowing the identification of class-I major histocompatibility complex (MHC) binding peptides based on single-nucleotide variation (SNV) mutations. To improve it, we here present an updated pipeline, i.e., ProGeo-neo v2.0, in which a one-stop software solution was proposed to identify neoantigens based on the paired tumor-normal whole genome sequencing (WGS)/whole exome sequencing (WES) data in FASTQ format. Preferably, in ProGeo-neo v2.0, several new features are provided. In addition to the identification of MHC-I neoantigens, the new version supports the prediction of MHC class II-restricted neoantigens, i.e., peptides up to 30-mer in length. Moreover, the source of neoantigens has been expanded, allowing more candidate neoantigens to be identified, such as in-frame insertion-deletion (indels) mutations, frameshift mutations, and gene fusion analysis. In addition, we propose two more efficient screening approaches, including an in-group authentic neoantigen peptides database and two more stringent thresholds. The range of candidate peptides was effectively narrowed down to those that are more likely to elicit an immune response, providing a more meaningful reference for subsequent experimental validation. Compared to ProGeo-neo, the ProGeo-neo v2.0 performed well based on the same dataset, including updated functionality and improved accuracy.

Single-Cell Sequencing Identifies the Heterogeneity of CD8+ T Cells and Novel Biomarker Genes in Hepatocellular Carcinoma

CD8+ T cells are required for the establishment of antitumor immunity, and their substantial infiltration is associated with a good prognosis. However, CD8+ T cell subsets in the tumor microenvironment may play distinct roles in tumor progression, prognosis, and immunotherapy. In this study, we used the scRNA-seq data of hepatocellular carcinoma (HCC) to reveal the heterogeneity of different CD8+ T cell subsets. The scRNA-seq data set GSE149614 was obtained from the GEO database, and the transcriptome and sample phenotypic data of TCGA-LIHC were obtained from the TCGA database. CD8+ T cell subtypes and metabolic gene sets were obtained from published reports. The data processing and analysis of CD8+ T cell groups was performed by R language. The PPI network was constructed to obtain the hub genes, and the KM survival curve of the hub genes was further plotted to determine the hub genes with differences in survival. CD8+ T cells in HCC were divided into 7 subsets, and the cytotoxic CD8 T cells 4 subset showed considerable differences between the TP53-mutant and nonmutant groups, as well as between different degrees of cirrhosis, HCC grades, stages, ages, and body weights. Cytotoxic CD8 T cells 4 differential genes were analyzed by TCGA-LIHC data and single-cell sequencing data set. 10 hub genes were found: FGA, ApoA1, ApoH, AHSG, FGB, HP, TTR, TF, HPX, and APOC3. Different subsets of CD8+ T cells were found to contribute to heterogeneous prognosis and pathway activity in HCC. Alterations in the cytotoxic and immune checkpoint gene expression during CD8+ T cell differentiation were also identified. We found that cytotoxic CD8 T cells 4 is closely associated with survival and prognosis of HCC and identified four differential genes that can be used as biological markers for survival, prognosis, and clinically relevant characteristics of HCC. Results of this study could help finding targets for immunotherapy of HCC and aid in the accelerated development of immunotherapy for HCC.

Computational cancer neoantigen prediction: current status and recent advances

Over the last few decades, immunotherapy has shown significant therapeutic efficacy in a broad range of cancer types. Antitumor immune responses are contingent on the recognition of tumor-specific antigens, which are termed neoantigens. Tumor neoantigens are ideal targets for immunotherapy since they can be recognized as non-self antigens by the host immune system and thus are able to elicit an antitumor T-cell response. There are an increasing number of studies that highlight the importance of tumor neoantigens in immunoediting and in the sensitivity to immune checkpoint blockade. Therefore, one of the most fundamental tasks in the field of immuno-oncology research is the identification of patient-specific neoantigens. To this end, a plethora of computational approaches have been developed in order to predict tumor-specific aberrant peptides and quantify their likelihood of binding to patients' human leukocyte antigen molecules in order to be recognized by T cells. In this review, we systematically summarize and present the most recent advances in computational neoantigen prediction, and discuss the challenges and novel methods that are being developed to resolve them.

•

Tumors have the ability to acquire immune escape mechanisms.

•

Tumor-specific aberrant peptides (neoantigens) can elicit an immune response by the host immune system.

•

The identification of neoantigens is one of the most fundamental tasks in the field of immuno-oncology research.

•

A plethora of computational approaches have been developed in order to predict patient-specificneoantigens.

Protein-gene Expression Nexus: Comprehensive characterization of human cancer cell lines with proteogenomic analysis

Researchers have gained new therapeutic insights using multi-omics platform approaches to study DNA, RNA, and proteins of comprehensively characterized human cancer cell lines. To improve our understanding of the molecular features associated with oncogenic modulation in cancer, we proposed a proteogenomic database for human cancer cell lines, called Protein-gene Expression Nexus (PEN). We have expanded the characterization of cancer cell lines to include genetic, mRNA, and protein data of 145 cancer cell lines from various public studies. PEN contains proteomic and phosphoproteomic data on 4,129,728 peptides, 13,862 proteins, 7,138 phosphorylation site-associated genomic variations, 117 studies, and 12 cancer. We analyzed functional characterizations along with the integrated datasets, such as cis/trans association for copy number alteration (CNA), single amino acid variation for coding genes, post-translation modification site variation for Single Amino Acid Variation, and novel peptide expression for noncoding regions and fusion genes. PEN provides a user-friendly interface for searching, browsing, and downloading data and also supports the visualization of genome-wide association between CNA and expression, novel peptide landscape, mRNA-protein abundance, and functional annotation. Together, this dataset and PEN data portal provide a resource to accelerate cancer research using model cancer cell lines. PEN is freely accessible at http://combio.snu.ac.kr/pen.

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.

Prospects and challenges of cancer systems medicine: from genes to disease networks

It is becoming evident that holistic perspectives toward cancer are crucial in deciphering the overwhelming complexity of tumors. Single-layer analysis of genome-wide data has greatly contributed to our understanding of cellular systems and their perturbations. However, fundamental gaps in our knowledge persist and hamper the design of effective interventions. It is becoming more apparent than ever, that cancer should not only be viewed as a disease of the genome but as a disease of the cellular system. Integrative multilayer approaches are emerging as vigorous assets in our endeavors to achieve systemic views on cancer biology. Herein, we provide a comprehensive review of the approaches, methods and technologies that can serve to achieve systemic perspectives of cancer. We start with genome-wide single-layer approaches of omics analyses of cellular systems and move on to multilayer integrative approaches in which in-depth descriptions of proteogenomics and network-based data analysis are provided. Proteogenomics is a remarkable example of how the integration of multiple levels of information can reduce our blind spots and increase the accuracy and reliability of our interpretations and network-based data analysis is a major approach for data interpretation and a robust scaffold for data integration and modeling. Overall, this review aims to increase cross-field awareness of the approaches and challenges regarding the omics-based study of cancer and to facilitate the necessary shift toward holistic approaches.

Deciphering Human Leukocyte Antigen Susceptibility Maps From Immunopeptidomics Characterization in Oncology and Infections

Genetic variability across the three major histocompatibility complex (MHC) class I genes (human leukocyte antigen [HLA] A, B, and C) may affect susceptibility to many diseases such as cancer, auto-immune or infectious diseases. Individual genetic variation may help to explain different immune responses to microorganisms across a population. HLA typing can be fast and inexpensive; however, deciphering peptides loaded on MHC-I and II which are presented to T cells, require the design and development of high-sensitivity methodological approaches and subsequently databases. Hence, these novel strategies and databases could help in the generation of vaccines using these potential immunogenic peptides and in identifying high-risk HLA types to be prioritized for vaccination programs. Herein, the recent developments and approaches, in this field, focusing on the identification of immunogenic peptides have been reviewed and the next steps to promote their translation into biomedical and clinical practice are discussed.

Main Strategies for the Identification of Neoantigens

Genetic instability of tumors leads to the appearance of numerous tumor-specific somatic mutations that could potentially result in the production of mutated peptides that are presented on the cell surface by the MHC molecules. Peptides of this kind are commonly called neoantigens. Their presence on the cell surface specifically distinguishes tumors from healthy tissues. This feature makes neoantigens a promising target for immunotherapy. The rapid evolution of high-throughput genomics and proteomics makes it possible to implement these techniques in clinical practice. In particular, they provide useful tools for the investigation of neoantigens. The most valuable genomic approach to this problem is whole-exome sequencing coupled with RNA-seq. High-throughput mass-spectrometry is another option for direct identification of MHC-bound peptides, which is capable of revealing the entire MHC-bound peptidome. Finally, structure-based predictions could significantly improve the understanding of physicochemical and structural features that affect the immunogenicity of peptides. The development of pipelines combining such tools could improve the accuracy of the peptide selection process and decrease the required time. Here we present a review of the main existing approaches to investigating the neoantigens and suggest a possible ideal pipeline that takes into account all modern trends in the context of neoantigen discovery.

The International Conference on Intelligent Biology and Medicine 2019 (ICIBM 2019): computational methods and applications in medical genomics

In this editorial, we briefly summarized the International Conference on Intelligent Biology and Medicine 2019 (ICIBM 2019) that was held on June 9-11, 2019 at Columbus, Ohio, USA. We further introduced the 19 research articles included in this supplement issue, covering four major areas, namely computational method development, genomics analysis, network-based analysis and biomarker prediction. The selected papers perform cutting edge computational research applied to a broad range of human diseases such as cancer, neural degenerative and chronic inflammatory disease. They also proposed solutions for fundamental medical genomics problems range from basic data processing and quality control to functional interpretation, biomarker and drug prediction, and database releasing.