Comprehensive Analysis Identifies Variability in PI3K Pathway Alterations in Triple-Negative Breast Cancer Subtypes.

PURPOSE: The PI3K pathway is frequently altered in triple-negative breast cancer (TNBC). Limited cell line and human data suggest that TNBC tumors characterized as mesenchymal (M) and luminal androgen receptor (LAR) subtypes have increased incidence of alterations in the PI3K pathway. The impact of PI3K pathway alterations across TNBC subtypes is poorly understood. METHODS: Pretreatment tumor was evaluated from operable TNBC patients enrolled on a clinical trial of neoadjuvant therapy (NAT; A Robust TNBC Evaluation fraMework to Improve Survival [ClinicalTrials.gov identifier: NCT02276443]). Tumors were characterized into seven TNBC subtypes per Pietenpol criteria (basal-like 1, basal-like 2, immunomodulatory, M, mesenchymal stem-like, LAR, and unstable). Using whole-exome sequencing, RNA sequencing, and immunohistochemistry for PTEN, alterations were identified in 32 genes known to activate the PI3K pathway. Alterations in each subtype were associated with pathologic response to NAT. RESULTS: In evaluated patients (N = 177), there was a significant difference in the incidence of PI3K pathway alterations across TNBC subtypes (P < .01). The highest incidence of alterations was seen in LAR (81%), BL2 (79%), and M (62%) subtypes. The odds ratio for pathologic complete response (pCR) in the presence of PIK3CA mutation, PTEN mutation, and/or PTEN loss was highest in the LAR subtype and lowest in the M subtype, but these findings did not reach statistical significance. Presence of PIK3CA mutation was associated with pCR in the LAR subtype (P = .02). CONCLUSION: PI3K pathway alteration can affect response to NAT in TNBC, and targeted agents may improve outcomes, particularly in patients with M and LAR TNBC.

PTEN in triple-negative breast carcinoma: protein expression and genomic alteration in pretreatment and posttreatment specimens.

Background: Recent advances have been made in targeting the phosphoinositide 3-kinase pathway in breast cancer. Phosphatase and tensin homolog (PTEN) is a key component of that pathway. Objective: To understand the changes in PTEN expression over the course of the disease in patients with triple-negative breast cancer (TNBC) and whether PTEN copy number variation (CNV) by next-generation sequencing (NGS) can serve as an alternative to immunohistochemistry (IHC) to identify PTEN loss. Methods: We compared PTEN expression by IHC between pretreatment tumors and residual tumors in the breast and lymph nodes after neoadjuvant chemotherapy in 96 patients enrolled in a TNBC clinical trial. A correlative analysis between PTEN protein expression and PTEN CNV by NGS was also performed. Results: With a stringent cutoff for PTEN IHC scoring, PTEN expression was discordant between pretreatment and posttreatment primary tumors in 5% of patients (n = 96) and between posttreatment primary tumors and lymph node metastases in 9% (n = 33). A less stringent cutoff yielded similar discordance rates. Intratumoral heterogeneity for PTEN loss was observed in 7% of the patients. Among pretreatment tumors, PTEN copy numbers by whole exome sequencing (n = 72) were significantly higher in the PTEN-positive tumors by IHC compared with the IHC PTEN-loss tumors (p < 0.0001). However, PTEN-positive and PTEN-loss tumors by IHC overlapped in copy numbers: 14 of 60 PTEN-positive samples showed decreased copy numbers in the range of those of the PTEN-loss tumors. Conclusion: Testing various specimens by IHC may generate different PTEN results in a small proportion of patients with TNBC; therefore, the decision of testing one versus multiple specimens in a clinical trial should be defined in the patient inclusion criteria. Although a distinct cutoff by which CNV differentiated PTEN-positive tumors from those with PTEN loss was not identified, higher copy number of PTEN may confer positive PTEN, whereas lower copy number of PTEN would necessitate additional testing by IHC to assess PTEN loss. Trial registration: NCT02276443.

Prognostic Impact of High Baseline Stromal Tumor-Infiltrating Lymphocytes in the Absence of Pathologic Complete Response in Early-Stage Triple-Negative Breast Cancer.

High stromal tumor-infiltrating lymphocytes (sTILs) are associated with an improved pathologic complete response (pCR) and survival in triple-negative breast cancer (TNBC). We hypothesized that high baseline sTILs would have a favorable prognostic impact in TNBC patients without a pCR after neoadjuvant chemotherapy (NACT). In this prospective NACT study, pretreatment biopsies from 318 patients with early-stage TNBC were evaluated for sTILs. Recursive partitioning analysis (RPA) was applied to search for the sTIL cutoff best associated with a pCR. With ≥20% sTILs identified as the optimal cutoff, 33% patients had high sTILs (pCR rate 64%) and 67% had low sTILs (pCR rate 29%). Patients were stratified according to the sTIL cutoff (low vs. high) and response to NACT (pCR vs. residual disease (RD)). The primary endpoint was event-free survival (EFS), with hazard ratios calculated using the Cox proportional hazards regression model and the 3-year restricted mean survival time (RMST) as primary measures. Within the high-sTIL group, EFS was better in patients with a pCR compared with those with RD (HR 0.05; 95% CI 0.01-0.39; p = 0.004). The difference in the 3-year RMST for EFS between the two groups was 5.6 months (95% CI 2.3-8.8; p = 0.001). However, among patients with RD, EFS was not significantly different between those with high sTILs and those with low sTILs (p = 0.7). RNA-seq analysis predicted more CD8+ T cells in the high-sTIL group with favorable EFS compared with the high-sTIL group with unfavorable EFS. This study did not demonstrate that high baseline sTILs confer a benefit in EFS in the absence of a pCR.

A functional genomic approach to actionable gene fusions for precision oncology.

Fusion genes represent a class of attractive therapeutic targets. Thousands of fusion genes have been identified in patients with cancer, but the functional consequences and therapeutic implications of most of these remain largely unknown. Here, we develop a functional genomic approach that consists of efficient fusion reconstruction and sensitive cell viability and drug response assays. Applying this approach, we characterize ~100 fusion genes detected in patient samples of The Cancer Genome Atlas, revealing a notable fraction of low-frequency fusions with activating effects on tumor growth. Focusing on those in the RTK-RAS pathway, we identify a number of activating fusions that can markedly affect sensitivity to relevant drugs. Last, we propose an integrated, level-of-evidence classification system to prioritize gene fusions systematically. Our study reiterates the urgent clinical need to incorporate similar functional genomic approaches to characterize gene fusions, thereby maximizing the utility of gene fusions for precision oncology.

The cytosolic iron-sulfur cluster assembly (CIA) pathway is required for replication stress tolerance of cancer cells to Chk1 and ATR inhibitors.

The relationship between ATR/Chk1 activity and replication stress, coupled with the development of potent and tolerable inhibitors of this pathway, has led to the clinical exploration of ATR and Chk1 inhibitors (ATRi/Chk1i) as anticancer therapies for single-agent or combinatorial application. The clinical efficacy of these therapies relies on the ability to ascertain which patient populations are most likely to benefit, so there is intense interest in identifying predictive biomarkers of response. To comprehensively evaluate the components that modulate cancer cell sensitivity to replication stress induced by Chk1i, we performed a synthetic-lethal drop-out screen in a cell line derived from a patient with triple-negative breast cancer (TNBC), using a pooled barcoded shRNA library targeting ~350 genes involved in DNA replication, DNA damage repair, and cycle progression. In addition, we sought to compare the relative requirement of these genes when DNA fidelity is challenged by clinically relevant anticancer breast cancer drugs, including cisplatin and PARP1/2 inhibitors, that have different mechanisms of action. This global comparison is critical for understanding not only which agents should be used together for combinatorial therapies in breast cancer patients, but also the genetic context in which these therapies will be most effective, and when a single-agent therapy will be sufficient to provide maximum therapeutic benefit to the patient. We identified unique potentiators of response to ATRi/Chk1i and describe a new role for components of the cytosolic iron-sulfur assembly (CIA) pathway, MMS19 and CIA2B-FAM96B, in replication stress tolerance of TNBC.

Inhibition of histone acetyltransferase function radiosensitizes CREBBP/EP300 mutants via repression of homologous recombination, potentially targeting a gain of function.

Despite radiation forming the curative backbone of over 50% of malignancies, there are no genomically-driven radiosensitizers for clinical use. Herein we perform in vivo shRNA screening to identify targets generally associated with radiation response as well as those exhibiting a genomic dependency. This identifies the histone acetyltransferases CREBBP/EP300 as a target for radiosensitization in combination with radiation in cognate mutant tumors. Further in vitro and in vivo studies confirm this phenomenon to be due to repression of homologous recombination following DNA damage and reproducible using chemical inhibition of histone acetyltransferase (HAT), but not bromodomain function. Selected mutations in CREBBP lead to a hyperacetylated state that increases CBP and BRCA1 acetylation, representing a gain of function targeted by HAT inhibition. Additionally, mutations in CREBBP/EP300 are associated with recurrence following radiation in squamous cell carcinoma cohorts. These findings provide both a mechanism of resistance and the potential for genomically-driven treatment.

Immune Phenotype and Response to Neoadjuvant Therapy in Triple-Negative Breast Cancer.

PURPOSE: Increasing tumor-infiltrating lymphocytes (TIL) is associated with higher rates of pathologic complete response (pCR) to neoadjuvant therapy (NAT) in patients with triple-negative breast cancer (TNBC). However, the presence of TILs does not consistently predict pCR, therefore, the current study was undertaken to more fully characterize the immune cell response and its association with pCR. EXPERIMENTAL DESIGN: We obtained pretreatment core-needle biopsies from 105 patients with stage I-III TNBC enrolled in ARTEMIS (NCT02276443) who received NAT from Oct 22, 2015 through July 24, 2018. The tumor-immune microenvironment was comprehensively profiled by performing T-cell receptor (TCR) sequencing, programmed death-ligand 1 (PD-L1) IHC, multiplex immunofluorescence, and RNA sequencing on pretreatment tumor samples. The primary endpoint was pathologic response to NAT. RESULTS: The pCR rate was 40% (42/105). Higher TCR clonality (median = 0.2 vs. 0.1, P = 0.03), PD-L1 positivity (OR: 2.91, P = 0.020), higher CD3+:CD68+ ratio (median = 14.70 vs. 8.20, P = 0.0128), and closer spatial proximity of T cells to tumor cells (median = 19.26 vs. 21.94 µm, P = 0.0169) were associated with pCR. In a multivariable model, closer spatial proximity of T cells to tumor cells and PD-L1 expression enhanced prediction of pCR when considered in conjunction with clinical stage. CONCLUSIONS: In patients receiving NAT for TNBC, deep immune profiling through detailed phenotypic characterization and spatial analysis can improve prediction of pCR in patients receiving NAT for TNBC when considered with traditional clinical parameters.

Medium-Chain Acyl-CoA Dehydrogenase Protects Mitochondria from Lipid Peroxidation in Glioblastoma.

Glioblastoma (GBM) is highly resistant to chemotherapies, immune-based therapies, and targeted inhibitors. To identify novel drug targets, we screened orthotopically implanted, patient-derived glioblastoma sphere-forming cells using an RNAi library to probe essential tumor cell metabolic programs. This identified high dependence on mitochondrial fatty acid metabolism. We focused on medium-chain acyl-CoA dehydrogenase (MCAD), which oxidizes medium-chain fatty acids (MCFA), due to its consistently high score and high expression among models and upregulation in GBM compared with normal brain. Beyond the expected energetics impairment, MCAD depletion in primary GBM models induced an irreversible cascade of detrimental metabolic effects characterized by accumulation of unmetabolized MCFAs, which induced lipid peroxidation and oxidative stress, irreversible mitochondrial damage, and apoptosis. Our data uncover a novel protective role for MCAD to clear lipid molecules that may cause lethal cell damage, suggesting that therapeutic targeting of MCFA catabolism may exploit a key metabolic feature of GBM. SIGNIFICANCE: MCAD exerts a protective role to prevent accumulation of toxic metabolic by-products in glioma cells, actively catabolizing lipid species that would otherwise affect mitochondrial integrity and induce cell death. This work represents a first demonstration of a nonenergetic role for dependence on fatty acid metabolism in cancer.This article is highlighted in the In This Issue feature, p. 2659.

Sequential Administration of XPO1 and ATR Inhibitors Enhances Therapeutic Response in TP53-mutated Colorectal Cancer.

BACKGROUND & AIMS: Understanding the mechanisms by which tumors adapt to therapy is critical for developing effective combination therapeutic approaches to improve clinical outcomes for patients with cancer. METHODS: To identify promising and clinically actionable targets for managing colorectal cancer (CRC), we conducted a patient-centered functional genomics platform that includes approximately 200 genes and paired this with a high-throughput drug screen that includes 262 compounds in four patient-derived xenografts (PDXs) from patients with CRC. RESULTS: Both screening methods identified exportin 1 (XPO1) inhibitors as drivers of DNA damage-induced lethality in CRC. Molecular characterization of the cellular response to XPO1 inhibition uncovered an adaptive mechanism that limited the duration of response in TP53-mutated, but not in TP53-wild-type CRC models. Comprehensive proteomic and transcriptomic characterization revealed that the ATM/ATR-CHK1/2 axes were selectively engaged in TP53-mutant CRC cells upon XPO1 inhibitor treatment and that this response was required for adapting to therapy and escaping cell death. Administration of KPT-8602, an XPO1 inhibitor, followed by AZD-6738, an ATR inhibitor, resulted in dramatic antitumor effects and prolonged survival in TP53-mutant models of CRC. CONCLUSIONS: Our findings anticipate tremendous therapeutic benefit and support the further evaluation of XPO1 inhibitors, especially in combination with DNA damage checkpoint inhibitors, to elicit an enduring clinical response in patients with CRC harboring TP53 mutations.

Somatic mutation distributions in cancer genomes vary with three-dimensional chromatin structure.

Somatic mutations in driver genes may ultimately lead to the development of cancer. Understanding how somatic mutations accumulate in cancer genomes and the underlying factors that generate somatic mutations is therefore crucial for developing novel therapeutic strategies. To understand the interplay between spatial genome organization and specific mutational processes, we studied 3,000 tumor-normal-pair whole-genome datasets from 42 different human cancer types. Our analyses reveal that the change in somatic mutational load in cancer genomes is co-localized with topologically-associating-domain boundaries. Domain boundaries constitute a better proxy to track mutational load change than replication timing measurements. We show that different mutational processes lead to distinct somatic mutation distributions where certain processes generate mutations in active domains, and others generate mutations in inactive domains. Overall, the interplay between three-dimensional genome organization and active mutational processes has a substantial influence on the large-scale mutation-rate variations observed in human cancers.

Allosteric SHP2 Inhibitor, IACS-13909, Overcomes EGFR-Dependent and EGFR-Independent Resistance Mechanisms toward Osimertinib.

Src homology 2 domain-containing phosphatase (SHP2) is a phosphatase that mediates signaling downstream of multiple receptor tyrosine kinases (RTK) and is required for full activation of the MAPK pathway. SHP2 inhibition has demonstrated tumor growth inhibition in RTK-activated cancers in preclinical studies. The long-term effectiveness of tyrosine kinase inhibitors such as the EGFR inhibitor (EGFRi), osimertinib, in non-small cell lung cancer (NSCLC) is limited by acquired resistance. Multiple clinically identified mechanisms underlie resistance to osimertinib, including mutations in EGFR that preclude drug binding as well as EGFR-independent activation of the MAPK pathway through alternate RTK (RTK-bypass). It has also been noted that frequently a tumor from a single patient harbors more than one resistance mechanism, and the plasticity between multiple resistance mechanisms could restrict the effectiveness of therapies targeting a single node of the oncogenic signaling network. Here, we report the discovery of IACS-13909, a specific and potent allosteric inhibitor of SHP2, that suppresses signaling through the MAPK pathway. IACS-13909 potently impeded proliferation of tumors harboring a broad spectrum of activated RTKs as the oncogenic driver. In EGFR-mutant osimertinib-resistant NSCLC models with EGFR-dependent and EGFR-independent resistance mechanisms, IACS-13909, administered as a single agent or in combination with osimertinib, potently suppressed tumor cell proliferation in vitro and caused tumor regression in vivo. Together, our findings provide preclinical evidence for using a SHP2 inhibitor as a therapeutic strategy in acquired EGFRi-resistant NSCLC. SIGNIFICANCE: These findings highlight the discovery of IACS-13909 as a potent, selective inhibitor of SHP2 with drug-like properties, and targeting SHP2 may serve as a therapeutic strategy to overcome tumor resistance to osimertinib.

Global impact of somatic structural variation on the DNA methylome of human cancers.

BACKGROUND: Genomic rearrangements exert a heavy influence on the molecular landscape of cancer. New analytical approaches integrating somatic structural variants (SSVs) with altered gene features represent a framework by which we can assign global significance to a core set of genes, analogous to established methods that identify genes non-randomly targeted by somatic mutation or copy number alteration. While recent studies have defined broad patterns of association involving gene transcription and nearby SSV breakpoints, global alterations in DNA methylation in the context of SSVs remain largely unexplored. RESULTS: By data integration of whole genome sequencing, RNA sequencing, and DNA methylation arrays from more than 1400 human cancers, we identify hundreds of genes and associated CpG islands (CGIs) for which the nearby presence of a somatic structural variant (SSV) breakpoint is recurrently associated with altered expression or DNA methylation, respectively, independently of copy number alterations. CGIs with SSV-associated increased methylation are predominantly promoter-associated, while CGIs with SSV-associated decreased methylation are enriched for gene body CGIs. Rearrangement of genomic regions normally having higher or lower methylation is often involved in SSV-associated CGI methylation alterations. Across cancers, the overall structural variation burden is associated with a global decrease in methylation, increased expression in methyltransferase genes and DNA damage response genes, and decreased immune cell infiltration. CONCLUSION: Genomic rearrangement appears to have a major role in shaping the cancer DNA methylome, to be considered alongside commonly accepted mechanisms including histone modifications and disruption of DNA methyltransferases.

Resistance to neoadjuvant chemotherapy in triple-negative breast cancer mediated by a reversible drug-tolerant state.

Eradicating triple-negative breast cancer (TNBC) resistant to neoadjuvant chemotherapy (NACT) is a critical unmet clinical need. In this study, patient-derived xenograft (PDX) models of treatment-naïve TNBC and serial biopsies from TNBC patients undergoing NACT were used to elucidate mechanisms of chemoresistance in the neoadjuvant setting. Barcode-mediated clonal tracking and genomic sequencing of PDX tumors revealed that residual tumors remaining after treatment with standard frontline chemotherapies, doxorubicin (Adriamycin) combined with cyclophosphamide (AC), maintained the subclonal architecture of untreated tumors, yet their transcriptomes, proteomes, and histologic features were distinct from those of untreated tumors. Once treatment was halted, residual tumors gave rise to AC-sensitive tumors with similar transcriptomes, proteomes, and histological features to those of untreated tumors. Together, these results demonstrated that tumors can adopt a reversible drug-tolerant state that does not involve clonal selection as an AC resistance mechanism. Serial biopsies obtained from patients with TNBC undergoing NACT revealed similar histologic changes and maintenance of stable subclonal architecture, demonstrating that AC-treated PDXs capture molecular features characteristic of human TNBC chemoresistance. Last, pharmacologic inhibition of oxidative phosphorylation using an inhibitor currently in phase 1 clinical development delayed residual tumor regrowth. Thus, AC resistance in treatment-naïve TNBC can be mediated by nonselective mechanisms that confer a reversible chemotherapy-tolerant state with targetable vulnerabilities.

Syndecan 1 is a critical mediator of macropinocytosis in pancreatic cancer.

Pancreatic ductal adenocarcinoma (PDAC) remains recalcitrant to all forms of cancer treatment and carries a five-year survival rate of only 8%1. Inhibition of oncogenic KRAS (hereafter KRAS*), the earliest lesion in disease development that is present in more than 90% of PDACs, and its signalling surrogates has yielded encouraging preclinical results with experimental agents2-4. However, KRAS*-independent disease recurrence following genetic extinction of Kras* in mouse models anticipates the need for co-extinction strategies5,6. Multiple oncogenic processes are initiated at the cell surface, where KRAS* physically and functionally interacts to direct signalling that is essential for malignant transformation and tumour maintenance. Insights into the complexity of the functional cell-surface-protein repertoire (surfaceome) have been technologically limited until recently and-in the case of PDAC-the genetic control of the function and composition of the PDAC surfaceome in the context of KRAS* signalling remains largely unknown. Here we develop an unbiased, functional target-discovery platform to query KRAS*-dependent changes of the PDAC surfaceome, which reveals syndecan 1 (SDC1, also known as CD138) as a protein that is upregulated at the cell surface by KRAS*. Localization of SDC1 at the cell surface-where it regulates macropinocytosis, an essential metabolic pathway that fuels PDAC cell growth-is essential for disease maintenance and progression. Thus, our study forges a mechanistic link between KRAS* signalling and a targetable molecule driving nutrient salvage pathways in PDAC and validates oncogene-driven surfaceome annotation as a strategy to identify cancer-specific vulnerabilities.

Pre-existing Functional Heterogeneity of Tumorigenic Compartment as the Origin of Chemoresistance in Pancreatic Tumors.

Adaptive drug-resistance mechanisms allow human tumors to evade treatment through selection and expansion of treatment-resistant clones. Here, studying clonal evolution of tumor cells derived from human pancreatic tumors, we demonstrate that in vitro cultures and in vivo tumors are maintained by a common set of tumorigenic cells that can be used to establish clonal replica tumors (CRTs), large cohorts of animals bearing human tumors with identical clonal composition. Using CRTs to conduct quantitative assessments of adaptive responses to therapeutics, we uncovered a multitude of functionally heterogeneous subpopulations of cells with differential degrees of drug sensitivity. High-throughput isolation and deep characterization of unique clonal lineages showed genetic and transcriptomic diversity underlying functionally diverse subpopulations. Molecular annotation of gemcitabine-naive clonal lineages with distinct responses to treatment in the context of CRTs generated signatures that can predict the response to chemotherapy, representing a potential biomarker to stratify patients with pancreatic cancer.

Functional Genomics Reveals Synthetic Lethality between Phosphogluconate Dehydrogenase and Oxidative Phosphorylation.

The plasticity of a preexisting regulatory circuit compromises the effectiveness of targeted therapies, and leveraging genetic vulnerabilities in cancer cells may overcome such adaptations. Hereditary leiomyomatosis renal cell carcinoma (HLRCC) is characterized by oxidative phosphorylation (OXPHOS) deficiency caused by fumarate hydratase (FH) nullizyogosity. To identify metabolic genes that are synthetically lethal with OXPHOS deficiency, we conducted a genetic loss-of-function screen and found that phosphogluconate dehydrogenase (PGD) inhibition robustly blocks the proliferation of FH mutant cancer cells both in vitro and in vivo. Mechanistically, PGD inhibition blocks glycolysis, suppresses reductive carboxylation of glutamine, and increases the NADP+/NADPH ratio to disrupt redox homeostasis. Furthermore, in the OXPHOS-proficient context, blocking OXPHOS using the small-molecule inhibitor IACS-010759 enhances sensitivity to PGD inhibition in vitro and in vivo. Together, our study reveals a dependency on PGD in OXPHOS-deficient tumors that might inform therapeutic intervention in specific patient populations.

High-resolution clonal mapping of multi-organ metastasis in triple negative breast cancer.

Most triple negative breast cancers (TNBCs) are aggressively metastatic with a high degree of intra-tumoral heterogeneity (ITH), but how ITH contributes to metastasis is unclear. Here, clonal dynamics during metastasis were studied in vivo using two patient-derived xenograft (PDX) models established from the treatment-naive primary breast tumors of TNBC patients diagnosed with synchronous metastasis. Genomic sequencing and high-complexity barcode-mediated clonal tracking reveal robust alterations in clonal architecture between primary tumors and corresponding metastases. Polyclonal seeding and maintenance of heterogeneous populations of low-abundance subclones is observed in each metastasis. However, lung, liver, and brain metastases are enriched for an identical population of high-abundance subclones, demonstrating that primary tumor clones harbor properties enabling them to seed and thrive in multiple organ sites. Further, clones that dominate multi-organ metastases share a genomic lineage. Thus, intrinsic properties of rare primary tumor subclones enable the seeding and colonization of metastases in secondary organs in these models.

A Pan-Cancer Compendium of Genes Deregulated by Somatic Genomic Rearrangement across More Than 1,400 Cases.

A systematic cataloging of genes affected by genomic rearrangement, using multiple patient cohorts and cancer types, can provide insight into cancer-relevant alterations outside of exomes. By integrative analysis of whole-genome sequencing (predominantly low pass) and gene expression data from 1,448 cancers involving 18 histopathological types in The Cancer Genome Atlas, we identified hundreds of genes for which the nearby presence (within 100 kb) of a somatic structural variant (SV) breakpoint is associated with altered expression. While genomic rearrangements are associated with widespread copy-number alteration (CNA) patterns, approximately 1,100 genes-including overexpressed cancer driver genes (e.g., TERT, ERBB2, CDK12, CDK4) and underexpressed tumor suppressors (e.g., TP53, RB1, PTEN, STK11)-show SV-associated deregulation independent of CNA. SVs associated with the disruption of topologically associated domains, enhancer hijacking, or fusion transcripts are implicated in gene upregulation. For cancer-relevant pathways, SVs considerably expand our understanding of how genes are affected beyond point mutation or CNA.

Discordant inheritance of chromosomal and extrachromosomal DNA elements contributes to dynamic disease evolution in glioblastoma.

To understand how genomic heterogeneity of glioblastoma (GBM) contributes to poor therapy response, we performed DNA and RNA sequencing on GBM samples and the neurospheres and orthotopic xenograft models derived from them. We used the resulting dataset to show that somatic driver alterations including single-nucleotide variants, focal DNA alterations and oncogene amplification on extrachromosomal DNA (ecDNA) elements were in majority propagated from tumor to model systems. In several instances, ecDNAs and chromosomal alterations demonstrated divergent inheritance patterns and clonal selection dynamics during cell culture and xenografting. We infer that ecDNA was unevenly inherited by offspring cells, a characteristic that affects the oncogenic potential of cells with more or fewer ecDNAs. Longitudinal patient tumor profiling found that oncogenic ecDNAs are frequently retained throughout the course of disease. Our analysis shows that extrachromosomal elements allow rapid increase of genomic heterogeneity during GBM evolution, independently of chromosomal DNA alterations.