ZygosityPredictor.

Summary: ZygosityPredictor provides functionality to evaluate how many copies of a gene are affected by mutations in next generation sequencing data. In cancer samples, the tool processes both somatic and germline mutations. In particular, ZygosityPredictor computes the number of affected copies for single nucleotide variants and small insertions and deletions (Indels). In addition, the tool integrates information at gene level via phasing of several variants and subsequent logic to derive how strongly a gene is affected by mutations and provides a measure of confidence. This information is of particular interest in precision oncology, e.g. when assessing whether unmutated copies of tumor-suppressor genes remain. Availability and implementation: ZygosityPredictor was implemented as an R-package and is available via Bioconductor at https://bioconductor.org/packages/ZygosityPredictor. Detailed documentation is provided in the vignette including application to an example genome.

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.

BRAFΔß3-αC in-frame deletion mutants differ in their dimerization propensity, HSP90 dependence, and druggability.

In-frame BRAF exon 12 deletions are increasingly identified in various tumor types. The resultant BRAFΔß3-αC oncoproteins usually lack five amino acids in the ß3-αC helix linker and sometimes contain de novo insertions. The dimerization status of BRAFΔß3-αC oncoproteins, their precise pathomechanism, and their direct druggability by RAF inhibitors (RAFi) has been under debate. Here, we functionally characterize BRAFΔLNVTAP>F and two novel mutants, BRAFdelinsFS and BRAFΔLNVT>F, and compare them with other BRAFΔß3-αC oncoproteins. We show that BRAFΔß3-αC oncoproteins not only form stable homodimers and large multiprotein complexes but also require dimerization. Nevertheless, details matter as aromatic amino acids at the deletion junction of some BRAFΔß3-αC oncoproteins, e.g., BRAFΔLNVTAP>F, increase their stability and dimerization propensity while conferring resistance to monomer-favoring RAFi such as dabrafenib or HSP 90/CDC37 inhibition. In contrast, dimer-favoring inhibitors such as naporafenib inhibit all BRAFΔß3-αC mutants in cell lines and patient-derived organoids, suggesting that tumors driven by such oncoproteins are vulnerable to these compounds.

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.

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.

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.

Successful BRAF/MEK inhibition in a patient with BRAFV600E-mutated extrapancreatic acinar cell carcinoma.

Pancreatic acinar cell carcinoma (PAC) is a rare disease with a poor prognosis. Treatment options for metastatic PAC are limited and often follow chemotherapeutic regimens for pancreatic ductal adenocarcinoma. Although recurrent genomic alterations, such as BRAF fusions and defects in genes involved in homologous recombination DNA repair, have been described in PAC, data on the clinical efficacy of molecularly guided, targeted treatment are scarce. Here we describe the case of a 27-yr-old patient with BRAFV600E-mutated PAC who was successfully treated with a combination of BRAF and MEK inhibitors. The patient presented to our clinic with abdominal pain and weight loss. Imaging showed extensive retroperitoneal disease as well as mediastinal lymphadenopathy. Because of elevated α-fetoprotein (AFP) levels and inconclusive histologic findings, a germ cell tumor was suspected; however, PEI chemotherapy was unsuccessful. A repeat biopsy yielded the diagnosis of PAC and treatment with FOLFIRINOX was initiated. Comprehensive molecular profiling within the MASTER (Molecularly Aided Stratification for Tumor Eradication Research) precision oncology program revealed a somatic BRAFV600E mutation and a germline PALB2 stop-gain mutation. Therapy was therefore switched to BRAF/MEK inhibition, resulting in almost complete remission and disease control for 12 mo and a remarkable improvement in the patient's general condition. These results indicate that BRAF alterations are a valid therapeutic target in PAC that should be routinely assessed in this patient population.

Germline SDHB-inactivating mutation in gastric spindle cell sarcoma.

Gastrointestinal stromal tumors (GISTs) are the most frequent mesenchymal tumors of the gastrointestinal tract. Inactivating mutations or epigenetic deregulation of succinate dehydrogenase complex (SDH) genes are considered defining features of a subset of GIST occurring in the stomach. Based on comprehensive molecular profiling and biochemical analysis within a precision oncology program, we identified hallmarks of SDH deficiency (germline SDHB-inactivating mutation accompanied by somatic loss of heterozygosity, lack of SDHB expression, global DNA hypermethylation, and elevated succinate/fumarate ratio) in a 40-year-old woman with undifferentiated gastric spindle cell sarcoma that did not meet the diagnostic criteria for other mesenchymal tumors of the stomach, including GIST. These data reveal that the loss of SDH function can be involved in the pathogenesis of non-GIST sarcoma of the gastrointestinal tract.

The landscape of chromothripsis across adult cancer types.

Chromothripsis is a recently identified mutational phenomenon, by which a presumably single catastrophic event generates extensive genomic rearrangements of one or a few chromosome(s). Considered as an early event in tumour development, this form of genome instability plays a prominent role in tumour onset. Chromothripsis prevalence might have been underestimated when using low-resolution methods, and pan-cancer studies based on sequencing are rare. Here we analyse chromothripsis in 28 tumour types covering all major adult cancers (634 tumours, 316 whole-genome and 318 whole-exome sequences). We show that chromothripsis affects a substantial proportion of human cancers, with a prevalence of 49% across all cases. Chromothripsis generates entity-specific genomic alterations driving tumour development, including clinically relevant druggable fusions. Chromothripsis is linked with specific telomere patterns and univocal mutational signatures in distinct tumour entities. Longitudinal analysis of chromothriptic patterns in 24 matched tumour pairs reveals insights in the clonal evolution of tumours with chromothripsis.

Identification and characterization of a BRAF fusion oncoprotein with retained autoinhibitory domains.

Fusion proteins involving the BRAF serine/threonine kinase occur in many cancers. The oncogenic potential of BRAF fusions has been attributed to the loss of critical N-terminal domains that mediate BRAF autoinhibition. We used whole-exome and RNA sequencing in a patient with glioblastoma multiforme to identify a rearrangement between TTYH3, encoding a membrane-resident, calcium-activated chloride channel, and BRAF intron 1, resulting in a TTYH3-BRAF fusion protein that retained all features essential for BRAF autoinhibition. Accordingly, the BRAF moiety of the fusion protein alone, which represents full-length BRAF without the amino acids encoded by exon 1 (BRAFΔE1), did not induce MEK/ERK phosphorylation or transformation. Likewise, neither the TTYH3 moiety of the fusion protein nor full-length TTYH3 provoked ERK pathway activity or transformation. In contrast, TTYH3-BRAF displayed increased MEK phosphorylation potential and transforming activity, which were caused by TTYH3-mediated tethering of near-full-length BRAF to the (endo)membrane system. Consistent with this mechanism, a synthetic approach, in which BRAFΔE1 was tethered to the membrane by fusing it to the cytoplasmic tail of CD8 also induced transformation. Furthermore, we demonstrate that TTYH3-BRAF signals largely independent of a functional RAS binding domain, but requires an intact BRAF dimer interface and activation loop phosphorylation sites. Cells expressing TTYH3-BRAF exhibited increased MEK/ERK signaling, which was blocked by clinically achievable concentrations of sorafenib, trametinib, and the paradox breaker PLX8394. These data provide the first example of a fully autoinhibited BRAF protein whose oncogenic potential is dictated by a distinct fusion partner and not by a structural change in BRAF itself.

Variant classification in precision oncology.

Next-generation sequencing has become a cornerstone of therapy guidance in cancer precision medicine and an indispensable research tool in translational oncology. Its rapidly increasing use during the last decade has expanded the options for targeted tumor therapies, and molecular tumor boards have grown accordingly. However, with increasing detection of genetic alterations, their interpretation has become more complex and error-prone, potentially introducing biases and reducing benefits in clinical practice. To facilitate interdisciplinary discussions of genetic alterations for treatment stratification between pathologists, oncologists, bioinformaticians, genetic counselors and medical scientists in specialized molecular tumor boards, several systems for the classification of variants detected by large-scale sequencing have been proposed. We review three recent and commonly applied classifications and discuss their individual strengths and weaknesses. Comparison of the classifications underlines the need for a clinically useful and universally applicable variant reporting system, which will be instrumental for efficient decision making based on sequencing analysis in oncology. Integrating these data, we propose a generalizable classification concept featuring a conservative and a more progressive scheme, which can be readily applied in a clinical setting.

Defective homologous recombination DNA repair as therapeutic target in advanced chordoma.

Chordomas are rare bone tumors with few therapeutic options. Here we show, using whole-exome and genome sequencing within a precision oncology program, that advanced chordomas (n = 11) may be characterized by genomic patterns indicative of defective homologous recombination (HR) DNA repair and alterations affecting HR-related genes, including, for example, deletions and pathogenic germline variants of BRCA2, NBN, and CHEK2. A mutational signature associated with HR deficiency was significantly enriched in 72.7% of samples and co-occurred with genomic instability. The poly(ADP-ribose) polymerase (PARP) inhibitor olaparib, which is preferentially toxic to HR-incompetent cells, led to prolonged clinical benefit in a patient with refractory chordoma, and whole-genome analysis at progression revealed a PARP1 p.T910A mutation predicted to disrupt the autoinhibitory PARP1 helical domain. These findings uncover a therapeutic opportunity in chordoma that warrants further exploration, and provide insight into the mechanisms underlying PARP inhibitor resistance.

NRG1 Fusions in KRAS Wild-Type Pancreatic Cancer.

We used whole-genome and transcriptome sequencing to identify clinically actionable genomic alterations in young adults with pancreatic ductal adenocarcinoma (PDAC). Molecular characterization of 17 patients with PDAC enrolled in a precision oncology program revealed gene fusions amenable to pharmacologic inhibition by small-molecule tyrosine kinase inhibitors in all patients with KRAS wild-type (KRASWT) tumors (4 of 17). These alterations included recurrent NRG1 rearrangements predicted to drive PDAC development through aberrant ERBB receptor-mediated signaling, and pharmacologic ERBB inhibition resulted in clinical improvement and remission of liver metastases in 2 patients with NRG1-rearranged tumors that had proved resistant to standard treatment. Our findings demonstrate that systematic screening of KRASWT tumors for oncogenic fusion genes will substantially improve the therapeutic prospects for a sizeable fraction of patients with PDAC.Significance: Advanced PDAC is a malignancy with few treatment options that lacks molecular mechanism-based therapies. Our study uncovers recurrent gene rearrangements such as NRG1 fusions as disease-driving events in KRASwt tumors, thereby providing novel insights into oncogenic signaling and new therapeutic options in this entity. Cancer Discov; 8(9); 1087-95. ©2018 AACR.This article is highlighted in the In This Issue feature, p. 1047.

From somatic variants towards precision oncology: Evidence-driven reporting of treatment options in molecular tumor boards.

BACKGROUND: A comprehensive understanding of cancer has been furthered with technological improvements and decreasing costs of next-generation sequencing (NGS). However, the complexity of interpreting genomic data is hindering the implementation of high-throughput technologies in the clinical context: increasing evidence on gene-drug interactions complicates the task of assigning clinical significance to genomic variants. METHODS: Here we present a method that automatically matches patient-specific genomic alterations to treatment options. The method relies entirely on public knowledge of somatic variants with predictive evidence on drug response. The output report is aimed at supporting clinicians in the task of finding the clinical meaning of genomic variants. We applied the method to 1) The Cancer Genome Atlas (TCGA) and Genomics Evidence Neoplasia Information Exchange (GENIE) cohorts and 2) 11 patients from the NCT MASTER trial whose treatment discussions included information on their genomic profiles. RESULTS: Our reporting strategy showed a substantial number of patients with actionable variants in the analyses of TCGA and GENIE samples. Notably, it was able to reproduce experts' treatment suggestions in a retrospective study of 11 patients from the NCT MASTER trial. Our results establish a proof of concept for comprehensive, evidence-based reports as a supporting tool for discussing treatment options in tumor boards. CONCLUSIONS: We believe that a standardized method to report actionable somatic variants will smooth the incorporation of NGS in the clinical context. We anticipate that tools like the one we present here will become essential in summarizing for clinicians the growing evidence in the field of precision medicine. The R code of the presented method is provided in Additional file 6 and available at https://github.com/jperera-bel/MTB-Report .

Validating Comprehensive Next-Generation Sequencing Results for Precision Oncology: The NCT/DKTK Molecularly Aided Stratification for Tumor Eradication Research Experience.

PURPOSE: Rapidly evolving genomics technologies, in particular comprehensive next-generation sequencing (NGS), have led to exponential growth in the understanding of cancer biology, shifting oncology toward personalized treatment strategies. However, comprehensive NGS approaches, such as whole-exome sequencing, have limitations that are related to the technology itself as well as to the input source. Hence, clinical implementation of comprehensive NGS in a quality-controlled diagnostic workflow requires both the standardization of sequencing procedures and continuous validation of sequencing results by orthogonal methods in an ongoing program to enable the determination of key test parameters and continuous improvement of NGS and bioinformatics pipelines. PATIENTS AND METHODS: We present validation data on 220 patients who were enrolled between 2013 and 2016 in a multi-institutional, genomics-guided precision oncology program (Molecularly Aided Stratification for Tumor Eradication Research) of the National Center for Tumor Diseases Heidelberg and the German Cancer Consortium. RESULTS: More than 90% of clinically actionable genomic alterations identified by combined whole-exome sequencing and transcriptome sequencing were successfully validated, with varying frequencies of discordant results across different types of alterations (fusions, 3.7%; single-nucleotide variants, 2.6%; amplifications, 1.1%; overexpression, 0.9%; deletions, 0.6%). The implementation of new computational methods for NGS data analysis led to a substantial improvement of gene fusion calling over time. CONCLUSION: Collectively, these data demonstrate the value of a rigorous validation program that partners with comprehensive NGS to successfully implement and continuously improve cancer precision medicine in a clinical setting.

Precision oncology based on omics data: The NCT Heidelberg experience.

Precision oncology implies the ability to predict which patients will likely respond to specific cancer therapies based on increasingly accurate, high-resolution molecular diagnostics as well as the functional and mechanistic understanding of individual tumors. While molecular stratification of patients can be achieved through different means, a promising approach is next-generation sequencing of tumor DNA and RNA, which can reveal genomic alterations that have immediate clinical implications. Furthermore, certain genetic alterations are shared across multiple histologic entities, raising the fundamental question of whether tumors should be treated by molecular profile and not tissue of origin. We here describe MASTER (Molecularly Aided Stratification for Tumor Eradication Research), a clinically applicable platform for prospective, biology-driven stratification of younger adults with advanced-stage cancer across all histologies and patients with rare tumors. We illustrate how a standardized workflow for selection and consenting of patients, sample processing, whole-exome/genome and RNA sequencing, bioinformatic analysis, rigorous validation of potentially actionable findings, and data evaluation by a dedicated molecular tumor board enables categorization of patients into different intervention baskets and formulation of evidence-based recommendations for clinical management. Critical next steps will be to increase the number of patients that can be offered comprehensive molecular analysis through collaborations and partnering, to explore ways in which additional technologies can aid in patient stratification and individualization of treatment, to stimulate clinically guided exploratory research projects, and to gradually move away from assessing the therapeutic activity of targeted interventions on a case-by-case basis toward controlled clinical trials of genomics-guided treatments.

Integration of genomics and histology revises diagnosis and enables effective therapy of refractory cancer of unknown primary with PDL1 amplification.

Identification of the tissue of origin in cancer of unknown primary (CUP) poses a diagnostic challenge and is critical for directing site-specific therapy. Currently, clinical decision-making in patients with CUP primarily relies on histopathology and clinical features. Comprehensive molecular profiling has the potential to contribute to diagnostic categorization and, most importantly, guide CUP therapy through identification of actionable lesions. We here report the case of an advanced-stage malignancy initially mimicking poorly differentiated soft-tissue sarcoma that did not respond to multiagent chemotherapy. Molecular profiling within a clinical whole-exome and transcriptome sequencing program revealed a heterozygous, highly amplified KRAS G12S mutation, compound-heterozygous TP53 mutation/deletion, high mutational load, and focal high-level amplification of Chromosomes 9p (including PDL1 [CD274] and JAK2) and 10p (including GATA3). Integrated analysis of molecular data and histopathology provided a rationale for immune checkpoint inhibitor (ICI) therapy with pembrolizumab, which resulted in rapid clinical improvement and a lasting partial remission. Histopathological analyses ruled out sarcoma and established the diagnosis of a poorly differentiated adenocarcinoma. Although neither histopathology nor molecular data were able to pinpoint the tissue of origin, our analyses established several differential diagnoses including triple-negative breast cancer (TNBC). We analyzed 157 TNBC samples from The Cancer Genome Atlas, revealing PDL1 copy number gains coinciding with excessive PDL1 mRNA expression in 24% of cases. Collectively, these results illustrate the impact of multidimensional tumor profiling in cases with nondescript histology and immunophenotype, show the predictive potential of PDL1 amplification for immune checkpoint inhibitors (ICIs), and suggest a targeted therapeutic strategy in Chromosome 9p24.1/PDL1-amplified cancers.

A Thermodynamic Model of Monovalent Cation Homeostasis in the Yeast Saccharomyces cerevisiae.

Cationic and heavy metal toxicity is involved in a substantial number of diseases in mammals and crop plants. Therefore, the understanding of tightly regulated transporter activities, as well as conceiving the interplay of regulatory mechanisms, is of substantial interest. A generalized thermodynamic description is developed for the complex interplay of the plasma membrane ion transporters, membrane potential and the consumption of energy for maintaining and restoring specific intracellular cation concentrations. This concept is applied to the homeostasis of cation concentrations in the yeast cells of S. cerevisiae. The thermodynamic approach allows to model passive ion fluxes driven by the electrochemical potential differences, but also primary or secondary active transport processes driven by the inter- play of different ions (symport, antiport) or by ATP consumption (ATPases). The model-confronted with experimental data-reproduces the experimentally observed potassium and proton fluxes induced by the external stimuli KCl and glucose. The estimated phenomenological constants combine kinetic parameters and transport coefficients. These are in good agreement with the biological understanding of the transporters thus providing a better understanding of the control exerted by the coupled fluxes. The model predicts the flux of additional ion species, like e.g. chloride, as a potential candidate for counterbalancing positive charges. Furthermore, the effect of a second KCl stimulus is simulated, predicting a reduced cellular response for cells that were first exposed to a high KCl stimulus compared to cells pretreated with a mild KCl stimulus. By describing the generalized forces that are responsible for a given flow, the model provides information and suggestions for new experiments. Furthermore, it can be extended to other systems such as e.g. Candida albicans, or selected plant cells.