Personalized RNA mutanome vaccines mobilize poly-specific therapeutic immunity against cancer.

T cells directed against mutant neo-epitopes drive cancer immunity. However, spontaneous immune recognition of mutations is inefficient. We recently introduced the concept of individualized mutanome vaccines and implemented an RNA-based poly-neo-epitope approach to mobilize immunity against a spectrum of cancer mutations. Here we report the first-in-human application of this concept in melanoma. We set up a process comprising comprehensive identification of individual mutations, computational prediction of neo-epitopes, and design and manufacturing of a vaccine unique for each patient. All patients developed T cell responses against multiple vaccine neo-epitopes at up to high single-digit percentages. Vaccine-induced T cell infiltration and neo-epitope-specific killing of autologous tumour cells were shown in post-vaccination resected metastases from two patients. The cumulative rate of metastatic events was highly significantly reduced after the start of vaccination, resulting in a sustained progression-free survival. Two of the five patients with metastatic disease experienced vaccine-related objective responses. One of these patients had a late relapse owing to outgrowth of ß2-microglobulin-deficient melanoma cells as an acquired resistance mechanism. A third patient developed a complete response to vaccination in combination with PD-1 blockade therapy. Our study demonstrates that individual mutations can be exploited, thereby opening a path to personalized immunotherapy for patients with cancer.

ArtiFuse-computational validation of fusion gene detection tools without relying on simulated reads.

MOTIVATION: Gene fusions are an important class of transcriptional variants that can influence cancer development and can be predicted from RNA sequencing (RNA-seq) data by multiple existing tools. However, the real-world performance of these tools is unclear due to the lack of known positive and negative events, especially with regard to fusion genes in individual samples. Often simulated reads are used, but these cannot account for all technical biases in RNA-seq data generated from real samples. RESULTS: Here, we present ArtiFuse, a novel approach that simulates fusion genes by sequence modification to the genomic reference, and therefore, can be applied to any RNA-seq dataset without the need for any simulated reads. We demonstrate our approach on eight RNA-seq datasets for three fusion gene prediction tools: average recall values peak for all three tools between 0.4 and 0.56 for high-quality and high-coverage datasets. As ArtiFuse affords total control over involved genes and breakpoint position, we also assessed performance with regard to gene-related properties, showing a drop-in recall value for low-expressed genes in high-coverage samples and genes with co-expressed paralogues. Overall tool performance assessed from ArtiFusions is lower compared to previously reported estimates on simulated reads. Due to the use of real RNA-seq datasets, we believe that ArtiFuse provides a more realistic benchmark that can be used to develop more accurate fusion gene prediction tools for application in clinical settings. AVAILABILITY AND IMPLEMENTATION: ArtiFuse is implemented in Python. The source code and documentation are available at https://github.com/TRON-Bioinformatics/ArtiFusion. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Immunomic, genomic and transcriptomic characterization of CT26 colorectal carcinoma.

BACKGROUND: Tumor models are critical for our understanding of cancer and the development of cancer therapeutics. Here, we present an integrated map of the genome, transcriptome and immunome of an epithelial mouse tumor, the CT26 colon carcinoma cell line. RESULTS: We found that Kras is homozygously mutated at p.G12D, Apc and Tp53 are not mutated, and Cdkn2a is homozygously deleted. Proliferation and stem-cell markers, including Top2a, Birc5 (Survivin), Cldn6 and Mki67, are highly expressed while differentiation and top-crypt markers Muc2, Ms4a8a (MS4A8B) and Epcam are not. Myc, Trp53 (tp53), Mdm2, Hif1a, and Nras are highly expressed while Egfr and Flt1 are not. MHC class I but not MHC class II is expressed. Several known cancer-testis antigens are expressed, including Atad2, Cep55, and Pbk. The highest expressed gene is a mutated form of the mouse tumor antigen gp70. Of the 1,688 non-synonymous point variations, 154 are both in expressed genes and in peptides predicted to bind MHC and thus potential targets for immunotherapy development. Based on its molecular signature, we predicted that CT26 is refractory to anti-EGFR mAbs and sensitive to MEK and MET inhibitors, as have been previously reported. CONCLUSIONS: CT26 cells share molecular features with aggressive, undifferentiated, refractory human colorectal carcinoma cells. As CT26 is one of the most extensively used syngeneic mouse tumor models, our data provide a map for the rationale design of mode-of-action studies for pre-clinical evaluation of targeted- and immunotherapies.

Galaxy LIMS for next-generation sequencing.

SUMMARY: We have developed a laboratory information management system (LIMS) for a next-generation sequencing (NGS) laboratory within the existing Galaxy platform. The system provides lab technicians standard and customizable sample information forms, barcoded submission forms, tracking of input sample quality, multiplex-capable automatic flow cell design and automatically generated sample sheets to aid physical flow cell preparation. In addition, the platform provides the researcher with a user-friendly interface to create a request, submit accompanying samples, upload sample quality measurements and access to the sequencing results. As the LIMS is within the Galaxy platform, the researcher has access to all Galaxy analysis tools and workflows. The system reports requests and associated information to a message queuing system, such that information can be posted and stored in external systems, such as a wiki. Through an API, raw sequencing results can be automatically pre-processed and uploaded to the appropriate request folder. Developed for the Illumina HiSeq 2000 instrument, many features are directly applicable to other instruments. AVAILABILITY AND IMPLEMENTATION: The code and documentation are available at http://tron-mainz.de/tron-facilities/computational-medicine/galaxy-lims/

Prediction of tumor-specific splicing from somatic mutations as a source of neoantigen candidates.

Motivation: Neoantigens are promising targets for cancer immunotherapies and might arise from alternative splicing. However, detecting tumor-specific splicing is challenging because many non-canonical splice junctions identified in tumors also appear in healthy tissues. To increase tumor-specificity, we focused on splicing caused by somatic mutations as a source for neoantigen candidates in individual patients. Results: We developed the tool splice2neo with multiple functionalities to integrate predicted splice effects from somatic mutations with splice junctions detected in tumor RNA-seq and to annotate the resulting transcript and peptide sequences. Additionally, we provide the tool EasyQuant for targeted RNA-seq read mapping to candidate splice junctions. Using a stringent detection rule, we predicted 1.7 splice junctions per patient as splice targets with a false discovery rate below 5% in a melanoma cohort. We confirmed tumor-specificity using independent, healthy tissue samples. Furthermore, using tumor-derived RNA, we confirmed individual exon-skipping events experimentally. Most target splice junctions encoded neoepitope candidates with predicted major histocompatibility complex (MHC)-I or MHC-II binding. Compared to neoepitope candidates from non-synonymous point mutations, the splicing-derived MHC-I neoepitope candidates had lower self-similarity to corresponding wild-type peptides. In conclusion, we demonstrate that identifying mutation-derived, tumor-specific splice junctions can lead to additional neoantigen candidates to expand the target repertoire for cancer immunotherapies. Availability and implementation: The R package splice2neo and the python package EasyQuant are available at https://github.com/TRON-Bioinformatics/splice2neo and https://github.com/TRON-Bioinformatics/easyquant, respectively.

Multiple instance learning to predict immune checkpoint blockade efficacy using neoantigen candidates.

Previous studies showed that the neoantigen candidate load is an imperfect predictor of immune checkpoint blockade (ICB) efficacy. Further studies provided evidence that the response to ICB is also affected by the qualitative properties of a few or even single candidates, limiting the predictive power based on candidate quantity alone. Here, we predict ICB efficacy based on neoantigen candidates and their neoantigen features in the context of the mutation type, using Multiple-Instance Learning via Embedded Instance Selection (MILES). Multiple instance learning is a type of supervised machine learning that classifies labeled bags that are formed by a set of unlabeled instances. MILES performed better compared with neoantigen candidate load alone for low-abundant fusion genes in renal cell carcinoma. Our findings suggest that MILES is an appropriate method to predict the efficacy of ICB therapy based on neoantigen candidates without requiring direct T cell response information.

CoVigator-A Knowledge Base for Navigating SARS-CoV-2 Genomic Variants.

BACKGROUND: The outbreak of the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) resulted in the global COVID-19 pandemic. The urgency for an effective SARS-CoV-2 vaccine has led to the development of the first series of vaccines at unprecedented speed. The discovery of SARS-CoV-2 spike-glycoprotein mutants, however, and consequentially the potential to escape vaccine-induced protection and increased infectivity, demonstrates the persisting importance of monitoring SARS-CoV-2 mutations to enable early detection and tracking of genomic variants of concern. RESULTS: We developed the CoVigator tool with three components: (1) a knowledge base that collects new SARS-CoV-2 genomic data, processes it and stores its results; (2) a comprehensive variant calling pipeline; (3) an interactive dashboard highlighting the most relevant findings. The knowledge base routinely downloads and processes virus genome assemblies or raw sequencing data from the COVID-19 Data Portal (C19DP) and the European Nucleotide Archive (ENA), respectively. The results of variant calling are visualized through the dashboard in the form of tables and customizable graphs, making it a versatile tool for tracking SARS-CoV-2 variants. We put a special emphasis on the identification of intrahost mutations and make available to the community what is, to the best of our knowledge, the largest dataset on SARS-CoV-2 intrahost mutations. In the spirit of open data, all CoVigator results are available for download. The CoVigator dashboard is accessible via covigator.tron-mainz.de. CONCLUSIONS: With increasing demand worldwide in genome surveillance for tracking the spread of SARS-CoV-2, CoVigator will be a valuable resource of an up-to-date list of mutations, which can be incorporated into global efforts.

Accurate detection of tumor-specific gene fusions reveals strongly immunogenic personal neo-antigens.

Cancer-associated gene fusions are a potential source for highly immunogenic neoantigens, but the lack of computational tools for accurate, sensitive identification of personal gene fusions has limited their targeting in personalized cancer immunotherapy. Here we present EasyFuse, a machine learning computational pipeline for detecting cancer-specific gene fusions in transcriptome data obtained from human cancer samples. EasyFuse predicts personal gene fusions with high precision and sensitivity, outperforming previously described tools. By testing immunogenicity with autologous blood lymphocytes from patients with cancer, we detected pre-established CD4+ and CD8+ T cell responses for 10 of 21 (48%) and for 1 of 30 (3%) identified gene fusions, respectively. The high frequency of T cell responses detected in patients with cancer supports the relevance of individual gene fusions as neoantigens that might be targeted in personalized immunotherapies, especially for tumors with low mutation burden.

Large-scale analysis of SARS-CoV-2 spike-glycoprotein mutants demonstrates the need for continuous screening of virus isolates.

Due to the widespread of the COVID-19 pandemic, the SARS-CoV-2 genome is evolving in diverse human populations. Several studies already reported different strains and an increase in the mutation rate. Particularly, mutations in SARS-CoV-2 spike-glycoprotein are of great interest as it mediates infection in human and recently approved mRNA vaccines are designed to induce immune responses against it. We analyzed 1,036,030 SARS-CoV-2 genome assemblies and 30,806 NGS datasets from GISAID and European Nucleotide Archive (ENA) focusing on non-synonymous mutations in the spike protein. Only around 2.5% of the samples contained the wild-type spike protein with no variation from the reference. Among the spike protein mutants, we confirmed a low mutation rate exhibiting less than 10 non-synonymous mutations in 99.6% of the analyzed sequences, but the mean and median number of spike protein mutations per sample increased over time. 5,472 distinct variants were found in total. The majority of the observed variants were recurrent, but only 21 and 14 recurrent variants were found in at least 1% of the mutant genome assemblies and NGS samples, respectively. Further, we found high-confidence subclonal variants in about 2.6% of the NGS data sets with mutant spike protein, which might indicate co-infection with various SARS-CoV-2 strains and/or intra-host evolution. Lastly, some variants might have an effect on antibody binding or T-cell recognition. These findings demonstrate the continuous importance of monitoring SARS-CoV-2 sequences for an early detection of variants that require adaptations in preventive and therapeutic strategies.

OXPHOS remodeling in high-grade prostate cancer involves mtDNA mutations and increased succinate oxidation.

Rewiring of energy metabolism and adaptation of mitochondria are considered to impact on prostate cancer development and progression. Here, we report on mitochondrial respiration, DNA mutations and gene expression in paired benign/malignant human prostate tissue samples. Results reveal reduced respiratory capacities with NADH-pathway substrates glutamate and malate in malignant tissue and a significant metabolic shift towards higher succinate oxidation, particularly in high-grade tumors. The load of potentially deleterious mitochondrial-DNA mutations is higher in tumors and associated with unfavorable risk factors. High levels of potentially deleterious mutations in mitochondrial Complex I-encoding genes are associated with a 70% reduction in NADH-pathway capacity and compensation by increased succinate-pathway capacity. Structural analyses of these mutations reveal amino acid alterations leading to potentially deleterious effects on Complex I, supporting a causal relationship. A metagene signature extracted from the transcriptome of tumor samples exhibiting a severe mitochondrial phenotype enables identification of tumors with shorter survival times.

HLA and proteasome expression body map.

BACKGROUND: The presentation of HLA peptide complexes to T cells is a highly regulated and tissue specific process involving multiple transcriptionally controlled cellular components. The extensive polymorphism of HLA genes and the complex composition of the proteasome make it difficult to map their expression profiles across tissues. METHODS: Here we applied a tailored gene quantification pipeline to 4323 publicly available RNA-Seq datasets representing 55 normal tissues and cell types to examine expression profiles of (classical and non-classical) HLA class I, class II and proteasomal genes. RESULTS: We generated the first comprehensive expression atlas of antigen presenting-related genes across 56 normal tissues and cell types, including immune cells, pancreatic islets, platelets and hematopoietic stem cells. We found a surprisingly heterogeneous HLA expression pattern with up to 100-fold difference in intra-tissue median HLA abundances. Cells of the immune system and lymphatic organs expressed the highest levels of classical HLA class I (HLA-A,-B,-C), class II (HLA-DQA1,-DQB1,-DPA1,-DPB1,-DRA,-DRB1) and non-classical HLA class I (HLA-E,-F) molecules, whereas retina, brain, muscle, megakaryocytes and erythroblasts showed the lowest abundance. In contrast, we identified a distinct and highly tissue-restricted expression pattern of the non-classical class I gene HLA-G in placenta, pancreatic islets, pituitary gland and testis. While the constitutive proteasome showed relatively constant expression across all tissues, we found the immunoproteasome to be enriched in lymphatic organs and almost absent in immune privileged tissues. CONCLUSIONS: Here, we not only provide a reference catalog of tissue and cell type specific HLA expression, but also highlight extremely variable expression of the basic components of antigen processing and presentation in different cell types. Our findings indicate that low expression of classical HLA class I molecules together with lack of immunoproteasome components as well as upregulation of HLA-G may be of key relevance to maintain tolerance in immune privileged tissues.

TCLP: an online cancer cell line catalogue integrating HLA type, predicted neo-epitopes, virus and gene expression.

Human cancer cell lines are an important resource for research and drug development. However, the available annotations of cell lines are sparse, incomplete, and distributed in multiple repositories. Re-analyzing publicly available raw RNA-Seq data, we determined the human leukocyte antigen (HLA) type and abundance, identified expressed viruses and calculated gene expression of 1,082 cancer cell lines. Using the determined HLA types, public databases of cell line mutations, and existing HLA binding prediction algorithms, we predicted antigenic mutations in each cell line. We integrated the results into a comprehensive knowledgebase. Using the Django web framework, we provide an interactive user interface with advanced search capabilities to find and explore cell lines and an application programming interface to extract cell line information. The portal is available at http://celllines.tron-mainz.de.