Search

Glucocorticoids (GCs) are a mainstay of contemporary, multidrug chemotherapy in the treatment of childhood acute lymphoblastic leukemia (ALL), and resistance to GCs remains a major clinical concern. Resistance to GCs is predictive of ALL relapse and poor clinical outcome, and therefore represents a major hurdle limiting further improvements in survival rates. While advances have been made in identifying genes implicated in GC resistance, there remains an insufficient understanding of the impact of cis-regulatory disruptions in resistance. To address this, we mapped the gene regulatory response to GCs in two ALL cell lines using functional genomics and high-throughput reporter assays and identified thousands of GC-responsive changes to chromatin state, including the formation of over 250 GC-responsive super-enhancers and a depletion of AP-1 bound cis-regulatory elements implicated in cell proliferation and anti-apoptotic processes. By integrating our GC response maps with genetic and epigenetic datasets in primary ALL cells from patients, we further uncovered cis-regulatory disruptions at GC-responsive genes that impact GC resistance in childhood ALL. Overall, these data indicate that GCs initiate pervasive effects on the leukemia epigenome, and that alterations to the GC gene regulatory network contribute to GC resistance.

Understanding the genomic and epigenetic mechanisms of drug resistance in pediatric acute lymphoblastic leukemia (ALL) is critical for further improvements in treatment outcomes. The role of transcriptomic response in conferring resistance to l-asparaginase (LASP) is poorly understood beyond asparagine synthetase (ASNS). We defined reproducible LASP response genes in LASP-resistant and LASP-sensitive ALL cell lines as well as primary leukemia samples from newly diagnosed patients. Defining target genes of the amino acid stress response-related transcription factor activating transcription factor 4 (ATF4) in ALL cell lines using chromatin immunoprecipitation sequencing (ChIP-seq) revealed 45% of genes that changed expression after LASP treatment were direct targets of the ATF4 transcription factor, and 34% of these genes harbored LASP-responsive ATF4 promoter binding events. SLC7A11 was found to be a response gene in cell lines and patient samples as well as a direct target of ATF4. SLC7A11 was also one of only 2.4% of LASP response genes with basal level gene expression that also correlated with LASP ex vivo resistance in primary leukemia cells. Experiments using chemical inhibition of SLC7A11 with sulfasalazine, gene overexpression, and partial gene knockout recapitulated LASP resistance or sensitivity in ALL cell lines. These findings show the importance of assessing changes in gene expression following treatment with an antileukemic agent for its association with drug resistance and highlight that many response genes may not differ in their basal expression in drug-resistant leukemia cells.

Background There is growing evidence for the inherited basis of susceptibility to childhood acute lymphoblastic leukemia (ALL). Genome-wide association studies have identified non-coding ALL risk variants at the ARID5B gene locus, but their exact functional effects and the molecular mechanism linking ARID5B to B-cell ALL leukemogenesis remain largely unknown. Methods We performed targeted sequencing of ARID5B in germline DNA of 5008 children with ALL. Variants were evaluated for association with ALL susceptibility using 3644 patients from the UK10K cohort as non-ALL controls, under an additive model. Cis-regulatory elements in ARID5B were systematically identified using dCas9-KRAB–mediated enhancer interference system enhancer screen in ALL cells. Disruption of transcription factor binding by ARID5B variant was predicted informatically and then confirmed using chromatin immunoprecipitation and coimmunoprecipitation. ARID5B variant association with hematological traits was examined using UK Biobank dataset. All statistical tests were 2-sided. Results We identified 54 common variants in ARID5B statistically significantly associated with leukemia risk, all of which were noncoding. Six cis-regulatory elements at the ARID5B locus were discovered using CRISPR-based high-throughput enhancer screening. Strikingly, the top ALL risk variant (rs7090445, P = 5.57 × 10–45) is located precisely within the strongest enhancer element, which is also distally tethered to the ARID5B promoter. The variant allele disrupts the MEF2C binding motif sequence, resulting in reduced MEF2C affinity and decreased local chromosome accessibility. MEF2C influences ARID5B expression in ALL, likely via a transcription factor complex with RUNX1. Using the UK Biobank dataset (n = 349 861), we showed that rs7090445 was also associated with lymphocyte percentage and count in the general population (P = 8.6 × 10–22 and 2.1 × 10–18, respectively). Conclusions Our results indicate that ALL risk variants in ARID5B function by modulating cis-regulatory elements at this locus.

Inhibition of members of the bromodomain and extraterminal (BET) family of proteins has proven a valid strategy for cancer chemotherapy. All BET identified to date contain two bromodomains (BD; BD1 and BD2) that are necessary for recognition of acetylated lysine residues in the N-terminal regions of histones. Chemical matter that targets BET (BETi) also interact via these domains. Molecular and cellular data indicate that BD1 and BD2 have different biological roles depending upon their cellular context, with BD2 particularly associated with cancer. We have therefore pursued the development of BD2-selective molecules both as chemical probes and as potential leads for drug development. Here we report the structure-based generation of a novel series of tetrahydroquinoline analogs that exhibit >50-fold selectivity for BD2 versus BD1. This selective targeting resulted in engagement with BD-containing proteins in cells, resulting in modulation of MYC proteins and downstream targets. These compounds were potent cytotoxins toward numerous pediatric cancer cell lines and were minimally toxic to nontumorigenic cells. In addition, unlike the pan BETi (+)-JQ1, these BD2-selective inhibitors demonstrated no rebound expression effects. Finally, we report a pharmacokinetic-optimized, metabolically stable derivative that induced growth delay in a neuroblastoma xenograft model with minimal toxicity. We conclude that BD2-selective agents are valid candidates for antitumor drug design for pediatric malignancies driven by the MYC oncogene.Significance:This study presents bromodomain-selective BET inhibitors that act as antitumor agents and demonstrates that these molecules have in vivo activity towards neuroblastoma, with essentially no toxicity.

SUMMARYB-cell lineage acute lymphoblastic leukemia (B-ALL) is comprised of diverse molecular subtypes and while transcriptional and DNA methylation profiling of B-ALL subtypes has been extensively examined, the accompanying chromatin landscape is not well characterized for many subtypes. We therefore mapped chromatin accessibility using ATAC-seq for 10 B-ALL molecular subtypes in primary ALL cells from 154 patients. Comparisons with B-cell progenitors identified candidate B-ALL cell-of-origin and AP-1-associatedcis-regulatory rewiring in B-ALL.Cis-regulatory rewiring promoted B-ALL-specific gene regulatory networks impacting oncogenic signaling pathways that perturb normal B-cell development. We also identified that over 20% of B-ALL accessible chromatin sites exhibit strong subtype enrichment, with transcription factor (TF) footprint profiling identifying candidate TFs that maintain subtype-specific chromatin architectures. Over 9000 inherited genetic variants were further uncovered that contribute to variability in chromatin accessibility among individual patient samples. Overall, our data suggest that distinct chromatin architectures are driven by diverse TFs and inherited genetic variants which promote unique gene regulatory networks that contribute to transcriptional differences among B-ALL subtypes.HIGHLIGHTSPro-B progenitor cells as the most common cell-of-origin for B-ALLAP-1 TF-associatedcis-regulatory rewiring in B-ALLSubtype-specific accessible chromatin signatures representing 20% of all B-ALL sitesRole for distinct TFs in promoting subtype-specific chromatin architecturesThousands of inherited genetic variants identified impacting chromatin state

Although acute lymphoblastic leukemia (ALL) is the most common childhood cancer, there is limited understanding of the contribution of inherited genetic variation on inter-individual differences in chemotherapy response. Defining genetic factors impacting therapy failure can help better predict response and identify drug resistance mechanisms. We therefore mapped inherited noncoding variants associated with chemotherapeutic drug resistance and/or treatment outcome to ALLcis-regulatory elements and investigated their gene regulatory potential and genomic connectivity using massively parallel reporter assays and promoter capture Hi-C, respectively. We identified 53 variants with reproducible allele-specific effects on transcription and high-confidence gene targets. Subsequent functional interrogation of the top variant (rs1247117) determined that it disrupted a PU.1 consensus motif and PU.1 binding affinity. Importantly, deletion of the genomic interval containing rs1247117 sensitized ALL cells to vincristine. Together, these data demonstrate that noncoding regulatory variation associated with diverse pharmacological traits harbor significant effects on allele-specific transcriptional activity and impact sensitivity to chemotherapeutic agents in ALL.

Glucocorticoids (GCs; i.e., steroids) are important chemotherapeutic agents in the treatment of B-cell precursor acute lymphoblastic leukemia (B-ALL) andde novoGC resistance predicts relapse and poor clinical outcome in patients. Glucocorticoids induce B-ALL cell apoptosis through activation of glucocorticoid receptor (GR), a ligand-induced nuclear receptor transcription factor (TF). We previously identified disruptions to glucocorticoid receptor (GR)-boundcis-regulatory elements controllingTLE1expression in GC-resistant primary B-ALL cells from patients.TLE1is a GC-response gene up-regulated by steroids and functions as a canonical Wnt signaling repressor. To better understand the mechanistic relationship between GC signaling and canonical Wnt signaling, we performed diverse functional analyses that identified extensive crosstalk and mutual antagonism between these two signaling pathways in B-ALL. We determined that crosstalk and antagonism was driven by the binding of GR and the canonical Wnt signaling TFs LEF1 and TCF7L2 to overlapping sets ofcis-regulatory elements associated with genes impacting cell death and cell proliferation, and was further accompanied by overlapping and opposing transcriptional programs. Our data additionally suggest thatcis-regulatory disruptions atTLE1are linked to GC resistance through a dampening of the GC response and GC-mediated apoptosis via enhanced canonical Wnt signaling. As a result of the extensive genomic and gene regulatory connectivity between these two signaling pathways, our data supports the importance of canonical Wnt signaling in mediating GC resistance in B-ALL.

Glucocorticoids (GCs) are a mainstay of contemporary, multi-drug chemotherapy in the treatment of childhood acute lymphoblastic leukemia (ALL), and resistance to GCs remains a major clinical concern. Resistance to GCs is predictive of ALL relapse and poor clinical outcome, and therefore represents a major hurdle limiting further improvements in survival rates. While advances have been made in identifying genes implicated in GC resistance, there remains an insufficient understanding of the impact of cis-regulatory disruptions in resistance. To address this, we mapped the gene regulatory response to GCs in two ALL cells using functional genomics and high-throughput reporter assays and identified thousands of GC-responsive changes to chromatin state, including the formation of over 250 GC-responsive super-enhancers and a depletion of AP-1 bound cis-regulatory elements implicated in cell proliferation and anti-apoptotic processes. By integrating our GC response maps with genetic and epigenetic datasets in primary ALL cells from patients, we further uncovered cis-regulatory disruptions at GC-responsive genes that impact GC resistance in childhood ALL. Overall, these data indicate that GCs initiate pervasive effects on the leukemia epigenome, and that alterations to the GC gene regulatory network contribute to GC resistance.

Glucocorticoids (GCs) are a mainstay of contemporary, multidrug chemotherapy in the treatment of childhood acute lymphoblastic leukemia (ALL), and resistance to GCs remains a major clinical concern. Resistance to GCs is predictive of ALL relapse and poor clinical outcome, and therefore represents a major hurdle limiting further improvements in survival rates. While advances have been made in identifying genes implicated in GC resistance, there remains an insufficient understanding of the impact of cis-regulatory disruptions in resistance. To address this, we mapped the gene regulatory response to GCs in two ALL cell lines using functional genomics and high-throughput reporter assays and identified thousands of GC-responsive changes to chromatin state, including the formation of over 250 GC-responsive super-enhancers and a depletion of AP-1 bound cis-regulatory elements implicated in cell proliferation and anti-apoptotic processes. By integrating our GC response maps with genetic and epigenetic datasets in primary ALL cells from patients, we further uncovered cis-regulatory disruptions at GC-responsive genes that impact GC resistance in childhood ALL. Overall, these data indicate that GCs initiate pervasive effects on the leukemia epigenome, and that alterations to the GC gene regulatory network contribute to GC resistance.

Inhibition of members of the bromodomain and extraterminal (BET) family of proteins has proven a valid strategy for cancer chemotherapy. All BET identified to date contain two bromodomains (BD; BD1 and BD2) that are necessary for recognition of acetylated lysine residues in the N-terminal regions of histones. Chemical matter that targets BET (BETi) also interact via these domains. Molecular and cellular data indicate that BD1 and BD2 have different biological roles depending upon their cellular context, with BD2 particularly associated with cancer. We have therefore pursued the development of BD2-selective molecules both as chemical probes and as potential leads for drug development. Here we report the structure-based generation of a novel series of tetrahydroquinoline analogs that exhibit >50-fold selectivity for BD2 versus BD1. This selective targeting resulted in engagement with BD-containing proteins in cells, resulting in modulation of MYC proteins and downstream targets. These compounds were potent cytotoxins toward numerous pediatric cancer cell lines and were minimally toxic to nontumorigenic cells. In addition, unlike the pan BETi (+)-JQ1, these BD2-selective inhibitors demonstrated no rebound expression effects. Finally, we report a pharmacokinetic-optimized, metabolically stable derivative that induced growth delay in a neuroblastoma xenograft model with minimal toxicity. We conclude that BD2-selective agents are valid candidates for antitumor drug design for pediatric malignancies driven by the MYC oncogene. Significance: This study presents bromodomain-selective BET inhibitors that act as antitumor agents and demonstrates that these molecules have in vivo activity towards neuroblastoma, with essentially no toxicity.

Transcription factors are DNA-binding proteins that have key roles in gene regulation1,2. Genome-wide occupancy maps of transcriptional regulators are important for understanding gene regulation and its effects on diverse biological processes3–6. However, only a minority of the more than 1,600 transcription factors encoded in the human genome has been assayed. Here we present, as part of the ENCODE (Encyclopedia of DNA Elements) project, data and analyses from chromatin immunoprecipitation followed by high-throughput sequencing (ChIP–seq) experiments using the human HepG2 cell line for 208 chromatin-associated proteins (CAPs). These comprise 171 transcription factors and 37 transcriptional cofactors and chromatin regulator proteins, and represent nearly one-quarter of CAPs expressed in HepG2 cells. The binding profiles of these CAPs form major groups associated predominantly with promoters or enhancers, or with both. We confirm and expand the current catalogue of DNA sequence motifs for transcription factors, and describe motifs that correspond to other transcription factors that are co-enriched with the primary ChIP target. For example, FOX family motifs are enriched in ChIP–seq peaks of 37 other CAPs. We show that motif content and occupancy patterns can distinguish between promoters and enhancers. This catalogue reveals high-occupancy target regions at which many CAPs associate, although each contains motifs for only a minority of the numerous associated transcription factors. These analyses provide a more complete overview of the gene regulatory networks that define this cell type, and demonstrate the usefulness of the large-scale production efforts of the ENCODE Consortium., ChIP–seq and CETCh–seq data are used to analyse binding maps for 208 transcription factors and other chromatin-associated proteins in a single human cell type, providing a comprehensive catalogue of the transcription factor landscape and gene regulatory networks in these cells.

FLT3 is an attractive therapeutic target in acute lymphoblastic leukemia (ALL) but the mechanism for its activation in this cancer is incompletely understood. Profiling global gene expression in large ALL cohorts, we identify over-expression of FLT3 in ZNF384-rearranged ALL, consistently across cases harboring different fusion partners with ZNF384. Mechanistically, we discover an intergenic enhancer element at the FLT3 locus that is exclusively activated in ZNF384-rearranged ALL, with the enhancer-promoter looping directly mediated by the fusion protein. There is also a global enrichment of active enhancers within ZNF384 binding sites across the genome in ZNF384-rearranged ALL cells. Downregulation of ZNF384 blunts FLT3 activation and decreases ALL cell sensitivity to FLT3 inhibitor gilteritinib in vitro. In patient-derived xenograft models of ZNF384-rearranged ALL, gilteritinib exhibits significant anti-leukemia efficacy as a monotherapy in vivo. Collectively, our results provide insights into FLT3 regulation in ALL and point to potential genomics-guided targeted therapy for this patient population.

Acute lymphoblastic leukemia (ALL) is a hematopoietic malignancy comprised of molecular subtypes largely characterized by aneuploidy or recurring chromosomal rearrangements. Despite extensive information on the ALL transcriptome and methylome, there is limited understanding of the ALL chromatin landscape. We therefore mapped accessible chromatin in 24 primary ALL cell biospecimens comprising three common molecular subtypes (DUX4/ERG, ETV6-RUNX1 and hyperdiploid) from patients treated at St. Jude Children's Research Hospital. Our findings highlight extensive chromatin reprogramming in ALL, including the identification ALL subtype-specific chromatin landscapes that are additionally modulated by genetic variation. Chromatin accessibility differences between ALL and normal B-cells implicate the activation of B-cell repressed chromatin domains and detail the disruption of normal B-cell development in ALL. Among ALL subtypes, we uncovered roles for basic helix-loop-helix, homeodomain and activator protein 1 transcription factors in promoting subtype-specific chromatin accessibility and distinct gene regulatory networks. In addition to chromatin subtype-specificity, we further identified over 3500 DNA sequence variants that alter the ALL chromatin landscape and contribute to inter-individual variability in chromatin accessibility. Collectively, our data suggest that subtype-specific chromatin landscapes and gene regulatory networks impact ALL biology and contribute to transcriptomic differences among ALL subtypes.

Acute lymphoblastic leukemia (ALL) is the most common cancer in children. TCF3-PBX1 fusion defines a common molecular subtype of ALL with unique clinical features, but the molecular basis of its inherited susceptibility is unknown. In a genome-wide association study of 1494 ALL cases and 2057 non-ALL controls, we identified a germline risk locus located in an intergenic region between BCL11A and PAPOLG: rs2665658, P = 1.88 × 10–8 for TCF3-PBX1 ALL vs non-ALL, and P = 1.70 × 10–8 for TCF3-PBX1 ALL vs other-ALL. The lead variant was validated in a replication cohort, and conditional analyses pointed to a single causal variant with subtype-specific effect. The risk variant is located in a regulatory DNA element uniquely activated in ALL cells with the TCF3-PBX1 fusion and may distally modulate the transcription of the adjacent gene REL. Our results expand the understanding of subtype-specific ALL susceptibility and highlight plausible interplay between germline variants and somatic genomic abnormalities in ALL pathogenesis.

Acute lymphoblastic leukemia (ALL) is the most prevalent childhood cancer and despite improved survival rates, relapsed ALL is still among the most common causes of cancer death in children. Although changes in the expression of specific genes have been linked to chemotherapeutic resistance, relatively little is understood of the pharmacogenomic impact of the noncoding, cis-regulatory landscape governing gene regulation. Glucocorticoids (GCs; i.e. steroids) are a mainstay of contemporary, multi-drug chemotherapy in ALL, and GC resistance is predictive of both relapse and poor clinical outcome in ALL. Because GCs function through activation of glucocorticoid receptor (GR), a nuclear receptor transcription factor that interacts directly with cis-regulatory elements, unveiling the glucocorticoid gene regulatory network (GC-GRN) in leukemia cells is crucial to understanding not only the biological mechanism of apoptosis, but also illuminating gene regulatory mechanisms contributing to GC resistance. To test the hypothesis that alterations to the GC-GRN are important contributors to steroid resistance in ALL, we comprehensively mapped cellular responses to GCs in human ALL cell lines using >100 independent functional genomic datasets. This comprehensive approach uncovered thousands of genes and cis-regulatory elements that were responsive to GCs, and further identified >38,000 high-confidence glucocorticoid response elements (GREs) in the ALL genome. A closer examination of these data revealed GR binding profiles that were consistent with the long-range flexible billboard model of gene regulation. By further integrating our results with genetic and epigenetic data in primary ALL cells from patients enrolled on St. Jude clinical trials, we identified 45 DNA sequence variants associated with ex vivo GC resistance that map to GREs and functionally validated an associated variant within the TLE1 gene locus. We also uncovered 1929 accessible chromatin sites (FDR10,000 sgRNAs targeting >1000 GR binding events at putative GC-resistance accessible chromatin sites identified a subset of GR binding sites implicated in GC resistance. Overall, these data indicate that GCs initiate pervasive, genome-wide effects on the leukemia epigenome and transcriptome, and that genetic and epigenetic alterations to GREs are mechanisms contributing to GC resistance in childhood ALL. Disclosures Pui: Adaptive Biotechnologies: Membership on an entity's Board of Directors or advisory committees; Novartis: Other: Data Monitoring Committee. Evans: Princess Máxima Center for Pediatric Oncology, Scientific Advisory Board, Chair: Membership on an entity's Board of Directors or advisory committees; BioSkryb, Inc.: Membership on an entity's Board of Directors or advisory committees; St. Jude Children's Research Hospital, Emeritus Member (began Jan 2021): Ended employment in the past 24 months.

Understanding the genomic and epigenetic mechanisms of drug resistance in pediatric acute lymphoblastic leukemia (ALL) is critical for further improvements in treatment outcome. The role of transcriptomic response in conferring resistance to l-asparaginase (LASP) is poorly understood, beyond asparagine synthetase (ASNS). We defined reproducible LASP response genes in LASP resistant and sensitive ALL cell lines (n = 7) as well as primary leukemia samples from newly diagnosed patients. We identified 2219 response genes (absolute log 2FC > 1.5, FDR p-value Disclosures Stock: Pfizer: Consultancy, Honoraria, Research Funding; amgen: Honoraria; agios: Honoraria; jazz: Honoraria; kura: Honoraria; kite: Honoraria; morphosys: Honoraria; servier: Honoraria; syndax: Consultancy, Honoraria; Pluristeem: Consultancy, Honoraria. Mullighan: Amgen: Current equity holder in publicly-traded company; Illumina: Membership on an entity's Board of Directors or advisory committees; AbbVie: Research Funding; Pfizer: Research Funding. Pui: Adaptive Biotechnologies: Membership on an entity's Board of Directors or advisory committees; Novartis: Other: Data Monitoring Committee. Evans: Princess Máxima Center for Pediatric Oncology, Scientific Advisory Board, Chair: Membership on an entity's Board of Directors or advisory committees; BioSkryb, Inc.: Membership on an entity's Board of Directors or advisory committees; St. Jude Children's Research Hospital, Emeritus Member (began Jan 2021): Ended employment in the past 24 months.

Acute lymphoblastic leukemia (ALL) is the most common malignancy in pediatric patients spanning both B- and T-cell lineages. ALL also occurs less commonly in adults but with cure rates of only around 30%. Past work has characterized ALL into molecular subtypes spanning a range of aberrant chromosomal rearrangements and oncogene chimeric fusions driving malignancy. While transcriptional profiling of these subtypes has been extensively examined, the accompanying chromatin accessibility landscape and corresponding gene regulatory repertoire is not well characterized for many subtypes. To better profile the ALL epigenomic and gene regulatory repertoire we examined chromatin accessibility of 12 distinct molecular subtypes (BCR-ABL1, ETV6-RUNX1, Hyperdiploid, Hypodiploid, KMT2A rearranged, Ph-Like, PAX5, DUX4/ERG, TCF3-PBX1, T-ALL, Early T Precursor and B-other) across 189 primary patient samples of pediatric ALL (n = 106) and adult ALL (n = 83) origin using ATAC-seq. To our knowledge, this represents the largest collection of chromatin accessibility data in primary ALL samples spanning multiple molecular subtypes to date. Collectively, we identified over 600,000 accessible chromatin sites in the ALL genome with over 50,000 regions of differentially accessible chromatin encompassing both common and subtype-specific modalities. Further, transcription factor (TF) footprint profiling of ATAC-seq yielded tens of thousands of candidate TF binding events and identified key TF drivers within distinct molecular subtypes. We additionally performed H3K27ac ChIP-seq in a subset of 12 primary ALL patient samples, with integration of these histone data for select patient samples allowing inferences about candidate super-enhancer drivers of ALL molecular subtypes. Overall, these analyses and data offer a window into the gene regulatory and epigenetic landscape of ALL, and further highlight the complexity and heterogeneity of accessible chromatin landscapes among distinct molecular subtypes of ALL. Disclosures Stock: Pfizer: Consultancy, Honoraria, Research Funding; amgen: Honoraria; agios: Honoraria; jazz: Honoraria; kura: Honoraria; kite: Honoraria; morphosys: Honoraria; servier: Honoraria; syndax: Consultancy, Honoraria; Pluristeem: Consultancy, Honoraria. Pui: Adaptive Biotechnologies: Membership on an entity's Board of Directors or advisory committees; Novartis: Other: Data Monitoring Committee. Evans: Princess Máxima Center for Pediatric Oncology, Scientific Advisory Board, Chair: Membership on an entity's Board of Directors or advisory committees; St. Jude Children's Research Hospital, Emeritus Member (began Jan 2021): Ended employment in the past 24 months; BioSkryb, Inc.: Membership on an entity's Board of Directors or advisory committees.

// Joy M. McDaniel 1, 2 , Katherine E. Varley 3 , Jason Gertz 3 , Daniel S. Savic 1 , Brian S. Roberts 1 , Sarah K. Bailey 4 , Lalita A. Shevde 4, 6 , Ryne C. Ramaker 1, 7 , Brittany N. Lasseigne 1 , Marie K. Kirby 1 , Kimberly M. Newberry 1 , E. Christopher Partridge 1 , Angela L. Jones 1 , Braden Boone 1 , Shawn E. Levy 1 , Patsy G. Oliver 5 , Katherine C. Sexton 6 , William E. Grizzle 6 , Andres Forero 6 , Donald J. Buchsbaum 5 , Sara J. Cooper 1 , Richard M. Myers 1 1 HudsonAlpha Institute for Biotechnology, Huntsville, AL 35806, USA 2 The University of Alabama in Huntsville, Huntsville, AL 35899, USA 3 Department of Oncological Sciences, Huntsman Cancer Institute, University of Utah, Salt Lake City, UT, 84112, USA 4 Department of Pathology, University of Alabama at Birmingham, Birmingham, AL, 35294, USA 5 Department of Radiation Oncology, University of Alabama at Birmingham, Birmingham, AL, 35294, USA 6 University of Alabama at Birmingham Comprehensive Cancer Center, Birmingham, AL 35294, USA 7 Department of Genetics, University of Alabama at Birmingham, Birmingham, AL, 35294, USA Correspondence to: Richard M. Myers, email: rmyers@hudsonalpha.org Keywords: TNBC, STAT3, ChIP-seq, RNA-seq, invasion Received: August 15, 2016 Accepted: November 14, 2016 Published: December 24, 2016 ABSTRACT Breast cancer is a heterogeneous disease comprised of four molecular subtypes defined by whether the tumor-originating cells are luminal or basal epithelial cells. Breast cancers arising from the luminal mammary duct often express estrogen receptor (ER), progesterone receptor (PR), and human epidermal growth receptor 2 (HER2). Tumors expressing ER and/or PR are treated with anti-hormonal therapies, while tumors overexpressing HER2 are targeted with monoclonal antibodies. Immunohistochemical detection of ER, PR, and HER2 receptors/proteins is a critical step in breast cancer diagnosis and guided treatment. Breast tumors that do not express these proteins are known as “triple negative breast cancer” (TNBC) and are typically basal-like. TNBCs are the most aggressive subtype, with the highest mortality rates and no targeted therapy, so there is a pressing need to identify important TNBC tumor regulators. The signal transducer and activator of transcription 3 (STAT3) transcription factor has been previously implicated as a constitutively active oncogene in TNBC. However, its direct regulatory gene targets and tumorigenic properties have not been well characterized. By integrating RNA-seq and ChIP-seq data from 2 TNBC tumors and 5 cell lines, we discovered novel gene signatures directly regulated by STAT3 that were enriched for processes involving inflammation, immunity, and invasion in TNBC. Functional analysis revealed that STAT3 has a key role regulating invasion and metastasis, a characteristic often associated with TNBC. Our findings suggest therapies targeting STAT3 may be important for preventing TNBC metastasis.

Deciphering the intricate molecular processes that orchestrate the spatial and temporal regulation of genes has become an increasingly major focus of biological research. The differential expression of genes by diverse cell types with a common genome is a hallmark of complex cellular functions, as well as the basis for multicellular life. Importantly, a more coherent understanding of gene regulation is critical for defining developmental processes, evolutionary principles and disease etiologies. Here we present our current understanding of gene regulation by focusing on the role of enhancer elements in these complex processes. Although functional genomic methods have provided considerable advances to our understanding of gene regulation, these assays, which are usually performed on a genome-wide scale, typically provide correlative observations that lack functional interpretation. Recent innovations in genome editing technologies have placed gene regulatory studies at an exciting crossroads, as systematic, functional evaluation of enhancers and other transcriptional regulatory elements can now be performed in a coordinated, high-throughput manner across the entire genome. This review provides insights on transcriptional enhancer function, their role in development and disease, and catalogues experimental tools commonly used to study these elements. Additionally, we discuss the crucial role of novel techniques in deciphering the complex gene regulatory landscape and how these studies will shape future research.

Summary Despite considerable research connecting cellular metabolism with differentiation decisions, the underlying mechanisms that translate metabolite-sensitive activities into unique gene programs are still unclear. We found that aspects of the interleukin-2 (IL-2)-sensitive effector gene program in CD4 + and CD8 + T cells in type 1 conditions (Th1) were regulated by glutamine and alpha-ketoglutarate (αKG)-induced events, in part through changes in DNA and histone methylation states. We further identified a mechanism by which IL-2- and αKG-sensitive metabolic changes regulated the association of CCCTC-binding factor (CTCF) with select genomic sites. αKG-sensitive CTCF sites were often associated with loci containing IL-2- and αKG-sensitive genome organization patterns and gene expression in T cells. IL-2- and αKG-sensitive CTCF sites in T cells were also associated with genes from developmental pathways that had αKG-sensitive expression in embryonic stem cells. The data collectively support a mechanism wherein CTCF serves to translate αKG-sensitive metabolic changes into context-dependent differentiation gene programs.

Large-scale efforts like the Encyclopedia of DNA Elements (ENCODE) Project have made tremendous progress in cataloging the genomic binding patterns of DNA-associated proteins (DAPs), such as transcription factors (TFs). However most chromatin immunoprecipitation-sequencing (ChIP-seq) analyses have focused on a few immortalized cell lines whose activities and physiology deviate in important ways from endogenous cells and tissues. Consequently, binding data from primary human tissue are essential to improving our understanding of in vivo gene regulation. Here we analyze ChIP-seq data for 20 DAPs assayed in two healthy human liver tissue samples, identifying more than 450,000 binding sites. We integrated binding data with transcriptome and phased whole genome data to investigate allelic DAP interactions and the impact of heterozygous sequence variation on the expression of neighboring genes. We find our tissue-based dataset demonstrates binding patterns more consistent with liver biology than cell lines, and describe uses of these data to better prioritize impactful non-coding variation. Collectively, our rich dataset offers novel insights into genome function in healthy liver tissue and provides a valuable research resource for assessing disease-related disruptions.

Non-coding variation within TCF7L2 remains the strongest genetic determinant of type 2 diabetes risk in humans. A considerable effort has been placed in understanding the functional roles of TCF7L2 in pancreatic beta cells, despite evidence of TCF7L2 expression in various peripheral tissues important in glucose homeostasis. Here, we use a humanized mouse model overexpressing Tcf7l2, resulting in glucose intolerance, to infer the contribution of Tcf7l2 overexpression in beta cells and in other tissues to the metabolic phenotypes displayed by these mice. Restoring Tcf7l2 expression specifically in beta cells to endogenous levels, in face of its overexpression elsewhere, results in impaired insulin secretion, reduced beta cell number and islet area, corroborating data obtained in humans showing similar phenotypes as a result of manipulations leading to Tcf7l2 loss of function. Interestingly, the persistent overexpression of Tcf7l2 in non-pancreatic tissues results in a significant worsening in glucose tolerance in vivo, indicating that Tcf7l2 overexpression in beta cells does not account for the glucose intolerance in the Tcf7l2 overexpression mouse model. Collectively, these data posit that Tcf7l2 plays key roles in glucose metabolism through actions beyond pancreatic beta cells, and further points to functionally opposing cell-type specific effects for Tcf7l2 on the maintenance of balanced glucose metabolism, thereby urging a careful examination of its role in non-pancreatic tissues as well as its composite metabolic effects across distinct tissues. Uncovering these roles may lead to new therapeutic targets for type 2 diabetes.

SUMMARY The basic leucine zipper transcription factor activating transcription factor-like (Batf) contributes to transcriptional programming of multiple effector T cells and is required for T helper 17 (Th17) and T follicular helper (Tfh) cell development. Here, we examine mechanisms by which Batf initiates gene transcription in developing effector CD4 T cells. We find that, in addition to its pioneering function, Batf controls developmentally regulated recruitment of the architectural factor Ctcf to promote chromatin looping that is associated with lineage-specific gene transcription. The chromatin-organizing actions of Batf are largely dependent on Ets1, which appears to be indispensable for the Batf-dependent recruitment of Ctcf. Moreover, most of the Batf-dependent sites to which Ctcf is recruited lie outside of activating protein-1-interferon regulatory factor (Ap-1-Irf) composite elements (AICEs), indicating that direct involvement of Batf-Irf complexes is not required. These results identify a cooperative role for Batf, Ets1, and Ctcf in chromatin reorganization that underpins the transcriptional programming of effector T cells., In Brief Pham et al. uncover mechanisms by which Batf restructures the chromatin landscape during CD4+ effector T cell differentiation. Batf controls Ctcf recruitment to lineage-specifying gene loci in an Ets1-dependent manner to promote chromatin looping and lineage-specific gene transcription, thereby identifying a heretofore unknown cooperativity of these factors in effector T cell development., Graphical abstract

Background The liver X receptors (LXRs, NR1H2 and NR1H3) and peroxisome proliferator-activated receptor gamma (PPARG, NR1C3) nuclear receptor transcription factors (TFs) are master regulators of energy homeostasis. Intriguingly, recent studies suggest that these metabolic regulators also impact tumor cell proliferation. However, a comprehensive temporal molecular characterization of the LXR and PPARG gene regulatory responses in tumor cells is still lacking. Methods To better define the underlying molecular processes governing the genetic control of cellular growth in response to extracellular metabolic signals, we performed a comprehensive, genome-wide characterization of the temporal regulatory cascades mediated by LXR and PPARG signaling in HT29 colorectal cancer cells. For this analysis, we applied a multi-tiered approach that incorporated cellular phenotypic assays, gene expression profiles, chromatin state dynamics, and nuclear receptor binding patterns. Results Our results illustrate that the activation of both nuclear receptors inhibited cell proliferation and further decreased glutathione levels, consistent with increased cellular oxidative stress. Despite a common metabolic reprogramming, the gene regulatory network programs initiated by these nuclear receptors were widely distinct. PPARG generated a rapid and short-term response while maintaining a gene activator role. By contrast, LXR signaling was prolonged, with initial, predominantly activating functions that transitioned to repressive gene regulatory activities at late time points. Conclusions Through the use of a multi-tiered strategy that integrated various genomic datasets, our data illustrate that distinct gene regulatory programs elicit common phenotypic effects, highlighting the complexity of the genome. These results further provide a detailed molecular map of metabolic reprogramming in cancer cells through LXR and PPARG activation. As ligand-inducible TFs, these nuclear receptors can potentially serve as attractive therapeutic targets for the treatment of various cancers. Electronic supplementary material The online version of this article (doi:10.1186/s13073-016-0328-6) contains supplementary material, which is available to authorized users.

Pioneer transcription factors are a unique class of proteins that access heterochromatin, facilitate chromatin accessibility, and recruit other transcription factors to regulate gene transcription. The basic leucine zipper transcription factor ATF-like (BATF) is a critical factor for the development of T helper cell lineages. While BATF has been proposed as a pioneer factor in both T and B cells, its molecular functions are incompletely understood. In this study, we have identified a novel aspect of the pioneer function of BATF in controlling the chromatin landscape in developing effector T cells. We find that BATF is indispensable for the recruitment of chromatin-organizing complexes containing the DNA binding factor CTCF and other co-factors such as Oct1 and Ets1 to a subset of enhancer and chromatin-looping nodes during effector T cell lineage specification. Accordingly, BATF has a unique function in modulating genomic spatial organization in a CTCF-dependent manner. This study establishes BATF as a bona fide pioneer factor and provides novel insight into its role as a spatial genome organizer that modulates chromatin accessibility and initiates chromatin looping for the coordinate and differential gene expression in CD4+ T cells.

Chromatin immunoprecipitation followed by next-generation DNA sequencing (ChIP-seq) is a widely used technique for identifying transcription factor (TF) binding events throughout an entire genome. However, ChIP-seq is limited by the availability of suitable ChIP-seq grade antibodies, and the vast majority of commercially available antibodies fail to generate usable data sets. To ameliorate these technical obstacles, we present a robust methodological approach for performing ChIP-seq through epitope tagging of endogenous TFs. We used clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9-based genome editing technology to develop CRISPR epitope tagging ChIP-seq (CETCh-seq) of DNA-binding proteins. We assessed the feasibility of CETCh-seq by tagging several DNA-binding proteins spanning a wide range of endogenous expression levels in the hepatocellular carcinoma cell line HepG2. Our data exhibit strong correlations between both replicate types as well as with standard ChIP-seq approaches that use TF antibodies. Notably, we also observed minimal changes to the cellular transcriptome and to the expression of the tagged TF. To examine the robustness of our technique, we further performed CETCh-seq in the breast adenocarcinoma cell line MCF7 as well as mouse embryonic stem cells and observed similarly high correlations. Collectively, these data highlight the applicability of CETCh-seq to accurately define the genome-wide binding profiles of DNA-binding proteins, allowing for a straightforward methodology to potentially assay the complete repertoire of TFs, including the large fraction for which ChIP-quality antibodies are not available.

Transcription factors (TFs) bind to thousands of DNA sequences in mammalian genomes, but most of these binding events appear to have no direct effect on gene expression. It is unclear why only a subset of TF bound sites are actively involved in transcriptional regulation. Moreover, the key genomic features that accurately discriminate between active and inactive TF binding events remain ambiguous. Recent studies have identified promoter-distal RNA polymerase II (RNAP2) binding at enhancer elements, suggesting that these interactions may serve as a marker for active regulatory sequences. Despite these correlative analyses, a thorough functional validation of these genomic co-occupancies is still lacking. To characterize the gene regulatory activity of DNA sequences underlying promoter-distal TF binding events that co-occur with RNAP2 and TF sites devoid of RNAP2 occupancy using a functional reporter assay, we performed cis-regulatory element sequencing (CRE-seq). We tested more than 1000 promoter-distal CCAAT/enhancer-binding protein beta (CEBPB)-bound sites in HepG2 and K562 cells, and found that CEBPB-bound sites co-occurring with RNAP2 were more likely to exhibit enhancer activity. CEBPB-bound sites further maintained substantial cell-type specificity, indicating that local DNA sequence can accurately convey cell-type–specific regulatory information. By comparing our CRE-seq results to a comprehensive set of genome annotations, we identified a variety of genomic features that are strong predictors of regulatory element activity and cell-type–specific activity. Collectively, our functional assay results indicate that RNAP2 occupancy can be used as a key genomic marker that can distinguish active from inactive TF bound sites.

Glucocorticoids are universally used in the treatment of acute lymphoblastic leukemia (ALL), and resistance to glucocorticoids in leukemia cells confers poor prognosis. To elucidate mechanisms of glucocorticoid resistance, we determined the prednisolone sensitivity of primary leukemia cells from 444 patients newly diagnosed with ALL and found significantly higher expression of CASP1 (encoding caspase 1) and its activator NLRP3 in glucocorticoid-resistant leukemia cells, resulting from significantly lower somatic methylation of the CASP1 and NLRP3 promoters. Overexpression of CASP1 resulted in cleavage of the glucocorticoid receptor, diminished the glucocorticoid-induced transcriptional response and increased glucocorticoid resistance. Knockdown or inhibition of CASP1 significantly increased glucocorticoid receptor levels and mitigated glucocorticoid resistance in CASP1-overexpressing ALL. Our findings establish a new mechanism by which the NLRP3-CASP1 inflammasome modulates cellular levels of the glucocorticoid receptor and diminishes cell sensitivity to glucocorticoids. The broad impact on the glucocorticoid transcriptional response suggests that this mechanism could also modify glucocorticoid effects in other diseases.

Most human transcription factors bind a small subset of potential genomic sites and often use different subsets in different cell types. To identify mechanisms that govern cell-type-specific transcription factor binding, we used an integrative approach to study estrogen receptor α (ER). We found that ER exhibits two distinct modes of binding. Shared sites, bound in multiple cell types, are characterized by high-affinity estrogen response elements (EREs), inaccessible chromatin, and a lack of DNA methylation, while cell-specific sites are characterized by a lack of EREs, co-occurrence with other transcription factors, and cell-type-specific chromatin accessibility and DNA methylation. These observations enabled accurate quantitative models of ER binding that suggest tethering of ER to one-third of cell-specific sites. The distinct properties of cell-specific binding were also observed with glucocorticoid receptor and for ER in primary mouse tissues, representing an elegant genomic encoding scheme for generating cell-type-specific gene regulation.

Recent functional characterisations of genome-wide association study (GWAS) loci suggest that cis-regulatory variation may be a common paradigm for complex disease susceptibility. Several studies point to a similar mechanism at the transcription factor 7-like 2 (TCF7L2) GWAS locus for type 2 diabetes. To address this possibility, we carried out an in vitro scan of this diabetes-associated locus to fine-map cis-regulatory sequences within this genomic interval.A systematic cell-based enhancer strategy was employed to interrogate all sequences within the 92 kb type-2-diabetes-association interval for cis-regulatory activity in a panel of cell lines (HCT-116, Neuro-2a, C2C12, U2OS, MIN6 and HepG2). We further evaluated chromatin state at a subset of these regions in HCT-116 and U2OS cells and examined allelic-specific enhancer properties at the type-2-diabetes-associated single nucleotide polymorphism (SNP) rs7903146.In total, we assigned cis-regulatory activity to approximately 30% (9/28) of constructs tested. Notably, a subset of enhancers was active across multiple cell lines and overlapped with key epigenetic markers suggestive of cis-regulatory sequences. We further replicated the allelic-specific properties for SNP rs7903146 in pancreatic beta cells and additionally demonstrate identical allelic-specific enhancer effects in other cell lines.These results provide a detailed map of cis-regulatory elements within the TCF7L2 GWAS locus and support the hypothesis of cis-regulatory variation leading to type 2 diabetes predisposition. The detection of allelic-specific effects for SNP rs7903146 in multiple cell lines further alludes to the likelihood of a peripheral defect in disease aetiology.

Genome-wide association studies (GWAS) have repeatedly shown an association between non-coding variants in the TCF7L2 locus and risk for type 2 diabetes (T2D), implicating a role for cis-regulatory variation within this locus in disease etiology. Supporting this hypothesis, we previously localized complex regulatory activity to the TCF7L2 T2D-associated interval using an in vivo bacterial artificial chromosome (BAC) enhancer-trapping reporter strategy. To follow-up on this broad initial survey of the TCF7L2 regulatory landscape, we performed a fine-mapping enhancer scan using in vivo mouse transgenic reporter assays. We functionally interrogated approximately 50% of the sequences within the T2D-associated interval, utilizing sequence conservation within this 92-kb interval to determine the regulatory potential of all evolutionary conserved sequences that exhibited conservation to the non-eutherian mammal opossum. Included in this study was a detailed functional interrogation of sequences spanning both protective and risk alleles of single nucleotide polymorphism (SNP) rs7903146, which has exhibited allele-specific enhancer function in pancreatic beta cells. Using these assays, we identified nine segments regulating various aspects of the TCF7L2 expression profile and that constitute nearly 70% of the sequences tested. These results highlight the regulatory complexity of this interval and support the notion that a TCF7L2 cis-regulatory disruption leads to T2D predisposition.

The comorbidity of type 2 diabetes (T2D) with several psychiatric diseases is well established. While environmental factors may partially account for these co-occurrences, common genetic susceptibilities could also be implicated in the confluence of these diseases. In support of shared genetic burdens, TCF7L2, the strongest genetic determinant for T2D risk in the human population, has been recently implicated in schizophrenia (SCZ) risk, suggesting that this may be one of many loci that pleiotropically influence both diseases. To investigate whether Tcf7l2 is involved in behavioral phenotypes in addition to its roles in glucose metabolism, we conducted several behavioral tests in mice with null alleles of Tcf7l2 or overexpressing Tcf7l2. We identified a role for Tcf7l2 in anxiety-like behavior and a dose-dependent effect of Tcf7l2 alleles on fear learning. None of the mutant mice showed differences in prepulse inhibition (PPI), which is a well-established endophenotype for SCZ. These results show that Tcf7l2 alters behavior in mice. Importantly, these differences are observed prior to the onset of detectable glucose metabolism abnormalities. Whether these differences are related to human anxiety-disorders or schizophrenia remains to be determined. These animal models have the potential to elucidate the molecular basis of psychiatric comorbidities in diabetes and should therefore be studied further.

Genome-wide association studies (GWAS) have consistently implicated noncoding variation within the TCF7L2 locus with type 2 diabetes (T2D) risk. While this locus represents the strongest genetic determinant for T2D risk in humans, it remains unclear how these noncoding variants affect disease etiology. To test the hypothesis that the T2D-associated interval harbors cis-regulatory elements controlling TCF7L2 expression, we conducted in vivo transgenic reporter assays to characterize the TCF7L2 regulatory landscape. We found that the 92-kb genomic interval associated with T2D harbors long-range enhancers regulating various aspects of the spatial–temporal expression patterns of TCF7L2, including expression in tissues involved in the control of glucose homeostasis. By selectively deleting this interval, we establish a critical role for these enhancers in robust TCF7L2 expression. To further determine whether variation in Tcf7l2 expression may lead to diabetes, we developed a Tcf7l2 copy-number allelic series in mice. We show that a null Tcf7l2 allele leads, in a dose-dependent manner, to lower glycemic profiles. Tcf7l2 null mice also display enhanced glucose tolerance coupled to significantly lowered insulin levels, suggesting that these mice are protected against T2D. Confirming these observations, transgenic mice harboring multiple Tcf7l2 copies and overexpressing this gene display reciprocal phenotypes, including glucose intolerance. These results directly demonstrate that Tcf7l2 plays a role in regulating glucose tolerance, suggesting that overexpression of this gene is associated with increased risk of T2D. These data highlight the role of enhancer elements as mediators of T2D risk in humans, strengthening the evidence that variation in cis-regulatory elements may be a paradigm for genetic predispositions to common disease.

Distal transcription enhancers are cis-regulatory elements that promote gene expression, enabling spatiotemporal control of genetic programs such as those required in metazoan developmental processes. Because of their importance, their disruption can lead to disease.