Genomic Patterns of Malignant Peripheral Nerve Sheath Tumor (MPNST) Evolution Correlate with Clinical Outcome and Are Detectable in Cell-Free DNA.

Malignant peripheral nerve sheath tumor (MPNST), an aggressive soft-tissue sarcoma, occurs in people with neurofibromatosis type 1 (NF1) and sporadically. Whole-genome and multiregional exome sequencing, transcriptomic, and methylation profiling of 95 tumor samples revealed the order of genomic events in tumor evolution. Following biallelic inactivation of NF1, loss of CDKN2A or TP53 with or without inactivation of polycomb repressive complex 2 (PRC2) leads to extensive somatic copy-number aberrations (SCNA). Distinct pathways of tumor evolution are associated with inactivation of PRC2 genes and H3K27 trimethylation (H3K27me3) status. Tumors with H3K27me3 loss evolve through extensive chromosomal losses followed by whole-genome doubling and chromosome 8 amplification, and show lower levels of immune cell infiltration. Retention of H3K27me3 leads to extensive genomic instability, but an immune cell-rich phenotype. Specific SCNAs detected in both tumor samples and cell-free DNA (cfDNA) act as a surrogate for H3K27me3 loss and immune infiltration, and predict prognosis. SIGNIFICANCE: MPNST is the most common cause of death and morbidity for individuals with NF1, a relatively common tumor predisposition syndrome. Our results suggest that somatic copy-number and methylation profiling of tumor or cfDNA could serve as a biomarker for early diagnosis and to stratify patients into prognostic and treatment-related subgroups. This article is highlighted in the In This Issue feature, p. 517.

Bayesian networks elucidate complex genomic landscapes in cancer.

Bayesian networks (BNs) are disciplined, explainable Artificial Intelligence models that can describe structured joint probability spaces. In the context of understanding complex relations between a number of variables in biological settings, they can be constructed from observed data and can provide a guiding, graphical tool in exploring such relations. Here we propose BNs for elucidating the relations between driver events in large cancer genomic datasets. We present a methodology that is specifically tailored to biologists and clinicians as they are the main producers of such datasets. We achieve this by using an optimal BN learning algorithm based on well established likelihood functions and by utilising just two tuning parameters, both of which are easy to set and have intuitive readings. To enhance value to clinicians, we introduce (a) the use of heatmaps for families in each network, and (b) visualising pairwise co-occurrence statistics on the network. For binary data, an optional step of fitting logic gates can be employed. We show how our methodology enhances pairwise testing and how biologists and clinicians can use BNs for discussing the main relations among driver events in large genomic cohorts. We demonstrate the utility of our methodology by applying it to 5 cancer datasets revealing complex genomic landscapes. Our networks identify central patterns in all datasets including a central 4-way mutual exclusivity between HDR, t(4,14), t(11,14) and t(14,16) in myeloma, and a 3-way mutual exclusivity of three major players: CALR, JAK2 and MPL, in myeloproliferative neoplasms. These analyses demonstrate that our methodology can play a central role in the study of large genomic cancer datasets.

Functional Genomic Identification of Predictors of Sensitivity and Mechanisms of Resistance to Multivalent Second-Generation TRAIL-R2 Agonists.

Multivalent second-generation TRAIL-R2 agonists are currently in late preclinical development and early clinical trials. Herein, we use a representative second-generation agent, MEDI3039, to address two major clinical challenges facing these agents: lack of predictive biomarkers to enable patient selection and emergence of resistance. Genome-wide CRISPR knockout screens were notable for the lack of resistance mechanisms beyond the canonical TRAIL-R2 pathway (caspase-8, FADD, BID) as well as p53 and BAX in TP53 wild-type models, whereas a CRISPR activatory screen identified cell death inhibitors MCL-1 and BCL-XL as mechanisms to suppress MEDI3039-induced cell death. High-throughput drug screening failed to identify genomic alterations associated with response to MEDI3039; however, transcriptomics analysis revealed striking association between MEDI3039 sensitivity and expression of core components of the extrinsic apoptotic pathway, most notably its main apoptotic effector caspase-8 in solid tumor cell lines. Further analyses of colorectal cell lines and patient-derived xenografts identified caspase-8 expression ratio to its endogenous regulator FLIP(L) as predictive of sensitivity to MEDI3039 in several major solid tumor types and a further subset indicated by caspase-8:MCL-1 ratio. Subsequent MEDI3039 combination screening of TRAIL-R2, caspase-8, FADD, and BID knockout models with 60 compounds with varying mechanisms of action identified two inhibitor of apoptosis proteins (IAP) that exhibited strong synergy with MEDI3039 that could reverse resistance only in BID-deleted models. In summary, we identify the ratios of caspase-8:FLIP(L) and caspase-8:MCL-1 as potential predictive biomarkers for second-generation TRAIL-R2 agonists and loss of key effectors such as FADD and caspase-8 as likely drivers of clinical resistance in solid tumors.

Patterns of genomic change in residual disease after neoadjuvant chemotherapy for estrogen receptor-positive and HER2-negative breast cancer.

BACKGROUND: Treatment of patients with residual disease after neoadjuvant chemotherapy for breast cancer is an unmet clinical need. We hypothesised that tumour subclones showing expansion in residual disease after chemotherapy would contain mutations conferring drug resistance. METHODS: We studied oestrogen receptor and/or progesterone receptor-positive, HER2-negative tumours from 42 patients in the EORTC 10994/BIG 00-01 trial who failed to achieve a pathological complete response. Genes commonly mutated in breast cancer were sequenced in pre and post-treatment samples. RESULTS: Oncogenic driver mutations were commonest in PIK3CA (38% of tumours), GATA3 (29%), CDH1 (17%), TP53 (17%) and CBFB (12%); and amplification was commonest for CCND1 (26% of tumours) and FGFR1 (26%). The variant allele fraction frequently changed after treatment, indicating that subclones had expanded and contracted, but there were changes in both directions for all of the commonly mutated genes. CONCLUSIONS: We found no evidence that expansion of clones containing recurrent oncogenic driver mutations is responsible for resistance to neoadjuvant chemotherapy. The persistence of classic oncogenic mutations in pathways for which targeted therapies are now available highlights their importance as drug targets in patients who have failed chemotherapy but provides no support for a direct role of driver oncogenes in resistance to chemotherapy. CLINICALTRIALS.GOV: EORTC 10994/BIG 1-00 Trial registration number NCT00017095.

Image-based consensus molecular subtype (imCMS) classification of colorectal cancer using deep learning.

OBJECTIVE: Complex phenotypes captured on histological slides represent the biological processes at play in individual cancers, but the link to underlying molecular classification has not been clarified or systematised. In colorectal cancer (CRC), histological grading is a poor predictor of disease progression, and consensus molecular subtypes (CMSs) cannot be distinguished without gene expression profiling. We hypothesise that image analysis is a cost-effective tool to associate complex features of tissue organisation with molecular and outcome data and to resolve unclassifiable or heterogeneous cases. In this study, we present an image-based approach to predict CRC CMS from standard H&E sections using deep learning. DESIGN: Training and evaluation of a neural network were performed using a total of n=1206 tissue sections with comprehensive multi-omic data from three independent datasets (training on FOCUS trial, n=278 patients; test on rectal cancer biopsies, GRAMPIAN cohort, n=144 patients; and The Cancer Genome Atlas (TCGA), n=430 patients). Ground truth CMS calls were ascertained by matching random forest and single sample predictions from CMS classifier. RESULTS: Image-based CMS (imCMS) accurately classified slides in unseen datasets from TCGA (n=431 slides, AUC)=0.84) and rectal cancer biopsies (n=265 slides, AUC=0.85). imCMS spatially resolved intratumoural heterogeneity and provided secondary calls correlating with bioinformatic prediction from molecular data. imCMS classified samples previously unclassifiable by RNA expression profiling, reproduced the expected correlations with genomic and epigenetic alterations and showed similar prognostic associations as transcriptomic CMS. CONCLUSION: This study shows that a prediction of RNA expression classifiers can be made from H&E images, opening the door to simple, cheap and reliable biological stratification within routine workflows.

APOBEC3-dependent kataegis and TREX1-driven chromothripsis during telomere crisis.

Chromothripsis and kataegis are frequently observed in cancer and may arise from telomere crisis, a period of genome instability during tumorigenesis when depletion of the telomere reserve generates unstable dicentric chromosomes1-5. Here we examine the mechanism underlying chromothripsis and kataegis by using an in vitro telomere crisis model. We show that the cytoplasmic exonuclease TREX1, which promotes the resolution of dicentric chromosomes4, plays a prominent role in chromothriptic fragmentation. In the absence of TREX1, the genome alterations induced by telomere crisis primarily involve breakage-fusion-bridge cycles and simple genome rearrangements rather than chromothripsis. Furthermore, we show that the kataegis observed at chromothriptic breakpoints is the consequence of cytosine deamination by APOBEC3B. These data reveal that chromothripsis and kataegis arise from a combination of nucleolytic processing by TREX1 and cytosine editing by APOBEC3B.

Molecular Evolution of IDH Wild-Type Glioblastomas Treated With Standard of Care Affects Survival and Design of Precision Medicine Trials: A Report From the EORTC 1542 Study.

PURPOSE: Precision medicine trials in glioblastoma (GBM) are often conducted at tumor recurrence. However, second surgeries for recurrent GBM are not routinely performed, and therefore, molecular data for trial inclusion are predominantly derived from the primary sample. This study aims to establish whether molecular targets change during tumor progression and, if so, whether this affects precision medicine trial design. MATERIALS AND METHODS: We collected 186 pairs of primary-recurrent GBM samples from patients receiving chemoradiotherapy with temozolomide and sequenced approximately 300 cancer genes. MGMT, TERT, and EGFRvIII status was individually determined. RESULTS: The molecular profile of our cohort was identical to that of other GBM cohorts (IDH wild-type [WT], 95%; EGFR amplified, approximately 50%), indicating that patients amenable to second surgery do not represent a specific molecular subtype. Molecular events in IDH WT GBMs were stable in approximately 80% of events, but changes in mutation status were observed for all examined genes (range, approximately 90% and 60% for TERT and EGFR mutations, respectively), and such changes strongly affected targeted trial size and design. A similar pattern of GBM driver instability was observed within MGMT promoter-methylated tumors. MGMT promoter methylation status remained prognostic at tumor recurrence. The observation that hypermutation at GBM recurrence was rare (8%) and not correlated with outcome was relevant for immunotherapy-based treatments. CONCLUSION: This large cohort of matched primary and recurrent IDH WT tumors establishes the frequency of GBM driver instability after chemoradiotherapy with temozolomide. This allows per gene or pathway calculation of trial size at tumor recurrence, using molecular data of the primary tumor only. We also identify genes for which repeat surgery is necessary because of low mutation retention rate.

PDX Finder: A portal for patient-derived tumor xenograft model discovery.

Patient-derived tumor xenograft (PDX) mouse models are a versatile oncology research platform for studying tumor biology and for testing chemotherapeutic approaches tailored to genomic characteristics of individual patients' tumors. PDX models are generated and distributed by a diverse group of academic labs, multi-institution consortia and contract research organizations. The distributed nature of PDX repositories and the use of different metadata standards for describing model characteristics presents a significant challenge to identifying PDX models relevant to specific cancer research questions. The Jackson Laboratory and EMBL-EBI are addressing these challenges by co-developing PDX Finder, a comprehensive open global catalog of PDX models and their associated datasets. Within PDX Finder, model attributes are harmonized and integrated using a previously developed community minimal information standard to support consistent searching across the originating resources. Links to repositories are provided from the PDX Finder search results to facilitate model acquisition and/or collaboration. The PDX Finder resource currently contains information for 1985 PDX models of diverse cancers including those from large resources such as the Patient-Derived Models Repository, PDXNet and EurOPDX. Individuals or organizations that generate and distribute PDXs are invited to increase the 'findability' of their models by participating in the PDX Finder initiative at www.pdxfinder.org.

Roles of cofactors and chromatin accessibility in Hox protein target specificity.

BACKGROUND: The regulation of specific target genes by transcription factors is central to our understanding of gene network control in developmental and physiological processes yet how target specificity is achieved is still poorly understood. This is well illustrated by the Hox family of transcription factors as their limited in vitro DNA-binding specificity contrasts with their clear in vivo functional specificity. RESULTS: We generated genome-wide binding profiles for three Hox proteins, Ubx, Abd-A and Abd-B, following transient expression in Drosophila Kc167 cells, revealing clear target specificity and a striking influence of chromatin accessibility. In the absence of the TALE class homeodomain cofactors Exd and Hth, Ubx and Abd-A bind at a very similar set of target sites in accessible chromatin, whereas Abd-B binds at an additional specific set of targets. Provision of Hox cofactors Exd and Hth considerably modifies the Ubx genome-wide binding profile enabling Ubx to bind at an additional novel set of targets. Both the Abd-B specific targets and the cofactor-dependent Ubx targets are in chromatin that is relatively DNase1 inaccessible prior to the expression of Hox proteins/Hox cofactors. CONCLUSIONS: Our experiments demonstrate a strong role for chromatin accessibility in Hox protein binding and suggest that Hox protein competition with nucleosomes has a major role in Hox protein target specificity in vivo.