Glioblastoma multiforme (GBM) is a lethal, therapy-resistant brain cancer consisting of numerous tumor cell subpopulations, including stem-like glioma-initiating cells (GICs), which contribute to tumor recurrence following initial response to therapy. Here, we identified miR-182 as a regulator of apoptosis, growth, and differentiation programs whose expression level is correlated with GBM patient survival. Repression of Bcl2-like12 (Bcl2L12), c-Met, and hypoxia-inducible factor 2α (HIF2A) is of central importance to miR-182 anti-tumor activity, as it results in enhanced therapy susceptibility, decreased GIC sphere size, expansion, and stemness in vitro. To evaluate the tumor-suppressive function of miR-182 in vivo, we synthesized miR-182-based spherical nucleic acids (182-SNAs); i.e., gold nanoparticles covalently functionalized with mature miR-182 duplexes. Intravenously administered 182-SNAs penetrated the blood-brain/blood-tumor barriers (BBB/BTB) in orthotopic GBM xenografts and selectively disseminated throughout extravascular glioma parenchyma, causing reduced tumor burden and increased animal survival. Our results indicate that harnessing the anti-tumor activities of miR-182 via safe and robust delivery of 182-SNAs represents a novel strategy for therapeutic intervention in GBM.

To assess telomerase as a cancer therapeutic target and determine adaptive mechanisms to telomerase inhibition, we modeled telomerase reactivation and subsequent extinction in T cell lymphomas arising in Atm(-/-) mice engineered with an inducible telomerase reverse transcriptase allele. Telomerase reactivation in the setting of telomere dysfunction enabled full malignant progression with alleviation of telomere dysfunction-induced checkpoints. These cancers possessed copy number alterations targeting key loci in human T cell lymphomagenesis. Upon telomerase extinction, tumor growth eventually slowed with reinstatement of telomere dysfunction-induced checkpoints, yet growth subsequently resumed as tumors acquired alternative lengthening of telomeres (ALT) and aberrant transcriptional networks centering on mitochondrial biology and oxidative defense. ALT+ tumors acquired amplification/overexpression of PGC-1β, a master regulator of mitochondrial biogenesis and function, and they showed marked sensitivity to PGC-1β or SOD2 knockdown. Genetic modeling of telomerase extinction reveals vulnerabilities that motivate coincidental inhibition of mitochondrial maintenance and oxidative defense mechanisms to enhance antitelomerase cancer therapy.

We describe the landscape of somatic genomic alterations of 66 chromophobe renal cell carcinomas (ChRCCs) on the basis of multidimensional and comprehensive characterization, including mtDNA and whole-genome sequencing. The result is consistent that ChRCC originates from the distal nephron compared with other kidney cancers with more proximal origins. Combined mtDNA and gene expression analysis implicates changes in mitochondrial function as a component of the disease biology, while suggesting alternative roles for mtDNA mutations in cancers relying on oxidative phosphorylation. Genomic rearrangements lead to recurrent structural breakpoints within TERT promoter region, which correlates with highly elevated TERT expression and manifestation of kataegis, representing a mechanism of TERT upregulation in cancer distinct from previously observed amplifications and point mutations.

Comprehensive multiplatform analysis of 80 uveal melanomas (UM) identifies four molecularly distinct, clinically relevant subtypes: two associated with poor-prognosis monosomy 3 (M3) and two with better-prognosis disomy 3 (D3). We show that BAP1 loss follows M3 occurrence and correlates with a global DNA methylation state that is distinct from D3-UM. Poor-prognosis M3-UM divide into subsets with divergent genomic aberrations, transcriptional features, and clinical outcomes. We report change-of-function SRSF2 mutations. Within D3-UM, EIF1AX- and SRSF2/SF3B1-mutant tumors have distinct somatic copy number alterations and DNA methylation profiles, providing insight into the biology of these low- versus intermediate-risk clinical mutation subtypes.

Aneuploidy, whole chromosome or chromosome arm imbalance, is a near-universal characteristic of human cancers. In 10,522 cancer genomes from The Cancer Genome Atlas, aneuploidy was correlated with TP53 mutation, somatic mutation rate, and expression of proliferation genes. Aneuploidy was anti-correlated with expression of immune signaling genes, due to decreased leukocyte infiltrates in high-aneuploidy samples. Chromosome arm-level alterations show cancer-specific patterns, including loss of chromosome arm 3p in squamous cancers. We applied genome engineering to delete 3p in lung cells, causing decreased proliferation rescued in part by chromosome 3 duplication. This study defines genomic and phenotypic correlates of cancer aneuploidy and provides an experimental approach to study chromosome arm aneuploidy.

We conducted comprehensive integrative molecular analyses of the complete set of tumors in The Cancer Genome Atlas (TCGA), consisting of approximately 10,000 specimens and representing 33 types of cancer. We performed molecular clustering using data on chromosome-arm-level aneuploidy, DNA hypermethylation, mRNA, and miRNA expression levels and reverse-phase protein arrays, of which all, except for aneuploidy, revealed clustering primarily organized by histology, tissue type, or anatomic origin. The influence of cell type was evident in DNA-methylation-based clustering, even after excluding sites with known preexisting tissue-type-specific methylation. Integrative clustering further emphasized the dominant role of cell-of-origin patterns. Molecular similarities among histologically or anatomically related cancer types provide a basis for focused pan-cancer analyses, such as pan-gastrointestinal, pan-gynecological, pan-kidney, and pan-squamous cancers, and those related by stemness features, which in turn may inform strategies for future therapeutic development.

Genetic alterations in signaling pathways that control cell-cycle progression, apoptosis, and cell growth are common hallmarks of cancer, but the extent, mechanisms, and co-occurrence of alterations in these pathways differ between individual tumors and tumor types. Using mutations, copy-number changes, mRNA expression, gene fusions and DNA methylation in 9,125 tumors profiled by The Cancer Genome Atlas (TCGA), we analyzed the mechanisms and patterns of somatic alterations in ten canonical pathways: cell cycle, Hippo, Myc, Notch, Nrf2, PI-3-Kinase/Akt, RTK-RAS, TGFβ signaling, p53 and β-catenin/Wnt. We charted the detailed landscape of pathway alterations in 33 cancer types, stratified into 64 subtypes, and identified patterns of co-occurrence and mutual exclusivity. Eighty-nine percent of tumors had at least one driver alteration in these pathways, and 57% percent of tumors had at least one alteration potentially targetable by currently available drugs. Thirty percent of tumors had multiple targetable alterations, indicating opportunities for combination therapy.

People of different ancestries vary in cancer risk and outcome, and their molecular differences may indicate sources of these variations. Determining the "local" ancestry composition at each genetic locus across ancestry-admixed populations can suggest causal associations. We present a protocol to identify local ancestry and detect the associated molecular changes, using data from the Cancer Genome Atlas. This workflow can be applied to cancer cohorts with matched tumor and normal data from admixed patients to examine germline contributions to cancer. For complete details on the use and execution of this protocol, please refer to Carrot-Zhang et al. (2020).

Eukaryotic messenger RNA (mRNA) typically contains a methylated guanosine (m7G) cap, which mediates major steps of mRNA metabolism. Recently, some RNAs in both prokaryotic and eukaryotic organisms have been found to carry a non-canonical cap such as the NAD cap. Here we report that Arabidopsis DXO family protein AtDXO1, which was previously known to be a decapping enzyme for NAD-capped RNAs (NAD-RNA), is an essential component for m7G capping. AtDXO1 associates with and activates RNA guanosine-7 methyltransferase (AtRNMT1) to catalyze conversion of the guanosine cap to the m7G cap. AtRNMT1 is an essential gene. Partial loss-of-function mutations of AtRNMT1 and knockout mutation of AtDXO1 reduce m7G-capped mRNA but increase G-capped mRNAs, leading to similar pleiotropic phenotypes, whereas overexpression of AtRNMT1 partially restores the atdxo1 phenotypes. This work reveals an important mechanism in m7G capping in plants by which the NAD-RNA decapping enzyme AtDXO1 is required for efficient guanosine cap methylation.

A high tumour mutational burden (hypermutation) is observed in some gliomas1-5; however, the mechanisms by which hypermutation develops and whether it predicts the response to immunotherapy are poorly understood. Here we comprehensively analyse the molecular determinants of mutational burden and signatures in 10,294 gliomas. We delineate two main pathways to hypermutation: a de novo pathway associated with constitutional defects in DNA polymerase and mismatch repair (MMR) genes, and a more common post-treatment pathway, associated with acquired resistance driven by MMR defects in chemotherapy-sensitive gliomas that recur after treatment with the chemotherapy drug temozolomide. Experimentally, the mutational signature of post-treatment hypermutated gliomas was recapitulated by temozolomide-induced damage in cells with MMR deficiency. MMR-deficient gliomas were characterized by a lack of prominent T cell infiltrates, extensive intratumoral heterogeneity, poor patient survival and a low rate of response to PD-1 blockade. Moreover, although bulk analyses did not detect microsatellite instability in MMR-deficient gliomas, single-cell whole-genome sequencing analysis of post-treatment hypermutated glioma cells identified microsatellite mutations. These results show that chemotherapy can drive the acquisition of hypermutated populations without promoting a response to PD-1 blockade and supports the diagnostic use of mutational burden and signatures in cancer.

Approximately 8% to 10% of pancreatic ductal adenocarcinomas (PDAC) do not harbor mutations in KRAS. Understanding the unique molecular and clinical features of this subset of pancreatic cancer is important to guide patient stratification for clinical trials of molecularly targeted agents.

The molecular underpinnings of HER2-low and HER2-0 (IHC 0) breast tumors remain poorly defined. Using genomic findings from 1039 patients with HER2-negative metastatic breast cancer undergoing next-generation sequencing from 7/2013-12/2020, we compare results between HER2-low (n = 487, 47%) and HER2-0 tumors (n = 552, 53%). A significantly higher number of ERBB2 alleles (median copy count: 2.05) are observed among HER2-low tumors compared to HER2-0 (median copy count: 1.79; P = 2.36e-6), with HER2-0 tumors harboring a higher rate of ERBB2 hemideletions (31.1% vs. 14.5%). No other genomic alteration reaches significance after accounting for multiple hypothesis testing, and no significant differences in tumor mutational burden are observed between HER2-low and HER2-0 tumors (median: 7.26 mutations/megabase vs. 7.60 mutations/megabase, p = 0.24). Here, we show that the genomic landscape of HER2-low and HER2-0 tumors does not differ significantly, apart from a higher ERBB2 copy count among HER2-low tumors, and a higher rate of ERBB2 hemideletions in HER2-0 tumors.

Obesity is linked to increased incidence of endometrioid endometrial cancer (EEC) and complex atypical hyperplasia (CAH). We here explore pattern and sequence of molecular alterations characterizing endometrial carcinogenesis in general and related to body mass index (BMI), to improve diagnostic stratification and treatment strategies. We performed molecular characterization of 729 prospectively collected EEC and CAH. Candidate biomarkers were identified in frozen samples by whole-exome and Sanger sequencing, oligonucleotide gene expression and Reverse Phase Protein Arrays (investigation cohort) and further explored in formalin fixed tissues by immunohistochemistry and Fluorescent in Situ Hybridization (validation cohort). We here demonstrate that PIK3CA mutations, PTEN loss, PI3K and KRAS activation are early events in endometrial carcinogenesis. Molecular changes related to KRAS activation and inflammation are more common in obese CAH patients, suggesting different prevention and systemic treatment strategies in obese and non-obese patients. We also found that oncoprotein Stathmin might improve preoperative diagnostic distinction between premalignant and malignant endometrial lesions.

Alterations in non-driver genes represent an emerging class of potential therapeutic targets in cancer. Hundreds to thousands of non-driver genes undergo loss of heterozygosity (LOH) events per tumor, generating discrete differences between tumor and normal cells. Here we interrogate LOH of polymorphisms in essential genes as a novel class of therapeutic targets. We hypothesized that monoallelic inactivation of the allele retained in tumors can selectively kill cancer cells but not somatic cells, which retain both alleles. We identified 5664 variants in 1278 essential genes that undergo LOH in cancer and evaluated the potential for each to be targeted using allele-specific gene-editing, RNAi, or small-molecule approaches. We further show that allele-specific inactivation of either of two essential genes (PRIM1 and EXOSC8) reduces growth of cells harboring that allele, while cells harboring the non-targeted allele remain intact. We conclude that LOH of essential genes represents a rich class of non-driver cancer vulnerabilities.

Cellular and molecular aberrations contribute to the disparity of human cancer incidence and etiology between ancestry groups. Multiomics profiling in The Cancer Genome Atlas (TCGA) allows for querying of the molecular underpinnings of ancestry-specific discrepancies in human cancer. Here, we provide a protocol for integrative associative analysis of ancestry with molecular correlates, including somatic mutations, DNA methylation, mRNA transcription, miRNA transcription, and pathway activity, using TCGA data. This protocol can be generalized to analyze other cancer cohorts and human diseases. For complete details on the use and execution of this protocol, please refer to Carrot-Zhang et al. (2020).

Several risk factors have been established for colorectal cancer, yet their direct mutagenic effects in patients' tumors remain to be elucidated. Here, we leveraged whole-exome sequencing data from 900 colorectal cancer cases that had occurred in three U.S.-wide prospective studies with extensive dietary and lifestyle information. We found an alkylating signature that was previously undescribed in colorectal cancer and then showed the existence of a similar mutational process in normal colonic crypts. This alkylating signature is associated with high intakes of processed and unprocessed red meat prior to diagnosis. In addition, this signature was more abundant in the distal colorectum, predicted to target cancer driver mutations KRAS p.G12D, KRAS p.G13D, and PIK3CA p.E545K, and associated with poor survival. Together, these results link for the first time a colorectal mutational signature to a component of diet and further implicate the role of red meat in colorectal cancer initiation and progression. SIGNIFICANCE: Colorectal cancer has several lifestyle risk factors, but the underlying mutations for most have not been observed directly in tumors. Analysis of 900 colorectal cancers with whole-exome sequencing and epidemiologic annotations revealed an alkylating mutational signature that was associated with red meat consumption and distal tumor location, as well as predicted to target KRAS p.G12D/p.G13D.This article is highlighted in the In This Issue feature, p. 2355.

Emerging data suggest that induction of viral mimicry responses through activation of double-stranded RNA (dsRNA) sensors in cancer cells is a promising therapeutic strategy. One approach to induce viral mimicry is to target molecular regulators of dsRNA sensing pathways. Here, we show that the exoribonuclease XRN1 is a negative regulator of the dsRNA sensor protein kinase R (PKR) in cancer cells with high interferon-stimulated gene (ISG) expression. XRN1 deletion causes PKR activation and consequent cancer cell lethality. Disruption of interferon signaling with the JAK1/2 inhibitor ruxolitinib can decrease cellular PKR levels and rescue sensitivity to XRN1 deletion. Conversely, interferon-β stimulation can increase PKR levels and induce sensitivity to XRN1 inactivation. Lastly, XRN1 deletion causes accumulation of endogenous complementary sense/anti-sense RNAs, which may represent candidate PKR ligands. Our data demonstrate how XRN1 regulates PKR and nominate XRN1 as a potential therapeutic target in cancer cells with an activated interferon cell state.

RTK/RAS/RAF pathway alterations (RPAs) are a hallmark of lung adenocarcinoma (LUAD). In this study, we use whole-genome sequencing (WGS) of 85 cases found to be RPA(-) by previous studies from The Cancer Genome Atlas (TCGA) to characterize the minority of LUADs lacking apparent alterations in this pathway. We show that WGS analysis uncovers RPA(+) in 28 (33%) of the 85 samples. Among the remaining 57 cases, we observe focal deletions targeting the promoter or transcription start site of STK11 (n = 7) or KEAP1 (n = 3), and promoter mutations associated with the increased expression of ILF2 (n = 6). We also identify complex structural variations associated with high-level copy number amplifications. Moreover, an enrichment of focal deletions is found in TP53 mutant cases. Our results indicate that RPA(-) cases demonstrate tumor suppressor deletions and genome instability, but lack unique or recurrent genetic lesions compensating for the lack of RPAs. Larger WGS studies of RPA(-) cases are required to understand this important LUAD subset.

We analyzed 921 adenocarcinomas of the esophagus, stomach, colon, and rectum to examine shared and distinguishing molecular characteristics of gastrointestinal tract adenocarcinomas (GIACs). Hypermutated tumors were distinct regardless of cancer type and comprised those enriched for insertions/deletions, representing microsatellite instability cases with epigenetic silencing of MLH1 in the context of CpG island methylator phenotype, plus tumors with elevated single-nucleotide variants associated with mutations in POLE. Tumors with chromosomal instability were diverse, with gastroesophageal adenocarcinomas harboring fragmented genomes associated with genomic doubling and distinct mutational signatures. We identified a group of tumors in the colon and rectum lacking hypermutation and aneuploidy termed genome stable and enriched in DNA hypermethylation and mutations in KRAS, SOX9, and PCBP1.

The Cancer Genome Atlas (TCGA) has catalyzed systematic characterization of diverse genomic alterations underlying human cancers. At this historic junction marking the completion of genomic characterization of over 11,000 tumors from 33 cancer types, we present our current understanding of the molecular processes governing oncogenesis. We illustrate our insights into cancer through synthesis of the findings of the TCGA PanCancer Atlas project on three facets of oncogenesis: (1) somatic driver mutations, germline pathogenic variants, and their interactions in the tumor; (2) the influence of the tumor genome and epigenome on transcriptome and proteome; and (3) the relationship between tumor and the microenvironment, including implications for drugs targeting driver events and immunotherapies. These results will anchor future characterization of rare and common tumor types, primary and relapsed tumors, and cancers across ancestry groups and will guide the deployment of clinical genomic sequencing.