Current DNA motif models can lead to incorrect hypotheses about the genomic recruitment of transcription factors
Ontology highlight
ABSTRACT:
ORGANISM(S): synthetic construct Homo sapiens
PROVIDER: GSE61854 | GEO | 2014/11/04
SECONDARY ACCESSION(S): PRJNA262558
REPOSITORIES: GEO
Similar Datasets
Project description:Recent genome-scale ChIP-chip studies of transcription factors have shown that a low percentage of experimentally determined binding sites contain the consensus motif for the immunoprecipitated factor. In most cases, differences between in vivo target sites that contain or lack a consensus motif have not been explored. We have previously shown that most sites to which E2F family members are bound in vivo do not contain E2F consensus motifs. The main purpose of this study was to develop an understanding of how E2F binding specificity is achieved in vivo. In particular, we have addressed how E2F family members are recruited to core promoter regions that lack a consensus motif and are excluded from other regions that contain a consensus motif. Using promoter and ENCODE arrays, we have shown that the predominant factors specifying whether E2F is recruited to an in vivo binding site are a) the site must be in a core promoter and b) the promoter region must be utilized as a promoter by the transcriptional machinery in that particular cell type. We have tested three models for recruitment of E2F to core promoters lacking a consensus site, including a) indirect recruitment, b) looping to the core promoter mediated by an E2F bound to a distal consensus motif, and c) assisted binding of E2F to a site that weakly resembles an E2F consensus motif within the core promoter. To test these models, we developed a new in vivo assay, termed eChIP, which allows analysis of transcription factor binding to isolated promoter fragments. Our findings suggest that in vivo a) the presence of a consensus motif is not sufficient to recruit E2Fs, b) E2Fs can bind to isolated regions that lack a consensus motif, and c) binding can require regions other than the best match to the E2F PWM in the core promoter. Keywords: E2F, ChIP-chip, transcription factor binding, consensus motifs 37 ChIP-chip arrays (of these, 14 array sets are biological duplicates). 22 samples are included in this series, the rest can be found in supplementary info to the following papers: Xu 2007, Jin 2006, Komashko 2008
2008-10-01 | E-GEOD-12126 | biostudies-arrayexpress
Project description:E2F in vivo binding specificity: comparison of consensus vs. non-consensus binding sites
| PRJNA113485 | ENA
Project description:Sequence-specific DNA-binding transcription factors are central to gene regulation. They are often associated with consensus binding sites that predict far more genomic sites than are bound in vivo. One explanation is that most sites are blocked by nucleosomes, such that only sites in nucleosome-depleted regulatory regions are bound. Alternatively, the consensus site may be poorly defined. We compared the binding of the yeast transcription factor Gcn4 in vivo using published ChIP-seq data (546 sites) and in vitro, using a modified SELEX method ("G-SELEX"), which utilizes short genomic DNA fragments to quantify binding at all sites. We confirm that Gcn4 binds strongly to an AP1-like sequence (TGACTCA) and weakly to half-sites. However, Gcn4 binds only some of the 1078 exact matches to this sequence, even in vitro. We show that the high-affinity site is RTGACTCAY (exact match), of which there are only 166 copies in the yeast genome, all occupied in vivo. Generally, RTGACTCAR/YTGACTCAY sites are bound much more weakly and YTGACTCAR sites are unbound, with biological implications for determining induction levels. We conclude that Gcn4 binding can be predicted using the genome sequence. Our study suggests that transcription factor consensus sequences should be defined more precisely using quantitative data.
2021-12-06 | GSE180114 | GEO
Project description:Recent genome-scale ChIP-chip studies of transcription factors have shown that a low percentage of experimentally determined binding sites contain the consensus motif for the immunoprecipitated factor. In most cases, differences between in vivo target sites that contain or lack a consensus motif have not been explored. We have previously shown that most sites to which E2F family members are bound in vivo do not contain E2F consensus motifs. The main purpose of this study was to develop an understanding of how E2F binding specificity is achieved in vivo. In particular, we have addressed how E2F family members are recruited to core promoter regions that lack a consensus motif and are excluded from other regions that contain a consensus motif. Using promoter and ENCODE arrays, we have shown that the predominant factors specifying whether E2F is recruited to an in vivo binding site are a) the site must be in a core promoter and b) the promoter region must be utilized as a promoter by the transcriptional machinery in that particular cell type. We have tested three models for recruitment of E2F to core promoters lacking a consensus site, including a) indirect recruitment, b) looping to the core promoter mediated by an E2F bound to a distal consensus motif, and c) assisted binding of E2F to a site that weakly resembles an E2F consensus motif within the core promoter. To test these models, we developed a new in vivo assay, termed eChIP, which allows analysis of transcription factor binding to isolated promoter fragments. Our findings suggest that in vivo a) the presence of a consensus motif is not sufficient to recruit E2Fs, b) E2Fs can bind to isolated regions that lack a consensus motif, and c) binding can require regions other than the best match to the E2F PWM in the core promoter. Keywords: E2F, ChIP-chip, transcription factor binding, consensus motifs
2008-10-01 | GSE12126 | GEO
Project description:The specificity of STAT1 transcription factor was determined quantitatively by measuring the relative binding affinity to hundreds of variants of the gonsensus site. The binding specities of eight mutants of STAT1 were also investigated. The mutants showed either overall reduced binding specificity or altered specificity accross the binding site. The effect of CpG methylaton within the binding site was also investigated, which demonstrated verry little effect on binding. A previously described method, Spec-Seq, was used to determine the binding specificititys of STAT1 and the mutants. Libreries of DNAs were enzymatically methylated and the same method was used for studying their binding affinities in comparison with the non-methylated counterparts. The proteins and DNAs were incubated for 30 min in in suitable buffer and the bound and unbound bands were resolved in 12% polyacrylamide gels. the DNAs were purified from the bound and unbound bands and sequenced.
2017-04-01 | GSE95526 | GEO
Project description:In this work, we used our recently developed method Spec-seq to characterize the binding specificity of Glucocorticoid receptor in vitro. This GR Spec-seq experiment has been run twice separately (Dec 2014 and Mar 2015) with different library compositions. The basic workflow is the same as our previous work for lac repressor published in Genetics 198.3 (2014): 1329-1343. Recombinant human GR protein was used to facilitate in vitro DNA-binding and separation experiments. Bound and Unbound DNA fragments were separated in EMSA gels, purified, barcoded for further Illumina sequencing. An initial experiment was performed using the putative consensus sequence: AGAACA GGG TGTTCT; randomized library 1: AGAACN NSN NGTTCT [diversity =512]; randomized library 2: AGAANN GGG NNNTCT [diversity = 1024]; randomized library 3: AGAACA GGG TGNNNN [diversity =256]; randomized library 4: AGAACA GGGC NNNTCT [diversity = 64]. The initial total library diversity was ~1853 with a library composition of 10% positive control sequence + randomized libraries 1-4 + 5% negative control sequence. A 5th library containing 2304 sequences based on DDNACW KKN KGTTCT, where D="not C", N="any base", W="A or T" and K="G or T" was subsequently prepared and analyzed using similar ratios of control sequences. Binding conditions for EMSA were 100ng FAM-labelled dsDNA+ 0/0.5/1/2/4uM GR DBD protein for each lane, 1X NEB buffer 4. GR DBD was prepared as previously described. The EMSA was performed using a 9% 33:1 acrylamide gel and TB buffer, and was run at 200V for 30mins @ 0 degrees C. The 2uM protein lane used for final sequencing. Bound/unbound fractions resulting from EMSA of these libraries and conditions were used to generate PWMs as described. The GR PWM that was generated through this analysis was used to define relative binding energies using the patser program (35), which can be accessed online at (http://stormo.wustl.edu/consensus/cgi-bin/Server/Interface/patser.cgi ). Derived binding affinities are proportional to the inverse of the natural log of the calculated energies.
2015-07-01 | E-GEOD-69386 | biostudies-arrayexpress
Project description:Eukaryotic cells express transcription factor (TF) paralogues that bind to nearly identical DNA sequences in vitro but bind at different genomic loci and perform different functions in vivo. Predicting how 2 paralogous TFs bind in vivo using DNA sequence alone is an important open problem. Here, we analyzed 2 yeast bHLH TFs, Cbf1p and Tye7p, which have highly similar binding preferences in vitro, yet bind at almost completely non-overlapping target loci in vivo. We dissected the determinants of specificity for these 2 proteins by making a number of chimeric TFs in which we swapped different domains of Cbf1p and Tye7p and determined the effects on in vivo binding and cellular function. From these experiments, we learned that the Cbf1p dimer achieves its specificity by binding cooperatively with other Cbf1p dimers bound nearby. In contrast, we found that Tye7p achieves its specificity by binding cooperatively with three other DNA-binding proteins, Gcr1p, Gcr2p, and Rap1p. Remarkably, most promoters (63%) that are bound by Tye7p do not contain a consensus Tye7p binding site. Using this information, we were able to build simple models to accurately discriminate bound and unbound genomic loci for both Cbf1p and Tye7p. We then successfully reprogrammed the human bHLH NPAS2 to bind Cbf1p in vivo targets and a Tye7p target intergenic region to be bound by Cbf1p. These results demonstrate that the genome-wide binding targets of paralogous TFs can be discriminated using sequence information, and provide lessons about TF specificity that can be applied across the phylogenetic tree.
2019-07-18 | GSE134427 | GEO
Project description:Spec-seq directly mesaures specificity by sequencing. We run an EMSA experiment on TF of interest with a library of DNA targets. Then we extract and sequence the bound and unbound portion of the DNA separately. Finally we count the occurences of each varient in the library and calculate the ratio and the relative binding energy.
2017-10-11 | GSE102219 | GEO
Project description:In this work we introduce a new high-throughput method, Spec-seq, that directly measures specificity by sequencing. It has several advantages over existing methods to quantify large collections (thousands) of binding site energies in one experiment. Using lac repressor as an example, we show that this method has excellent reproducibility giving energy measurements generally consistent within about 0.1kT. While not measuring in parallel as many different sequences as some of the higher throughput methods, we obtain high accuracies because we measure exactly what is necessary, the relative affinity to a large collection of sequences, without any mathematical fitting procedures or approximations required. It is similar to MITOMI in the number of different sites that can be analyzed in parallel, but it is much simpler to perform, requiring only a means to separate bound and unbound fractions which are then sequenced using standard high throughput, short read sequencing machines that are now readily available. When applied to the lac repressor we obtain, in a single experiment, data covering a large fraction of all the previous studies, plus thousands of additional variants, allowing us to compile a much more comprehensive profile of its specificity. We learn that the lac repressor can bind to sites of different lengths but that the preferred sequence, and the mode of binding, depends on the length. We also apply the method to two other members of the LacI/GalR protein family, PurR and YcjW, to obtain extensive models of their specificity and test the generality of the lac repressor's ability to bind to operators of variable length. We find that the lac repressor is apparently unique in its ability to bind with high affinity to sites of different lengths and with different modes of binding. 4 independent experiments.
2014-10-01 | E-GEOD-61223 | biostudies-arrayexpress
Project description:Motivation: The DNA binding specificity of a transcription factor (TF) is typically represented using a position weight matrix (PWM) model, which implicitly assumes that individual bases in a TF binding site contribute independently to the binding affinity, an assumption that does not always hold. For this reason, more complex models of binding specificity have been developed. However, these models have their own caveats: they typically have a large number of parameters, which makes them hard to learn and interpret. Results: We propose novel regression-based models of TF-DNA binding specificity, trained using high resolution in vitro data from custom protein binding microarray (PBM) experiments. Our PBMs are specifically designed to cover a large number of putative DNA binding sites for the TFs of interest (yeast TFs Cbf1 and Tye7, and human TFs c-Myc, Max, and Mad2) in their native genomic context. These high-throughput, quantitative data are well suited for training complex models that take into account not only independent contributions from individual bases, but also contributions from di- and trinucleotides at various positions within or near the binding sites. To ensure that our models remain interpretable, we use feature selection to identify a small number of sequence features that accurately predict TF-DNA binding specificity. To further illustrate the accuracy of our regression models, we show that even in the case of paralogous TF with highly similar PWMs, our new models can distinguish the specificities of individual factors. Thus, our work represents an important step towards better sequence-based models of individual TF-DNA binding specificity. Four protein binding microarray (PBM) experiments of human transcription factors were performed. Briefly, the PBMs involved binding GST-tagged transcription factors c-Myc, Max, and Mad2(Mxi1) to double-stranded 180K Agilent microarrays in order to determine their binding specificity for putative DNA binding sites in native genomic context. Briefly, we represent three categories of 36-bp sequences: 1) bound probes, 2) unbound probes (or negative controls), and 3) test probes. Bound probes corresponded to genomic regions bound in vivo by c-Myc, Max, or Mad2 (ChIP-seq P < 10^(-10) in HeLaS3 or K562 celld (ENCODE)) that contain at least two consecutive 8-mers with universal PBM E-score > 0.4 (Munteanu and Gordan, LNCS 2013). All putative binding sites occurr at the same position within the probes on the array. M-bM-^@M-^\UnboundM-bM-^@M-^] probes corresponded to genomic regions with ChIP-seq P < 10^(-10) and a maximum 8-mer E-score < 0.2. We also designed test probes that contain, within constant flanking regions, all nnCACGTGnn 10-mers and 18 nnnCACGTGnnn 12-mers (where n = A, C, G, or T). Each DNA sequence represented on the array is present in 6 replicate spots. We report the PBM signal intensity for each spot. The PBM protocol is described in Berger et al., Nature Biotechnology 2006 (PMID 16998473).
2013-05-17 | E-GEOD-47026 | biostudies-arrayexpress