Trackplot: a fast and lightweight R script for epigenomic enrichment plots

Motivation

BigWig files serve as essential inputs in epigenomic data visualization. However, current R packages for visualizing these files are limited, slow, and burdened by numerous dependencies.

Results

We introduce trackplot, a minimal R script designed for the rapid generation of integrative genomics viewer (IGV) style track plots, profile plots, and heatmaps from bigWig files. This script offers speed, owing to its reliance on bwtool, resulting in performance gains of several magnitudes compared to equivalent packages. The script is lightweight, requiring only the data.table and bwtool packages as primary dependencies. Notably, the plots are generated in base R graphics, eliminating the need for additional packages. trackplot queries the University of California Santa Cruz (UCSC) genome browser for gene models thereby enhancing the reproducibility of analyses. The script extends its support to general transfer format (GTF) further enhancing its versatility. This tool addresses the gaps in existing bigWig visualization approaches by offering speed, simplicity, and minimal dependencies, thereby presenting a valuable asset to researchers in the fields of epigenomics.

Availability and implementation

trackplot is implemented in R is made available under MIT license at https://github.com/PoisonAlien/trackplot.

The GAPDH redox switch safeguards reductive capacity and enables survival of stressed tumour cells

Glyceraldehyde 3-phosphate dehydrogenase (GAPDH) is known to contain an active-site cysteine residue undergoing oxidation in response to hydrogen peroxide, leading to rapid inactivation of the enzyme. Here we show that human and mouse cells expressing a GAPDH mutant lacking this redox switch retain catalytic activity but are unable to stimulate the oxidative pentose phosphate pathway and enhance their reductive capacity. Specifically, we find that anchorage-independent growth of cells and spheroids is limited by an elevation of endogenous peroxide levels and is largely dependent on a functional GAPDH redox switch. Likewise, tumour growth in vivo is limited by peroxide stress and suppressed when the GAPDH redox switch is disabled in tumour cells. The induction of additional intratumoural oxidative stress by chemo- or radiotherapy synergized with the deactivation of the GAPDH redox switch. Mice lacking the GAPDH redox switch exhibit altered fatty acid metabolism in kidney and heart, apparently in compensation for the lack of the redox switch. Together, our findings demonstrate the physiological and pathophysiological relevance of oxidative GAPDH inactivation in mammals.

TAL1 activation in T-Cell acute lymphoblastic leukemia: A novel oncogenic 3' neoenhancer

T-cell acute lymphocytic leukemia protein 1 (TAL1) is one of the most frequently deregulated oncogenes in T-cell acute lymphoblastic leukemia (T-ALL). Its deregulation can occur through diverse in cis-alterations, including SIL-TAL1 microdeletions, translocations with Tcell Receptor (TCR) loci and, more recently described upstream intergenic non-coding mutations. These mutations consist of recurrent focal microinsertions that create an oncogenic neo-enhancer accompanied with activating epigenetic marks. This observation laid the groundwork for an innovative paradigm, the activation of proto-oncogenes via genomic alterations of non-coding intergenic regions. However, for the majority of TAL1 expressing (TAL1+) T-ALLs, the deregulation mechanism remains "unresolved". We took advantage of H3K27ac and H3K4me3 ChIP-seq data of eight T-ALLs, including five TAL1+ cases and identified a putative novel oncogenic neo-enhancer downstream of TAL1 in an "unresolved" monoallelic TAL1+ case. A rare but recurrent somatic heterozygous microinsertion within this region creates a de novo binding site for MYB transcription factor (TF). Here, we demonstrate that this mutation leads to increased enhancer activity, gain of active epigenetic marks and TAL1 activation via recruitment of MYB. These results highlight the diversity of non-coding mutations that can drive oncogene activation.

Genetic and methylation profiles distinguish benign, malignant and spitzoid melanocytic tumors

Accurate classification of melanocytic tumors is important for prognostic evaluation, treatment and follow-up protocols of patients. The majority of melanocytic proliferations can be classified solely based on clinical and pathological criteria, however in select cases a definitive diagnostic assessment remains challenging and additional diagnostic biomarkers would be advantageous. We analyzed melanomas, nevi, Spitz nevi and atypical spitzoid tumors using parallel sequencing (exons of 611 genes and 507 gene translocation analysis) and methylation arrays (850k Illumina EPIC). By combining detailed genetic and epigenetic analysis with reference-based and reference-free DNA methylome deconvolution we compared Spitz nevi to nevi and melanoma and assessed the potential for these methods in classifying challenging spitzoid tumors. Results were correlated with clinical and histologic features. Spitz nevi were found to cluster independently of nevi and melanoma and demonstrated a different mutation profile. Multiple copy number alterations and TERT promoter mutations were identified only in melanomas. Genome-wide methylation in Spitz nevi was comparable to benign nevi while the Leukocytes UnMethylation for Purity (LUMP) algorithm in Spitz nevi was comparable to melanoma. Histologically difficult to classify Spitz tumor cases were assessed which, based on methylation arrays, clustered between Spitz nevi and melanoma and in terms of genetic profile or copy number variations demonstrated worrisome features suggesting a malignant neoplasm. Comprehensive sequencing and methylation analysis verify Spitz nevi as an independent melanocytic entity distinct from both nevi and melanoma. Combined genetic and methylation assays can offer additional insights in diagnosing difficult to classify Spitzoid tumors.

Transcriptome analysis identifies TODL as a novel lncRNA associated with proliferation, differentiation, and tumorigenesis in liposarcoma through FOXM1 Running Title: TODL lncRNA as a potential therapeutic target for liposarcoma

Liposarcoma, the most common soft tissue sarcoma, is a group of fat cell mesenchymal tumors with different histological subtypes. The dysregulation of long non-coding RNAs (lncRNAs) has been observed in human cancers including a few studies in sarcoma. However, the global transcriptome analysis and potential role of lncRNAs remain unexplored in liposarcoma. The present investigation uncovers the transcriptomic profile of liposarcoma by RNA sequencing to gain insight into the global transcriptional changes in liposarcoma. Our RNA sequencing analysis has identified that many oncogenic lncRNAs are differentially expressed in different subtypes of liposarcoma including MALAT1, PVT1, SNHG15, LINC00152, and MIR210HG. Importantly, we identified a highly overexpressed, unannotated, and novel lncRNA in dedifferentiated liposarcomas. We have named it TODL, transcript overexpressed in dedifferentiated liposarcoma. TODL lncRNA displayed significantly higher expression in dedifferentiated liposarcoma cell lines and patient samples. Interestingly, functional studies revealed that TODL lncRNA has an oncogenic function in liposarcoma cells by regulating proliferation, cell cycle, apoptosis, differentiation, and tumorigenesis in the murine model. Silencing of TODL lncRNA highlighted the enrichment of several key oncogenic signaling pathways including cell cycle, transcriptional misregulation, FOXM1 network, p53 signaling, PLK1 signaling, FoxO, and signaling Aurora signaling pathways. RNA pull-down assay revealed the binding of TODL lncRNA with FOXM1, an oncogenic transcription factor, and the key regulator of the cell cycle. Silencing of TODL lncRNA also induces adipogenesis in dedifferentiated liposarcomas. Altogether, our finding indicates that TODL could be utilized as a novel, specific diagnostic biomarker, and a pharmacological target for therapeutic development in controlling aggressive and metastatic dedifferentiated liposarcomas.

Histone H3K27 demethylase KDM6A is an epigenetic gatekeeper of mTORC1 signalling in cancer

OBJECTIVE

Large-scale genome sequencing efforts of human tumours identified epigenetic modifiers as one of the most frequently mutated gene class in human cancer. However, how these mutations drive tumour development and tumour progression are largely unknown. Here, we investigated the function of the histone demethylase KDM6A in gastrointestinal cancers, such as liver cancer and pancreatic cancer.

DESIGN

Genetic alterations as well as expression analyses of KDM6A were performed in patients with liver cancer. Genetic mouse models of liver and pancreatic cancer coupled with Kdm6a-deficiency were investigated, transcriptomic and epigenetic profiling was performed, and in vivo and in vitro drug treatments were conducted.

RESULTS

KDM6A expression was lost in 30% of patients with liver cancer. Kdm6a deletion significantly accelerated tumour development in murine liver and pancreatic cancer models. Kdm6a-deficient tumours showed hyperactivation of mTORC1 signalling, whereas endogenous Kdm6a re-expression by inducible RNA-interference in established Kdm6a-deficient tumours diminished mTORC1 activity resulting in attenuated tumour progression. Genome-wide transcriptional and epigenetic profiling revealed direct binding of Kdm6a to crucial negative regulators of mTORC1, such as Deptor, and subsequent transcriptional activation by epigenetic remodelling. Moreover, in vitro and in vivo genetic epistasis experiments illustrated a crucial function of Deptor and mTORC1 in Kdm6a-dependent tumour suppression. Importantly, KDM6A expression in human tumours correlates with mTORC1 activity and KDM6A-deficient tumours exhibit increased sensitivity to mTORC1 inhibition.

CONCLUSION

KDM6A is an important tumour suppressor in gastrointestinal cancers and acts as an epigenetic toggle for mTORC1 signalling. Patients with KDM6A-deficient tumours could benefit of targeted therapy focusing on mTORC1 inhibition.

Epigenetic analysis of patients with T-ALL identifies poor outcomes and a hypomethylating agent-responsive subgroup

Adult "T cell" acute lymphoblastic leukemia (T-ALL) is an aggressive hematological malignancy that is associated with poor outcomes, requiring additional therapeutic options. The DNA methylation landscapes of adult T-ALL remain undercharacterized. Here, we systematically analyzed the DNA methylation profiles of normal thymic-sorted T cell subpopulations and 143 primary adult T-ALLs as part of the French GRAALL 2003-2005 trial. Our results indicated that T-ALL is epigenetically heterogeneous consisting of five subtypes (C₁-C₅), which were either associated with co-occurring DNA methyltransferase 3 alpha (DNMT3A)/isocitrate dehydrogenase [NADP(+)] 2 (IDH2) mutations (C₁), TAL bHLH transcription factor 1, erythroid differentiation factor (TAL1) deregulation (C₂), T cell leukemia homeobox 3 (TLX3) (C₃), TLX1/in cis-homeobox A9 (HOXA9) (C₄), or in trans-HOXA9 overexpression (C₅). Integrative analysis of DNA methylation and gene expression identified potential cluster-specific oncogenes and tumor suppressor genes. In addition to an aggressive hypomethylated subgroup (C₁), our data identified an unexpected subset of hypermethylated T-ALL (C₅) associated with poor outcome and primary therapeutic response. Using mouse xenografts, we demonstrated that hypermethylated T-ALL samples exhibited therapeutic responses to the DNA hypomethylating agent 5-azacytidine, which significantly (survival probability; P = 0.001 for C₃, 0.01 for C₄, and 0.0253 for C₅) delayed tumor progression. These findings suggest that epigenetic-based therapies may provide an alternative treatment option in hypermethylated T-ALL.

Abstract 2123: Chromatin profiling of glioblastoma tissues identifies core oncogenic dependency and therapeutic opportunities

Background: Glioblastoma (GBM) is the most aggressive and therapy-refractory brain tumor in adults. Molecular profiling of GBM on the basis of gene expression, DNA methylation and genomic variations has advanced both cancer research and clinical diagnosis. However, the enhancer architectures and regulatory circuitries governing tumor-intrinsic transcriptional diversity and subtype identity are still elusive.

Methods: Chromatin immunoprecipitation followed by sequencing analysis was applied to examine H3K27ac deposition and to map the active regulatory landscapes across 95 GBM biopsies and 12 normal brain tissues. RNA-sequencing analysis was performed to measure transcriptome and to classify transcriptional subtypes. Super-enhancer associated genes and master transcriptional factors were identified. Function of novel cancer associated genes was explored using both patient-derived GBM propagating cells and orthotopic xenograft models.

Results: Analyses of differentially regulated enhancers and super-enhancers uncovered previously unrecognized layers of inter-tumor heterogeneity. Integrative analysis of variant enhancer loci and transcriptome identified topographies of transcriptional enhancers and core regulatory circuitries in four molecular subtypes of primary tumors: AC1-mesenchymal, AC1-classical, AC2-proneural and AC3-proneural. Subtype-specific enhancer domains contributed to transcriptional diversity and shaped subtype identity. Moreover, this study reveals novel oncogenic dependency of GBM on various super-enhancer-driven transcriptional factors and druggable targets (e.g., BRD4).

Conclusion: Through profiling of transcriptional enhancers, we provide clinically relevant insights into the molecular classification, pathogenesis and therapeutic intervention of GBM.

Citation Format: Liang Xu, Ye Chen, Yulun Huang, Edwin Sandanaraj, John S. Yu, Ruby Yu-Tong Lin, Pushkar Dakle, Xin-Yu Ke, Yuk Kien Chong, Lynnette Koh, Anand Mayakonda, Kassoum Nacro, Jeffrey Hill, Mo-Li Huang, Sigal Gery, See Wee Lim, Zhengyun Huang, Ying Xu, Jianxiang Chen, Longchuan Bai, Shaomeng Wang, Hiroaki Wakimoto, Tseng Tsai Yeo, Beng Ti Ang, Markus Müschen, Carol Tang, Tuan Zea Tan, Phillip H. Koeffler. Chromatin profiling of glioblastoma tissues identifies core oncogenic dependency and therapeutic opportunities [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 2123.

Topography of transcriptionally active chromatin in glioblastoma

Molecular profiling of the most aggressive brain tumor glioblastoma (GBM) on the basis of gene expression, DNA methylation, and genomic variations advances both cancer research and clinical diagnosis. The enhancer architectures and regulatory circuitries governing tumor-intrinsic transcriptional diversity and subtype identity are still elusive. Here, by mapping H3K27ac deposition, we analyze the active regulatory landscapes across 95 GBM biopsies, 12 normal brain tissues, and 38 cell line counterparts. Analyses of differentially regulated enhancers and super-enhancers uncovered previously unrecognized layers of intertumor heterogeneity. Integrative analysis of variant enhancer loci and transcriptome identified topographies of transcriptional enhancers and core regulatory circuitries in four molecular subtypes of primary tumors: AC1-mesenchymal, AC1-classical, AC2-proneural, and AC3-proneural. Moreover, this study reveals core oncogenic dependency on super-enhancer-driven transcriptional factors, long noncoding RNAs, and druggable targets in GBM. Through profiling of transcriptional enhancers, we provide clinically relevant insights into molecular classification, pathogenesis, and therapeutic intervention of GBM.

Methrix: an R/bioconductor package for systematic aggregation and analysis of bisulfite sequencing data

MOTIVATION

Whole genome bisulfite sequencing (WGBS), measures DNA methylation at base pair resolution resulting in large bedGraph like coverage files. Current options for processing such files are hindered by discrepancies in file format specification, speed and memory requirements.

RESULTS

We developed methrix, an R package, which provides a toolset for systematic analysis of large datasets. Core functionality of the package includes a comprehensive bedGraph or similar tab-separated text file reader - which summarizes methylation calls based on annotated reference indices, infers and collapses strands, and handles uncovered reference CpG sites while facilitating a flexible input file format specification. Additional optimized functions for quality control filtering, sub-setting, and visualization allow user-friendly and effective processing of WGBS results. Easy integration with tools for differentially methylated region (DMR) calling and annotation further eases the analysis of genome-wide methylation data. Overall, methrix enriches established WGBS workflows by bringing together computational efficiency and versatile functionality.

AVAILABILITY AND IMPLEMENTATION

Methrix is implemented as an R package, made available under MIT license at https://github.com/CompEpigen/methrix and can be installed from the Bioconductor repository.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Globally altered epigenetic landscape and delayed osteogenic differentiation in H3.3-G34W-mutant giant cell tumor of bone

The neoplastic stromal cells of giant cell tumor of bone (GCTB) carry a mutation in H3F3A, leading to a mutant histone variant, H3.3-G34W, as a sole recurrent genetic alteration. We show that in patient-derived stromal cells H3.3-G34W is incorporated into the chromatin and associates with massive epigenetic alterations on the DNA methylation, chromatin accessibility and histone modification level, that can be partially recapitulated in an orthogonal cell line system by the introduction of H3.3-G34W. These epigenetic alterations affect mainly heterochromatic and bivalent regions and provide possible explanations for the genomic instability, as well as the osteolytic phenotype of GCTB. The mutation occurs in differentiating mesenchymal stem cells and associates with an impaired osteogenic differentiation. We propose that the observed epigenetic alterations reflect distinct differentiation stages of H3.3 WT and H3.3 MUT stromal cells and add to H3.3-G34W-associated changes.

The histone variant mutation H3.3-G34W occurs in the majority of giant cell tumor of bone (GCTB). By profiling patient-derived GCTB tumor cells, the authors show that this mutation associates with epigenetic alterations in heterochromatic and bivalent regions that contribute to an impaired osteogenic differentiation and the osteolytic phenotype of GCTB.

TP63, SOX2, and KLF5 Establish a Core Regulatory Circuitry That Controls Epigenetic and Transcription Patterns in Esophageal Squamous Cell Carcinoma Cell Lines

BACKGROUND & AIMS

We investigated the transcriptome of esophageal squamous cell carcinoma (ESCC) cells, activity of gene regulatory (enhancer and promoter regions), and the effects of blocking epigenetic regulatory proteins.

METHODS

We performed chromatin immunoprecipitation sequencing with antibodies against H3K4me1, H3K4me3, and H3K27ac and an assay for transposase-accessible chromatin to map the enhancer regions and accessible chromatin in 8 ESCC cell lines. We used the CRC_Mapper algorithm to identify core regulatory circuitry transcription factors in ESCC cell lines, and determined genome occupancy profiles for 3 of these factors. In ESCC cell lines, expression of transcription factors was knocked down with small hairpin RNAs, promoter and enhancer regions were disrupted by CRISPR/Cas9 genome editing, or bromodomains and extraterminal (BET) family proteins and histone deacetylases (HDACs) were inhibited with ARV-771 and romidepsin, respectively. ESCC cell lines were then analyzed by whole-transcriptome sequencing, immunoprecipitation, immunoblots, immunohistochemistry, and viability assays. Interactions between distal enhancers and promoters were identified and verified with circular chromosome conformation capture sequencing. NOD-SCID mice were given injections of modified ESCC cells, some mice where given injections of HDAC or BET inhibitors, and growth of xenograft tumors was measured.

RESULTS

We identified super-enhancer-regulated circuits and transcription factors TP63, SOX2, and KLF5 as core regulatory factors in ESCC cells. Super-enhancer regulation of ALDH3A1 mediated by core regulatory factors was required for ESCC viability. We observed direct interactions between the promoter region of TP63 and functional enhancers, mediated by the core regulatory circuitry transcription factors. Deletion of enhancer regions from ESCC cells decreased expression of the core regulatory circuitry transcription factors and reduced cell viability; these same results were observed with knockdown of each core regulatory circuitry transcription factor. Incubation of ESCC cells with BET and HDAC disrupted the core regulatory circuitry program and the epigenetic modifications observed in these cells; mice given injections of HDAC or BET inhibitors developed smaller xenograft tumors from the ESCC cell lines. Xenograft tumors grew more slowly in mice given the combination of ARV-771 and romidepsin than mice given either agent alone.

CONCLUSIONS

In epigenetic and transcriptional analyses of ESCC cell lines, we found the transcription factors TP63, SOX2, and KLF5 to be part of a core regulatory network that determines chromatin accessibility, epigenetic modifications, and gene expression patterns in these cells. A combination of epigenetic inhibitors slowed growth of xenograft tumors derived from ESCC cells in mice.

Identification of somatic alterations in lipoma using whole exome sequencing

Lipomas are benign fatty tumors with a high prevalence rate, mostly found in adults but have a good prognosis. Until now, reason for lipoma occurrence not been identified. We performed whole exome sequencing to define the mutational spectrum in ten lipoma patients along with their matching control samples. We presented genomic insight into the development of lipomas, the most common benign tumor of soft tissue. Our analysis identified 412 somatic variants including missense mutations, splice site variants, frameshift indels, and stop gain/lost. Copy number variation analysis highlighted minor aberrations in patients. Kinase genes and transcriptions factors were among the validated mutated genes critical for cell proliferation and survival. Pathway analysis revealed enrichment of calcium, Wnt and phospholipase D signaling in patients. In conclusion, whole exome sequencing in lipomas identified mutations in genes with a possible role in development and progression of lipomas.

RNA-Binding Protein ZFP36L1 Suppresses Hypoxia and Cell-Cycle Signaling

ZFP36L1 is a tandem zinc-finger RNA-binding protein that recognizes conserved adenylate-uridylate-rich elements (ARE) located in 3'untranslated regions (UTR) to mediate mRNA decay. We hypothesized that ZFP36L1 is a negative regulator of a posttranscriptional hub involved in mRNA half-life regulation of cancer-related transcripts. Analysis of in silico data revealed that ZFP36L1 was significantly mutated, epigenetically silenced, and downregulated in a variety of cancers. Forced expression of ZFP36L1 in cancer cells markedly reduced cell proliferation in vitro and in vivo, whereas silencing of ZFP36L1 enhanced tumor cell growth. To identify direct downstream targets of ZFP36L1, systematic screening using RNA pull-down of wild-type and mutant ZFP36L1 as well as whole transcriptome sequencing of bladder cancer cells {plus minus} tet-on ZFP36L1 was performed. A network of 1,410 genes was identified as potential direct targets of ZFP36L1. These targets included a number of key oncogenic transcripts such as HIF1A, CCND1, and E2F1. ZFP36L1 specifically bound to the 3'UTRs of these targets for mRNA degradation, thus suppressing their expression. Dual luciferase reporter assays and RNA electrophoretic mobility shift assays showed that wild-type, but not zinc-finger mutant ZFP36L1, bound to HIF1A 3'UTR and mediated HIF1A mRNA degradation, leading to reduced expression of HIF1A and its downstream targets. Collectively, our findings reveal an indispensable role of ZFP36L1 as a posttranscriptional safeguard against aberrant hypoxic signaling and abnormal cell-cycle progression. SIGNIFICANCE: RNA-binding protein ZFP36L1 functions as a tumor suppressor by regulating the mRNA stability of a number of mRNAs involved in hypoxia and cell-cycle signaling.

Super-enhancer-associated MEIS1 promotes transcriptional dysregulation in Ewing sarcoma in co-operation with EWS-FLI1

As the second most common malignant bone tumor in children and adolescents, Ewing sarcoma is initiated and exacerbated by a chimeric oncoprotein, most commonly, EWS-FLI1. In this study, we apply epigenomic analysis to characterize the transcription dysregulation in this cancer, focusing on the investigation of super-enhancer and its associated transcriptional regulatory mechanisms. We demonstrate that super-enhancer-associated transcripts are significantly enriched in EWS-FLI1 target genes, contribute to the aberrant transcriptional network of the disease, and mediate the exceptional sensitivity of Ewing sarcoma to transcriptional inhibition. Through integrative analysis, we identify MEIS1 as a super-enhancer-driven oncogene, which co-operates with EWS-FLI1 in transcriptional regulation, and plays a key pro-survival role in Ewing sarcoma. Moreover, APCDD1, another super-enhancer-associated gene, acting as a downstream target of both MEIS1 and EWS-FLI1, is also characterized as a novel tumor-promoting factor in this malignancy. These data delineate super-enhancer-mediated transcriptional deregulation in Ewing sarcoma, and uncover numerous candidate oncogenes which can be exploited for further understanding of the molecular pathogenesis for this disease.

Chromatin remodeling mediated by ARID1A is indispensable for normal hematopoiesis in mice

Precise regulation of chromatin architecture is vital to physiological processes including hematopoiesis. ARID1A is a core component of the mammalian SWI/SNF complex, which is one of the ATP-dependent chromatin remodeling complexes. To uncover the role of ARID1A in hematopoietic development, we utilized hematopoietic cell-specific deletion of Arid1a in mice. We demonstrate that ARID1A is essential for maintaining the frequency and function of hematopoietic stem cells and its loss impairs the differentiation of both myeloid and lymphoid lineages. ARID1A deficiency led to a global reduction in open chromatin and ensuing transcriptional changes affected key genes involved in hematopoietic development. We also observed that silencing of ARID1A affected ATRA-induced differentiation of NB4 cells, suggesting its role in granulocytic differentiation of human leukemic cells. Overall, our study provides a comprehensive elucidation of the function of ARID1A in hematopoiesis and highlights the central role of ARID1A-containing SWI/SNF complex in maintaining chromatin dynamics in hematopoietic cells.

Maftools: efficient and comprehensive analysis of somatic variants in cancer

Numerous large-scale genomic studies of matched tumor-normal samples have established the somatic landscapes of most cancer types. However, the downstream analysis of data from somatic mutations entails a number of computational and statistical approaches, requiring usage of independent software and numerous tools. Here, we describe an R Bioconductor package, Maftools, which offers a multitude of analysis and visualization modules that are commonly used in cancer genomic studies, including driver gene identification, pathway, signature, enrichment, and association analyses. Maftools only requires somatic variants in Mutation Annotation Format (MAF) and is independent of larger alignment files. With the implementation of well-established statistical and computational methods, Maftools facilitates data-driven research and comparative analysis to discover novel results from publicly available data sets. In the present study, using three of the well-annotated cohorts from The Cancer Genome Atlas (TCGA), we describe the application of Maftools to reproduce known results. More importantly, we show that Maftools can also be used to uncover novel findings through integrative analysis.

Profiling the B/T cell receptor repertoire of lymphocyte derived cell lines

BACKGROUND

Clonal VDJ rearrangement of B/T cell receptors (B/TCRs) occurring during B/T lymphocyte development has been used as a marker to track the clonality of B/T cell populations.

METHODS

We systematically profiled the B/T cell receptor repertoire of 936 cancer cell lines across a variety of cancer types as well as 462 Epstein-Barr Virus (EBV) transformed normal B lymphocyte lines using RNA sequencing data.

RESULTS

Rearranged B/TCRs were readily detected in cell lines derived from lymphocytes, and subclonality or potential biclonality were found in a number of blood cancer cell lines. Clonal BCR/TCR rearrangements were detected in several blast phase CML lines and unexpectedly, one gastric cancer cell line (KE-97), reflecting a lymphoid origin of these cells. Notably, clonality was highly prevalent in EBV transformed B lymphocytes, suggesting either transformation only occurred in a few B cells or those with a growth advantage dominated the transformed population through clonal evolution.

CONCLUSIONS

Our analysis reveals the complexity and heterogeneity of the BCR/TCR rearrangement repertoire and provides a unique insight into the clonality of lymphocyte derived cell lines.

Identification of distinct mutational patterns and new driver genes in oesophageal squamous cell carcinomas and adenocarcinomas

OBJECTIVES

Oesophageal squamous cell carcinoma (OSCC) and adenocarcinoma (OAC) are distinct cancers in terms of a number of clinical and epidemiological characteristics, complicating the design of clinical trials and biomarker developments. We analysed 1048 oesophageal tumour-germline pairs from both subtypes, to characterise their genomic features, and biological and clinical significance.

DESIGN

Previously exome-sequenced samples were re-analysed to identify significantly mutated genes (SMGs) and mutational signatures. The biological functions of novel SMGs were investigated using cell line and xenograft models. We further performed whole-genome bisulfite sequencing and chromatin immunoprecipitation (ChIP)-seq to characterise epigenetic alterations.

RESULTS

OSCC and OAC displayed nearly mutually exclusive sets of driver genes, indicating that they follow independent developmental paths. The combined sample size allowed the statistical identification of a number of novel subtype-specific SMGs, mutational signatures and prognostic biomarkers. Particularly, we identified a novel mutational signature similar to Catalogue Of Somatic Mutations In Cancer (COSMIC)signature 16, which has prognostic value in OSCC. Two newly discovered SMGs, CUL3 and ZFP36L2, were validated as important tumour-suppressors specific to the OSCC subtype. We further identified their additional loss-of-function mechanisms. CUL3 was homozygously deleted specifically in OSCC and other squamous cell cancers (SCCs). Notably, ZFP36L2 is associated with super-enhancer in healthy oesophageal mucosa; DNA hypermethylation in its super-enhancer reduced active histone markers in squamous cancer cells, suggesting an epigenetic inactivation of a super-enhancer-associated SCC suppressor.

CONCLUSIONS

These data comprehensively contrast differences between OSCC and OAC at both genomic and epigenomic levels, and reveal novel molecular features for further delineating the pathophysiological mechanisms and treatment strategies for these cancers.

ASXL2 regulates hematopoiesis in mice and its deficiency promotes myeloid expansion

Chromosomal translocation t(8;21)(q22;q22) which leads to the generation of oncogenic RUNX1-RUNX1T1 (AML1-ETO) fusion is observed in approximately 10% of acute myelogenous leukemia (AML). To identify somatic mutations that co-operate with t(8;21)-driven leukemia, we performed whole and targeted exome sequencing of an Asian cohort at diagnosis and relapse. We identified high frequency of truncating alterations in ASXL2 along with recurrent mutations of KIT, TET2, MGA, FLT3, and DHX15 in this subtype of AML. To investigate in depth the role of ASXL2 in normal hematopoiesis, we utilized a mouse model of ASXL2 deficiency. Loss of ASXL2 caused progressive hematopoietic defects characterized by myeloid hyperplasia, splenomegaly, extramedullary hematopoiesis, and poor reconstitution ability in transplantation models. Parallel analyses of young and >1-year old Asxl2-deficient mice revealed age-dependent perturbations affecting, not only myeloid and erythroid differentiation, but also maturation of lymphoid cells. Overall, these findings establish a critical role for ASXL2 in maintaining steady state hematopoiesis, and provide insights into how its loss primes the expansion of myeloid cells.

The c-MYC-BMI1 axis is essential for SETDB1-mediated breast tumourigenesis

Characterising the activated oncogenic signalling that leads to advanced breast cancer is of clinical importance. Here, we showed that SET domain, bifurcated 1 (SETDB1), a histone H3 lysine 9 methyltransferase, is aberrantly expressed and behaves as an oncogenic driver in breast cancer. SETDB1 enhances c-MYC and cyclin D1 expression by promoting the internal ribosome entry site (IRES)-mediated translation of MYC/CCND1 mRNA, resulting in prominent signalling of c-MYC to promote cell cycle progression, and provides a growth/self-renewal advantage to breast cancer cells. The activated c-MYC-BMI1 axis is essential for SETDB1-mediated breast tumourigenesis, because silencing of either c-MYC or BMI1 profoundly impairs the enhanced growth/colony formation conferred by SETDB1. Furthermore, c-MYC directly binds to the SETDB1 promoter region and enhances its transcription, suggesting a positive regulatory interplay between SETDB1 and c-MYC. In this study, we identified SETDB1 as a prominent oncogene and characterised the underlying mechanism whereby SETDB1 drives breast cancer, providing a therapeutic rationale for targeting SETDB1-BMI1 signalling in breast cancer. Copyright © 2018 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Super-Enhancer-Driven Long Non-Coding RNA LINC01503, Regulated by TP63, Is Over-Expressed and Oncogenic in Squamous Cell Carcinoma

BACKGROUND & AIMS

Long non-coding RNAs (lncRNAs) are expressed in tissue-specific pattern, but it is not clear how these are regulated. We aimed to identify squamous cell carcinoma (SCC)-specific lncRNAs and investigate mechanisms that control their expression and function.

METHODS

We studied expression patterns and functions of 4 SCC-specific lncRNAs. We obtained 113 esophageal SCC (ESCC) and matched non-tumor esophageal tissues from a hospital in Shantou City, China, and performed quantitative reverse transcription polymerase chain reaction assays to measure expression levels of LINC01503. We collected clinical data from patients and compared expression levels with survival times. LINC01503 was knocked down using small interfering RNAs and oligonucleotides in TE7, TE5, and KYSE510 cell lines and overexpressed in KYSE30 cells. Cells were analyzed by chromatin immunoprecipitation sequencing, luciferase reporter assays, colony formation, migration and invasion, and mass spectrometry analyses. Cells were injected into nude mice and growth of xenograft tumors was measured. LINC01503 interaction with proteins was studied using fluorescence in situ hybridization, RNA pulldown, and RNA immunoprecipitation analyses.

RESULTS

We identified a lncRNA, LINC01503, which is regulated by a super enhancer and is expressed at significantly higher levels in esophageal and head and neck SCCs than in non-tumor tissues. High levels in SCCs correlated with shorter survival times of patients. The transcription factor TP63 bound to the super enhancer at the LINC01503 locus and activated its transcription. Expression of LINC01503 in ESCC cell lines increased their proliferation, colony formation, migration, and invasion. Knockdown of LINC01503 in SCC cells reduced their proliferation, colony formation, migration, and invasion, and the growth of xenograft tumors in nude mice. Expression of LINC01503 in ESCC cell lines reduced ERK2 dephosphorylation by DUSP6, leading to activation of ERK signaling via MAPK. LINC01503 disrupted the interaction between EBP1 and the p85 subunit of PI3K, increasing AKT signaling.

CONCLUSIONS

We identified an lncRNA, LINC01503, which is increased in SCC cells compared with non-tumor cells. Increased expression of LINC01503 promotes ESCC cell proliferation, migration, invasion, and growth of xenograft tumors. It might be developed as a biomarker of aggressive SCCs in patients.

CARD10, a CEBPE target involved in granulocytic differentiation

Maturation of granulocytes is dependent on controlled gene expression by myeloid lineage restricted transcription factors. CEBPE is one of the essential transcription factors required for granulocytic differentiation. Identification of downstream targets of CEBPE is vital to understand better its role in terminal granulopoiesis. In this study, we have identified Card10 as a novel target of CEBPE. We show that CEBPE binds to regulatory elements upstream of the murine Card10 locus, and expression of CARD10 is significantly reduced in Cebpe knock-out mice. Silencing Card10 in a human cell line and in murine primary cells impaired granulopoiesis, affecting expression of genes involved in myeloid cell development and function. Taken together, our data demonstrate for the first time that Card10 is expressed in granulocytes and is a direct target of CEBPE with functions extending to myeloid differentiation.

Super-Enhancers Promote Transcriptional Dysregulation in Nasopharyngeal Carcinoma

Nasopharyngeal carcinoma (NPC) is an invasive cancer with particularly high incidence in Southeast Asia and Southern China. The pathogenic mechanisms of NPC, particularly those involving epigenetic dysregulation, remain largely elusive, hampering clinical management of this malignancy. To identify novel druggable targets, we carried out an unbiased high-throughput chemical screening and observed that NPC cells were highly sensitive to inhibitors of cyclin-dependent kinases (CDK), especially THZ1, a covalent inhibitor of CDK7. THZ1 demonstrated pronounced antineoplastic activities both in vitro and in vivo An integrative analysis using both whole-transcriptome sequencing and chromatin immunoprecipitation sequencing pinpointed oncogenic transcriptional amplification mediated by super-enhancers (SE) as a key mechanism underlying the vulnerability of NPC cells to THZ1 treatment. Further characterization of SE-mediated networks identified many novel SE-associated oncogenic transcripts, such as BCAR1, F3, LDLR, TBC1D2, and the long noncoding RNA TP53TG1. These transcripts were highly and specifically expressed in NPC and functionally promoted NPC malignant phenotypes. Moreover, DNA-binding motif analysis within the SE segments suggest that several transcription factors (including ETS2, MAFK, and TEAD1) may help establish and maintain SE activity across the genome. Taken together, our data establish the landscape of SE-associated oncogenic transcriptional network in NPC, which can be exploited for the development of more effective therapeutic regimens for this disease. Cancer Res; 77(23); 6614-26. ©2017 AACR.

Selinexor (KPT-330) has antitumor activity against anaplastic thyroid carcinoma in vitro and in vivo and enhances sensitivity to doxorubicin

Anaplastic thyroid carcinoma (ATC) is one of the most lethal malignancies having no effective treatment. Exportin-1 (XPO1) is the key mediator of nuclear export of many tumor suppressor proteins and is overexpressed in human cancers. In this study, we examined the therapeutic potential of selinexor (XPO1 inhibitor) against human ATC cells both in vitro and in vivo. Here, we showed that XPO1 is robustly expressed in primary ATC samples and human ATC cell lines. Silencing of XPO1 by either shRNA or selinexor significantly reduced cellular growth and induced cell cycle arrest, apoptosis of ATC cells by altering the protein expression of cancer-related genes. Moreover, selinexor significantly inhibited tumor growth of ATC xenografts. Microarray analysis showed enrichment of DNA replication, cell cycle, cell cycle checkpoint and TNF pathways in selinexor treated ATC cells. Importantly, selinexor decreased AXL and GAS6 levels in CAL62 and HTH83 cells and suppressed the phosphorylation of downstream targets of AXL signaling such as AKT and P70S6K. Finally, a combination of selinexor with doxorubicin demonstrated a synergistic decrease in the cellular proliferation of several ATC cells. These results provide a rationale for investigating the efficacy of combining selinexor and doxorubicin therapy to improve the outcome of ATC patients.

Targeting super-enhancer-associated oncogenes in oesophageal squamous cell carcinoma

OBJECTIVES

Oesophageal squamous cell carcinoma (OSCC) is an aggressive malignancy and the major histological subtype of oesophageal cancer. Although recent large-scale genomic analysis has improved the description of the genetic abnormalities of OSCC, few targetable genomic lesions have been identified, and no molecular therapy is available. This study aims to identify druggable candidates in this tumour.

DESIGN

High-throughput small-molecule inhibitor screening was performed to identify potent anti-OSCC compounds. Whole-transcriptome sequencing (RNA-Seq) and chromatin immunoprecipitation sequencing (ChIP-Seq) were conducted to decipher the mechanisms of action of CDK7 inhibition in OSCC. A variety of in vitro and in vivo cellular assays were performed to determine the effects of candidate genes on OSCC malignant phenotypes.

RESULTS

The unbiased high-throughput small-molecule inhibitor screening led us to discover a highly potent anti-OSCC compound, THZ1, a specific CDK7 inhibitor. RNA-Seq revealed that low-dose THZ1 treatment caused selective inhibition of a number of oncogenic transcripts. Notably, further characterisation of the genomic features of these THZ1-sensitive transcripts demonstrated that they were frequently associated with super-enhancer (SE). Moreover, SE analysis alone uncovered many OSCC lineage-specific master regulators. Finally, integrative analysis of both THZ1-sensitive and SE-associated transcripts identified a number of novel OSCC oncogenes, including PAK4, RUNX1, DNAJB1, SREBF2 and YAP1, with PAK4 being a potential druggable kinase.

CONCLUSIONS

Our integrative approaches led to a catalogue of SE-associated master regulators and oncogenic transcripts, which may significantly promote both the understanding of OSCC biology and the development of more innovative therapies.

Abstract 2450: Mutational profiling of MLL-PTD acute myeloid leukemia

In this study, we performed whole-exome and targeted sequencing on 85 MLL-PTD AML patients. These AMLs have oncogenic tandem duplication of the MLL gene. At least one well-known oncogenic driver mutation was identified in over 90% of the MLL-PTD patients. In line with earlier sequencing studies of other AML subtypes and the TCGA-AML-sequencing project, DNMT3A was the most often mutated epigenetic regulator (25%); IDH1/2 hotspot mutations were identified in 31% of patients. TET family was the third most prominently mutated epigenetic regulator (TET1 (5%), TET2 (16.3%). Mutations of epigenetic regulators also occurred in polycomb-associated proteins (EZH2, ASXL family members), chromatin remodelers (ARID2, ARID1A), genes associated with histone acetylation (CREBBP, EP300, KAT6A, KAT6B) and histone methylation (MLL2, MLL3).

Proliferation-related pathway was extensively mutated, with 54 of 80 MLL-PTD patients (67.5%) carrying at least one mutation of proliferative genes. Specifically, FLT3 mutations were found in 46% of patient samples. Notably, some FLT3-ITD patients had more than one type of internal tandem duplication (ITD) insertion, probably reflecting existence of multiple subclones in these leukemias.

We found highly prevalent mutations of cohesin genes: STAG2 (16%), SMC1A (6%), SMC3 (1%), RAD21 (1%) and CTCF (6%). Cohesin pathway is more frequently mutated in MLL-PTD patients (26%) than the AML samples from either TCGA (13%) or a meta-analysis of 1000 AML (9.1%). Remarkably, an extremely high proportion of the mutations had a strong tendency to disrupt the coding sequence in STAG2, emphasizing their crucial tumor-suppressor role in this AML subtype (16% in MLL-PTD vs 3% in TCGA-AML.

RNA processing pathway was also strikingly altered in MLL-PTD patients. The most prominently mutated genes within this category were the splicing factors. They included U2AF1 (13%, S34F/Y), SRSF2 (3%), SF3A1 (5%), ZRSR2 (3%), DHX15 (1%) and CWC22 (1%).

Multiple mutations co-occur with MLL-PTD which are usually acquired in a sequential manner. A potential ordering for acquisition of many mutations include IDH2/DNMT3A/U2AF1/TET2→MLL-PTD→RAS-receptor tyrosine kinase based on the following reasons: #1, real-time-PCR showed that MLL-PTD was absent in remission while mutations of IDH2, DNMT3A, TET2 and U2AF1 were still retained with a high VAF. This suggests that MLL-PTD was acquired after mutations of IDH2, DNMT3A, TET2 and U2AF1; #2, MLL-PTD is highly stable during disease progression as compared with mutations of the RAS-RTK. On the other hand, RAS-RTK mutations frequently exist as subclonal mutations and tend to be unstable during disease progression. These observations support a notion that MLL-PTD was acquired prior to RAS-RTK. Taken together, MLL-PTD is acquired after those remission-persisting, initiating mutations (IDH2, DNMT3A, TET2 and U2AF1), but prior to lesions of the proliferation-related drivers.

Citation Format: Lingwen Ding, Qiaoyang Sun, Kar-Tong Tan, Wenwen Chien, Anand Mayakonda, Dechen Lin, Xinyi Loh, Jinfen Xiao, Manja Meggendorfer, Tamara Alpermann, Manoj Garg, Su-Lin Lim, Vikas Madan, Norimichi Hattori, Yasunobu Nagata, Satoru Miyano, Allen Yeoh Eng Juh, Hsin-An Hou, Yan-Yi Jiang, Yan-Yi Jiang, Sumiko Takao, Li-Zhen Liu, Siew-Zhuan Tan, Siew-Zhuan Tan, Michael Lill, Mutsumi Hayashi, Akitoshi Kinoshita, Hagop M. Kantarjian, Steven M. Kornblau, Seishi Ogawa, Torsten Haferlach, Henry Yang, H. Phillip Koeffler. Mutational profiling of MLL-PTD acute myeloid leukemia [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 2450. doi:10.1158/1538-7445.AM2017-2450

Abstract 2569: BCL6 promotes glioma and serves as a novel therapeutic target

ZBTB transcription factors orchestrate gene transcription during tissue development. However, their roles in glioblastoma multiforme (GBM) remain unexplored. Here, through a functional screening of ZBTB genes, we identify that BCL6 is required for GBM cell proliferation and is overexpressed in GBM samples. Using a somatic transgenic mouse model, Bcl6 confers the proliferative and invasive stimuli to glioma cells in vivo. Moreover, we discover AXL as a novel downstream target of BCL6. Through AXL, BCL6 enhances both MEK-ERK and S6K-RPS6 axes. Pharmacological inhibition of BCL6 activity effectively blocks GBM growth and inhibits AXL expression. Together, these findings uncover a novel glioma-promoting role of BCL6, and provide the rationale of targeting BCL6 as a potential therapeutic approach.

Citation Format: Liang Xu, Ye CHEN, Marina Dutra-Clarke, Anand Mayakonda, Masaharu Hazawa, Steve E. Savinoff, Ngan Doan, Jonathan W. Said, William H. Yong, Henry Yang, Ling-Wen Ding, Yan-Yi Jiang, Jeffrey W. Tyner, Jianhong Ching, Jean-Paul Kovalik, Markus Müschen, Joshua J. Breunig, De-Chen Lin, Phillip Koeffler. BCL6 promotes glioma and serves as a novel therapeutic target [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 2569. doi:10.1158/1538-7445.AM2017-2569

Abstract 1524: BCL6 modulates the TP53 and STAT pathways in glioma

Glioblastoma multiforme (GBM) remains the most aggressive brain malignancy with little improvement in prognosis or therapy for decades. Recently, we identified BCL6, also known as ZBTB27, to be a novel oncogene in GBM. In this study, we performed IHC analysis of 153 primary human glioma specimens and 8 normal brain samples. BCL6 expression is robustly elevated in tumor samples and positively correlated with glioma pathological grade. High BCL6 expression strongly predicts a worse prognosis of GBM patients. Depletion of BCL6 in human GBM cells reduced the incorporation of BrdU, promoted the cellular senescence and inhibited the growth of human GBM cells in vivo. Next, genome-wide occupancy of BCL6 in GBM cells was characterized by ChIP-seq assay. Genomic regions centered on BCL6 peaks are co-enriched with RNA-Pol II and flanked with strong H3K27ac and H3K4me3 modifications. MYC and two long non-coding RNAs MALAT1 and NEAT1 were identified as novel BCL6 targets in GBM. Moreover, pathway enrichment analysis of BCL6 peak-associated genes reveals a significant enrichment of JAK-STAT, TP53, ERBB and MAPK pathways. We demostrated further that BCL6 represses the TP53 pathway and promotes the JAK-STAT pathway activation in GBM cells. Together, our findings uncover potential downstream targets and provide a better understanding of BCL6 function in GBM.

Citation Format: Ye Chen, Liang Xu, Marina Dutra-Clarke, Anand Mayakonda, De-Chen Lin, Lynnette Koh, Yuk Kien Chong, Edwin Sandanaraj, Vikas Madan, Henry Yang, Ngan Doan, Jonathan W. Said, William H. Yong, Markus Müschen, Beng Ti Ang, Carol Tang, Joshua J. Breunig, Phillip Koeffler. BCL6 modulates the TP53 and STAT pathways in glioma [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 1524. doi:10.1158/1538-7445.AM2017-1524

Abstract 806: MGA is a potential tumor suppressor in acute myeloid leukemia

MGA is an incompletely studied gene with a high mutation frequency in MLL-PTD AML (9%) and in core bind factor AML (8%). This gene encodes a MAX-interacting protein and is believed to act as a transcription factor that suppresses MYC binding to its target. By in silico analysis, we found that MGA is expressed in normal myeloid hematopoietic cells and AML, and the expression level is comparable with TET2 or DNMT3A. Further data mining of TCGA revealed a high frequency of inactivating mutations of the MGA gene in a variety of cancers such as various adenocarcinomas. To interrogate functionally its role in leukemogenesis, lentiviral constructs containing either shRNA or CRISPR-sgRNA targeted to different regions of the MGA gene were generated. MGA expressing AML cell line EOL-1 was silenced by shRNA or CRISPER system. Silencing was confirmed by western blot (shRNA) and Sanger Sequencing (sgRNA). An increase of methylcellulose colony number (~30%) was observed in MGA silenced cell lines. Control EOL-1 cells or EOL-1 cells silenced with MGA CRISPR sgRNAs were injected into both flanks of NSG mice, and tumor masses were harvested 21 days after injection. Silencing of MGA by CRISPR-sgRNA consistently enhanced in vivo xenograft cell growth. In addition, western blot analyses revealed silencing of MGA in EOL-1 cells increased protein levels of Cyclin E1 and phos-RB (S807 phosphorylation inhibits the ability of RB to target protein allowing cell cycle progression), indicative of a proliferative advantage conferred by the silencing of MGA.

MGA may be a potential regulator of the MYC pathway. We, therefore, examined whether silencing of MGA alters MYC transcriptional activity. Luciferase reporter assay was carried out in 293FT cells stabilized with either scramble or shRNA- targeting MGA. Luciferase activities were measured 48 h after transfection of cells with MYC activity reporter pMyc4ElbLuc and normalized to the corresponding co-transfected Renilla luciferase activity. A fourfold increase in luciferase activity was observed in MGA silenced cells when compared with non- targeting shRNA controls. Furthermore, Kaplan–Meier survival analysis was performed in the TCGA-AML patients by comparison of cases with highest versus lowest expression of MGA. P-values were calculated by log-rank test. MGA expression data and patient survival data were retrieved from TCGA-AML patients RNA seq, or microarray (70 AML patients). The MGA expression ‘high’ and ‘low’ groups were defined by 15% higher than the median or 15% lower than the median, respectively. AML patients with lower levels of MGA in their leukemic samples had a worse outcome compared with those whose leukemic cells expressed higher levels of MGA. Collectively, our results suggest that MGA may function as a potential tumor-suppressor in AML.

Citation Format: Qiaoyang Sun, Lingwen Ding, Kar-Tong Tan, Wenwen Chien, Xinyi Loh, Jinfen Xiao, Anand Mayakonda, Dechen Lin, Yanyi Jiang, Henry Yang, Sigal Gery, H. Phillip Koeffler. MGA is a potential tumor suppressor in acute myeloid leukemia [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 806. doi:10.1158/1538-7445.AM2017-806

Mutational profiling of acute lymphoblastic leukemia with testicular relapse

Relapsed acute lymphoblastic leukemia (ALL) is the leading cause of deaths of childhood cancer. Although relapse usually happens in the bone marrow, extramedullary relapse occasionally occurs including either the central nervous system or testis (<1-2%). We selected two pediatric ALL patients who experienced testicular relapse and interrogated their leukemic cells with exome sequencing. The sequencing results and clonality analyses suggest that relapse of patient D483 directly evolved from the leukemic clone at diagnosis which survived chemotherapy. In contrast, relapse leukemia cells (both bone marrow and testis) of patient D727 were likely derived from a common ancestral clone, and testicular relapse likely arose independently from the bone marrow relapsed leukemia. Our findings decipher the mutational spectra and shed light on the clonal evolution of two cases of pediatric ALL with testicular relapse. Presence of CREBBP/NT5C2 mutations suggests that a personalized therapeutic approach should be applied to these two patients.

Genomic and Epigenomic Heterogeneity of Hepatocellular Carcinoma

Understanding the intratumoral heterogeneity of hepatocellular carcinoma is instructive for developing personalized therapy and identifying molecular biomarkers. Here we applied whole-exome sequencing to 69 samples from 11 patients to resolve the genetic architecture of subclonal diversification. Spatial genomic diversity was found in all 11 hepatocellular carcinoma cases, with 29% of driver mutations being heterogeneous, including TERT, ARID1A, NOTCH2, and STAG2. Similar with other cancer types, TP53 mutations were always shared between all tumor regions, that is, located on the "trunk" of the evolutionary tree. In addition, we found that variants within several drug targets such as KIT, SYK, and PIK3CA were mutated in a fully clonal manner, indicating their therapeutic potentials for hepatocellular carcinoma. Temporal dissection of mutational signatures suggested that mutagenic processes associated with exposure to aristolochic acid and aflatoxin might play a more important role in early, as opposed to late, stages of hepatocellular carcinoma development. Moreover, we observed extensive intratumoral epigenetic heterogeneity in hepatocellular carcinoma based on multiple independent analytical methods and showed that intratumoral methylation heterogeneity might play important roles in the biology of hepatocellular carcinoma cells. Our results also demonstrated prominent heterogeneity of intratumoral methylation even in a stable hepatocellular carcinoma genome. Together, these findings highlight widespread intratumoral heterogeneity at both the genomic and epigenomic levels in hepatocellular carcinoma and provide an important molecular foundation for better understanding the pathogenesis of this malignancy. Cancer Res; 77(9); 2255-65. ©2017 AACR.

Molecular mechanism and therapeutic implications of selinexor (KPT-330) in liposarcoma

Exportin-1 mediates nuclear export of multiple tumor suppressor and growth regulatory proteins. Aberrant expression of exportin-1 is noted in human malignancies, resulting in cytoplasmic mislocalization of its target proteins. We investigated the efficacy of selinexor against liposarcoma cells both in vitro and in vivo. Exportin-1 was highly expressed in liposarcoma samples and cell lines as determined by immunohistochemistry, western blot, and immunofluorescence assay. Knockdown of endogenous exportin-1 inhibited proliferation of liposarcoma cells. Selinexor also significantly decreased cell proliferation as well as induced cell cycle arrest and apoptosis of liposarcoma cells. The drug also significantly decreased tumor volumes and weights of liposarcoma xenografts. Importantly, selinexor inhibited insulin-like growth factor 1 (IGF1) activation of IGF-1R/AKT pathway through upregulation of insulin-like growth factor binding protein 5 (IGFBP5). Further, overexpression and knockdown experiments showed that IGFBP5 acts as a tumor suppressor and its expression was restored upon selinexor treatment of liposarcoma cells. Selinexor decreased aurora kinase A and B levels in these cells and inhibitors of these kinases suppressed the growth of the liposarcoma cells. Overall, our study showed that selinexor treatment restored tumor suppressive function of IGFBP5 and inhibited aurora kinase A and B in liposarcoma cells supporting the usefulness of selinexor as a potential therapeutic strategy for the treatment of this cancer.

Mutational Landscape of Pediatric Acute Lymphoblastic Leukemia

Current standard of care for patients with pediatric acute lymphoblastic leukemia (ALL) is mainly effective, with high remission rates after treatment. However, the genetic perturbations that give rise to this disease remain largely undefined, limiting the ability to address resistant tumors or develop less toxic targeted therapies. Here, we report the use of next-generation sequencing to interrogate the genetic and pathogenic mechanisms of 240 pediatric ALL cases with their matched remission samples. Commonly mutated genes fell into several categories, including RAS/receptor tyrosine kinases, epigenetic regulators, transcription factors involved in lineage commitment, and the p53/cell-cycle pathway. Unique recurrent mutational hotspots were observed in epigenetic regulators CREBBP (R1446C/H), WHSC1 (E1099K), and the tyrosine kinase FLT3 (K663R, N676K). The mutant WHSC1 was established as a gain-of-function oncogene, while the epigenetic regulator ARID1A and transcription factor CTCF were functionally identified as potential tumor suppressors. Analysis of 28 diagnosis/relapse trio patients plus 10 relapse cases revealed four evolutionary paths and uncovered the ordering of acquisition of mutations in these patients. This study provides a detailed mutational portrait of pediatric ALL and gives insights into the molecular pathogenesis of this disease. Cancer Res; 77(2); 390-400. ©2016 AACR.

ZNF750 is a lineage-specific tumour suppressor in squamous cell carcinoma

ZNF750 controls epithelial homeostasis by regulating epidermal-differentiation genes, a role underscored by its pathogenic mutations in esophageal squamous cell cancers (SCCs). However, the precise role of ZNF750 in SCC cell biology remains unclear. In this study, we report that ZNF750 is exclusively deleted, mutated and underexpressed in human SCCs, and low ZNF750 expression is associated with poor survival. Restoration of wildtype, but not mutant ZNF750 protein uniquely inhibited the malignant phenotypes of SCC cells both in vitro and in vivo. Notably, ZNF750 promoted the expression of a long non-coding RNA (TINCR), which mediated both cancer-inhibition and differentiation-induction effects of ZNF750. In addition, ZNF750 potently suppressed cell migration by directly inhibiting the transactivation of LAMC2. Together, our findings characterize ZNF750 as a crucial SCC-specific suppressor and uncover its novel anticancer-associated functions.

Genomic landscape of liposarcoma

Liposarcoma (LPS) is the most common type of soft tissue sarcoma accounting for 20% of all adult sarcomas. Due to absence of clinically effective treatment options in inoperable situations and resistance to chemotherapeutics, a critical need exists to identify novel therapeutic targets. We analyzed LPS genomic landscape using SNP arrays, whole exome sequencing and targeted exome sequencing to uncover the genomic information for development of specific anti-cancer targets. SNP array analysis indicated known amplified genes (MDM2, CDK4, HMGA2) and important novel genes (UAP1, MIR557, LAMA4, CPM, IGF2, ERBB3, IGF1R). Carboxypeptidase M (CPM), recurrently amplified gene in well-differentiated/de-differentiated LPS was noted as a putative oncogene involved in the EGFR pathway. Notable deletions were found at chromosome 1p (RUNX3, ARID1A), chromosome 11q (ATM, CHEK1) and chromosome 13q14.2 (MIR15A, MIR16-1). Significantly and recurrently mutated genes (false discovery rate < 0.05) included PLEC (27%), MXRA5 (21%), FAT3 (24%), NF1 (20%), MDC1 (10%), TP53 (7%) and CHEK2 (6%). Further, in vitro and in vivo functional studies provided evidence for the tumor suppressor role for Neurofibromin 1 (NF1) gene in different subtypes of LPS. Pathway analysis of recurrent mutations demonstrated signaling through MAPK, JAK-STAT, Wnt, ErbB, axon guidance, apoptosis, DNA damage repair and cell cycle pathways were involved in liposarcomagenesis. Interestingly, we also found mutational and copy number heterogeneity within a primary LPS tumor signifying the importance of multi-region sequencing for cancer-genome guided therapy. In summary, these findings provide insight into the genomic complexity of LPS and highlight potential druggable pathways for targeted therapeutic approach.

CRM1 Inhibition Promotes Cytotoxicity in Ewing Sarcoma Cells by Repressing EWS-FLI1-Dependent IGF-1 Signaling

Ewing sarcoma (EWS) is an aggressive bone malignancy that mainly affects children and young adults. The mechanisms by which EWS (EWSR1) fusion genes drive the disease are not fully understood. CRM1 (XPO1) traffics proteins from the nucleus, including tumor suppressors and growth factors, and is overexpressed in many cancers. A small-molecule inhibitor of CRM1, KPT-330, has shown therapeutic promise, but has yet to be investigated in the context of EWS. In this study, we demonstrate that CRM1 is also highly expressed in EWS. shRNA-mediated or pharmacologic inhibition of CRM1 in EWS cells dramatically decreased cell growth while inducing apoptosis, cell-cycle arrest, and protein expression alterations to several cancer-related factors. Interestingly, silencing of CRM1 markedly reduced EWS-FLI1 fusion protein expression at the posttranscriptional level and upregulated the expression of the well-established EWS-FLI1 target gene, insulin-like growth factor binding protein 3 (IGFBP3), which inhibits IGF-1. Accordingly, KPT-330 treatment attenuated IGF-1-induced activation of the IGF-1R/AKT pathway. Furthermore, knockdown of IGFBP3 increased cell growth and rescued the inhibitory effects on IGF-1 signaling triggered by CRM1 inhibition. Finally, treatment of EWS cells with a combination of KPT-330 and the IGF-1R inhibitor, linsitinib, synergistically decreased cell proliferation both in vitro and in vivo Taken together, these findings provide a strong rationale for investigating the efficacy of combinatorial inhibition of CRM1 and IGF-1R for the treatment of EWS. Cancer Res; 76(9); 2687-97. ©2016 AACR.

Comprehensive mutational analysis of primary and relapse acute promyelocytic leukemia

Acute promyelocytic leukemia (APL) is a subtype of myeloid leukemia characterized by differentiation block at the promyelocyte stage. Besides the presence of chromosomal rearrangement t(15;17), leading to the formation of PML-RARA (promyelocytic leukemia-retinoic acid receptor alpha) fusion, other genetic alterations have also been implicated in APL. Here, we performed comprehensive mutational analysis of primary and relapse APL to identify somatic alterations, which cooperate with PML-RARA in the pathogenesis of APL. We explored the mutational landscape using whole-exome (n=12) and subsequent targeted sequencing of 398 genes in 153 primary and 69 relapse APL. Both primary and relapse APL harbored an average of eight non-silent somatic mutations per exome. We observed recurrent alterations of FLT3, WT1, NRAS and KRAS in the newly diagnosed APL, whereas mutations in other genes commonly mutated in myeloid leukemia were rarely detected. The molecular signature of APL relapse was characterized by emergence of frequent mutations in PML and RARA genes. Our sequencing data also demonstrates incidence of loss-of-function mutations in previously unidentified genes, ARID1B and ARID1A, both of which encode for key components of the SWI/SNF complex. We show that knockdown of ARID1B in APL cell line, NB4, results in large-scale activation of gene expression and reduced in vitro differentiation potential.

Genomic landscape of liposarcoma

Liposarcoma (LPS) is the most common type of soft tissue sarcoma accounting for 20% of all adult sarcomas. Due to absence of clinically effective treatment options in inoperable situations and resistance to chemotherapeutics, a critical need exists to identify novel therapeutic targets. We analyzed LPS genomic landscape using SNP arrays, whole exome sequencing and targeted exome sequencing to uncover the genomic information for development of specific anti-cancer targets. SNP array analysis indicated known amplified genes (MDM2, CDK4, HMGA2) and important novel genes (UAP1, MIR557, LAMA4, CPM, IGF2, ERBB3, IGF1R). Carboxypeptidase M (CPM), recurrently amplified gene in well-differentiated/de-differentiated LPS was noted as a putative oncogene involved in the EGFR pathway. Notable deletions were found at chromosome 1p (RUNX3, ARID1A), chromosome 11q (ATM, CHEK1) and chromosome 13q14.2 (MIR15A, MIR16-1). Significantly and recurrently mutated genes (false discovery rate < 0.05) included PLEC (27%), MXRA5 (21%), FAT3 (24%), NF1 (20%), MDC1 (10%), TP53 (7%) and CHEK2 (6%). Further, in vitro and in vivo functional studies provided evidence for the tumor suppressor role for Neurofibromin 1 (NF1) gene in different subtypes of LPS. Pathway analysis of recurrent mutations demonstrated signaling through MAPK, JAK-STAT, Wnt, ErbB, axon guidance, apoptosis, DNA damage repair and cell cycle pathways were involved in liposarcomagenesis. Interestingly, we also found mutational and copy number heterogeneity within a primary LPS tumor signifying the importance of multi-region sequencing for cancer-genome guided therapy. In summary, these findings provide insight into the genomic complexity of LPS and highlight potential druggable pathways for targeted therapeutic approach.