Loss of chromosome Y in primary tumors

Certain cancer types afflict female and male patients disproportionately. The reasons include differences in male/female physiology, effect of sex hormones, risk behavior, environmental exposures, and genetics of the sex chromosomes X and Y. Loss of Y (LOY) is common in peripheral blood cells in aging men, and this phenomenon is associated with several diseases. However, the frequency and role of LOY in tumors is little understood. Here, we present a comprehensive catalog of LOY in >5,000 primary tumors from male patients in the TCGA. We show that LOY rates vary by tumor type and provide evidence for LOY being either a passenger or driver event depending on context. LOY in uveal melanoma specifically is associated with age and survival and is an independent predictor of poor outcome. LOY creates common dependencies on DDX3X and EIF1AX in male cell lines, suggesting that LOY generates unique vulnerabilities that could be therapeutically exploited.

False-Negative Testing for FIP1L1::PDGFRA by Fluorescence in situ Hybridization Is a Frequent Cause of Diagnostic Delay

The imatinib-sensitive fusion gene FIP1L1::PDGFRA is the most frequent molecular abnormality identified in patients with eosinophilic myeloid neoplasms. Rapid recognition of this mutation is essential given the poor prognosis of PDGFRA-associated myeloid neoplasms prior to the availability of imatinib therapy. We report a case of a patient in whom delayed diagnosis resulted in cardiac transplantation for eosinophilic endomyocardial fibrosis. The delay in diagnosis was due, in part, to a false-negative result in fluorescence in situ hybridization (FISH) testing for FIP1L1::PDGFRA. To explore this further, we examined our cohort of patients presenting with confirmed or suspected eosinophilic myeloid neoplasms and found 8 additional patients with negative FISH results despite a positive reverse-transcriptase polymerase chain reaction test for FIP1L1::PDGFRA. More importantly, false-negative FISH results delayed the median time to imatinib treatment by 257 days. These data emphasize the importance of empiric imatinib therapy in patients with clinical features suggestive of PDGFRA-associated disease.

Abstract 3133: Loss of the Y chromosome as a candidate driver in solid tumors

It has long been recognized that certain cancer types afflict female and male patients disproportionately. The reasons for this disparity include differences in male and female physiology, effect of sex hormones, risk behavior and environmental exposures, and different copy number of the sex chromosomes X and Y. Loss of Y (LOY) is common in peripheral blood cells in aging men, and this phenomenon is associated with several diseases. However, the frequency and role of LOY in tumors is little understood. Here, we present a comprehensive catalog of LOY in about 5,000 primary tumors from male patients in the TCGA. We show that LOY rates vary by tumor type, and provide evidence for LOY being either a passenger or driver event depending on context. LOY in uveal melanoma specifically is associated with age, survival and is an independent predictor of poor outcome. LOY creates common dependencies on DDX3X and EIF1AX in male cell lines, suggesting that LOY generates unique vulnerabilities that could be therapeutically exploited.

Citation Format: Meifang Qi, Esther Rheinbay, Andrew A. Lane, Jiali Pang, Irene Mitsiades. Loss of the Y chromosome as a candidate driver in solid tumors [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 3133.

Abstract P3-08-03: HOXB13 is coamplified with HER2 in a subset of breast cancers

Background: The Breast Cancer Index (BCI) is an algorithmic gene expression-based signature comprised of two functional biomarker panels, the Molecular Grade Index (MGI) and the two-gene ratio, HOXB13/IL17BR. The HOXB13/IL17BR ratio has been shown to predict endocrine response across several different treatment scenarios, and the homeobox factor, HOXB13, is the primary determinant of endocrine benefit and response. Previously, we found that HOXB13 gene expression was positively correlated with HER2 in estrogen receptor-positive (ER+) but not ER-negative (ER-) tumors. To validate these findings, we examined HOXB13 gene expression in The Cancer Genome Atlas (TCGA) breast cancer cohort. In addition, considering that the HOXB13 gene locus resides 9 megabases from the HER2 locus at chromosome 17q21-22, we examined whether HOXB13 was coamplified with HER2 in HER2-positive breast cancers. Methods: Messenger RNA (mRNA) sequencing (mRNAseq) data from the TCGA breast cancer dataset were used to correlate HOXB13 gene expression with immunohistochemical (IHC)-based HER2 expression in ER+ (ER+PR+, ER+PR-) and ER- (ER-PR+ and ER-PR-) breast cancers. Thresholded copy number data from GISTIC output and mRNAseq data from the TCGA datasets were used to correlate HER2 and HOXB13 gene amplification (at least 4 inferred copies) with HOXB13 gene expression in IHC-based ER+ (n= 562) and ER- (n=163) breast cancer subsets. Formalin-fixed paraffin embedded (FFPE) samples from 80 consecutive HER2-amplified invasive breast carcinomas (64 ER+, 16 ER-) diagnosed between 2010 and 2013 were retrospectively collected from the MGH clinical archives and processed for fluorescence in situ hybridization (FISH) for HER2 (spectrum orange) HOXB13 (spectrum green) and CEP17 (spectrum blue). A ratio of > 2.0 HER2 and HOXB13 to CEP17 signals in at least 60 interphase tumor cell nuclei was considered as amplification of HER2 and HOXB13.Results: In the TCGA breast tumors, HOXB13 gene expression was significantly higher in HER2 positive (by IHC status) tumors in both ER+ (p= 3.7 x10-5) and ER- (p= 4.7 x10-4) tumors compared to HER2-negative samples. In HER2+ TCGA breast tumors, genomic analysis revealed co-amplification of HOXB13 with HER2 in 27% of all tumors, and 25.3% and 31.7% in the ER+ and ER- subsets, respectively. Amplification of the HOXB13 locus alone was observed in 2.4% of HER2-non amplified ER+ and ER- TCGA tumors. In the TCGA cohort, median HOXB13 gene expression was not significantly different when comparing HER2-amplified and HER2/HOXB13 coamplified ER+ (p= 0.33) and ER- (p=0.052) tumors, HER2-amplified and HER2-nonamplified ER+ (p=0.42) and ER- (p=0.67) tumors, and HER2/HOXB13 coamplified and HER2-nonamplified ER+ (p=0.50) and ER- (p=0.11) tumors. Lastly, in the consecutive case series of 80 HER2-amplified breast cancers from MGH, HOXB13 was coamplified in 23% of all tumors, and 23.4% and 18.75% in the ER+ and ER- subsets, respectively. Conclusion: HOXB13 was coamplified with HER2 in a substantial subset of TCGA and MGH primary breast cancers irrespective of ER status.

Citation Format: Dennis Sgroi, Maristela Onozato, Irene Mitsiades, Ryan McMullin, Daniel Hicks, Angela Volorio, Esther Rheinbay. HOXB13 is coamplified with HER2 in a subset of breast cancers [abstract]. In: Proceedings of the 2021 San Antonio Breast Cancer Symposium; 2021 Dec 7-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2022;82(4 Suppl):Abstract nr P3-08-03.

cDNA-detector: detection and removal of cDNA contamination in DNA sequencing libraries

BACKGROUND

Exogenous cDNA introduced into an experimental system, either intentionally or accidentally, can appear as added read coverage over that gene in next-generation sequencing libraries derived from this system. If not properly recognized and managed, this cross-contamination with exogenous signal can lead to incorrect interpretation of research results. Yet, this problem is not routinely addressed in current sequence processing pipelines.

RESULTS

We present cDNA-detector, a computational tool to identify and remove exogenous cDNA contamination in DNA sequencing experiments. We demonstrate that cDNA-detector can identify cDNAs quickly and accurately from alignment files. A source inference step attempts to separate endogenous cDNAs (retrocopied genes) from potential cloned, exogenous cDNAs. cDNA-detector provides a mechanism to decontaminate the alignment from detected cDNAs. Simulation studies show that cDNA-detector is highly sensitive and specific, outperforming existing tools. We apply cDNA-detector to several highly-cited public databases (TCGA, ENCODE, NCBI SRA) and show that contaminant genes appear in sequencing experiments where they lead to incorrect coverage peak calls.

CONCLUSIONS

cDNA-detector is a user-friendly and accurate tool to detect and remove cDNA detection in NGS libraries. This two-step design reduces the risk of true variant removal since it allows for manual review of candidates. We find that contamination with intentionally and accidentally introduced cDNAs is an underappreciated problem even in widely-used consortium datasets, where it can lead to spurious results. Our findings highlight the importance of sensitive detection and removal of contaminant cDNA from NGS libraries before downstream analysis.

CXCR6 positions cytotoxic T cells to receive critical survival signals in the tumor microenvironment

Cytotoxic T lymphocyte (CTL) responses against tumors are maintained by stem-like memory cells that self-renew but also give rise to effector-like cells. The latter gradually lose their anti-tumor activity and acquire an epigenetically fixed, hypofunctional state, leading to tumor tolerance. Here, we show that the conversion of stem-like into effector-like CTLs involves a major chemotactic reprogramming that includes the upregulation of chemokine receptor CXCR6. This receptor positions effector-like CTLs in a discrete perivascular niche of the tumor stroma that is densely occupied by CCR7+ dendritic cells (DCs) expressing the CXCR6 ligand CXCL16. CCR7+ DCs also express and trans-present the survival cytokine interleukin-15 (IL-15). CXCR6 expression and IL-15 trans-presentation are critical for the survival and local expansion of effector-like CTLs in the tumor microenvironment to maximize their anti-tumor activity before progressing to irreversible dysfunction. These observations reveal a cellular and molecular checkpoint that determines the magnitude and outcome of anti-tumor immune responses.

RNA transcripts stimulate homologous recombination by forming DR-loops

Homologous recombination (HR) repairs DNA double-strand breaks (DSBs) in the S and G2 phases of the cell cycle^1-3. Several HR proteins are preferentially recruited to DSBs at transcriptionally active loci^4-10, but how transcription promotes HR is poorly understood. Here we develop an assay to assess the effect of local transcription on HR. Using this assay, we find that transcription stimulates HR to a substantial extent. Tethering RNA transcripts to the vicinity of DSBs recapitulates the effects of local transcription, which suggests that transcription enhances HR through RNA transcripts. Tethered RNA transcripts stimulate HR in a sequence- and orientation-dependent manner, indicating that they function by forming DNA-RNA hybrids. In contrast to most HR proteins, RAD51-associated protein 1 (RAD51AP1) only promotes HR when local transcription is active. RAD51AP1 drives the formation of R-loops in vitro and is required for tethered RNAs to stimulate HR in cells. Notably, RAD51AP1 is necessary for the DSB-induced formation of DNA-RNA hybrids in donor DNA, linking R-loops to D-loops. In vitro, RAD51AP1-generated R-loops enhance the RAD51-mediated formation of D-loops locally and give rise to intermediates that we term 'DR-loops', which contain both DNA-DNA and DNA-RNA hybrids and favour RAD51 function. Thus, at DSBs in transcribed regions, RAD51AP1 promotes the invasion of RNA transcripts into donor DNA, and stimulates HR through the formation of DR-loops.

Whole-genome characterization of lung adenocarcinomas lacking the RTK/RAS/RAF pathway

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.

Framework for quality assessment of whole genome cancer sequences

Bringing together cancer genomes from different projects increases power and allows the investigation of pan-cancer, molecular mechanisms. However, working with whole genomes sequenced over several years in different sequencing centres requires a framework to compare the quality of these sequences. We used the Pan-Cancer Analysis of Whole Genomes cohort as a test case to construct such a framework. This cohort contains whole cancer genomes of 2832 donors from 18 sequencing centres. We developed a non-redundant set of five quality control (QC) measurements to establish a star rating system. These QC measures reflect known differences in sequencing protocol and provide a guide to downstream analyses and allow for exclusion of samples of poor quality. We have found that this is an effective framework of quality measures. The implementation of the framework is available at: https://dockstore.org/containers/quay.io/jwerner_dkfz/pancanqc:1.2.2.

Working with cancer genomes from multiple projects can increase investigative power, but quality of sequences can vary. Here, the authors present a framework for comparing whole genome sequencing quality to help researchers guide downstream analyses and exclude poor quality samples.

Abstract 5895: Whole-genome characterization of lung adenocarcinomas lacking alterations in RTK/RAS/RAF/MAPK pathway

Lung adenocarcinomas (LUAD) are typically characterized by genetic activation of the receptor tyrosine kinase (RTK)/RAS/RAF/MAP kinase (MAPK) pathway. A minority of LUAD cases (20-25%) lack apparent genetic alterations in this pathway, and thus are ineligible for most targeted therapies. These candidate “oncogene negative” LUADs may harbor novel classes of oncogenic drivers or represent a biologically distinct class of tumors. To characterize the genomic landscape of oncogene-negative LUADs, we nominated 98 cases that were found to lack an activating RTK/RAS/RAF/MAPK pathway alteration in a TCGA study utilizing whole exome sequencing, microarray, and transcriptome data. We profiled these tumors with high-depth whole genome sequencing (WGS), with the goal of identifying noncoding and structural variant driver DNA alterations in both known and novel loci. Of the 98 cases, 20 harbored somatic KRAS mutations that had been missed in the prior WES and transcriptome studies because of insufficient coverage, including 8 cases with the recently targetable p.G12C mutation. 16 samples harbored oncogenic or loss-of-function structural variants in FGFR1, MAPK1, EGFR, NF1, RASA1, ARAF, NTRK2 and NRG1. 5 other samples with SNV or indels in EGFR, ERBB2 and SOS1 were reclassified as oncogene positive. Thus via comprehensive genomic analysis, we confirmed that 57 of the 98 WGS cases did not harbor any detectable alterations in genes encoding any known RTK/RAS/RAF/MAPK members, representing 13% cases chosen as “lung adenocarcinomas” for the TCGA study. Among the 57 confirmed oncogene-negative LUADs, we identified focal deletions targeting the promoter and transcription start site of tumor suppressor genes STK11, KEAP1 and SMARCA4 in 10 samples. Expression and methylation profiling suggested an enrichment of the TP53-deficient phenotype, including cell cycle and FOXM1 deregulation, among the oncogene-negative samples. Moreover, novel promoter mutations associated with increased expression were identified in ILF2, which regulates DNA damage response pathways. Finally, a subset of confirmed oncogene-negative LUADs harbored increased expression of neuroendocrine markers, suggesting that these oncogene-negative samples may either be mis-diagnosed as LUAD or represent LUAD with mixed features of other subtypes of lung cancer; indeed, 14 of the 57 confirmed oncogene-negative cases show histological features of large cell neuroendocrine lung carcinoma. This would suggest that 10% of the cases in this study are both lung adenocarcinoma and “oncogene-negative” to date. Our results provide some of the first comprehensive genomic characterization of oncogene-negative LUADs, implicating TP53 and structural variants in the pathogenesis of this common and difficult to treat entity.

Citation Format: Jian Carrot-Zhang, Siddhartha Devarakonda, Nicolas Robine, Xiaotong Yao, Tiago C. Silva, Jeff Damrauer, Aditya Deshpande, Ming-Sound Tsao, Christina Yao, Chris Wong, Lisui Bao, Hyo Young Choi, Ina Felau, Jean C. Zenklusen, Gordon Robertson, Tuan Trieua, Wei-Wei Liang, Meng Zhou, Esther Rheinbay, Neil Hayes, Ekta Khurana, Li Ding, Peter Laird, Olivier Elemento, John Weinstein, David Kwiatkowski, Chris Benz, Josh Stuart, Lixing Yang, Mauro Castro, William Travis, Katherine Hoadley, Ben Berman, TCGA Analysis Network, Matthew Meyerson, Ramaswamy Govindan, Marcin Imielinski. Whole-genome characterization of lung adenocarcinomas lacking alterations in RTK/RAS/RAF/MAPK pathway [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 5895.

Pan-cancer analysis of whole genomes

Cancer is driven by genetic change, and the advent of massively parallel sequencing has enabled systematic documentation of this variation at the whole-genome scale1-3. Here we report the integrative analysis of 2,658 whole-cancer genomes and their matching normal tissues across 38 tumour types from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). We describe the generation of the PCAWG resource, facilitated by international data sharing using compute clouds. On average, cancer genomes contained 4-5 driver mutations when combining coding and non-coding genomic elements; however, in around 5% of cases no drivers were identified, suggesting that cancer driver discovery is not yet complete. Chromothripsis, in which many clustered structural variants arise in a single catastrophic event, is frequently an early event in tumour evolution; in acral melanoma, for example, these events precede most somatic point mutations and affect several cancer-associated genes simultaneously. Cancers with abnormal telomere maintenance often originate from tissues with low replicative activity and show several mechanisms of preventing telomere attrition to critical levels. Common and rare germline variants affect patterns of somatic mutation, including point mutations, structural variants and somatic retrotransposition. A collection of papers from the PCAWG Consortium describes non-coding mutations that drive cancer beyond those in the TERT promoter4; identifies new signatures of mutational processes that cause base substitutions, small insertions and deletions and structural variation5,6; analyses timings and patterns of tumour evolution7; describes the diverse transcriptional consequences of somatic mutation on splicing, expression levels, fusion genes and promoter activity8,9; and evaluates a range of more-specialized features of cancer genomes8,10-18.

Analyses of non-coding somatic drivers in 2,658 cancer whole genomes

The discovery of drivers of cancer has traditionally focused on protein-coding genes1-4. Here we present analyses of driver point mutations and structural variants in non-coding regions across 2,658 genomes from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium5 of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). For point mutations, we developed a statistically rigorous strategy for combining significance levels from multiple methods of driver discovery that overcomes the limitations of individual methods. For structural variants, we present two methods of driver discovery, and identify regions that are significantly affected by recurrent breakpoints and recurrent somatic juxtapositions. Our analyses confirm previously reported drivers6,7, raise doubts about others and identify novel candidates, including point mutations in the 5' region of TP53, in the 3' untranslated regions of NFKBIZ and TOB1, focal deletions in BRD4 and rearrangements in the loci of AKR1C genes. We show that although point mutations and structural variants that drive cancer are less frequent in non-coding genes and regulatory sequences than in protein-coding genes, additional examples of these drivers will be found as more cancer genomes become available.

DeTiN: overcoming tumor-in-normal contamination

A key step in achieving accurate detection of somatic mutations is comparison of sequencing data from a tumor sample to its matched germline control. Sensitivity to detect somatic variants is greatly reduced when the matched normal sample is contaminated with tumor cells. To overcome this limitation, we developed deTiN, a method that first estimates the tumor-in-normal contamination (TiN) level, and then, in contaminated cases, improves sensitivity by reclassifying initially discarded variants as somatic.

Abstract LB-280: The landscape of somatic microsatellite indels across cancer: detection and identification of driver events

Microsatellites (MSs) are tracts of variable-length repeats of short DNA motifs that are abundant in the human genome and exhibit high rates of mutations in the form of insertions or deletions of the repeated motif (MS indels). Despite their prevalence, the contribution of somatic MS indels to cancer is largely unexplored due to difficulties in detecting them and assessing their significance. Here, we present a comprehensive analysis of MS indels across 20 tumor types. We characterize the overall MS indel landscape and detect genes with candidate driver MS indel events. We present two novel tools: MSMuTect for accurate detection of somatic MS indels and MSMutSig for identifying candidate cancer genes containing events at higher frequency than expected by chance. We observe high variability of the frequency of MS indels across tumors and demonstrate that the number and pattern of MS indels can accurately distinguish microsatellite stable (MSS) tumors from tumors with microsatellite instability (MSI). Applying MSMutSig across 6,788 tumors from 20 different tumor types identified 7 genes with significant MS indel hotspots: ACVR2A, RNF43, DOCK3, MSH3, ESRP1, PRDM2 and JAK1. In the four genes that have been previously implicated in cancer (ACVR2A, RNF43, JAK1 and MSH3), we identified previously unreported MS indels events. Three of the genes with significant loci - DOCK3, PRDM2 and ESRP1- had not been previously listed as cancer genes. MS indels in DOCK3, a negative regulator of the WNT pathway, were mutually exclusive with mutations in CTNNB1. MS indels in ESRP1, an RNA processing gene, correlated with alternative splicing of FGFR2, an event associated with the epithelial-to-mesenchymal transition. Overall, our comprehensive analysis of somatic MS indels across cancer highlights their importance, particularly in

MSI tumors, significantly contributes to the ongoing global efforts to detect cancer genes, and may improve classification of patients into clinically-relevant subgroups.

Citation Format: Yosef E. Maruvka, Kent W. Mouw, Rosa Karlic, Rasanna Parasuraman, Atanas Kamburov, Paz Polak, Nicholas J. Haradhvala, Julian M. Hess, Esther Rheinbay, Yehuda Brody, Lior Z. Braunstein, Alan D’Andrea, Michael S. Lawrence, Adam Bass, Andre Bernards, Franziska Michor, Gad Getz. The landscape of somatic microsatellite indels across cancer: detection and identification of driver events [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr LB-280. doi:10.1158/1538-7445.AM2017-LB-280

Recurrent and functional regulatory mutations in breast cancer

Genomic analysis of tumours has led to the identification of hundreds of cancer genes on the basis of the presence of mutations in protein-coding regions. By contrast, much less is known about cancer-causing mutations in non-coding regions. Here we perform deep sequencing in 360 primary breast cancers and develop computational methods to identify significantly mutated promoters. Clear signals are found in the promoters of three genes. FOXA1, a known driver of hormone-receptor positive breast cancer, harbours a mutational hotspot in its promoter leading to overexpression through increased E2F binding. RMRP and NEAT1, two non-coding RNA genes, carry mutations that affect protein binding to their promoters and alter expression levels. Our study shows that promoter regions harbour recurrent mutations in cancer with functional consequences and that the mutations occur at similar frequencies as in coding regions. Power analyses indicate that more such regions remain to be discovered through deep sequencing of adequately sized cohorts of patients.

Mutational Strand Asymmetries in Cancer Genomes Reveal Mechanisms of DNA Damage and Repair

Mutational processes constantly shape the somatic genome, leading to immunity, aging, cancer, and other diseases. When cancer is the outcome, we are afforded a glimpse into these processes by the clonal expansion of the malignant cell. Here, we characterize a less explored layer of the mutational landscape of cancer: mutational asymmetries between the two DNA strands. Analyzing whole-genome sequences of 590 tumors from 14 different cancer types, we reveal widespread asymmetries across mutagenic processes, with transcriptional ("T-class") asymmetry dominating UV-, smoking-, and liver-cancer-associated mutations and replicative ("R-class") asymmetry dominating POLE-, APOBEC-, and MSI-associated mutations. We report a striking phenomenon of transcription-coupled damage (TCD) on the non-transcribed DNA strand and provide evidence that APOBEC mutagenesis occurs on the lagging-strand template during DNA replication. As more genomes are sequenced, studying and classifying their asymmetries will illuminate the underlying biological mechanisms of DNA damage and repair.

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

EWS-FLI1 utilizes divergent chromatin remodeling mechanisms to directly activate or repress enhancer elements in Ewing sarcoma

The aberrant transcription factor EWS-FLI1 drives Ewing sarcoma, but its molecular function is not completely understood. We find that EWS-FLI1 reprograms gene regulatory circuits in Ewing sarcoma by directly inducing or repressing enhancers. At GGAA repeat elements, which lack evolutionary conservation and regulatory potential in other cell types, EWS-FLI1 multimers induce chromatin opening and create de novo enhancers that physically interact with target promoters. Conversely, EWS-FLI1 inactivates conserved enhancers containing canonical ETS motifs by displacing wild-type ETS transcription factors. These divergent chromatin-remodeling patterns repress tumor suppressors and mesenchymal lineage regulators while activating oncogenes and potential therapeutic targets, such as the kinase VRK1. Our findings demonstrate how EWS-FLI1 establishes an oncogenic regulatory program governing both tumor survival and differentiation.

Reconstructing and reprogramming the tumor-propagating potential of glioblastoma stem-like cells

Developmental fate decisions are dictated by master transcription factors (TFs) that interact with cis-regulatory elements to direct transcriptional programs. Certain malignant tumors may also depend on cellular hierarchies reminiscent of normal development but superimposed on underlying genetic aberrations. In glioblastoma (GBM), a subset of stem-like tumor-propagating cells (TPCs) appears to drive tumor progression and underlie therapeutic resistance yet remain poorly understood. Here, we identify a core set of neurodevelopmental TFs (POU3F2, SOX2, SALL2, and OLIG2) essential for GBM propagation. These TFs coordinately bind and activate TPC-specific regulatory elements and are sufficient to fully reprogram differentiated GBM cells to "induced" TPCs, recapitulating the epigenetic landscape and phenotype of native TPCs. We reconstruct a network model that highlights critical interactions and identifies candidate therapeutic targets for eliminating TPCs. Our study establishes the epigenetic basis of a developmental hierarchy in GBM, provides detailed insight into underlying gene regulatory programs, and suggests attendant therapeutic strategies. PAPERCLIP:

An aberrant transcription factor network essential for Wnt signaling and stem cell maintenance in glioblastoma

Glioblastoma (GBM) is thought to be driven by a subpopulation of cancer stem cells (CSCs) that self-renew and recapitulate tumor heterogeneity yet remain poorly understood. Here, we present a comparative analysis of chromatin state in GBM CSCs that reveals widespread activation of genes normally held in check by Polycomb repressors. These activated targets include a large set of developmental transcription factors (TFs) whose coordinated activation is unique to the CSCs. We demonstrate that a critical factor in the set, ASCL1, activates Wnt signaling by repressing the negative regulator DKK1. We show that ASCL1 is essential for the maintenance and in vivo tumorigenicity of GBM CSCs. Genome-wide binding profiles for ASCL1 and the Wnt effector LEF-1 provide mechanistic insight and suggest widespread interactions between the TF module and the signaling pathway. Our findings demonstrate regulatory connections among ASCL1, Wnt signaling, and collaborating TFs that are essential for the maintenance and tumorigenicity of GBM CSCs.

H2A.Z landscapes and dual modifications in pluripotent and multipotent stem cells underlie complex genome regulatory functions

Background

The histone variant H2A.Z has been implicated in nucleosome exchange, transcriptional activation and Polycomb repression. However, the relationships among these seemingly disparate functions remain obscure.

Results

We mapped H2A.Z genome-wide in mammalian ES cells and neural progenitors. H2A.Z is deposited promiscuously at promoters and enhancers, and correlates strongly with H3K4 methylation. Accordingly, H2A.Z is present at poised promoters with bivalent chromatin and at active promoters with H3K4 methylation, but is absent from stably repressed promoters that are specifically enriched for H3K27 trimethylation. We also characterized post-translational modification states of H2A.Z, including a novel species dually-modified by ubiquitination and acetylation that is enriched at bivalent chromatin.

Conclusions

Our findings associate H2A.Z with functionally distinct genomic elements, and suggest that post-translational modifications may reconcile its contrasting locations and roles.

Imp2 controls oxidative phosphorylation and is crucial for preserving glioblastoma cancer stem cells

Growth of numerous cancer types is believed to be driven by a subpopulation of poorly differentiated cells, often referred to as cancer stem cells (CSCs), that have the capacity for self-renewal, tumor initiation, and generation of nontumorigenic progeny. Despite their potentially key role in tumor establishment and maintenance, the energy requirements of these cells and the mechanisms that regulate their energy production are unknown. Here, we show that the oncofetal insulin-like growth factor 2 mRNA-binding protein 2 (IMP2, IGF2BP2) regulates oxidative phosphorylation (OXPHOS) in primary glioblastoma (GBM) sphere cultures (gliomaspheres), an established in vitro model for CSC expansion. We demonstrate that IMP2 binds several mRNAs that encode mitochondrial respiratory chain complex subunits and that it interacts with complex I (NADH:ubiquinone oxidoreductase) proteins. Depletion of IMP2 in gliomaspheres decreases their oxygen consumption rate and both complex I and complex IV activity that results in impaired clonogenicity in vitro and tumorigenicity in vivo. Importantly, inhibition of OXPHOS but not of glycolysis abolishes GBM cell clonogenicity. Our observations suggest that gliomaspheres depend on OXPHOS for their energy production and survival and that IMP2 expression provides a key mechanism to ensure OXPHOS maintenance by delivering respiratory chain subunit-encoding mRNAs to mitochondria and contributing to complex I and complex IV assembly.

A tell-tail sign of chromatin: histone mutations drive pediatric glioblastoma

Recent genomic analyses of pediatric glioblastoma, a poorly understood tumor with dismal outcome, have identified mutations in histone H3 variants that affect critical amino acids in the tail. The findings extend discoveries of chromatin regulator inactivation and gain-of-function mutations by documenting alteration of a modifiable histone residue in human cancer.

Genomewide analysis of PRC1 and PRC2 occupancy identifies two classes of bivalent domains

In embryonic stem (ES) cells, bivalent chromatin domains with overlapping repressive (H3 lysine 27 tri-methylation) and activating (H3 lysine 4 tri-methylation) histone modifications mark the promoters of more than 2,000 genes. To gain insight into the structure and function of bivalent domains, we mapped key histone modifications and subunits of Polycomb-repressive complexes 1 and 2 (PRC1 and PRC2) genomewide in human and mouse ES cells by chromatin immunoprecipitation, followed by ultra high-throughput sequencing. We find that bivalent domains can be segregated into two classes—the first occupied by both PRC2 and PRC1 (PRC1-positive) and the second specifically bound by PRC2 (PRC2-only). PRC1-positive bivalent domains appear functionally distinct as they more efficiently retain lysine 27 tri-methylation upon differentiation, show stringent conservation of chromatin state, and associate with an overwhelming number of developmental regulator gene promoters. We also used computational genomics to search for sequence determinants of Polycomb binding. This analysis revealed that the genomewide locations of PRC2 and PRC1 can be largely predicted from the locations, sizes, and underlying motif contents of CpG islands. We propose that large CpG islands depleted of activating motifs confer epigenetic memory by recruiting the full repertoire of Polycomb complexes in pluripotent cells.

Polycomb-group (PcG) proteins play essential roles in the epigenetic regulation of gene expression during development. PcG proteins are repressors that catalyze lysine 27 tri-methylation on histone H3. They are antagonized by trithorax-group proteins that catalyze lysine 4 tri-methylation. Recent studies of ES cells revealed a novel chromatin pattern consisting of overlapping lysine 27 and lysine 4 tri-methylation. Genomic regions with these opposing modifications were termed “bivalent domains” and proposed to silence developmental regulators while keeping them “poised” for alternate fates. However, our understanding of PcG regulation and bivalent domains remains limited. For instance, bivalent domains affect over 2,000 promoters with diverse functions, which suggests that they may function in diverse cellular processes. Moreover, the mechanisms that underlie the targeting of PcG complexes to specific genomic regions remain completely unknown. To gain insight into these issues, we used ultra high-throughput sequencing to map PcG complexes and related modifications genomewide in human and mouse ES cells. The data identify two classes of bivalent domains with distinct regulatory properties. They also reveal striking relationships between genome sequence and chromatin state that suggest a prominent role for the DNA sequence in dictating the genomewide localization of PcG complexes and, consequently, bivalent domains in ES cells.

Phylogenetic and structural analysis of centromeric DNA and kinetochore proteins

Analysis of centromeric DNA and kinetochore proteins suggests that critical structural features of kinetochores have been well conserved from yeast to man.

Background

Kinetochores are large multi-protein structures that assemble on centromeric DNA (CEN DNA) and mediate the binding of chromosomes to microtubules. Comprising 125 base-pairs of CEN DNA and 70 or more protein components, Saccharomyces cerevisiae kinetochores are among the best understood. In contrast, most fungal, plant and animal cells assemble kinetochores on CENs that are longer and more complex, raising the question of whether kinetochore architecture has been conserved through evolution, despite considerable divergence in CEN sequence.

Results

Using computational approaches, ranging from sequence similarity searches to hidden Markov model-based modeling, we show that organisms with CENs resembling those in S. cerevisiae (point CENs) are very closely related and that all contain a set of 11 kinetochore proteins not found in organisms with complex CENs. Conversely, organisms with complex CENs (regional CENs) contain proteins seemingly absent from point-CEN organisms. However, at least three quarters of known kinetochore proteins are present in all fungi regardless of CEN organization. At least six of these proteins have previously unidentified human orthologs. When fungi and metazoa are compared, almost all have kinetochores constructed around Spc105 and three conserved multi-protein linker complexes (MIND, COMA, and the NDC80 complex).

Conclusion

Our data suggest that critical structural features of kinetochores have been well conserved from yeast to man. Surprisingly, phylogenetic analysis reveals that human kinetochore proteins are as similar in sequence to their yeast counterparts as to presumptive Drosophila melanogaster or Caenorhabditis elegans orthologs. This finding is consistent with evidence that kinetochore proteins have evolved very rapidly relative to components of other complex cellular structures.