Multi-omic and functional analysis for classification and treatment of sarcomas with FUS-TFCP2 or EWSR1-TFCP2 fusions.

Linking clinical multi-omics with mechanistic studies may improve the understanding of rare cancers. We leverage two precision oncology programs to investigate rhabdomyosarcoma with FUS/EWSR1-TFCP2 fusions, an orphan malignancy without effective therapies. All tumors exhibit outlier ALK expression, partly accompanied by intragenic deletions and aberrant splicing resulting in ALK variants that are oncogenic and sensitive to ALK inhibitors. Additionally, recurrent CKDN2A/MTAP co-deletions provide a rationale for PRMT5-targeted therapies. Functional studies show that FUS-TFCP2 blocks myogenic differentiation, induces transcription of ALK and truncated TERT, and inhibits DNA repair. Unlike other fusion-driven sarcomas, TFCP2-rearranged tumors exhibit genomic instability and signs of defective homologous recombination. DNA methylation profiling demonstrates a close relationship with undifferentiated sarcomas. In two patients, sarcoma was preceded by benign lesions carrying FUS-TFCP2, indicating stepwise sarcomagenesis. This study illustrates the potential of linking precision oncology with preclinical research to gain insight into the classification, pathogenesis, and therapeutic vulnerabilities of rare cancers.

NCT/DKFZ MASTER handbook of interpreting whole-genome, transcriptome, and methylome data for precision oncology.

Analysis of selected cancer genes has become an important tool in precision oncology but cannot fully capture the molecular features and, most importantly, vulnerabilities of individual tumors. Observational and interventional studies have shown that decision-making based on comprehensive molecular characterization adds significant clinical value. However, the complexity and heterogeneity of the resulting data are major challenges for disciplines involved in interpretation and recommendations for individualized care, and limited information exists on how to approach multilayered tumor profiles in clinical routine. We report our experience with the practical use of data from whole-genome or exome and RNA sequencing and DNA methylation profiling within the MASTER (Molecularly Aided Stratification for Tumor Eradication Research) program of the National Center for Tumor Diseases (NCT) Heidelberg and Dresden and the German Cancer Research Center (DKFZ). We cover all relevant steps of an end-to-end precision oncology workflow, from sample collection, molecular analysis, and variant prioritization to assigning treatment recommendations and discussion in the molecular tumor board. To provide insight into our approach to multidimensional tumor profiles and guidance on interpreting their biological impact and diagnostic and therapeutic implications, we present case studies from the NCT/DKFZ molecular tumor board that illustrate our daily practice. This manual is intended to be useful for physicians, biologists, and bioinformaticians involved in the clinical interpretation of genome-wide molecular information.

Signaling-induced systematic repression of miRNAs uncovers cancer vulnerabilities and targeted therapy sensitivity.

Targeted therapies are effective in treating cancer, but success depends on identifying cancer vulnerabilities. In our study, we utilize small RNA sequencing to examine the impact of pathway activation on microRNA (miRNA) expression patterns. Interestingly, we discover that miRNAs capable of inhibiting key members of activated pathways are frequently diminished. Building on this observation, we develop an approach that integrates a low-miRNA-expression signature to identify druggable target genes in cancer. We train and validate our approach in colorectal cancer cells and extend it to diverse cancer models using patient-derived in vitro and in vivo systems. Finally, we demonstrate its additional value to support genomic and transcriptomic-based drug prediction strategies in a pan-cancer patient cohort from the National Center for Tumor Diseases (NCT)/German Cancer Consortium (DKTK) Molecularly Aided Stratification for Tumor Eradication (MASTER) precision oncology trial. In conclusion, our strategy can predict cancer vulnerabilities with high sensitivity and accuracy and might be suitable for future therapy recommendations in a variety of cancer subtypes.

Proteogenomic analysis reveals RNA as a source for tumor-agnostic neoantigen identification.

Systemic pan-tumor analyses may reveal the significance of common features implicated in cancer immunogenicity and patient survival. Here, we provide a comprehensive multi-omics data set for 32 patients across 25 tumor types for proteogenomic-based discovery of neoantigens. By using an optimized computational approach, we discover a large number of tumor-specific and tumor-associated antigens. To create a pipeline for the identification of neoantigens in our cohort, we combine DNA and RNA sequencing with MS-based immunopeptidomics of tumor specimens, followed by the assessment of their immunogenicity and an in-depth validation process. We detect a broad variety of non-canonical HLA-binding peptides in the majority of patients demonstrating partially immunogenicity. Our validation process allows for the selection of 32 potential neoantigen candidates. The majority of neoantigen candidates originates from variants identified in the RNA data set, illustrating the relevance of RNA as a still understudied source of cancer antigens. This study underlines the importance of RNA-centered variant detection for the identification of shared biomarkers and potentially relevant neoantigen candidates.

aRgus: Multilevel visualization of non-synonymous single nucleotide variants & advanced pathogenicity score modeling for genetic vulnerability assessment.

The widespread use of high-throughput sequencing techniques is leading to a rapidly increasing number of disease-associated variants of unknown significance and candidate genes. Integration of knowledge concerning their genetic, protein as well as functional and conservational aspects is necessary for an exhaustive assessment of their relevance and for prioritization of further clinical and functional studies investigating their role in human disease. To collect the necessary information, a multitude of different databases has to be accessed and data extraction from the original sources commonly is not user-friendly and requires advanced bioinformatics skills. This leads to a decreased data accessibility for a relevant number of potential users such as clinicians, geneticist, and clinical researchers. Here, we present aRgus (https://argus.urz.uni-heidelberg.de/), a standalone webtool for simple extraction and intuitive visualization of multi-layered gene, protein, variant, and variant effect prediction data. aRgus provides interactive exploitation of these data within seconds for any known gene of the human genome. In contrast to existing online platforms for compilation of variant data, aRgus complements visualization of chromosomal exon-intron structure and protein domain annotation with ClinVar and gnomAD variant distributions as well as position-specific variant effect prediction score modeling. aRgus thereby enables timely assessment of protein regions vulnerable to variation with single amino acid resolution and provides numerous applications in variant and protein domain interpretation as well as in the design of in vitro experiments.

IgE type multiple myeloma exhibits hypermutated phenotype and tumor reactive T cells.

Multiple myeloma (MM) is a hematological malignancy originating from malignant and clonally expanding plasma cells. MM can be molecularly stratified, and its clonal evolution deciphered based on the Ig heavy and light chains of the respective malignant plasma cell clone. Of all MM subtypes, IgE type MM accounts for only <0.1% of cases and is associated with an aggressive clinical course and consequentially dismal prognosis. In such malignancies, adoptive transfer of autologous lymphocytes specifically targeting presented (neo)epitopes encoded by either somatically mutated or specifically overexpressed genes has resulted in substantial objective clinical regressions even in relapsed/refractory disease. However, there are no data on the genetic and immunological characteristics of this rare and aggressive entity. Here, we comprehensively profiled IgE type kappa MM on a genomic and immune repertoire level by integrating DNA- and single-cell RNA sequencing and comparative profiling against non-IgE type MM samples. We demonstrate distinct pathophysiological mechanisms as well as novel opportunities for targeting IgE type MM. Our data further provides the rationale for patient-individualized neoepitope-targeting cell therapy in high tumor mutation burden MM.

Bcl-xL as prognostic marker and potential therapeutic target in cholangiocarcinoma.

Intrahepatic, perihilar, and distal cholangiocarcinoma (iCCA, pCCA, dCCA) are highly malignant tumours with increasing mortality rates due to therapy resistances. Among the mechanisms mediating resistance, overexpression of anti-apoptotic Bcl-2 proteins (Bcl-2, Bcl-xL , Mcl-1) is particularly important. In this study, we investigated whether antiapoptotic protein patterns are prognostically relevant and potential therapeutic targets in CCA. Bcl-2 proteins were analysed in a pan-cancer cohort from the NCT/DKFZ/DKTK MASTER registry trial (n = 1140, CCA n = 72) via RNA-sequencing and transcriptome-based protein activity interference revealing high ranks of CCA for Bcl-xL and Mcl-1. Expression of Bcl-xL , Mcl-1, and Bcl-2 was assessed in human CCA tissue and cell lines compared with cholangiocytes by immunohistochemistry, immunoblotting, and quantitative-RT-PCR. Immunohistochemistry confirmed the upregulation of Bcl-xL and Mcl-1 in iCCA tissues. Cell death of CCA cell lines upon treatment with specific small molecule inhibitors of Bcl-xL (Wehi-539), of Mcl-1 (S63845), and Bcl-2 (ABT-199), either alone, in combination with each other or together with chemotherapeutics was assessed by flow cytometry. Targeting Bcl-xL induced cell death and augmented the effect of chemotherapy in CCA cells. Combined inhibition of Bcl-xL and Mcl-1 led to a synergistic increase in cell death in CCA cell lines. Correlation between Bcl-2 protein expression and survival was analysed within three independent patient cohorts from cancer centers in Germany comprising 656 CCA cases indicating a prognostic value of Bcl-xL in CCA depending on the CCA subtype. Collectively, these observations identify Bcl-xL as a key protein in cell death resistance of CCA and may pave the way for clinical application.

Comprehensive genomic and epigenomic analysis in cancer of unknown primary guides molecularly-informed therapies despite heterogeneity.

The benefit of molecularly-informed therapies in cancer of unknown primary (CUP) is unclear. Here, we use comprehensive molecular characterization by whole genome/exome, transcriptome and methylome analysis in 70 CUP patients to reveal substantial mutational heterogeneity with TP53, MUC16, KRAS, LRP1B and CSMD3 being the most frequently mutated known cancer-related genes. The most common fusion partner is FGFR2, the most common focal homozygous deletion affects CDKN2A. 56/70 (80%) patients receive genomics-based treatment recommendations which are applied in 20/56 (36%) cases. Transcriptome and methylome data provide evidence for the underlying entity in 62/70 (89%) cases. Germline analysis reveals five (likely) pathogenic mutations in five patients. Recommended off-label therapies translate into a mean PFS ratio of 3.6 with a median PFS1 of 2.9 months (17 patients) and a median PFS2 of 7.8 months (20 patients). Our data emphasize the clinical value of molecular analysis and underline the need for innovative, mechanism-based clinical trials.

Genomic architecture of FGFR2 fusions in cholangiocarcinoma and its implication for molecular testing.

BACKGROUND: Cholangiocarcinoma (CCA) is a primary malignancy of the biliary tract with a dismal prognosis. Recently, several actionable genetic aberrations were identified with significant enrichment in intrahepatic CCA, including FGFR2 gene fusions with a prevalence of 10-15%. Recent clinical data demonstrate that these fusions are druggable in a second-line setting in advanced/metastatic disease and the efficacy in earlier lines of therapy is being evaluated in ongoing clinical trials. This scenario warrants standardised molecular profiling of these tumours. METHODS: A detailed analysis of the original genetic data from the FIGHT-202 trial, on which the approval of Pemigatinib was based, was conducted. RESULTS: Comparing different detection approaches and displaying representative cases, we described the genetic landscape and architecture of FGFR2 fusions in iCCA and show biological and technical aspects to be considered for their detection. We elaborated parameters, including a suggestion for annotation, that should be stated in a molecular diagnostic FGFR2 report to allow a complete understanding of the analysis performed and the information provided. CONCLUSION: This study provides a detailed presentation and dissection of the technical and biological aspects regarding FGFR2 fusion detection, which aims to support molecular pathologists, pathologists and clinicians in diagnostics, reporting of the results and decision-making.

MCL1 as putative target in pancreatoblastoma.

Pancreatoblastoma (PB) is a rare tumor of the pancreas. In case of metastases, the treatment options are sparse and targeted approaches are not developed. We here evaluate MCL1 amplification as a putative target in PB.Thirteen samples from adult (10/13) and pediatric patients (3/13) were collected. Three of these samples had been previously subjected to whole-exome sequencing (2 cases) or whole-genome sequencing (1 case) within a precision oncology program (NCT/DKTK MASTER), and this analysis had shown copy number gains of MCL1 gene. We established a fluorescence in situ hybridization (FISH) test to assess the copy number alterations of MCL1 gene in 13 formalin-fixed paraffin-embedded PBs, including the 3 cases assessed by genome sequencing. FISH analysis showed the amplification of MCL1 in 2 cases (both were adult PB), one of which was a case with the highest copy number gain at genomic analysis. In both cases, the average gene copy number per cell was ≥ 5.7 and the MCL1/1p12 ratio was ≥ 2.4. Our data support MCL1 as a putative target in PB. Patients with MCL1-amplified PB might benefit from MCL1 inhibition. Sequencing data is useful to screen for amplification; however, the established FISH for MCL1 can help to determine the level and cellular heterogeneity of MCL1 amplification more accurately.

The genomic and transcriptional landscape of primary central nervous system lymphoma.

Primary lymphomas of the central nervous system (PCNSL) are mainly diffuse large B-cell lymphomas (DLBCLs) confined to the central nervous system (CNS). Molecular drivers of PCNSL have not been fully elucidated. Here, we profile and compare the whole-genome and transcriptome landscape of 51 CNS lymphomas (CNSL) to 39 follicular lymphoma and 36 DLBCL cases outside the CNS. We find recurrent mutations in JAK-STAT, NFkB, and B-cell receptor signaling pathways, including hallmark mutations in MYD88 L265P (67%) and CD79B (63%), and CDKN2A deletions (83%). PCNSLs exhibit significantly more focal deletions of HLA-D (6p21) locus as a potential mechanism of immune evasion. Mutational signatures correlating with DNA replication and mitosis are significantly enriched in PCNSL. TERT gene expression is significantly higher in PCNSL compared to activated B-cell (ABC)-DLBCL. Transcriptome analysis clearly distinguishes PCNSL and systemic DLBCL into distinct molecular subtypes. Epstein-Barr virus (EBV)+ CNSL cases lack recurrent mutational hotspots apart from IG and HLA-DRB loci. We show that PCNSL can be clearly distinguished from DLBCL, having distinct expression profiles, IG expression and translocation patterns, as well as specific combinations of genetic alterations.

Comprehensive genomic analysis of refractory multiple myeloma reveals a complex mutational landscape associated with drug resistance and novel therapeutic vulnerabilities.

The outcomes of patients with multiple myeloma (MM) refractory to immunomodulatory agents (IMiDs) and proteasome inhibitors (PIs) remain poor. In this study, we performed whole genome and transcriptome sequencing of 39 heavily pretreated relapsed/refractory MM (RRMM) patients to identify mechanisms of resistance and potential therapeutic targets. We observed a high mutational load and indications of increased genomic instability. Recurrently mutated genes in RRMM, which had not been previously reported or only observed at a lower frequency in newly diagnosed MM, included NRAS, BRAF, TP53, SLC4A7, MLLT4, EWSR1, HCFC2, and COPS3. We found multiple genomic regions with bi-allelic events affecting tumor suppressor genes and demonstrated a significant adverse impact of bi-allelic TP53 alterations on survival. With regard to potentially resistance conferring mutations, recurrently mutated gene networks included genes with relevance for PI and IMiD activity; the latter particularly affecting members of the Cereblon and the COP9 signalosome complex. We observed a major impact of signatures associated with exposure to melphalan or impaired DNA double-strand break homologous recombination repair in RRMM. The latter coincided with mutations in genes associated with PARP inhibitor sensitivity in 49% of RRMM patients; a finding with potential therapeutic implications. In conclusion, this comprehensive genomic characterization revealed a complex mutational and structural landscape in RRMM and highlights potential implications for therapeutic strategies.

Comprehensive Genomic and Transcriptomic Analysis for Guiding Therapeutic Decisions in Patients with Rare Cancers.

The clinical relevance of comprehensive molecular analysis in rare cancers is not established. We analyzed the molecular profiles and clinical outcomes of 1,310 patients (rare cancers, 75.5%) enrolled in a prospective observational study by the German Cancer Consortium that applies whole-genome/exome and RNA sequencing to inform the care of adults with incurable cancers. On the basis of 472 single and six composite biomarkers, a cross-institutional molecular tumor board provided evidence-based management recommendations, including diagnostic reevaluation, genetic counseling, and experimental treatment, in 88% of cases. Recommended therapies were administered in 362 of 1,138 patients (31.8%) and resulted in significantly improved overall response and disease control rates (23.9% and 55.3%) compared with previous therapies, translating into a progression-free survival ratio >1.3 in 35.7% of patients. These data demonstrate the benefit of molecular stratification in rare cancers and represent a resource that may promote clinical trial access and drug approvals in this underserved patient population. SIGNIFICANCE: Rare cancers are difficult to treat; in particular, molecular pathogenesis-oriented medical therapies are often lacking. This study shows that whole-genome/exome and RNA sequencing enables molecularly informed treatments that lead to clinical benefit in a substantial proportion of patients with advanced rare cancers and paves the way for future clinical trials.See related commentary by Eggermont et al., p. 2677.This article is highlighted in the In This Issue feature, p. 2659.

A community challenge to evaluate RNA-seq, fusion detection, and isoform quantification methods for cancer discovery.

The accurate identification and quantitation of RNA isoforms present in the cancer transcriptome is key for analyses ranging from the inference of the impacts of somatic variants to pathway analysis to biomarker development and subtype discovery. The ICGC-TCGA DREAM Somatic Mutation Calling in RNA (SMC-RNA) challenge was a crowd-sourced effort to benchmark methods for RNA isoform quantification and fusion detection from bulk cancer RNA sequencing (RNA-seq) data. It concluded in 2018 with a comparison of 77 fusion detection entries and 65 isoform quantification entries on 51 synthetic tumors and 32 cell lines with spiked-in fusion constructs. We report the entries used to build this benchmark, the leaderboard results, and the experimental features associated with the accurate prediction of RNA species. This challenge required submissions to be in the form of containerized workflows, meaning each of the entries described is easily reusable through CWL and Docker containers at https://github.com/SMC-RNA-challenge. A record of this paper's transparent peer review process is included in the supplemental information.

Deconvolution of sarcoma methylomes reveals varying degrees of immune cell infiltrates with association to genomic aberrations.

BACKGROUND: Soft-tissue sarcomas (STS) are a heterogeneous group of mesenchymal tumors for which response to immunotherapies is not well established. Therefore, it is important to risk-stratify and identify STS patients who will most likely benefit from these treatments. RESULTS: To reveal shared and distinct methylation signatures present in STS, we performed unsupervised deconvolution of DNA methylation data from the TCGA sarcoma and an independent validation cohort. We showed that leiomyosarcoma can be subclassified into three distinct methylation groups. More importantly, we identified a component associated with tumor-infiltrating leukocytes, which suggests varying degrees of immune cell infiltration in STS subtypes and an association with prognosis. We further investigated the genomic alterations that may influence tumor infiltration by leukocytes including RB1 loss in undifferentiated pleomorphic sarcomas and ELK3 amplification in dedifferentiated liposarcomas. CONCLUSIONS: In summary, we have leveraged unsupervised methylation-based deconvolution to characterize the immune compartment and molecularly stratify subtypes in STS, which may benefit precision medicine in the future.

Accurate and efficient detection of gene fusions from RNA sequencing data.

The identification of gene fusions from RNA sequencing data is a routine task in cancer research and precision oncology. However, despite the availability of many computational tools, fusion detection remains challenging. Existing methods suffer from poor prediction accuracy and are computationally demanding. We developed Arriba, a novel fusion detection algorithm with high sensitivity and short runtime. When applied to a large collection of published pancreatic cancer samples (n = 803), Arriba identified a variety of driver fusions, many of which affected druggable proteins, including ALK, BRAF, FGFR2, NRG1, NTRK1, NTRK3, RET, and ROS1. The fusions were significantly associated with KRAS wild-type tumors and involved proteins stimulating the MAPK signaling pathway, suggesting that they substitute for activating mutations in KRAS In addition, we confirmed the transforming potential of two novel fusions, RRBP1-RAF1 and RASGRP1-ATP1A1, in cellular assays. These results show Arriba's utility in both basic cancer research and clinical translation.