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.

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.

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.

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.

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.

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.

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.

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.

Integrative genomic analyses reveal mechanisms of glucocorticoid resistance in acute lymphoblastic leukemia.

Identification of genomic and epigenomic determinants of drug resistance provides important insights for improving cancer treatment. Using agnostic genome-wide interrogation of mRNA and miRNA expression, DNA methylation, SNPs, CNAs and SNVs/Indels in primary human acute lymphoblastic leukemia cells, we identified 463 genomic features associated with glucocorticoid resistance. Gene-level aggregation identified 118 overlapping genes, 15 of which were confirmed by genome-wide CRISPR screen. Collectively, this identified 30 of 38 (79%) known glucocorticoid-resistance genes/miRNAs and all 38 known resistance pathways, while revealing 14 genes not previously associated with glucocorticoid-resistance. Single cell RNAseq and network-based transcriptomic modelling corroborated the top previously undiscovered gene, CELSR2. Manipulation of CELSR2 recapitulated glucocorticoid resistance in human leukemia cell lines and revealed a synergistic drug combination (prednisolone and venetoclax) that mitigated resistance in mouse xenograft models. These findings illustrate the power of an integrative genomic strategy for elucidating genes and pathways conferring drug resistance in cancer cells.

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.

Bromodomain-Selective BET Inhibitors Are Potent Antitumor Agents against MYC-Driven Pediatric Cancer.

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.

Batf Pioneers the Reorganization of Chromatin in Developing Effector T Cells via Ets1-Dependent Recruitment of Ctcf.

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.

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.

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.

Genomic regulation of invasion by STAT3 in triple negative breast cancer.

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.