ABSTRACT
Advanced computational methods exploit gene expression and epigenetic datasets to predict gene regulatory networks controlled by transcription factors (TFs). These methods have identified cell fate determining TFs but require large amounts of reference data and experimental expertise. Here, we present an easy to use network-based computational framework that exploits enhancers defined by bidirectional transcription, using as sole input CAGE sequencing data to correctly predict TFs key to various human cell types. Next, we applied this Analysis Algorithm for Networks Specified by Enhancers based on CAGE (ANANSE-CAGE) to predict TFs driving red and white blood cell development, and THP-1 leukemia cell immortalization. Further, we predicted TFs that are differentially important to either cell line- or primary- associated MLL-AF9-driven gene programs, and in primary MLL-AF9 acute leukemia. Our approach identified experimentally validated as well as thus far unexplored TFs in these processes. ANANSE-CAGE will be useful to identify transcription factors that are key to any cell fate change using only CAGE-seq data as input.
Subject(s)
Gene Regulatory Networks , Leukemia, Myeloid, Acute , Humans , Transcription Factors/genetics , Transcription Factors/metabolism , Leukemia, Myeloid, Acute/genetics , Algorithms , Blood Cells/metabolism , Computational BiologyABSTRACT
Overexpression of the BRE (brain and reproductive organ-expressed) gene defines a distinct pediatric and adult acute myeloid leukemia (AML) subgroup. Here we identify a promoter enriched for active chromatin marks in BRE intron 4 causing strong biallelic expression of a previously unknown C-terminal BRE transcript. This transcript starts with BRE intron 4 sequences spliced to exon 5 and downstream sequences, and if translated might code for an N terminally truncated BRE protein. Remarkably, the new BRE transcript was highly expressed in over 50% of 11q23/KMT2A (lysine methyl transferase 2A)-rearranged and t(8;16)/KAT6A-CREBBP cases, while it was virtually absent from other AML subsets and normal tissues. In gene reporter assays, the leukemia-specific fusion protein KMT2A-MLLT3 transactivated the intragenic BRE promoter. Further epigenome analyses revealed 97 additional intragenic promoter marks frequently bound by KMT2A in AML with C-terminal BRE expression. The corresponding genes may be part of a context-dependent KMT2A-MLLT3-driven oncogenic program, because they were higher expressed in this AML subtype compared with other groups. C-terminal BRE might be an important contributor to this program because in a case with relapsed AML, we observed an ins(11;2) fusing CHORDC1 to BRE at the region where intragenic transcription starts in KMT2A-rearranged and KAT6A-CREBBP AML.
Subject(s)
Gene Rearrangement , Leukemia, Myeloid, Acute/genetics , Nerve Tissue Proteins/genetics , Protein Interaction Domains and Motifs/genetics , Transcriptional Activation , Translocation, Genetic , Cell Line , Chromosomes, Human, Pair 11 , Chromosomes, Human, Pair 16 , Epigenesis, Genetic , Exons , Gene Expression Regulation, Leukemic , Histone-Lysine N-Methyltransferase/genetics , Histones/metabolism , Humans , Introns , Myeloid-Lymphoid Leukemia Protein/genetics , Nuclear Proteins/genetics , Promoter Regions, GeneticABSTRACT
Colorectal cancer (CRC) is one of the most common malignancies in the Western world showing an increasing incidence, and has been associated with genetic and lifestyle factors. Individual susceptibility to CRC may be due partly to variations in detoxification capacity in the gastrointestinal tract. Genetic polymorphisms in detoxification enzymes may result in variations in detoxification activities, which subsequently might influence the levels of toxic/carcinogenic compounds, and this may influence the risk for CRC. To determine whether genetic polymorphisms in detoxification enzymes predispose to the development of CRC, 371 patients with sporadic CRC and 415 healthy controls were genotyped for polymorphisms in the important detoxification enzymes UDP-glucuronosyltransferase UGT1A1, UGT1A6, UGT1A7 and UGT1A8, and glutathione S-transferase GSTA1, GSTM1, GSTP1 and GSTT1. Patients and controls were all of Caucasian origin. DNA was isolated from either blood or tissue and tested by polymerase chain reaction followed by restriction fragment length polymorphism analyses. Logistic regression analyses showed significant age- and gender-adjusted risks for CRC associated with variant genotypes of UGT1A6 [OR 1.5, 95% (confidence interval) CI 1.03-2.3] and UGT1A7 (OR 2.4, 95% CI 1.3-4.6), whereas no associations were found between CRC and the other polymorphic genes as mentioned above. In conclusion, the data suggest that the presence of variant UGT1A6 and UGT1A7 genotypes with expected reduced enzyme activities, might enhance susceptibility to CRC.
