TAL1 activation in T-cell acute lymphoblastic leukemia: a novel oncogenic 3' neo-enhancer.

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 cis-alterations, including SIL-TAL1 microdeletions, translocations with T-cell Receptor loci, and more recently described upstream intergenic non-coding mutations. These mutations consist of recurrent focal microinsertions that create an oncogenic neo-enhancer accompanied by activating epigenetic marks. This observation laid the groundwork for an innovative paradigm concerning the activation of proto-oncogenes via genomic alterations of non-coding intergenic regions. However, for the majority of T-ALL expressing TAL1 (TAL1+), the deregulation mechanism remains 'unresolved'. We took advantage of H3K27ac and H3K4me3 chromatin immunoprecipitation sequencing data of eight cases of T-ALL, including five TAL1+ cases. We 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. 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.

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.

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.

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, subsetting 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.

Transcriptome analysis identifies TODL as a novel lncRNA associated with proliferation, differentiation, and tumorigenesis in liposarcoma through FOXM1.

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.

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.

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.

Identification of a novel enhancer of CEBPE essential for granulocytic differentiation.

CCAAT/enhancer binding protein ε (CEBPE) is an essential transcription factor for granulocytic differentiation. Mutations of CEBPE occur in individuals with neutrophil-specific granule deficiency (SGD), which is characterized by defects in neutrophil maturation. Cebpe-knockout mice also exhibit defects in terminal differentiation of granulocytes, a phenotype reminiscent of SGD. Analysis of DNase I hypersensitive sites sequencing data revealed an open chromatin region 6 kb downstream of the transcriptional start site of Cebpe in murine myeloid cells. We identified an interaction between this +6-kb region and the core promoter of Cebpe using circular chromosome conformation capture sequencing (4C-seq). To understand the role of this putative enhancer in transcriptional regulation of Cebpe, we targeted it using catalytically inactive Cas9 fused to Krüppel-associated box (KRAB) domain and observed a significant downregulation of transcript and protein levels of CEBPE in cells expressing guide RNA targeting the +6-kb region. To further investigate the role of this novel enhancer further in myelopoiesis, we generated mice with deletion of this region using CRISPR/Cas9 technology. Germline deletion of the +6-kb enhancer resulted in reduced levels of CEBPE and its target genes and caused a severe block in granulocytic differentiation. We also identified binding of CEBPA and CEBPE to the +6-kb enhancer, which suggests their role in regulating the expression of Cebpe In summary, we have identified a novel enhancer crucial for regulating expression of Cebpe and required for normal granulocytic differentiation.

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.

Targetable BET proteins- and E2F1-dependent transcriptional program maintains the malignancy of glioblastoma.

Competitive BET bromodomain inhibitors (BBIs) targeting BET proteins (BRD2, BRD3, BRD4, and BRDT) show promising preclinical activities against brain cancers. However, the BET protein-dependent glioblastoma (GBM)-promoting transcriptional network remains elusive. Here, with mechanistic exploration of a next-generation chemical degrader of BET proteins (dBET6), we reveal a profound and consistent impact of BET proteins on E2F1- dependent transcriptional program in both differentiated GBM cells and brain tumor-initiating cells. dBET6 treatment drastically reduces BET protein genomic occupancy, RNA-Pol2 activity, and permissive chromatin marks. Subsequently, dBET6 represses the proliferation, self-renewal, and tumorigenic ability of GBM cells. Moreover, dBET6-induced degradation of BET proteins exerts superior antiproliferation effects compared to conventional BBIs and overcomes both intrinsic and acquired resistance to BBIs in GBM cells. Our study reveals crucial functions of BET proteins and provides the rationale and therapeutic merits of targeted degradation of BET proteins in GBM.

BCL6 promotes glioma and serves as a therapeutic target.

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.

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.

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.

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.

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.

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.

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.

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.

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.