Glioma progression is shaped by genetic evolution and microenvironment interactions.

The factors driving therapy resistance in diffuse glioma remain poorly understood. To identify treatment-associated cellular and genetic changes, we analyzed RNA and/or DNA sequencing data from the temporally separated tumor pairs of 304 adult patients with isocitrate dehydrogenase (IDH)-wild-type and IDH-mutant glioma. Tumors recurred in distinct manners that were dependent on IDH mutation status and attributable to changes in histological feature composition, somatic alterations, and microenvironment interactions. Hypermutation and acquired CDKN2A deletions were associated with an increase in proliferating neoplastic cells at recurrence in both glioma subtypes, reflecting active tumor growth. IDH-wild-type tumors were more invasive at recurrence, and their neoplastic cells exhibited increased expression of neuronal signaling programs that reflected a possible role for neuronal interactions in promoting glioma progression. Mesenchymal transition was associated with the presence of a myeloid cell state defined by specific ligand-receptor interactions with neoplastic cells. Collectively, these recurrence-associated phenotypes represent potential targets to alter disease progression.

Molecular Profiling Reveals Biologically Discrete Subsets and Pathways of Progression in Diffuse Glioma.

Therapy development for adult diffuse glioma is hindered by incomplete knowledge of somatic glioma driving alterations and suboptimal disease classification. We defined the complete set of genes associated with 1,122 diffuse grade II-III-IV gliomas from The Cancer Genome Atlas and used molecular profiles to improve disease classification, identify molecular correlations, and provide insights into the progression from low- to high-grade disease. Whole-genome sequencing data analysis determined that ATRX but not TERT promoter mutations are associated with increased telomere length. Recent advances in glioma classification based on IDH mutation and 1p/19q co-deletion status were recapitulated through analysis of DNA methylation profiles, which identified clinically relevant molecular subsets. A subtype of IDH mutant glioma was associated with DNA demethylation and poor outcome; a group of IDH-wild-type diffuse glioma showed molecular similarity to pilocytic astrocytoma and relatively favorable survival. Understanding of cohesive disease groups may aid improved clinical outcomes.

Silent mutations make some noise.

Synonymous mutations do not alter amino acids and are generally considered nonfunctional in cancer. Supek et al. now present a compelling analysis suggesting that such silent mutations can be oncogenic by altering transcript splicing and thereby affecting protein function.

Extrachromosomal DNA in the cancerous transformation of Barrett's oesophagus.

Oncogene amplification on extrachromosomal DNA (ecDNA) drives the evolution of tumours and their resistance to treatment, and is associated with poor outcomes for patients with cancer1-6. At present, it is unclear whether ecDNA is a later manifestation of genomic instability, or whether it can be an early event in the transition from dysplasia to cancer. Here, to better understand the development of ecDNA, we analysed whole-genome sequencing (WGS) data from patients with oesophageal adenocarcinoma (EAC) or Barrett's oesophagus. These data included 206 biopsies in Barrett's oesophagus surveillance and EAC cohorts from Cambridge University. We also analysed WGS and histology data from biopsies that were collected across multiple regions at 2 time points from 80 patients in a case-control study at the Fred Hutchinson Cancer Center. In the Cambridge cohorts, the frequency of ecDNA increased between Barrett's-oesophagus-associated early-stage (24%) and late-stage (43%) EAC, suggesting that ecDNA is formed during cancer progression. In the cohort from the Fred Hutchinson Cancer Center, 33% of patients who developed EAC had at least one oesophageal biopsy with ecDNA before or at the diagnosis of EAC. In biopsies that were collected before cancer diagnosis, higher levels of ecDNA were present in samples from patients who later developed EAC than in samples from those who did not. We found that ecDNAs contained diverse collections of oncogenes and immunomodulatory genes. Furthermore, ecDNAs showed increases in copy number and structural complexity at more advanced stages of disease. Our findings show that ecDNA can develop early in the transition from high-grade dysplasia to cancer, and that ecDNAs progressively form and evolve under positive selection.

Extrachromosomal DNA amplifications in cancer.

Extrachromosomal DNA (ecDNA) amplification is an important driver alteration in cancer. It has been observed in most cancer types and is associated with worse patient outcome. The functional impact of ecDNA has been linked to its unique properties, such as its circular structure that is associated with altered chromatinization and epigenetic regulatory landscape, as well as its ability to randomly segregate during cell division, which fuels intercellular copy number heterogeneity. Recent investigations suggest that ecDNA is structurally more complex than previously anticipated and that it localizes to specialized nuclear bodies (hubs) and can act in trans as an enhancer for genes on other ecDNAs or chromosomes. In this Review, we synthesize what is currently known about how ecDNA is generated and how its genetic and epigenetic architecture affects proto-oncogene deregulation in cancer. We discuss how recently identified ecDNA functions may impact oncogenesis but also serve as new therapeutic vulnerabilities in cancer.

The somatic genomic landscape of glioblastoma.

We describe the landscape of somatic genomic alterations based on multidimensional and comprehensive characterization of more than 500 glioblastoma tumors (GBMs). We identify several novel mutated genes as well as complex rearrangements of signature receptors, including EGFR and PDGFRA. TERT promoter mutations are shown to correlate with elevated mRNA expression, supporting a role in telomerase reactivation. Correlative analyses confirm that the survival advantage of the proneural subtype is conferred by the G-CIMP phenotype, and MGMT DNA methylation may be a predictive biomarker for treatment response only in classical subtype GBM. Integrative analysis of genomic and proteomic profiles challenges the notion of therapeutic inhibition of a pathway as an alternative to inhibition of the target itself. These data will facilitate the discovery of therapeutic and diagnostic target candidates, the validation of research and clinical observations and the generation of unanticipated hypotheses that can advance our molecular understanding of this lethal cancer.

Mosaic analysis with double markers reveals tumor cell of origin in glioma.

Cancer cell of origin is difficult to identify by analyzing cells within terminal stage tumors, whose identity could be concealed by the acquired plasticity. Thus, an ideal approach to identify the cell of origin is to analyze proliferative abnormalities in distinct lineages prior to malignancy. Here, we use mosaic analysis with double markers (MADM) in mice to model gliomagenesis by initiating concurrent p53/Nf1 mutations sporadically in neural stem cells (NSCs). Surprisingly, MADM-based lineage tracing revealed significant aberrant growth prior to malignancy only in oligodendrocyte precursor cells (OPCs), but not in any other NSC-derived lineages or NSCs themselves. Upon tumor formation, phenotypic and transcriptome analyses of tumor cells revealed salient OPC features. Finally, introducing the same p53/Nf1 mutations directly into OPCs consistently led to gliomagenesis. Our findings suggest OPCs as the cell of origin in this model, even when initial mutations occur in NSCs, and highlight the importance of analyzing premalignant stages to identify the cancer cell of origin.

Circular ecDNA promotes accessible chromatin and high oncogene expression.

Oncogenes are commonly amplified on particles of extrachromosomal DNA (ecDNA) in cancer1,2, but our understanding of the structure of ecDNA and its effect on gene regulation is limited. Here, by integrating ultrastructural imaging, long-range optical mapping and computational analysis of whole-genome sequencing, we demonstrate the structure of circular ecDNA. Pan-cancer analyses reveal that oncogenes encoded on ecDNA are among the most highly expressed genes in the transcriptome of the tumours, linking increased copy number with high transcription levels. Quantitative assessment of the chromatin state reveals that although ecDNA is packaged into chromatin with intact domain structure, it lacks higher-order compaction that is typical of chromosomes and displays significantly enhanced chromatin accessibility. Furthermore, ecDNA is shown to have a significantly greater number of ultra-long-range interactions with active chromatin, which provides insight into how the structure of circular ecDNA affects oncogene function, and connects ecDNA biology with modern cancer genomics and epigenetics.

Quality of life after surgery for lower grade gliomas.

BACKGROUND: Few large studies have investigated quality of life (QOL) for adults diagnosed with lower grade glioma (LGG). METHODS: QOL was assessed for 320 adults with LGG (World Health Organization grade 2/3) enrolled in the International Low Grade Glioma Registry by using the Medical Outcomes Study 36-Item Short Form health survey. Data on symptoms were also collected. QOL outcomes were examined by treatment group and also compared to those from a population-based case-control study of meningioma (the Meningioma Consortium), in which 1722 meningioma cases diagnosed among residents of Connecticut, Massachusetts, California, Texas, and North Carolina from May 1, 2006 through March 14, 2013 were enrolled and frequency matched to 1622 controls by age, sex, and geography. RESULTS: The LGG sample average age is 45 years at the time of interview and 53.1% male. Almost 55% of patients had received radiation and chemotherapy (primarily temozolomide); 32.4% had received neither treatment. Two-thirds of participants with LGG report difficulty with speaking, memory, or thinking, and over one of three reports personality change or difficulty driving. After controlling for age and other comorbidities, individuals with LGG report levels of physical, emotional, and mental health functioning below those reported in a meningioma as well as a general healthy population. CONCLUSIONS: Despite being relatively young, persons with LGG report significantly reduced QOL compared to persons with nonmalignant brain tumors and to a control population, which highlights the need to better acknowledge and manage these symptoms for this group of patients diagnosed in the prime of life.

Extrachromosomal DNA (ecDNA): an origin of tumor heterogeneity, genomic remodeling, and drug resistance.

The genome of cancer cells contains circular extrachromosomal DNA (ecDNA) elements not found in normal cells. Analysis of clinical samples reveal they are common in most cancers and their presence indicates poor prognosis. They often contain enhancers and driver oncogenes that are highly expressed. The circular ecDNA topology leads to an open chromatin conformation and generates new gene regulatory interactions, including with distal enhancers. The absence of centromeres leads to random distribution of ecDNAs during cell division and genes encoded on them are transmitted in a non-mendelian manner. ecDNA can integrate into and exit from chromosomal DNA. The numbers of specific ecDNAs can change in response to treatment. This dynamic ability to remodel the cancer genome challenges long-standing fundamentals, providing new insights into tumor heterogeneity, cancer genome remodeling, and drug resistance.

Integrated cistromic and expression analysis of amplified NKX2-1 in lung adenocarcinoma identifies LMO3 as a functional transcriptional target.

The NKX2-1 transcription factor, a regulator of normal lung development, is the most significantly amplified gene in human lung adenocarcinoma. To study the transcriptional impact of NKX2-1 amplification, we generated an expression signature associated with NKX2-1 amplification in human lung adenocarcinoma and analyzed DNA-binding sites of NKX2-1 by genome-wide chromatin immunoprecipitation. Integration of these expression and cistromic analyses identified LMO3, itself encoding a transcription regulator, as a candidate direct transcriptional target of NKX2-1. Further cistromic and overexpression analyses indicated that NKX2-1 can cooperate with the forkhead box transcription factor FOXA1 to regulate LMO3 gene expression. RNAi analysis of NKX2-1-amplified cells compared with nonamplified cells demonstrated that LMO3 mediates cell survival downstream from NKX2-1. Our findings provide new insight into the transcriptional regulatory network of NKX2-1 and suggest that LMO3 is a transcriptional signal transducer in NKX2-1-amplified lung adenocarcinomas.

A survey of intragenic breakpoints in glioblastoma identifies a distinct subset associated with poor survival.

With the advent of high-throughput sequencing technologies, much progress has been made in the identification of somatic structural rearrangements in cancer genomes. However, characterization of the complex alterations and their associated mechanisms remains inadequate. Here, we report a comprehensive analysis of whole-genome sequencing and DNA copy number data sets from The Cancer Genome Atlas to relate chromosomal alterations to imbalances in DNA dosage and describe the landscape of intragenic breakpoints in glioblastoma multiforme (GBM). Gene length, guanine-cytosine (GC) content, and local presence of a copy number alteration were closely associated with breakpoint susceptibility. A dense pattern of repeated focal amplifications involving the murine double minute 2 (MDM2)/cyclin-dependent kinase 4 (CDK4) oncogenes and associated with poor survival was identified in 5% of GBMs. Gene fusions and rearrangements were detected concomitant within the breakpoint-enriched region. At the gene level, we noted recurrent breakpoints in genes such as apoptosis regulator FAF1. Structural alterations of the FAF1 gene disrupted expression and led to protein depletion. Restoration of the FAF1 protein in glioma cell lines significantly increased the FAS-mediated apoptosis response. Our study uncovered a previously underappreciated genomic mechanism of gene deregulation that can confer growth advantages on tumor cells and may generate cancer-specific vulnerabilities in subsets of GBM.

TumorFusions: an integrative resource for cancer-associated transcript fusions.

Gene fusion represents a class of molecular aberrations in cancer and has been exploited for therapeutic purposes. In this paper we describe TumorFusions, a data portal that catalogues 20 731 gene fusions detected in 9966 well characterized cancer samples and 648 normal specimens from The Cancer Genome Atlas (TCGA). The portal spans 33 cancer types in TCGA. Fusion transcripts were identified via a uniform pipeline, including filtering against a list of 3838 transcript fusions detected in a panel of 648 non-neoplastic samples. Fusions were mapped to somatic DNA rearrangements identified using whole genome sequencing data from 561 cancer samples as a means of validation. We observed that 65% of transcript fusions were associated with a chromosomal alteration, which is annotated in the portal. Other features of the portal include links to SNP array-based copy number levels and mutational patterns, exon and transcript level expressions of the partner genes, and a network-based centrality score for prioritizing functional fusions. Our portal aims to be a broadly applicable and user friendly resource for cancer gene annotation and is publicly available at http://www.tumorfusions.org.

PAF promotes stemness and radioresistance of glioma stem cells.

An integrated genomic and functional analysis to elucidate DNA damage signaling factors promoting self-renewal of glioma stem cells (GSCs) identified proliferating cell nuclear antigen (PCNA)-associated factor (PAF) up-regulation in glioblastoma. PAF is preferentially overexpressed in GSCs. Its depletion impairs maintenance of self-renewal without promoting differentiation and reduces tumor-initiating cell frequency. Combined transcriptomic and metabolomic analyses revealed that PAF supports GSC maintenance, in part, by influencing DNA replication and pyrimidine metabolism pathways. PAF interacts with PCNA and regulates PCNA-associated DNA translesion synthesis (TLS); consequently, PAF depletion in combination with radiation generated fewer tumorspheres compared with radiation alone. Correspondingly, pharmacological impairment of DNA replication and TLS phenocopied the effect of PAF depletion in compromising GSC self-renewal and radioresistance, providing preclinical proof of principle that combined TLS inhibition and radiation therapy may be a viable therapeutic option in the treatment of glioblastoma multiforme (GBM).

Comprehensive, Integrative Genomic Analysis of Diffuse Lower-Grade Gliomas.

BACKGROUND: Diffuse low-grade and intermediate-grade gliomas (which together make up the lower-grade gliomas, World Health Organization grades II and III) have highly variable clinical behavior that is not adequately predicted on the basis of histologic class. Some are indolent; others quickly progress to glioblastoma. The uncertainty is compounded by interobserver variability in histologic diagnosis. Mutations in IDH, TP53, and ATRX and codeletion of chromosome arms 1p and 19q (1p/19q codeletion) have been implicated as clinically relevant markers of lower-grade gliomas. METHODS: We performed genomewide analyses of 293 lower-grade gliomas from adults, incorporating exome sequence, DNA copy number, DNA methylation, messenger RNA expression, microRNA expression, and targeted protein expression. These data were integrated and tested for correlation with clinical outcomes. RESULTS: Unsupervised clustering of mutations and data from RNA, DNA-copy-number, and DNA-methylation platforms uncovered concordant classification of three robust, nonoverlapping, prognostically significant subtypes of lower-grade glioma that were captured more accurately by IDH, 1p/19q, and TP53 status than by histologic class. Patients who had lower-grade gliomas with an IDH mutation and 1p/19q codeletion had the most favorable clinical outcomes. Their gliomas harbored mutations in CIC, FUBP1, NOTCH1, and the TERT promoter. Nearly all lower-grade gliomas with IDH mutations and no 1p/19q codeletion had mutations in TP53 (94%) and ATRX inactivation (86%). The large majority of lower-grade gliomas without an IDH mutation had genomic aberrations and clinical behavior strikingly similar to those found in primary glioblastoma. CONCLUSIONS: The integration of genomewide data from multiple platforms delineated three molecular classes of lower-grade gliomas that were more concordant with IDH, 1p/19q, and TP53 status than with histologic class. Lower-grade gliomas with an IDH mutation either had 1p/19q codeletion or carried a TP53 mutation. Most lower-grade gliomas without an IDH mutation were molecularly and clinically similar to glioblastoma. (Funded by the National Institutes of Health.).

Whole-genome and multisector exome sequencing of primary and post-treatment glioblastoma reveals patterns of tumor evolution.

Glioblastoma (GBM) is a prototypical heterogeneous brain tumor refractory to conventional therapy. A small residual population of cells escapes surgery and chemoradiation, resulting in a typically fatal tumor recurrence ∼ 7 mo after diagnosis. Understanding the molecular architecture of this residual population is critical for the development of successful therapies. We used whole-genome sequencing and whole-exome sequencing of multiple sectors from primary and paired recurrent GBM tumors to reconstruct the genomic profile of residual, therapy resistant tumor initiating cells. We found that genetic alteration of the p53 pathway is a primary molecular event predictive of a high number of subclonal mutations in glioblastoma. The genomic road leading to recurrence is highly idiosyncratic but can be broadly classified into linear recurrences that share extensive genetic similarity with the primary tumor and can be directly traced to one of its specific sectors, and divergent recurrences that share few genetic alterations with the primary tumor and originate from cells that branched off early during tumorigenesis. Our study provides mechanistic insights into how genetic alterations in primary tumors impact the ensuing evolution of tumor cells and the emergence of subclonal heterogeneity.

Reconstructing the molecular life history of gliomas.

At the time of their clinical manifestation, the heterogeneous group of adult and pediatric gliomas carries a wide range of diverse somatic genomic alterations, ranging from somatic single-nucleotide variants to structural chromosomal rearrangements. Somatic abnormalities may have functional consequences, such as a decrease, increase or change in mRNA transcripts, and cells pay a penalty for maintaining them. These abnormalities, therefore, must provide cells with a competitive advantage to become engrained into the glioma genome. Here, we propose a model of gliomagenesis consisting of the following five consecutive phases that glioma cells have traversed prior to clinical manifestation: (I) initial growth; (II) oncogene-induced senescence; (III) stressed growth; (IV) replicative senescence/crisis; (V) immortal growth. We have integrated the findings from a large number of studies in biology and (neuro)oncology and relate somatic alterations and other results discussed in these papers to each of these five phases. Understanding the story that each glioma tells at presentation may ultimately facilitate the design of novel, more effective therapeutic approaches.

Transformation by the (R)-enantiomer of 2-hydroxyglutarate linked to EGLN activation.

The identification of succinate dehydrogenase (SDH), fumarate hydratase (FH) and isocitrate dehydrogenase (IDH) mutations in human cancers has rekindled the idea that altered cellular metabolism can transform cells. Inactivating SDH and FH mutations cause the accumulation of succinate and fumarate, respectively, which can inhibit 2-oxoglutarate (2-OG)-dependent enzymes, including the EGLN prolyl 4-hydroxylases that mark the hypoxia inducible factor (HIF) transcription factor for polyubiquitylation and proteasomal degradation. Inappropriate HIF activation is suspected of contributing to the pathogenesis of SDH-defective and FH-defective tumours but can suppress tumour growth in some other contexts. IDH1 and IDH2, which catalyse the interconversion of isocitrate and 2-OG, are frequently mutated in human brain tumours and leukaemias. The resulting mutants have the neomorphic ability to convert 2-OG to the (R)-enantiomer of 2-hydroxyglutarate ((R)-2HG). Here we show that (R)-2HG, but not (S)-2HG, stimulates EGLN activity, leading to diminished HIF levels, which enhances the proliferation and soft agar growth of human astrocytes. These findings define an enantiomer-specific mechanism by which the (R)-2HG that accumulates in IDH mutant brain tumours promotes transformation and provide a justification for exploring EGLN inhibition as a potential treatment strategy.

Suppression of lung adenocarcinoma progression by Nkx2-1.

Despite the high prevalence and poor outcome of patients with metastatic lung cancer the mechanisms of tumour progression and metastasis remain largely uncharacterized. Here we modelled human lung adenocarcinoma, which frequently harbours activating point mutations in KRAS and inactivation of the p53 pathway, using conditional alleles in mice. Lentiviral-mediated somatic activation of oncogenic Kras and deletion of p53 in the lung epithelial cells of Kras(LSL-G12D/+);p53(flox/flox) mice initiates lung adenocarcinoma development. Although tumours are initiated synchronously by defined genetic alterations, only a subset becomes malignant, indicating that disease progression requires additional alterations. Identification of the lentiviral integration sites allowed us to distinguish metastatic from non-metastatic tumours and determine the gene expression alterations that distinguish these tumour types. Cross-species analysis identified the NK2-related homeobox transcription factor Nkx2-1 (also called Ttf-1 or Titf1) as a candidate suppressor of malignant progression. In this mouse model, Nkx2-1 negativity is pathognomonic of high-grade poorly differentiated tumours. Gain- and loss-of-function experiments in cells derived from metastatic and non-metastatic tumours demonstrated that Nkx2-1 controls tumour differentiation and limits metastatic potential in vivo. Interrogation of Nkx2-1-regulated genes, analysis of tumours at defined developmental stages, and functional complementation experiments indicate that Nkx2-1 constrains tumours in part by repressing the embryonically restricted chromatin regulator Hmga2. Whereas focal amplification of NKX2-1 in a fraction of human lung adenocarcinomas has focused attention on its oncogenic function, our data specifically link Nkx2-1 downregulation to loss of differentiation, enhanced tumour seeding ability and increased metastatic proclivity. Thus, the oncogenic and suppressive functions of Nkx2-1 in the same tumour type substantiate its role as a dual function lineage factor.