Neodb: a comprehensive neoantigen database and discovery platform for cancer immunotherapy.

Neoantigens derived from somatic deoxyribonucleic acid alterations are ideal cancer-specific targets. However, integrated platform for neoantigen discovery is urgently needed. Recently, many scattered experimental evidences suggest that some neoantigens are immunogenic, and comprehensive collection of these experimentally validated neoantigens is still lacking. Here, we have integrated the commonly used tools in the current neoantigen discovery process to form a comprehensive web-based analysis platform. To identify experimental evidences supporting the immunogenicity of neoantigens, we performed comprehensive literature search and constructed the database. The collection of public neoantigens was obtained by using comprehensive features to filter the potential neoantigens from recurrent driver mutations. Importantly, we constructed a graph neural network (GNN) model (Immuno-GNN) using an attention mechanism to consider the spatial interactions between human leukocyte antigen and antigenic peptides for neoantigen immunogenicity prediction. The new easy-to-use R/Shiny web-based neoantigen database and discovery platform, Neodb, contains currently the largest number of experimentally validated neoantigens. In addition to validated neoantigen, Neodb also includes three additional modules for facilitating neoantigen prediction and analysis, including 'Tools' module (comprehensive neoantigen prediction tools); 'Driver-Neo' module (collection of public neoantigens derived from recurrent mutations) and 'Immuno-GNN' module (a novel immunogenicity prediction tool based on a GNN). Immuno-GNN shows improved performance compared with known methods and also represents the first application of GNN model in neoantigen immunogenicity prediction. The construction of Neodb will facilitate the study of neoantigen immunogenicity and the clinical application of neoantigen-based cancer immunotherapy. Database URL https://liuxslab.com/Neodb/.

Copy number alteration features in pan-cancer homologous recombination deficiency prediction and biology.

Homologous recombination deficiency (HRD) renders cancer cells vulnerable to unrepaired double-strand breaks and is an important therapeutic target as exemplified by the clinical efficacy of poly ADP-ribose polymerase (PARP) inhibitors as well as the platinum chemotherapy drugs applied to HRD patients. However, it remains a challenge to predict HRD status precisely and economically. Copy number alteration (CNA), as a pervasive trait of human cancers, can be extracted from a variety of data sources, including whole genome sequencing (WGS), SNP array, and panel sequencing, and thus can be easily applied clinically. Here we systematically evaluate the predictive performance of various CNA features and signatures in HRD prediction and build a gradient boosting machine model (HRDCNA) for pan-cancer HRD prediction based on these CNA features. CNA features BP10MB[1] (The number of breakpoints per 10MB of DNA is 1) and SS[ > 7 & <=8] (The log10-based size of segments is greater than 7 and less than or equal to 8) are identified as the most important features in HRD prediction. HRDCNA suggests the biallelic inactivation of BRCA1, BRCA2, PALB2, RAD51C, RAD51D, and BARD1 as the major genetic basis for human HRD, and may also be applied to effectively validate the pathogenicity of BRCA1/2 variants of uncertain significance (VUS). Together, this study provides a robust tool for cost-effective HRD prediction and also demonstrates the applicability of CNA features and signatures in cancer precision medicine.

Accurate prediction of pan-cancer types using machine learning with minimal number of DNA methylation sites.

DNA methylation analysis has been applied to determine the primary site of cancer; however, robust and accurate prediction of cancer types with a minimum number of sites is still a significant scientific challenge. To build an accurate and robust cancer type prediction tool with a minimum number of DNA methylation sites, we internally benchmarked different DNA methylation site selection and ranking procedures, as well as different classification models. We used The Cancer Genome Atlas dataset (26 cancer types with 8296 samples) to train and test models and used an independent dataset (17 cancer types with 2738 samples) for model validation. A deep neural network model using a combined feature selection procedure (named MethyDeep) can predict 26 cancer types using 30 methylation sites with superior performance compared with the known methods for both primary and metastatic cancers in independent validation datasets. In conclusion, MethyDeep is an accurate and robust cancer type predictor with the minimum number of DNA methylation sites; it could help the cost-effective clarification of cancer of unknown primary patients and the liquid biopsy-based early screening of cancers.

TLimmuno2: predicting MHC class II antigen immunogenicity through transfer learning.

Major histocompatibility complex (MHC) class II molecules play a pivotal role in antigen presentation and CD4+ T cell response. Accurate prediction of the immunogenicity of MHC class II-associated antigens is critical for vaccine design and cancer immunotherapies. However, current computational methods are limited by insufficient training data and algorithmic constraints, and the rules that govern which peptides are truly recognized by existing T cell receptors remain poorly understood. Here, we build a transfer learning-based, long short-term memory model named 'TLimmuno2' to predict whether epitope-MHC class II complex can elicit T cell response. Through leveraging binding affinity data, TLimmuno2 shows superior performance compared with existing models on independent validation datasets. TLimmuno2 can find real immunogenic neoantigen in real-world cancer immunotherapy data. The identification of significant MHC class II neoantigen-mediated immunoediting signal in the cancer genome atlas pan-cancer dataset further suggests the robustness of TLimmuno2 in identifying really immunogenic neoantigens that are undergoing negative selection during cancer evolution. Overall, TLimmuno2 is a powerful tool for the immunogenicity prediction of MHC class II presented epitopes and could promote the development of personalized immunotherapies.

The repertoire of copy number alteration signatures in human cancer.

Copy number alterations (CNAs) are a predominant source of genetic alterations in human cancer and play an important role in cancer progression. However comprehensive understanding of the mutational processes and signatures of CNA is still lacking. Here we developed a mechanism-agnostic method to categorize CNA based on various fragment properties, which reflect the consequences of mutagenic processes and can be extracted from different types of data, including whole genome sequencing (WGS) and single nucleotide polymorphism (SNP) array. The 14 signatures of CNA have been extracted from 2778 pan-cancer analysis of whole genomes WGS samples, and further validated with 10 851 the cancer genome atlas SNP array dataset. Novel patterns of CNA have been revealed through this study. The activities of some CNA signatures consistently predict cancer patients' prognosis. This study provides a repertoire for understanding the signatures of CNA in cancer, with potential implications for cancer prognosis, evolution and etiology.

Seq2Neo: A Comprehensive Pipeline for Cancer Neoantigen Immunogenicity Prediction.

Neoantigens derived from somatic DNA alterations are ideal cancer-specific targets. In recent years, the combination therapy of PD-1/PD-L1 blockers and neoantigen vaccines has shown clinical efficacy in original PD-1/PD-L1 blocker non-responders. However, not all somatic DNA mutations result in immunogenicity among cancer cells and efficient tools to predict the immunogenicity of neoepitopes are still urgently needed. Here, we present the Seq2Neo pipeline, which provides a one-stop solution for neoepitope feature prediction using raw sequencing data. Neoantigens derived from different types of genome DNA alterations, including point mutations, insertion deletions and gene fusions, are all supported. Importantly, a convolutional neural network (CNN)-based model was trained to predict the immunogenicity of neoepitopes and this model showed an improved performance compared to the currently available tools in immunogenicity prediction using independent datasets. We anticipate that the Seq2Neo pipeline could become a useful tool in the prediction of neoantigen immunogenicity and cancer immunotherapy. Seq2Neo is open-source software under an academic free license (AFL) v3.0 and is freely available at Github.

Oncogenic EFNA4 Amplification Promotes Lung Adenocarcinoma Lymph Node Metastasis.

Lymph nodes metastases are common in patients with lung cancer. Additionally, those patients are often at a higher risk for death from lung tumor than those with tumor-free lymph nodes. Somatic DNA alterations are key drivers of cancer, and copy number alterations (CNAs) are major types of DNA alteration that promote lung cancer progression. Here, we performed genome-wide DNA copy number analysis, and identified a novel lung-cancer-metastasis-related gene, EFNA4. The EFNA4 genome locus was significantly amplified, and EFNA4 mRNA expression was significantly up-regulated in lung cancer compared with normal lung tissue, and also in lung cancer with lymph node metastases compared with lung cancer without metastasis. EFNA4 encodes Ephrin A4, which is the ligand for Eph receptors. The function of EFNA4 in human lung cancer remains largely unknown. Through cell line experiments we showed that EFNA4 overexpression contributes to lung tumor cells growth, migration and adhesion. Conversely, EFNA4 knockdown or knockout led to the growth suppression of cells and tumor xenografts in mice. Lung cancer patients with EFNA4 overexpression have poor prognosis. Together, by elucidating a new layer of the role of EFNA4 in tumor proliferation and migration, our study demonstrates a better understanding of the function of the significantly amplified and overexpressed gene EFNA4 in lung tumor metastasis, and suggests EFNA4 as a potential target in metastatic lung cancer therapy.