We have profiled promoter DNA methylation alterations in 272 glioblastoma tumors in the context of The Cancer Genome Atlas (TCGA). We found that a distinct subset of samples displays concerted hypermethylation at a large number of loci, indicating the existence of a glioma-CpG island methylator phenotype (G-CIMP). We validated G-CIMP in a set of non-TCGA glioblastomas and low-grade gliomas. G-CIMP tumors belong to the proneural subgroup, are more prevalent among lower-grade gliomas, display distinct copy-number alterations, and are tightly associated with IDH1 somatic mutations. Patients with G-CIMP tumors are younger at the time of diagnosis and experience significantly improved outcome. These findings identify G-CIMP as a distinct subset of human gliomas on molecular and clinical grounds.

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.

The molecular biology involving neoadjuvant chemotherapy (NAC) response is poorly understood. To elucidate the impact of NAC on the breast cancer transcriptome and its association with clinical outcome, we analyzed gene expression data derived from serial tumor samples of patients with breast cancer who received NAC in the I-SPY 1 TRIAL.

Partial loss of large ribosomal subunit protein 24 (RPL24) function is known to protect mice against Akt or Myc-driven cancers, in part via translational inhibition of a subset of cap(eIF4E)-dependently translated mRNAs. The role of RPL24 in human malignancies is unknown. By analyzing a public dataset of matched human breast cancers and normal mammary tissue, we found that breast cancers express significantly more RPL24 than matched normal breast samples. Depletion of RPL24 in breast cancer cells by >70% reduced cell viability by 80% and decreased protein expression of the eIF4E-dependently translated proteins cyclin D1 (75%), survivin (46%) and NBS1 (30%) without altering GAPDH or beta-tubulin levels. RPL24 knockdown also reduced 80S subunit levels relative to 40S and 60S levels. These effects on expression of eIF4E-dependent proteins and ribosome assembly were mimicked by 2-24 h treatment with the pan-HDACi, trichostatin A (TSA), which induced acetylation of 15 different polysome-associated proteins including RPL24. Furthermore, HDAC6-selective inhibition or HDAC6 knockdown induced ribosomal protein acetylation. Via mass spectrometry, we found that 60S-associated, but not, polysome-associated, RPL24 undergoes HDACi-induced acetylation on K27. Thus, RPL24 K27 acetylation may play a role in ribosome assembly. These findings point toward a novel acetylation-dependent polysome assembly mechanism regulating tumorigenesis.

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.

In human prostate cancer, the microRNA biogenesis machinery increases with prostate cancer progression. Here, we show that deletion of the Dgcr8 gene, a critical component of this complex, inhibits tumor progression in a Pten-knockout mouse model of prostate cancer. Early stages of tumor development were unaffected, but progression to advanced prostatic intraepithelial neoplasia was severely inhibited. Dgcr8 loss blocked Pten null-induced expansion of the basal-like, but not luminal, cellular compartment. Furthermore, while late-stage Pten knockout tumors exhibit decreased senescence-associated beta-galactosidase activity and increased proliferation, the simultaneous deletion of Dgcr8 blocked these changes resulting in levels similar to wild type. Sequencing of small RNAs in isolated epithelial cells uncovered numerous miRNA changes associated with PTEN loss. Consistent with a Pten-Dgcr8 association, analysis of a large cohort of human prostate tumors shows a strong correlation between Akt activation and increased Dgcr8 mRNA levels. Together, these findings uncover a critical role for microRNAs in enhancing proliferation and enabling the expansion of the basal cell compartment associated with tumor progression following Pten loss.

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.

Age is one of the most important risk factors for human malignancies, including breast cancer; in addition, age at diagnosis has been shown to be an independent indicator of breast cancer prognosis. Except for inherited forms of breast cancer, however, there is little genetic or epigenetic understanding of the biological basis linking aging with sporadic breast cancer incidence and its clinical behavior.

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.

Cancer progression involves the gradual loss of a differentiated phenotype and acquisition of progenitor and stem-cell-like features. Here, we provide novel stemness indices for assessing the degree of oncogenic dedifferentiation. We used an innovative one-class logistic regression (OCLR) machine-learning algorithm to extract transcriptomic and epigenetic feature sets derived from non-transformed pluripotent stem cells and their differentiated progeny. Using OCLR, we were able to identify previously undiscovered biological mechanisms associated with the dedifferentiated oncogenic state. Analyses of the tumor microenvironment revealed unanticipated correlation of cancer stemness with immune checkpoint expression and infiltrating immune cells. We found that the dedifferentiated oncogenic phenotype was generally most prominent in metastatic tumors. Application of our stemness indices to single-cell data revealed patterns of intra-tumor molecular heterogeneity. Finally, the indices allowed for the identification of novel targets and possible targeted therapies aimed at tumor differentiation.

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.

A small fraction of cancer patients with advanced disease survive significantly longer than patients with clinically comparable tumors. Molecular mechanisms for exceptional responses to therapy have been identified by genomic analysis of tumor biopsies from individual patients. Here, we analyzed tumor biopsies from an unbiased cohort of 111 exceptional responder patients using multiple platforms to profile genetic and epigenetic aberrations as well as the tumor microenvironment. Integrative analysis uncovered plausible mechanisms for the therapeutic response in nearly a quarter of the patients. The mechanisms were assigned to four broad categories-DNA damage response, intracellular signaling, immune engagement, and genetic alterations characteristic of favorable prognosis-with many tumors falling into multiple categories. These analyses revealed synthetic lethal relationships that may be exploited therapeutically and rare genetic lesions that favor therapeutic success, while also providing a wealth of testable hypotheses regarding oncogenic mechanisms that may influence the response to cancer therapy.

Muscle-invasive bladder cancer (MIBC) is a molecularly diverse disease with heterogeneous clinical outcomes. Several molecular classifications have been proposed, but the diversity of their subtype sets impedes their clinical application.

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.

Transcript alterations often result from somatic changes in cancer genomes1. Various forms of RNA alterations have been described in cancer, including overexpression2, altered splicing3 and gene fusions4; however, it is difficult to attribute these to underlying genomic changes owing to heterogeneity among patients and tumour types, and the relatively small cohorts of patients for whom samples have been analysed by both transcriptome and whole-genome sequencing. Here we present, to our knowledge, the most comprehensive catalogue of cancer-associated gene alterations to date, obtained by characterizing tumour transcriptomes from 1,188 donors of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA)5. Using matched whole-genome sequencing data, we associated several categories of RNA alterations with germline and somatic DNA alterations, and identified probable genetic mechanisms. Somatic copy-number alterations were the major drivers of variations in total gene and allele-specific expression. We identified 649 associations of somatic single-nucleotide variants with gene expression in cis, of which 68.4% involved associations with flanking non-coding regions of the gene. We found 1,900 splicing alterations associated with somatic mutations, including the formation of exons within introns in proximity to Alu elements. In addition, 82% of gene fusions were associated with structural variants, including 75 of a new class, termed 'bridged' fusions, in which a third genomic location bridges two genes. We observed transcriptomic alteration signatures that differ between cancer types and have associations with variations in DNA mutational signatures. This compendium of RNA alterations in the genomic context provides a rich resource for identifying genes and mechanisms that are functionally implicated in cancer.

Improved risk stratification and predictive biomarkers of treatment response are needed for non-muscle-invasive bladder cancer (NMIBC). Here we assessed the clinical utility of targeted RNA and DNA molecular profiling in NMIBC.

HER2-targeted therapy dramatically improves outcomes in early breast cancer. Here we report the results of two HER2-targeted combinations in the neoadjuvant I-SPY2 phase 2 adaptive platform trial for early breast cancer at high risk of recurrence: ado-trastuzumab emtansine plus pertuzumab (T-DM1/P) and paclitaxel, trastuzumab and pertuzumab (THP). Eligible women have >2.5 cm clinical stage II/III HER2+ breast cancer, adaptively randomized to T-DM1/P, THP, or a common control arm of paclitaxel/trastuzumab (TH), followed by doxorubicin/cyclophosphamide, then surgery. Both T-DM1/P and THP arms 'graduate' in all subtypes: predicted pCR rates are 63%, 72% and 33% for T-DM1/P (n = 52), THP (n = 45) and TH (n = 31) respectively. Toxicity burden is similar between arms. Degree of HER2 pathway signaling and phosphorylation in pretreatment biopsy specimens are associated with response to both T-DM1/P and THP and can further identify highly responsive HER2+ tumors to HER2-directed therapy. This may help identify patients who can safely de-escalate cytotoxic chemotherapy without compromising excellent outcome.

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).