CCCTC-Binding Factor Translates Interleukin 2- and α-Ketoglutarate-Sensitive Metabolic Changes in T Cells into Context-Dependent Gene Programs.

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.

Occupancy maps of 208 chromatin-associated proteins in one human cell type.

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.

Distinct properties of cell-type-specific and shared transcription factor binding sites.

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.

A genome-wide interactome of DNA-associated proteins in the human liver.

Large-scale efforts like the 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 differ 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 identify and analyze more than 440,000 binding sites using ChIP-seq data for 20 DAPs in two human liver tissue samples. We integrated binding data with transcriptome and phased WGS data to investigate allelic DAP interactions and the impact of heterozygous sequence variation on the expression of neighboring genes. Our tissue-based data set exhibits binding patterns more consistent with liver biology than cell lines, and we describe uses of these data to better prioritize impactful noncoding variation. Collectively, our rich data set offers novel insights into genome function in human liver tissue and provides a valuable resource for assessing disease-related disruptions.

Decoding transcriptional enhancers: Evolving from annotation to functional interpretation.

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.

CETCh-seq: CRISPR epitope tagging ChIP-seq of DNA-binding proteins.

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.

Promoter-distal RNA polymerase II binding discriminates active from inactive CCAAT/ enhancer-binding protein beta binding sites.

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.

Evidence of non-pancreatic beta cell-dependent roles of Tcf7l2 in the regulation of glucose metabolism in mice.

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.

Alterations in TCF7L2 expression define its role as a key regulator of glucose metabolism.

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.

Investigation of inherited noncoding genetic variation impacting the pharmacogenomics of childhood acute lymphoblastic leukemia treatment.

Defining genetic factors impacting chemotherapy failure can help to better predict response and identify drug resistance mechanisms. However, there is limited understanding of the contribution of inherited noncoding genetic variation on inter-individual differences in chemotherapy response in childhood acute lymphoblastic leukemia (ALL). Here we map inherited noncoding variants associated with treatment outcome and/or chemotherapeutic drug resistance to ALL cis-regulatory elements and investigate their gene regulatory potential and target gene connectivity using massively parallel reporter assays and three-dimensional chromatin looping assays, respectively. We identify 54 variants with transcriptional effects and high-confidence gene connectivity. Additionally, functional interrogation of the top variant, rs1247117, reveals changes in chromatin accessibility, PU.1 binding affinity and gene expression, and deletion of the genomic interval containing rs1247117 sensitizes cells to vincristine. Together, these data demonstrate that noncoding regulatory variants associated with diverse pharmacological traits harbor significant effects on allele-specific transcriptional activity and impact sensitivity to antileukemic agents.

Mutual antagonism between glucocorticoid and canonical Wnt signaling pathways in B-cell acute lymphoblastic leukemia.

Glucocorticoids (GCs; i.e., steroids) are important chemotherapeutic agents in the treatment of B-cell precursor acute lymphoblastic leukemia (B-ALL) and de novo GC 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)-bound cis -regulatory elements controlling TLE1 expression in GC-resistant primary B-ALL cells from patients. TLE1 is 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 of cis -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 that cis -regulatory disruptions at TLE1 are 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.

Epigenomic mapping in B-cell acute lymphoblastic leukemia identifies transcriptional regulators and noncoding variants promoting distinct chromatin architectures.

B-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-associated cis-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.

Functional investigation of inherited noncoding genetic variation impacting the pharmacogenomics of childhood acute lymphoblastic leukemia treatment.

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 ALL cis-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.

Epigenomic mapping reveals distinct B cell acute lymphoblastic leukemia chromatin architectures and regulators.

B cell lineage acute lymphoblastic leukemia (B-ALL) is composed of diverse molecular subtypes, and while transcriptional and DNA methylation profiling has been extensively examined, the chromatin landscape is not well characterized for many subtypes. We therefore mapped chromatin accessibility using ATAC-seq in primary B-ALL cells from 156 patients spanning ten molecular subtypes and present this dataset as a resource. Differential chromatin accessibility and transcription factor (TF) footprint profiling were employed and identified B-ALL cell of origin, TF-target gene interactions enriched in B-ALL, and key TFs associated with accessible chromatin sites preferentially active in B-ALL. We further identified over 20% of accessible chromatin sites exhibiting strong subtype enrichment and candidate TFs that maintain subtype-specific chromatin architectures. Over 9,000 genetic variants were uncovered, contributing to variability in chromatin accessibility among patient samples. Our data suggest that distinct chromatin architectures are driven by diverse TFs and inherited genetic variants that promote unique gene-regulatory networks.

Molecular Mechanisms of ARID5B-Mediated Genetic Susceptibility to Acute Lymphoblastic Leukemia.

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.

Epigenetic activation of the FLT3 gene by ZNF384 fusion confers a therapeutic susceptibility in acute lymphoblastic leukemia.

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.

Epigenomic profiling of glucocorticoid responses identifies cis-regulatory disruptions impacting steroid resistance in childhood acute lymphoblastic leukemia.

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.

Amino acid stress response genes promote L-asparaginase resistance in pediatric acute lymphoblastic leukemia.

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.

Genome-Wide Association Study of Susceptibility Loci for TCF3-PBX1 Acute Lymphoblastic Leukemia in Children.

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.

Profiling chromatin accessibility in pediatric acute lymphoblastic leukemia identifies subtype-specific chromatin landscapes and gene regulatory networks.

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.