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.

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.

MC38 colorectal tumor cell lines from two different sources display substantial differences in transcriptome, mutanome and neoantigen expression.

Introduction: The cell line MC38 is a commonly used murine model for colorectal carcinoma. It has a high mutational burden, is sensitive to immune checkpoint immunotherapy and endogenous CD8+ T cell responses against neoantigens have been reported. Methods: Here, we re-sequenced exomes and transcriptomes of MC38 cells from two different sources, namely Kerafast (originating from NCI/NIH, MC38-K) and the Leiden University Medical Center cell line collection (MC38-L), comparing the cell lines on the genomic and transcriptomic level and analyzing their recognition by CD8+ T cells with known neo-epitope specificity. Results: The data reveals a distinct structural composition of MC38-K and MC38-L cell line genomes and different ploidies. Further, the MC38-L cell line harbored about 1.3-fold more single nucleotide variations and small insertions and deletions than the MC38-K cell line. In addition, the observed mutational signatures differed; only 35.3% of the non-synonymous variants and 5.4% of the fusion gene events were shared. Transcript expression values of both cell lines correlated strongly (p = 0.919), but we found different pathways enriched in the genes that were differentially upregulated in the MC38-L or MC38-K cells, respectively. Our data show that previously described neoantigens in the MC38 model such as Rpl18mut and Adpgkmut were absent in the MC38-K cell line resulting that such neoantigen-specific CD8+ T cells recognizing and killing MC38-L cells did not recognize or kill MC38-K cells. Conclusion: This strongly indicates that at least two sub-cell lines of MC38 exist in the field and underlines the importance of meticulous tracking of investigated cell lines to obtain reproducible results, and for correct interpretation of the immunological data without artifacts. We present our analyses as a reference for researchers to select the appropriate sub-cell line for their own studies.

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.

Bioinformatics for Cancer Immunotherapy.

Our immune system plays a key role in health and disease as it is capable of responding to foreign antigens as well as acquired antigens from cancer cells. Latter are caused by somatic mutations, the so-called neoepitopes, and might be recognized by T cells if they are presented by HLA molecules on the surface of cancer cells. Personalized mutanome vaccines are a class of customized immunotherapies, which is dependent on the detection of individual cancer-specific tumor mutations and neoepitope (i.e., prediction, followed by a rational vaccine design, before on-demand production. The development of next generation sequencing (NGS) technologies and bioinformatic tools allows a large-scale analysis of each parameter involved in this process. Here, we provide an overview of the bioinformatic aspects involved in the design of personalized, neoantigen-based vaccines, including the detection of mutations and the subsequent prediction of potential epitopes, as well as methods for associated biomarker research, such as high-throughput sequencing of T-cell receptors (TCRs), followed by data analysis and the bioinformatics quantification of immune cell infiltration in cancer samples.

Multi-Omics Characterization of the 4T1 Murine Mammary Gland Tumor Model.

Background: Tumor models are critical for our understanding of cancer and the development of cancer therapeutics. The 4T1 murine mammary cancer cell line is one of the most widely used breast cancer models. Here, we present an integrated map of the genome, transcriptome, and immunome of 4T1. Results: We found Trp53 (Tp53) and Pik3g to be mutated. Other frequently mutated genes in breast cancer, including Brca1 and Brca2, are not mutated. For cancer related genes, Nav3, Cenpf, Muc5Ac, Mpp7, Gas1, MageD2, Dusp1, Ros, Polr2a, Rragd, Ros1, and Hoxa9 are mutated. Markers for cell proliferation like Top2a, Birc5, and Mki67 are highly expressed, so are markers for metastasis like Msln, Ect2, and Plk1, which are known to be overexpressed in triple-negative breast cancer (TNBC). TNBC markers are, compared to a mammary gland control sample, lower (Esr1), comparably low (Erbb2), or not expressed at all (Pgr). We also found testis cancer antigen Pbk as well as colon/gastrointestinal cancer antigens Gpa33 and Epcam to be highly expressed. Major histocompatibility complex (MHC) class I is expressed, while MHC class II is not. We identified 505 single nucleotide variations (SNVs) and 20 insertions and deletions (indels). Neoantigens derived from 22 SNVs and one deletion elicited CD8+ or CD4+ T cell responses in IFNγ-ELISpot assays. Twelve high-confidence fusion genes were observed. We did not observe significant downregulation of mismatch repair (MMR) genes or SNVs/indels impairing their function, providing evidence for 6-thioguanine resistance. Effects of the integration of the murine mammary tumor virus were observed at the genome and transcriptome level. Conclusions: 4T1 cells share substantial molecular features with human TNBC. As 4T1 is a common model for metastatic tumors, our data supports the rational design of mode-of-action studies for pre-clinical evaluation of targeted immunotherapies.