Transcription Factor Occupancy Research Articles

Transcription Factor Occupancy Can Mediate Active Turnover of DNA Methylation at Regulatory Regions

Distal regulatory elements, including enhancers, play a critical role in regulating gene activity. Transcription factor binding to these elements correlates with Low Methylated Regions (LMRs) in a process that is poorly understood. Here we ask whether and how actual occupancy of DNA-binding factors is linked to DNA methylation at the level of individual molecules. Using CTCF as an example, we observe that frequency of binding correlates with the likelihood of a demethylated state and sites of low occupancy display heterogeneous DNA methylation within the CTCF motif. In line with a dynamic model of binding and DNA methylation turnover, we find that 5-hydroxymethylcytosine (5hmC), formed as an intermediate state of active demethylation, is enriched at LMRs in stem and somatic cells. Moreover, a significant fraction of changes in 5hmC during differentiation occurs at these regions, suggesting that transcription factor activity could be a key driver for active demethylation. Since deletion of CTCF is lethal for embryonic stem cells, we used genetic deletion of REST as another DNA-binding factor implicated in LMR formation to test this hypothesis. The absence of REST leads to a decrease of hydroxymethylation and a concomitant increase of DNA methylation at its binding sites. These data support a model where DNA-binding factors can mediate turnover of DNA methylation as an integral part of maintenance and reprogramming of regulatory regions.

Exonic Transcription Factor Binding Directs Codon Choice and Affects Protein Evolution

Genomes contain both a genetic code specifying amino acids and a regulatory code specifying transcription factor (TF) recognition sequences. We used genomic deoxyribonuclease I footprinting to map nucleotide resolution TF occupancy across the human exome in 81 diverse cell types. We found that ~15% of human codons are dual-use codons ("duons") that simultaneously specify both amino acids and TF recognition sites. Duons are highly conserved and have shaped protein evolution, and TF-imposed constraint appears to be a major driver of codon usage bias. Conversely, the regulatory code has been selectively depleted of TFs that recognize stop codons. More than 17% of single-nucleotide variants within duons directly alter TF binding. Pervasive dual encoding of amino acid and regulatory information appears to be a fundamental feature of genome evolution.

An integrative analysis reveals functional targets of GATA6 transcriptional regulation in gastric cancer

Lineage-restricted transcription factors (TFs) are frequently mutated or overexpressed in cancer and contribute toward malignant behaviors, but the molecular bases of their oncogenic properties are largely unknown. Because TF activities are difficult to inhibit directly with small molecules, the genes and pathways they regulate might represent more tractable targets for drug therapy. We studied GATA6, a TF gene that is frequently amplified or overexpressed in gastric, esophageal, and pancreatic adenocarcinomas. GATA6-overexpressing gastric cancer cell lines cluster in gene expression space, separate from non-overexpressing lines. This expression clustering signifies a shared pathogenic group of genes that GATA6 may regulate through direct cis-element binding. We used chromatin immunoprecipation and sequencing (ChIP-seq) to identify GATA6-bound genes and considered TF occupancy in relation to genes that respond to GATA6 depletion in cell lines and track with GATA6 mRNA (synexpression groups) in primary gastric cancers. Among other cellular functions, GATA6-occupied genes control apoptosis and govern M-phase of the cell cycle. Depletion of GATA6 reduced levels of the latter transcripts and arrested cells in G2 and M phases of the cell cycle. Synexpression in human tumor samples identified likely direct transcriptional targets substantially better than consideration only of transcripts that respond to GATA6 loss in cultured cells. Candidate target genes responded to loss of GATA6 or its homolog GATA4 and even more to depletion of both proteins. Many GATA6-dependent genes lacked nearby binding sites but several strongly dependent, synexpressed, and GATA6-bound genes encode TFs such as MYC, HES1, RARB, and CDX2. Thus, many downstream effects occur indirectly through other TFs and GATA6 activity in gastric cancer is partially redundant with GATA4. This integrative analysis of locus occupancy, gene dependency, and synexpression provides a functional signature of GATA6-overexpressing gastric cancers, revealing both limits and new therapeutic directions for a challenging and frequently fatal disease.

Refined DNase-seq protocol and data analysis reveals intrinsic bias in transcription factor footprint identification.

DNase-seq is a powerful technique for identifying cis-regulatory elements across the genome. We studied the key experimental parameters to optimize the performance of DNase-seq. We found that sequencing short 50-100bp fragments that accumulate in long inter-nucleosome linker regions is more efficient for identifying transcription factor binding sites than using longer fragments. We also assessed the potential of DNase-seq to predict transcription factor occupancy through the generation of nucleotide-resolution transcription factor footprints. In modeling the sequence-specific DNaseI cutting bias we found a surprisingly strong effect that varied over more than two orders of magnitude. This confounds DNaseI footprint analysis to the extent that the nucleotide resolution cleavage patterns at most transcription factor binding sites are derived from intrinsic DNaseI cleavage bias rather than from specific protein-DNA interactions. In contrast, quantitative comparison of DNaseI hypersensitivity between states can predict transcription factor occupancy associated with particular biological perturbations.

Discrimination between thermodynamic models ofcis-regulation using transcription factor occupancy data

Many studies have identified binding preferences for transcription factors (TFs), but few have yielded predictive models of how combinations of transcription factor binding sites generate specific levels of gene expression. Synthetic promoters have emerged as powerful tools for generating quantitative data to parameterize models of combinatorial cis-regulation. We sought to improve the accuracy of such models by quantifying the occupancy of TFs on synthetic promoters in vivo and incorporating these data into statistical thermodynamic models of cis-regulation. Using chromatin immunoprecipitation-seq, we measured the occupancy of Gcn4 and Cbf1 in synthetic promoter libraries composed of binding sites for Gcn4, Cbf1, Met31/Met32 and Nrg1. We measured the occupancy of these two TFs and the expression levels of all promoters in two growth conditions. Models parameterized using only expression data predicted expression but failed to identify several interactions between TFs. In contrast, models parameterized with occupancy and expression data predicted expression data, and also revealed Gcn4 self-cooperativity and a negative interaction between Gcn4 and Nrg1. Occupancy data also allowed us to distinguish between competing regulatory mechanisms for the factor Gcn4. Our framework for combining occupancy and expression data produces predictive models that better reflect the mechanisms underlying combinatorial cis-regulation of gene expression.

Role of SWI/SNF in acute leukemia maintenance and enhancer-mediated Myc regulation

Cancer cells frequently depend on chromatin regulatory activities to maintain a malignant phenotype. Here, we show that leukemia cells require the mammalian SWI/SNF chromatin remodeling complex for their survival and aberrant self-renewal potential. While Brg1, an ATPase subunit of SWI/SNF, is known to suppress tumor formation in several cell types, we found that leukemia cells instead rely on Brg1 to support their oncogenic transcriptional program, which includes Myc as one of its key targets. To account for this context-specific function, we identify a cluster of lineage-specific enhancers located 1.7 Mb downstream from Myc that are occupied by SWI/SNF as well as the BET protein Brd4. Brg1 is required at these distal elements to maintain transcription factor occupancy and for long-range chromatin looping interactions with the Myc promoter. Notably, these distal Myc enhancers coincide with a region that is focally amplified in ∼3% of acute myeloid leukemias. Together, these findings define a leukemia maintenance function for SWI/SNF that is linked to enhancer-mediated gene regulation, providing general insights into how cancer cells exploit transcriptional coactivators to maintain oncogenic gene expression programs.

Differences In The Genome-Wide Epigenetic Signatures Of mRNA and Long Non-Coding RNA Genes In Mouse Erythroblasts and Megakaryocytes

The ENCODE project has demonstrated that epigenetic signatures, including DNA methylation and transcription factor (TF) occupancy, define gene expression. However, ENCODE was constructed using static cells that were not capable of further differentiation. We hypothesize that specific epigenetic profiles are associated with erythroid and megakaryocytic differentiation. To test this hypothesis, we isolated primary erythroblasts (EBs) and megakaryocytes (Megs) from mouse bone marrow by flow cytometry and prepared: 1) DNA for genome-wide methylation analysis using MBD2 Methyl-Seq; 2) RNA for both RNA-Seq analysis and microarray analysis of novel and annotated lncRNA expression levels; and 3) chromatin for genome-wide chromatin immunoprecipitation (ChIP-Seq) analysis of occupancy by the TFs GATA1 and NFE2. We developed the web-based high-throughput sequencing tool suite SigSeeker (http://sigseeker.org) to predict regions of methylation and TF occupancy across the genome.High-confidence methylation, GATA1, and NFE2 profiles, represented by the intersection of two independent EB and Meg biological replicates, are shown in Table 1. Of the approximate 100,000 methylated regions in EBs and Megs, 45% were shared between the two cell types. Unlike methylation, GATA1 and NFE2 occupancy showed strong cell type-specific profiles, with most GATA1 occupied sites (79%) being EB-specific, and most NFE2 occupied sites (72%) being Meg-specific. While 26% of EB-specific GATA1 peaks were co-occupied by NFE2, co-occupancy by GATA1 and NFE2 in Megs was rare (0.6%).Table 1Methylation, GATA1 and NFE2 Peaks in EBs and MegsEB-specificCommonMeg-specificMethylation61,10645,18052,822GATA1 Only4,6691,23699NFE2 Only1,3501,2483,147GATA1 and NFE23,3821,03420We developed a second web-based tool called SBR (http://sbrblood.msseeker.org) to correlate ChIP-Seq and Methyl-Seq profiles with RNA-Seq and lncRNA transcriptional data sets. Almost 95% of RefSeq (coding) genes with methylation in the promoter regions were not expressed. Fewer than 5% of methylated RefSeq promoters also had TF occupancy. Unlike promoters, the bodies of RefSeq genes had a high degree of overlap between methylation and TF binding. In EBs, 38% of RefSeq genes with GATA1 occupancy in the body were methylated. Of genes with this profile, 81% were transcriptionally silent. In Megs, 42% of RefSeq genes with NFE2 occupancy in the body were methylated. However, 80% of these Meg-specific genes were transcriptionally active. In contrast to RefSeq genes, lncRNA genes have a different signature. More than 85% of EB-expressed lncRNA promoters are methylated and ∼25% of these promoters are occupied by GATA1 and/or NFE2. Over 90% of EB-expressed lncRNA gene bodies occupied by GATA1 and/or NFE2 are methylated. In contrast less than 10% of Meg-expressed lncRNA promoters are methylated and ∼only 2% of these promoters are occupied by NFE2. However, 90% of Meg-expressed lncRNA gene bodies occupied by GATA1 and/or NFE2 are methylated.Ingenuity IPA analysis of the transcriptional profiles associated with different epigenetic signatures revealed differentially regulated cellular pathways. In EBs, GATA1 occupancy in the promoters of silent genes was associated with cardiovascular (p≤ 10-7) and nervous system development (p≤ 10-7). In Megs, NFE2 promoter occupancy was associated with active genes involved in nucleic acid metabolism (p≤ 10-3) and nervous system development (p≤ 10-5). Genes expressed in both EBs and Megs were involved in transcription (p≤ 10-6), cell cycle progression (p≤ 10-5), and decreased hypoplasia (p≤ 10-20). GATA1 and NFE2 co-occupied genes expressed in both EBs and Megs were associated with suppression of bone-related (p≤ 10-6) and neuron-related transcripts (p≤ 10-8).In summary, we show that TF occupancy and methylation significantly overlap in RefSeq gene bodies, but not in promoters. These profiles have revealed that GATA1 occupancy, independent of NFE2 co-occupancy, correlates with EB-specific transcriptional silencing, whereas Meg-specific transcriptional activation is associated with NFE2 occupancy. In contrast, over 90% active and TF-occupied lncRNA (both novel and previously annotated) gene bodies are methylated. These epigenetic correlations will be important for future studies assessing the regulation of mRNAs and lncRNAs. Disclosures:No relevant conflicts of interest to declare.

HBS1L-MYB intergenic Variants Modulate Fetal Hemoglobin Via Long-Range MYB Enhancers

Three quantitative trait loci (QTLs) modifying fetal hemoglobin (HbF) levels have been identified, and these have been shown to have a predictive value of disease severity in β thalassemia and sickle cell disease in diverse ethnic groups. One of the HbF QTLs which consists of a set of common intergenic single nucleotide polymorphisms (SNPs) in the HBS1L-MYB intergenic region on chromosome 6q23, has also been consistently identified as highly associated with clinically important human erythroid traits. Despite extensive genetic evidence, a clear mechanistic basis for the association between the intergenic SNPs and erythroid biology has remained elusive, although the two flanking genes (HBS1L and MYB) are candidate target genes.Here, we set out to characterize the regulatory potential of the human HBS1L-MYB intergenic region in detail and to investigate its functional impact on the erythroid phenotype-associated variants.We profiled chromatin occupancy of the key erythroid LDB1 transcription factor (TF) complex in primary human erythroid progenitors (HEPs) using chromatin immunoprecipitation coupled to high-throughput sequencing (ChIP-Seq) and quantitative PCR (ChIP-qPCR). We detected an intergenic cluster containing 7 binding sites for the LDB1–complex, characterized by strong binding and co-occupancy of core complex proteins LDB1, GATA1, TAL1 and ETO2. One of these sites was co-occupied by the erythroid-specific TF KLF1, a protein which was also found to bind the murine intergenic region. Depletion of LDB1, TAL1 and KLF1 in K562 cells using RNA interference resulted in a specific downregulation of MYB expression while leaving HBS1L levels unaffected, demonstrating that the erythroid TFs occupying the intergenic enhancers are required for MYB expression.Using chromosome conformation capture (3C) coupled to high-throughput sequencing (3C-Seq), we profiled higher order chromatin structure within the locus, and detected several strong chromatin co-associations between the MYB promoter and intergenic sequences, almost all of which correlated with TF binding. The binding activities were distributed over a conserved core region identical to a 24-kb interval containing genetic variants in strong genetic association with erythroid traits in human populations. This block of SNPs is referred to as HBS1L-MYB intergenic polymorphisms block-2 (HMIP-2). The SNPs within HMIP-2 clustered in 2 regions positioned directly under the 2 LDB1-complex ChIP-Seq peaks, at -84 and -71 kb from the MYB transcription start site. Using allele-specific ChIP in K562 cells, we observed diminished (25-50%) binding of LDB1, GATA1, TAL1 and KLF1 to the minor rs66650371 allele, showing that rs66650371 affects local TF binding. An allele-specific 3C analysis also showed reduced interactions between the minor rs66650371 allele at -84 and MYB. We validated and expanded observations made in K562 cells using primary human cells: 1. HEPs from high HbF individuals homozygous for all minor alleles of the phenotype-associated -84kb and -71kb intergenic variants in the conserved core (‘SNP/SNP') showed consistently lower MYB levels throughout phase II of the culture as compared to wildtype control cells (‘WT/WT'); 37% lower MYB on average; 2. ChIP experiments on SNP/SNP and WT/WT HEPs showed reduced binding of GATA1 and KLF1 to the -84 and -71 regulatory elements (containing the associated variants), the results were further confirmed by allele-specific ChIP assays in SNP/WT HEPs; 3. 3C-qPCR assays on cultured SNP/SNP and WT/WT cells demonstrated diminished looping between the -84 element and the MYB promoter in SNP/SNP individuals; 4. SNP/WT HEPs also showed allelic imbalance of MYB but not HBS1L transcripts when compared to controls (WT/WT and SNP/SNP HEPs).In conclusion, we show that the HBS1L-MYB gene-free interval contains distal enhancer elements that interact with MYB, a critical regulator of erythroid development and HbF levels. Key variants in the intergenic interval affect MYB expression by reducing TF binding to its regulatory elements and disrupting long-range enhancer-gene interaction. Our study identifies the first causal link between the 6q23 HbF QTL and MYB regulation, provides novel insights into the molecular control of clinically important haematological traits and adds another layer of complexity to the regulation of MYB, suggesting potential targets for therapeutic intervention. Disclosures:Thein: Sangart: Consultancy; Shire: Research Funding; Novartis: Speakers Bureau.

Coupling transcription factor occupancy to nucleosome architecture with DNase-FLASH

It is currently not possible to resolve the genome-wide relationship of transcription factors (TFs) and nucleosomes at the level of individual chromatin templates despite rapidly increasing data on TF and nucleosome occupancy in the human genome. Here we describe DNase I-released fragment-length analysis of hypersensitivity (DNase-FLASH), an approach that directly couples mapping of TF occupancy, via quantification of DNA microfragments released from individual TF recognition sites in regulatory DNA, to the surrounding nucleosome architecture, via analysis of larger DNA fragments, in a single assay. DNase-FLASH enables coupling of individual TF footprints to nucleosome occupancy, identifying TFs that precisely demarcate the regulatory DNA-nucleosome interface.

Integrating and mining the chromatin landscape of cell-type specificity using self-organizing maps

We tested whether self-organizing maps (SOMs) could be used to effectively integrate, visualize, and mine diverse genomics data types, including complex chromatin signatures. A fine-grained SOM was trained on 72 ChIP-seq histone modifications and DNase-seq data sets from six biologically diverse cell lines studied by The ENCODE Project Consortium. We mined the resulting SOM to identify chromatin signatures related to sequence-specific transcription factor occupancy, sequence motif enrichment, and biological functions. To highlight clusters enriched for specific functions such as transcriptional promoters or enhancers, we overlaid onto the map additional data sets not used during training, such as ChIP-seq, RNA-seq, CAGE, and information on cis-acting regulatory modules from the literature. We used the SOM to parse known transcriptional enhancers according to the cell-type-specific chromatin signature, and we further corroborated this pattern on the map by EP300 (also known as p300) occupancy. New candidate cell-type-specific enhancers were identified for multiple ENCODE cell types in this way, along with new candidates for ubiquitous enhancer activity. An interactive web interface was developed to allow users to visualize and custom-mine the ENCODE SOM. We conclude that large SOMs trained on chromatin data from multiple cell types provide a powerful way to identify complex relationships in genomic data at user-selected levels of granularity.

The next generation of transcription factor binding site prediction.

Finding where transcription factors (TFs) bind to the DNA is of key importance to decipher gene regulation at a transcriptional level. Classically, computational prediction of TF binding sites (TFBSs) is based on basic position weight matrices (PWMs) which quantitatively score binding motifs based on the observed nucleotide patterns in a set of TFBSs for the corresponding TF. Such models make the strong assumption that each nucleotide participates independently in the corresponding DNA-protein interaction and do not account for flexible length motifs. We introduce transcription factor flexible models (TFFMs) to represent TF binding properties. Based on hidden Markov models, TFFMs are flexible, and can model both position interdependence within TFBSs and variable length motifs within a single dedicated framework. The availability of thousands of experimentally validated DNA-TF interaction sequences from ChIP-seq allows for the generation of models that perform as well as PWMs for stereotypical TFs and can improve performance for TFs with flexible binding characteristics. We present a new graphical representation of the motifs that convey properties of position interdependence. TFFMs have been assessed on ChIP-seq data sets coming from the ENCODE project, revealing that they can perform better than both PWMs and the dinucleotide weight matrix extension in discriminating ChIP-seq from background sequences. Under the assumption that ChIP-seq signal values are correlated with the affinity of the TF-DNA binding, we find that TFFM scores correlate with ChIP-seq peak signals. Moreover, using available TF-DNA affinity measurements for the Max TF, we demonstrate that TFFMs constructed from ChIP-seq data correlate with published experimentally measured DNA-binding affinities. Finally, TFFMs allow for the straightforward computation of an integrated TF occupancy score across a sequence. These results demonstrate the capacity of TFFMs to accurately model DNA-protein interactions, while providing a single unified framework suitable for the next generation of TFBS prediction.

ERNAs Promote Transcription by Establishing Chromatin Accessibility at Defined Genomic Loci

Transcription factors and DNA regulatory binding motifs are fundamental components of the gene regulatory network. Here, by using genome-wide binding profiling, we show extensive occupancy of transcription factors of myogenesis (MyoD and Myogenin) at extragenic enhancer regions coinciding with RNA synthesis (i.e., eRNA). In particular, multiple regions were transcribed to eRNA within the regulatory region of MYOD1, including previously characterized distal regulatory regions (DRR) and core enhancer (CE). While (CE)RNA enhanced RNA polymerase II (Pol II) occupancy and transcription at MYOD1, (DRR)RNA acted to activate the downstream myogenic genes. The deployment of transcriptional machinery to appropriate loci is contingent on chromatin accessibility, a rate-limiting step preceding Pol II assembly. By nuclease sensitivity assay, we found that eRNAs regulate genomic access of the transcriptional complex to defined regulatory regions. In conclusion, our data suggest that eRNAs contribute to establishing acell-type-specific transcriptional circuitry by directing chromatin-remodeling events.

Stepwise acetyltransferase association and histone acetylation at the Myod1 locus during myogenic differentiation

While chromatin modifications can offer a useful readout for enhancer activities, it is less clear whether these modification marks are a cause or consequence of transcription factor occupancy and enhancer activation. We have examined in details the temporal events of acetyltransferase associations and histone acetylations at different regulatory regions of the Myod1 locus. Our studies demonstrate that the histone acetyltransferase (HAT) p300 is stepwise enriched at distinct Myod1 regulatory regions during myogenic differentiation. This enrichment of p300 is associated with increased histone acetylation in a discrete pattern. Inhibition of p300 HAT activity impedes myogenic differentiation, which is coupled with decreased histone acetylation at specific Myod1 regulatory regions. We show for the first time that p300 is directly involved in the early regulation of Myod1 enhancer, and provide molecular insights into how p300 HAT activity and histone acetylation are related to enhancer activation and, consequently, gene transcription.

Computational Identification of Diverse Mechanisms Underlying Transcription Factor-DNA Occupancy

ChIP-based genome-wide assays of transcription factor (TF) occupancy have emerged as a powerful, high-throughput method to understand transcriptional regulation, especially on a global scale. This has led to great interest in the underlying biochemical mechanisms that direct TF-DNA binding, with the ultimate goal of computationally predicting a TF's occupancy profile in any cellular condition. In this study, we examined the influence of various potential determinants of TF-DNA binding on a much larger scale than previously undertaken. We used a thermodynamics-based model of TF-DNA binding, called “STAP,” to analyze 45 TF-ChIP data sets from Drosophila embryonic development. We built a cross-validation framework that compares a baseline model, based on the ChIP'ed (“primary”) TF's motif, to more complex models where binding by secondary TFs is hypothesized to influence the primary TF's occupancy. Candidates interacting TFs were chosen based on RNA-SEQ expression data from the time point of the ChIP experiment. We found widespread evidence of both cooperative and antagonistic effects by secondary TFs, and explicitly quantified these effects. We were able to identify multiple classes of interactions, including (1) long-range interactions between primary and secondary motifs (separated by ≤150 bp), suggestive of indirect effects such as chromatin remodeling, (2) short-range interactions with specific inter-site spacing biases, suggestive of direct physical interactions, and (3) overlapping binding sites suggesting competitive binding. Furthermore, by factoring out the previously reported strong correlation between TF occupancy and DNA accessibility, we were able to categorize the effects into those that are likely to be mediated by the secondary TF's effect on local accessibility and those that utilize accessibility-independent mechanisms. Finally, we conducted in vitro pull-down assays to test model-based predictions of short-range cooperative interactions, and found that seven of the eight TF pairs tested physically interact and that some of these interactions mediate cooperative binding to DNA.

Abstract 2314: The chromatin structure regulates selective glucocorticoid receptor recruitment to the genome in intact tissue.

Abstract The dissection of regulatory pathways in physiological and pathophysiological states and during pharmacological intervention necessitates the application of biochemical and genomic approaches in intact target tissues. The hypothalamus-pituitary-adrenal axis is a highly conserved neuroendocrine pathway that results in the synthesis of glucocorticoid hormones from the adrenal gland to control reactions to stress and inflammatory processes. Exploiting a mouse model in which animals have been adrenalectomized to remove endogenous glucocorticoids, we monitor changes of the chromatin landscape in liver, following a single glucocorticoid injection, using ChIP-seq and DNase-seq methodologies. We observe that GR exclusively occupies accessible chromatin, demonstrating that accessible chromatin can be highly predictive of transcription factor occupancy in tissue. Importantly, changes in RNA polymerase II recruitment to GR regulated genes are associated with proximal regions of induced and repressed chromatin accessibility. Underscoring the highly tissue-specific GR-dependent transcriptional response, GR binding sites in other cell types are occluded in liver through highly positioned nucleosomes and inaccessible chromatin. Remarkably, we demonstrate that sixty two percent of GR binding sites in liver are associated with C/EBPbeta occupancy. At the majority of these sites chromatin is pre-programmed, suggesting a priming function of C/EBPbeta in GR recruitment. However, at a subset of GR binding sites C/EBPbeta loading is assisted by GR recruitment in conjunction with de novo chromatin remodeling. Disruption of C/EBPbeta binding to chromatin using a dominant negative C/EBP targeted to the liver results in an attenuation of pre-programmed chromatin accessibility, GR recruitment and GR induced chromatin remodeling only at sites occupied by GR and C/EBP. Our results demonstrate a highly cooperative mechanism by which C/EBP regulates selective GR binding to the genome in liver tissue and provides a mechanism for phenotypic similarities described for GR and C/EBP mouse models. We suggest that the selective targeting of GR in other tissues is likely mediated by the combined action of cell-specific priming proteins and chromatin remodelers. Citation Format: Lars Grontved, Sam John, Songjoon Baek, Ying Liu, John Buckley, Charles Vinson, Greti Aguilera, Gordon L. Hager. The chromatin structure regulates selective glucocorticoid receptor recruitment to the genome in intact tissue. [abstract]. In: Proceedings of the 104th Annual Meeting of the American Association for Cancer Research; 2013 Apr 6-10; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2013;73(8 Suppl):Abstract nr 2314. doi:10.1158/1538-7445.AM2013-2314

HCFC1 is a common component of active human CpG-island promoters and coincides with ZNF143, THAP11, YY1, and GABP transcription factor occupancy

In human transcriptional regulation, DNA-sequence-specific factors can associate with intermediaries that orchestrate interactions with a diverse set of chromatin-modifying enzymes. One such intermediary is HCFC1 (also known as HCF-1). HCFC1, first identified in herpes simplex virus transcription, has a poorly defined role in cellular transcriptional regulation. We show here that, in HeLa cells, HCFC1 is observed bound to 5400 generally active CpG-island promoters. Examination of the DNA sequences underlying the HCFC1-binding sites revealed three sequence motifs associated with the binding of (1) ZNF143 and THAP11 (also known as Ronin), (2) GABP, and (3) YY1 sequence-specific transcription factors. Subsequent analysis revealed colocalization of HCFC1 with these four transcription factors at ∼90% of the 5400 HCFC1-bound promoters. These studies suggest that a relatively small number of transcription factors play a major role in HeLa-cell transcriptional regulation in association with HCFC1.

Taking into account nucleosomes for predicting gene expression

The eukaryotic genome is organized in a chain of nucleosomes that consist of 145–147bp of DNA wrapped around a histone octamer protein core. Binding of transcription factors (TF) to nucleosomal DNA is frequently impeded, which makes it a challenging task to calculate TF occupancy at a given regulatory genomic site for predicting gene expression. Here, we review methods to calculate TF binding to DNA in the presence of nucleosomes. The main theoretical problems are (i) the computation speed that is becoming a bottleneck when partial unwrapping of DNA from the nucleosome is considered, (ii) the perturbation of the binding equilibrium by the activity of ATP-dependent chromatin remodelers, which translocate nucleosomes along the DNA, and (iii) the model parameterization from high-throughput sequencing data and fluorescence microscopy experiments in living cells. We discuss strategies that address these issues to efficiently compute transcription factor binding in chromatin.

Temporal mapping of CEBPA and CEBPB binding during liver regeneration reveals dynamic occupancy and specific regulatory codes for homeostatic and cell cycle gene batteries

Dynamic shifts in transcription factor binding are central to the regulation of biological processes by allowing rapid changes in gene transcription. However, very few genome-wide studies have examined how transcription factor occupancy is coordinated temporally in vivo in higher animals. Here, we quantified the genome-wide binding patterns of two key hepatocyte transcription factors, CEBPA and CEBPB (also known as C/EBPalpha and C/EBPbeta), at multiple time points during the highly dynamic process of liver regeneration elicited by partial hepatectomy in mouse. Combining these profiles with RNA polymerase II binding data, we find three temporal classes of transcription factor binding to be associated with distinct sets of regulated genes involved in the acute phase response, metabolic/homeostatic functions, or cell cycle progression. Moreover, we demonstrate a previously unrecognized early phase of homeostatic gene expression prior to S-phase entry. By analyzing the three classes of CEBP bound regions, we uncovered mutually exclusive sets of sequence motifs, suggesting temporal codes of CEBP recruitment by differential cobinding with other factors. These findings were validated by sequential ChIP experiments involving a panel of central transcription factors and/or by comparison to external ChIP-seq data. Our quantitative investigation not only provides in vivo evidence for the involvement of many new factors in liver regeneration but also points to similarities in the circuitries regulating self-renewal of differentiated cells. Taken together, our work emphasizes the power of global temporal analyses of transcription factor occupancy to elucidate mechanisms regulating dynamic biological processes in complex higher organisms.

Transcription-factor occupancy at HOT regions quantitatively predicts RNA polymerase recruitment in five human cell lines.

BackgroundHigh-occupancy target (HOT) regions are compact genome loci occupied by many different transcription factors (TFs). HOT regions were initially defined in invertebrate model organisms, and we here show that they are a ubiquitous feature of the human gene-regulation landscape.ResultsWe identified HOT regions by a comprehensive analysis of ChIP-seq data from 96 DNA-associated proteins in 5 human cell lines. Most HOT regions co-localize with RNA polymerase II binding sites, but many are not near the promoters of annotated genes. At HOT promoters, TF occupancy is strongly predictive of transcription preinitiation complex recruitment and moderately predictive of initiating Pol II recruitment, but only weakly predictive of elongating Pol II and RNA transcript abundance. TF occupancy varies quantitatively within human HOT regions; we used this variation to discover novel associations between TFs. The sequence motif associated with any given TF’s direct DNA binding is somewhat predictive of its empirical occupancy, but a great deal of occupancy occurs at sites without the TF’s motif, implying indirect recruitment by another TF whose motif is present.ConclusionsMammalian HOT regions are regulatory hubs that integrate the signals from diverse regulatory pathways to quantitatively tune the promoter for RNA polymerase II recruitment.

The B cell transcription program mediates hypomethylation and overexpression of key genes in Epstein-Barr virus-associated proliferative conversion

BackgroundEpstein-Barr virus (EBV) infection is a well characterized etiopathogenic factor for a variety of immune-related conditions, including lymphomas, lymphoproliferative disorders and autoimmune diseases. EBV-mediated transformation of resting B cells to proliferating lymphoblastoid cells occurs in early stages of infection and is an excellent model for investigating the mechanisms associated with acquisition of unlimited growth.ResultsWe investigated the effects of experimental EBV infection of B cells on DNA methylation profiles by using high-throughput analysis. Remarkably, we observed hypomethylation of around 250 genes, but no hypermethylation. Hypomethylation did not occur at repetitive sequences, consistent with the absence of genomic instability in lymphoproliferative cells. Changes in methylation only occurred after cell divisions started, without the participation of the active demethylation machinery, and were concomitant with acquisition by B cells of the ability to proliferate. Gene Ontology analysis, expression profiling, and high-throughput analysis of the presence of transcription factor binding motifs and occupancy revealed that most genes undergoing hypomethylation are active and display the presence of NF-κB p65 and other B cell-specific transcription factors. Promoter hypomethylation was associated with upregulation of genes relevant for the phenotype of proliferating lymphoblasts. Interestingly, pharmacologically induced demethylation increased the efficiency of transformation of resting B cells to lymphoblastoid cells, consistent with productive cooperation between hypomethylation and lymphocyte proliferation.ConclusionsOur data provide novel clues on the role of the B cell transcription program leading to DNA methylation changes, which we find to be key to the EBV-associated conversion of resting B cells to proliferating lymphoblasts.