Anand Mayakonda

Profiling of somatic mutations in acute myeloid leukemia with FLT3-ITD at diagnosis and relapse

Acute myeloid leukemia (AML) with an FLT3 internal tandem duplication (FLT3-ITD) mutation is an aggressive hematologic malignancy with a grave prognosis. To identify the mutational spectrum associated with relapse, whole-exome sequencing was performed on 13 matched diagnosis, relapse, and remission trios followed by targeted sequencing of 299 genes in 67 FLT3-ITD patients. The FLT3-ITD genome has an average of 13 mutations per sample, similar to other AML subtypes, which is a low mutation rate compared with that in solid tumors. Recurrent mutations occur in genes related to DNA methylation, chromatin, histone methylation, myeloid transcription factors, signaling, adhesion, cohesin complex, and the spliceosome. Their pattern of mutual exclusivity and cooperation among mutated genes suggests that these genes have a strong biological relationship. In addition, we identified mutations in previously unappreciated genes such as MLL3, NSD1, FAT1, FAT4, and IDH3B. Mutations in 9 genes were observed in the relapse-specific phase. DNMT3A mutations are the most stable mutations, and this DNMT3A-transformed clone can be present even in morphologic complete remissions. Of note, all AML matched trio samples shared at least 1 genomic alteration at diagnosis and relapse, suggesting common ancestral clones. Two types of clonal evolution occur at relapse: either the founder clone recurs or a subclone of the founder clone escapes from induction chemotherapy and expands at relapse by acquiring new mutations. Relapse-specific mutations displayed an increase in transversions. Functional assays demonstrated that both MLL3 and FAT1 exert tumor-suppressor activity in the FLT3-ITD subtype. An inhibitor of XPO1 synergized with standard AML induction chemotherapy to inhibit FLT3-ITD growth. This study clearly shows that FLT3-ITD AML requires additional driver genetic alterations in addition to FLT3-ITD alone.

Spatial intratumoral heterogeneity and temporal clonal evolution in esophageal squamous cell carcinoma.

Esophageal squamous cell carcinoma (ESCC) is among the most common malignancies, but little is known about its spatial intratumoral heterogeneity (ITH) and temporal clonal evolutionary processes. To address this, we performed multiregion whole-exome sequencing on 51 tumor regions from 13 ESCC cases and multiregion global methylation profiling for 3 of these 13 cases. We found an average of 35.8% heterogeneous somatic mutations with strong evidence of ITH. Half of the driver mutations located on the branches of tumor phylogenetic trees targeted oncogenes, including PIK3CA, NFE2L2 and MTOR, among others. By contrast, the majority of truncal and clonal driver mutations occurred in tumor-suppressor genes, including TP53, KMT2D and ZNF750, among others. Interestingly, phyloepigenetic trees robustly recapitulated the topological structures of the phylogenetic trees, indicating a possible relationship between genetic and epigenetic alterations. Our integrated investigations of spatial ITH and clonal evolution provide an important molecular foundation for enhanced understanding of tumorigenesis and progression in ESCC.

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.

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.

Ordering of mutations in acute myeloid leukemia with partial tandem duplication of MLL (MLL-PTD)

Partial tandem duplication of MLL (MLL-PTD) characterizes acute myeloid leukemia (AML) patients often with a poor prognosis. To understand the order of occurrence of MLL-PTD in relation to other major AML mutations and to identify novel mutations that may be present in this unique AML molecular subtype, exome and targeted sequencing was performed on 85 MLL-PTD AML samples using HiSeq-2000. Genes involved in the cohesin complex (STAG2), a splicing factor (U2AF1) and a poorly studied gene, MGA were recurrently mutated, whereas NPM1, one of the most frequently mutated AML gene, was not mutated in MLL-PTD patients. Interestingly, clonality analysis suggests that IDH2/1, DNMT3A, U2AF1 and TET2 mutations are clonal and occur early, and MLL-PTD likely arises after these initial mutations. Conversely, proliferative mutations (FLT3, RAS), typically appear later, are largely subclonal and tend to be unstable. This study provides important insights for understanding the relative importance of different mutations for defining a targeted therapeutic strategy for MLL-PTD AML patients.

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.

Maftools: Efficient analysis, visualization and summarization of MAF files from large-scale cohort based cancer studies.

Mutation Annotation Format (MAF) has become a standard file format for storing somatic/germline variants derived from sequencing of large cohort of cancer samples. MAF files contain a list of all variants detected in a sample along with various annotations associated with the putative variant. MAF file forms the basis for many downstream analyses and provides complete landscape of the cohort. Here we introduce maftools–an R package that provides rich source of functions for performing various analyses, visualizations and summarization of MAF files. Maftools uses data.table library for faster processing/summarization and ggplot2 for generating rich and publication quality visualizations. Maftools also takes advantages of S4 class system for better data representation, with easy to use and flexible functions. Availability and Implementation maftools is implemented as an R package available at https://github.com/PoisonAlien/maftools Contact csiamt@nus.edu.sg

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.

BCL6 promotes glioma and serves as a therapeutic target

Significance Glioblastoma (GBM) is the most lethal brain malignancy lacking effective treatment. In this study, we demonstrate that BCL6 is a prognostic marker and a targetable GBM-promoting factor. Silencing of BCL6 inhibits the malignant phenotype of GBM cells and triggers cellular senescence. We also identify AXL as an important BCL6 transcriptional target, the expression of which is also regulated positively by NCoR, a BCL6 cofactor. Either silencing of BCL6 or targeted disruption of the BCL6/NCoR complex diminishes AXL expression and inhibits GBM growth. This study elucidates a crucial BCL6-mediated signaling pathway in GBM biology. More importantly, our results highlight the promise and merit of targeting BCL6 for treating this deadly disease. ZBTB transcription factors orchestrate gene transcription during tissue development. However, their roles in glioblastoma (GBM) remain unexplored. Here, through a functional screening of ZBTB genes, we identify that BCL6 is required for GBM cell viability and that BCL6 overexpression is associated with worse prognosis. In a somatic transgenic mouse model, depletion of Bcl6 inhibits the progression of KrasG12V-driven high-grade glioma. Transcriptome analysis demonstrates the involvement of BCL6 in tumor protein p53 (TP53), erythroblastic leukemia viral oncogene homolog (ErbB), and MAPK signaling pathways. Indeed, BCL6 represses the expression of wild-type p53 and its target genes in GBM cells. Knockdown of BCL6 augments the activation of TP53 pathway in response to radiation. Importantly, we discover that receptor tyrosine kinase AXL is a transcriptional target of BCL6 in GBM and mediates partially the regulatory effects of BCL6 on both MEK-ERK (mitogen-activated protein/extracellular signal-regulated kinase kinase–extracellular signal-regulated kinase) and S6K-RPS6 (ribosomal protein S6 kinase–ribosomal protein S6) axes. Similar to BCL6 silencing, depletion of AXL profoundly attenuates GBM proliferation both in vitro and in vivo. Moreover, targeted inhibition of BCL6/nuclear receptor corepressor 1 (NCoR) complex by peptidomimetic inhibitor not only significantly decreases AXL expression and the activity of MEK-ERK and S6K-RPS6 cascades but also displays a potent antiproliferative effect against GBM cells. Together, these findings uncover a glioma-promoting role of BCL6 and provide the rationale of targeting BCL6 as a potential therapeutic approach.