Skip to main content
Download PDF

Sequence specificity incompletely defines the genome-wide occupancy of Myc

Genome Biology volume 15, Article number: 482 (2014) Cite this article

Abstract

Background

The Myc-Max heterodimer is a transcription factor that regulates expression of a large number of genes. Genome occupancy of Myc-Max is thought to be driven by Enhancer box (E-box) DNA elements, CACGTG or variants, to which the heterodimer binds in vitro.

Results

By analyzing ChIP-Seq datasets, we demonstrate that the positions occupied by Myc-Max across the human genome correlate with the RNA polymerase II, Pol II, transcription machinery significantly better than with E-boxes. Metagene analyses show that in promoter regions, Myc is uniformly positioned about 100 bp upstream of essentially all promoter proximal paused polymerases with Max about 15 bp upstream of Myc. We re-evaluate the DNA binding properties of full length Myc-Max proteins. Electrophoretic mobility shift assay results demonstrate Myc-Max heterodimers display significant sequence preference, but have high affinity for any DNA. Quantification of the relative affinities of Myc-Max for all possible 8-mers using universal protein-binding microarray assays shows that sequences surrounding core 6-mers significantly affect binding. Compared to the in vitro sequence preferences, Myc-Max genomic occupancy measured by ChIP-Seq is largely, although not completely, independent of sequence specificity.

Conclusions

We quantified the affinity of Myc-Max to all possible 8-mers and compared this with the sites of Myc binding across the human genome. Our results indicate that the genomic occupancy of Myc cannot be explained by its intrinsic DNA specificity and suggest that the transcription machinery and associated promoter accessibility play a predominant role in Myc recruitment.

Background

c-Myc was initially identified as a proto-oncoprotein and subsequently demonstrated to be a global regulator of transcription [1]-[6]. In vitro, the basic helix-loop-helix leucine zipper (bHLHZip) domain of Myc binds preferentially, albeit very weakly, to double stranded DNA containing the palindrome CACGTG and this is considered the canonical E-box [7],[8]. Myc pairs with Max and the heterodimer binds to CACGTG with higher affinity [9]-[11]. A crystal structure of the bHLHZip domains of Myc-Max bound to DNA revealed that the two proteins interact through each protein’s bHLHZip domain and each make specific contacts with four bases [12]. These initial observations and a large number of studies on the effects of Myc on specific genes led to what is now the prevailing model that Myc, in conjunction with Max, binds to E-box sequences and subsequently regulates transcription by Pol II [1]-[4]. However, this model does not explain how specificity for Myc is achieved as there are a number of other bHLHZip protein family members that can bind to the same sequence [13],[14].

Many studies attempting to identify Myc target genes found that the sets of genes regulated by Myc displayed great variation depending on the cell types and conditions used [15]. A significant advancement in understanding Myc function was achieved by two comprehensive studies from the Young and Levens labs [5],[6]. Both used inducible systems to show that Myc, when switched on, utilized existing expression programs and globally amplified transcription leading to an increase in the majority of expressed mRNAs. One study concluded that induction of Myc in P493 cells led to increased binding by Myc-Max heterodimers at the E-box containing core promoter sequences of actively transcribed genes [5]. The other study using primary B cells treated with lipopolysaccharide to induce Myc expression found a only a `loose association’ of Myc with E-boxes due to the high frequency of random occurrence of degenerate E-box sequences [6].

Regulation of gene expression is controlled predominantly through the action of DNA-binding transcription factors that affect both initiation and elongation. A prominent feature of metazoan genomes is the promoter proximal paused Pol II that is found engaged in transcription about 30 to 80 bp downstream of the transcription start site (TSS) on most expressed genes [16]-[19]. The transition of these paused polymerases into productive elongation requires the kinase activity of the Positive Transcription Elongation Factor b, P-TEFb [16],[20]. MYC was the first gene shown to be regulated by elongation [21] and Myc itself associates with P-TEFb [22]-[24] and causes an increase in productive elongation on targeted genes [25]-[27].

The mechanism of Myc regulation of transcription has been assumed to involve Myc-Max heterodimers binding to E-boxes near TSSs and then influencing the function of the transcription machinery. The discovery that Myc globally regulates essentially all expressed genes [5],[6] hints that Myc recruitment may be more general. An earlier study showed that Myc occupancy primarily correlated with chromosomal loci with an `open conformation’ [28],[29]. These regions are often occupied by the transcription machinery [19],[30]. To test the basic assumption that Myc-Max heterodimers predominately occupy high affinity DNA elements in cells, we performed detailed analyses of available human ChIP-Seq datasets for Myc, Max, and Pol II and determined the relationship between sites of occupancy of Myc and Max and locations of high affinity DNA elements. Surprisingly, the global occupancy of Myc and Max strongly correlated with Pol II transcription machinery rather than with sequences that the heterodimer prefers in vitro.

Results

Genome occupancy of Myc and Max correlates with Pol II

Using the UCSC Genome Browser [31] and ChIP-Seq datasets generated from HeLa cells [32] with antibodies stringently validated by the ENCODE project [33], occupancies of Myc and Max visually correlate with Pol II better than with the E-box element CACGTG. For example, a broad view of 10 genes across a 200 kb region shows almost identical patterns for Myc and Max and a high level of visual correlation with promoter proximal paused polymerases on each of the genes (Figure 1A). Many genes exhibit divergent transcription as indicated by GRO-Seq [34] that can result in paused Pol II in both orientations. A closer view of one such gene demonstrates that Myc and Max reside in a position between the two peaks of Pol II (Figure 1B). It is important to remember that the position of the immunoprecipitated factor is not indicated by the envelope of mapped DNA fragments, but rather by the peak of that envelope. Visual analysis of highly expressed genes, exemplified by MYC, provides further evidence that Myc and Max occupancy is tied to Pol II, including polymerases within the transcribed regions and downstream of the Poly(A) addition site (Figure 1C). For the three regions shown there is almost no correlation of Myc or Max with the canonical CACGTG E-box (Figure 1). In comparison, distributions of CTCF [35] and a number of other DNA-binding transcription factors (Additional file 1: Figure S1) are distinct from Myc, Max, and Pol II. When entire datasets were analyzed, genomic regions occupied by Myc exhibited a much more significant overlap with Pol II ChIP-Seq peaks than with the E-box element CACGTG (Fisher’s exact test: P value < 10-300 vs. 4.5 × 10-7).

Figure 1
figure 1

Examples of Pol II, Myc, and Max occupancy. Genome browser tracks show occupancy determined by ChIP-Seq for Pol II, Myc, Max, and CTCF over the indicated gene regions in HeLa cells. The positions of the canonical CACGTG E-boxes are indicated. Regions around (A) chromosome 19 containing 10 genes, (B) PSMB2, and (C) MYC are shown. GRO-Seq data are for IMR90 cells from GSE13518 [34].

Full size image

Several straightforward bioinformatic tools were used to obtain a global view of the correlation of Myc and Max compared to Pol II and CTCF. The average occupancy around the TSS of 20,886 genes in HeLa cells was calculated and plotted. Promoter proximal paused Pol II peaked on average 83 bp downstream of the TSS. Myc and Max on average peaked upstream of the TSS at -20 and -35, respectively (Figure 2A). Myc and Max also exhibited a slope transition at around +300 which has been previously noted for Pol II, the Med1 subunit of Mediator, and other transcription factors [36],[37]. High resolution heatmaps were generated to assess the uniformity of these distributions in the 4 kb region centered on the TSSs across the same gene set (Figure 2B). Genes were ranked by the amount of Pol II in all four heatmaps. The patterns for Myc and Max occupancy are essentially identical and they closely match the occupancy pattern for Pol II, but not CTCF. These results indicate that Myc and Max are found about 100 bp upstream of the promoter proximal paused Pol II on most of the genes occupied by Pol II. In addition, Myc and Max were also positioned very closely with Pol II in enhancer regions (Additional file 1: Figure S1C).

Figure 2
figure 2

Correlation of Myc and Max with Pol II occupancy. (A) Metagene analysis showing the average of 20,886 genes. (B) High resolution heatmaps of the same genes rank-ordered by Pol II occupancy. The region shown is from -2 kb to +2 kb around the TSS. (C) Correlation of the occupancy of the indicated proteins. (D) Metagene analyses of Myc, Max, and Pol II ChIP-Seq datasets from eight different cell lines (HeLa, GM12878, K562, H128, H2171, MM1S, P493, and U87). (E) Metagene analysis of Myc, Max, Med1, and Pol II ChIP-Seq datasets from four different cell lines (H2171, MM1S, P493, and U87). Average occupancies of regions from -1,000 to +1,000 bp around the TSS are shown.

Full size image

These ChIP-Seq datasets were also compared using an algorithm that measures the similarity of peak positions and heights in any two datasets (Figure 2C). A value of 0 means there is no overlap of the signals at any position and 1 indicates the datasets are identical. Myc and Max most closely correlate with each other, as expected. Importantly, the second highest genome-wide correlation for both Myc and Max was Pol II. The correlation of Myc with Pol II would not be expected to be as high as its correlation with Max because of the approximately 100 bp offset of Myc (and Max) from the peaks of promoter proximal paused Pol II. As expected, CTCF was the least well correlated with all datasets because it is bound by its CTC-containing motif mainly in intragenic regions [35]. The correlation analysis was extended to include Fos, Jun, and E2F1 and none of these factors correlated as well with Pol II as Myc and Max (Additional file 1: Figure S2).

We extended our analyses to eight human cell lines with Myc, Max, and Pol II ChIP-Seq datasets. All eight datasets were combined into a multi-genome metagene analysis and the results clearly indicated that on average, as was found in HeLa cells, Myc and Max were about 100 bp upstream of the promoter proximal paused Pol II and Myc is shifted downstream from Max (Figure 2D). Datasets for the Med1 subunit of Mediator were available for four of these cell lines and the multi-genome analysis displayed a similar distribution for Myc and Med1 including a downstream bulge over the promoter proximal paused Pol II (Figure 2E). These analyses strongly suggest that the Myc might be recruited to these genomic loci by the transcription machinery, with Mediator as a reasonable candidate.

Under stoichiometric conditions with high concentrations of proteins and DNA, Myc-Max heterodimers display relaxed sequence specificity

Because of the low correlation between Myc-Max genome occupancy and CACGTG sequences, we re-examined the DNA binding properties of the Myc and Max proteins. Full length versions of Myc and two isoforms of Max, MaxS and MaxL, were expressed in E. coli and purified to homogeneity (Figure 3A). The two Max isoforms were also individually mixed with Myc under denaturing conditions, allowed to refold using a step dialysis protocol, and then purified to obtain native, homogeneous heterodimers of Myc-MaxS and Myc-MaxL (Figure 3A). Electrophoretic mobility shift assays were carried out using three 26 bp dsDNA oligos that were identical except for the center 6 bps that contained the canonical CACGTG E-box, GTGGTG, or a completely unrelated sequence ATCTAG (Figure 3B). Native gels were silver stained to examine the shift in the position of 200 ng of protein. As expected, both homodimeric Max isoforms bound stoichiometrically to the CACGTG containing probe producing protein/DNA complexes that migrated further than the free proteins. MaxS displayed only very weak, transient binding to the other two probes while MaxL had reduced, but significant affinity for GTGGTG and low affinity for the ATCTAG probe (Figure 3B). Both Myc-Max complexes, regardless of Max isoform, produced a discrete protein DNA complex with the CACGTG probe. Surprisingly, both heterodimers bound stoichiometrically to the other two non-E-box probes (Figure 3B). Two individual studies assaying DNA binding with the same full-length proteins yielded identical shifting patterns [5],[38]. The differences in the relative levels of staining of free and DNA-bound forms of Max versus Myc-Max was caused by differences in the staining (development time) of the four representative gels shown. When MaxL and Myc-MaxL were analyzed on the same gel they displayed similar staining levels and comparable increases in staining when bound to DNA (Figure 3C). It is important to understand that these EMSAs (Figure 3B and C) were carried out under stoichiometric conditions with high concentrations proteins and DNA. These conditions do not allow the determination of dissociation constants and, especially for Myc-Max, do not display the sequence specific differences in binding that are known to exist. Instead they show that Myc-Max can bind to any DNA sequence at the high, but not unreasonable concentration tested (125 nM). The Myc-Max-DNA complexes showed only a small change in mobility comparing to the free proteins. This could be due to a change in conformation of Myc-Max that leads to a lowering of the mobility like that seen for HEXIM1 bound to 7SK RNA [39].

Figure 3
figure 3

Biochemical analysis of Myc and Max. (A) SDS-PAGE of the indicated recombinant proteins that were expressed in E. coli and purified as described in Methods. (B) EMSA using native polyacrylamide gel electrophoresis with 200 ng of the indicated proteins (250 nM Max dimer and 125 nM Myc-Max) with 0, 0.1, 0.3, 1, or 3-fold molar excess of the indicated dsDNA. The gels were silver stained to show the mobility of the proteins. The arrows indicate protein-DNA complexes. (C) EMSA with simultaneous staining of MaxL and Myc-MaxL. A total of 2.5 pmole of each protein (125 nM) per lane with two levels of the indicated DNA probes. Complexes containing indicated proteins are indicated with arrows. Note that in the Myc-Max prep some dissociation of Max has occurred leading to a low level of Max and Max-DNA species. (D, E, F, G, and H) EMSAs using 0.01 nM of the indicated radiolabeled probe (blue) with the indicated concentration of proteins and competitor DNAs.

Full size image

Dissociation constants of the protein-DNA complexes were determined under the required non-stoichiometric conditions using 0.01 nM radiolabeled probe. MaxL and Myc-MaxL displayed tight binding to CACGTG (Kds of 0.4 nM and 0.1 nM, respectively) (Figure 3D). MaxL did not form a discrete complex with the ATCTAG probe with the concentrations of protein tested (Kd >1 ?M), but instead gave only a smeary band below the position of a tightly bound complex (arrow) (Figure 3E). This is due to initial binding followed by release of the probe during the running of the gel. Myc-MaxL displayed significant affinity for the ATCTAG probe (Kd = 20 nM) (Figure 3F). Competition binding assays under these non-stoichiometric conditions demonstrated that CACGTG containing DNA was able to compete with the binding of MaxL and Myc-MaxL to the CACGTG probe (Figure 3G and H). At 1,000-fold higher concentration, the ATCTAG containing DNA was also able to compete for binding of both Max and Myc-Max to the CACGTG probe (Figure 3G and H). These results indicate that both Max and Myc-Max prefer to bind to the probe containing CACGTG as expected. In the stoichiometric assay described above, 125 nM Myc-Max but not 250 nM Max dimer formed discrete complexes with ATCTAG DNA. In the non-stoichiometric assay, Myc-Max displayed significantly higher affinity for the ATCTAG probe than Max and this difference was seen at 10 and 100 nM protein (Figure 3F). In the competition assay (1 nM protein) the difference between Myc-Max and Max was not seen. The concentration dependent change in the relative binding of Myc-Max and Max to non-specific DNA we observed could be related to the different on and off rates for the two proteins [40]. From all of the in vitro binding studies shown so far, we conclude that Myc-Max demonstrates a sequence preference, but that it also has significant affinity for DNA lacking a canonical E-box.

Determination of the complete sequence preference for Myc-Max and comparison with occupancy in cells

In our first attempts at trying to compare the in vivo occupancy of Myc and Max to the location of E-boxes, we ran into difficulty because of the existence of a large number of reported non-canonical E-boxes. Without quantification of the relative affinity of Myc-Max for all these sites it was difficult to correlate them with in vivo occupancy. Because of this, protein-binding microarray (PBM) assays using `all 10-mer’ universal array designs [41],[42] were used to quantify the relative occupancies of the Myc-MaxL heterodimer and the MaxL homodimer across all possible 8 bp sequences (that is, 8-mers). After normalization, relative Myc-Max occupancy for each of the 32,896 8-mers exhibited a 56-fold range, from 0.018 to 1 (Figure 4A, inset). Although the method is very different from the EMSA assay described above, the PBM results also reflect the relaxed sequence preferences of Myc-Max. Most of the sequences containing CACGTG had high occupancy, but flanking bases had a significant influence (Figure 4A). In addition, we found several E-box variants and other core 6-mers with relatively high Myc-Max occupancy. The top 12 core 6-mers and the effect of the flanking bases are shown in Figure 4A. Like the canonical CACGTG core, Myc-Max occupancy of the other core 6-mers was significantly affected by flanking bases.

Figure 4
figure 4

Binding of Myc to all possible 8-mers and comparison with genomic occupancy. (A) Fluorescent signal generated by Myc in vitro binding with an array containing all possible 8-mers was normalized. Twelve core 6-mer sequences with the highest in vitro occupancy are shown. The relative affinity of all 8-mers for each 6-mer is shown (10 points if the 6-mer is a palindrome or 16 if it is not). The inset shows the sorted in vitro binding signal for all possible 8-mers. (B) Genome browser view of a region on chromosome 19 comparing Myc, Max, and Pol II occupancy with the distribution of the top 12 6-mers (from A). The height of each 6-mer peak is equal to its relative in vitro occupancy (shown as percent). (C, D) Zoomed in views of two regions shown in (B) that demonstrate the lack of correlation of Myc and Max occupancy with the intrinsic affinity for the underlying DNA determined in vitro.

Full size image

The problem of not knowing the relative affinity of Myc-Max for the previously proposed non-canonical E-boxes was resolved by the PBM assays so we used that information to examine the role intrinsic DNA affinity plays in the occupancy of the heterodimer in cells. A genome browser track comprising the location and relative in vitro occupancy (percent of the top binding site) of each of the top 12 6-mers was generated that graphically displays the range of intrinsic affinities across the genome (Figure 4B). This is an improvement compared to just marking canonical and non-canonical E-boxes without regard to relative affinities of the different sites. Visual comparison of the occupancy of Myc, Max, and Pol II in HeLa cells to the accurate distribution of intrinsic affinities does not provide evidence for a strong correlation between intrinsic affinity and occupancy in cells (Figure 4B). Closer inspection revealed that strong binding sites were not occupied and Myc and Max were found in regions that did not have any of the top 12 6-mer sites (Figure 4C and D).

Several analyses were performed to compute the correlation between the 8-mer sequence preferences determined by PBM and the actual genomic occupancy of Myc, as measured by ChIP-Seq. The ChIP-Seq Peak algorithm [36] was used to determine the genomic location of each of the top 30,000 Myc peaks in HeLa cells. A 100 bp interval surrounding each peak was scanned to find the 8-mer with the highest possible in vitro occupancy and this score was assigned to each ChIP-Seq peak. These in vitro occupancy scores were normalized to 1, rank-ordered from highest to lowest values, and then plotted for all 30,000 peaks (Figure 5A, blue plot). Seventy-four percent of these Myc peaks were associated with low affinity 8-mers with in vitro occupancies below 0.2. To determine if the distribution of 8-mers around sites of Myc occupancy is different from what occurs by chance, the same analysis was performed on 30,000 100 bp regions randomly chosen from accessible DNA (DNase I sensitive regions [43]) (Figure 5A, black plot). The choice of DNase I sensitive regions as control sequences for this analysis is justified by the fact that 95% of the Myc peaks fall within such regions. Comparison of the two plots indicated that, as expected, genomic loci occupied by Myc contain more sites with high in vitro Myc occupancy compared to random accessible DNA regions (Wilcoxon rank-sum test: P value < 2.2 × 10-16). This enrichment is further shown by means of a receiver operating characteristic (ROC) curve (Figure 5A, inset). ROCs are commonly used in genomic analyses to assess whether a specific quantitative feature (here, in vitro Myc occupancy) can distinguish between two classes of sequences (here, ChIP-Seq peaks versus random accessible regions). Although the area under the ROC curve is better than expected by chance (0.637 vs. 0.5), the ROC analysis shows that the in vitro 8-mer occupancies cannot be used to accurately predict whether an accessible genomic region will be bound by Myc in cells. Here, the ROC plot shows that at a false positive rate of 0.1, the true positive rate is only 0.25. To make only 10% false positive predictions of Myc in vivo binding using the in vitro 8-mer scores, we would only be able to capture 25% of the true Myc ChIP-Seq peaks. This means that the vast majority of sites occupied by Myc are associated with low scoring 8-mers, as graphically indicated in Figure 5A.

Figure 5
figure 5

Comparison of Myc ChIP-Seq occupancy with in vitro binding affinities. (A) The top 30,000 sites occupied by Myc (blue) were rank-ordered and scored by the in vitro occupancy of the best 8-mer in a 100 bp window (y-axis). This was repeated at 30,000 random locations of DNase I-sensitivity (black) and the results were directly compared by ROC analysis (inset). (B) The top 30,000 sites occupied by Myc were rank-ordered by ChIP-Seq signal and scored logarithmically by either normalized ChIP-Seq signal (blue line) or the in vitro occupancy of the best 8-mer in a 100 bp window (black dots). (C) The data in (B) are presented using a default R boxplot (box: 1st to 3rd quartile, line: median, whiskers: 1.5 × interquartile range beyond the box, outliers are stacked) with ChIP-Seq signal in blue and in vitro 8-mers in grey.

Full size image

To further assess whether the intrinsic binding specificity of Myc-Max determines its level of genomic occupancy in the cell, the same Myc sites were rank-ordered by their ChIP-Seq occupancy and compared to the signal of the best 8-mer within a 100 bp window around each peak. The Myc ChIP-Seq signal of the top 30,000 peaks varies about 30-fold (Figure 5B, blue line showing decreasing occupancy from left to right). Using the same x-axis, a second plot was generated that displays the relative affinity of the best 8-mer associated with each of these Myc peaks (Figure 5B, black dots). A slight preference for high affinity 8-mers is visible over the top 5,000 Myc peaks, but the overwhelming conclusion is that 8-mers with a wide range of in vitro occupancies are found around Myc peaks irrespective of the level of in vivo occupancy (Figure 5B). While a statistically significant correlation can be observed between Myc ChIP-Seq occupancy and in vitro 8-mer binding strength, this relationship is weak (Spearman correlation coefficient: ρ = 0.22, P value < 2.2 × 10-16). Had the cellular occupancy correlated well with the affinity for the underlying DNA sequences, there would have been a cloud of black dots clustered around the blue curve in Figure 5B and the Spearman correlation coefficient would have been close to 1. A plot of the same data after ChIP-Seq peaks were grouped into log-scaled bins provides a more detailed view of the high occupancy sites in cells that might be expected to correlate better with intrinsic DNA affinities. However, the huge range of in vitro occupancy scores is clearly found even for the highest occupancy sites (Figure 5C). All these analyses suggest that Myc occupancy is driven only to a small extent by its intrinsic sequence preference, and additional mechanisms are required to recruit Myc to its genomic binding locations in the cell.

Genomic sites with higher relative levels of Max

Apart from associating with Myc, Max can form Max-Max homodimers or bind with Mad proteins to form Mad-Max heterodimers [44] and these can also bind E-box DNA sites [45]. We reasoned that such sites might have more Max than Myc. To identify these sites the HeLa Myc and Max datasets were normalized and a new track was generated in which the ChIP-Seq signal for Myc was subtracted from the signal for Max. Several thousand peaks with significant levels of extra Max were found. A representative region of chromosome 17, covering about 1 million bps that contains more than a dozen genes occupied by Pol II, Myc, and Max, is shown in Figure 6. The region contains about 20 peaks of Myc and Max and two of these sites have significant levels of extra Max. Both peaks of extra Max are on top of high scoring CACGTG sites (Figure 6B and C). Interestingly, the top 5,000 sites with extra Max (difference values greater than 0.5) were more tightly associated with high scoring 8-mers than were Myc sites (Additional file 1: Figure S3A) and had a more significant overlap with CACGTG than the Myc sites (Fisher’s exact test: P value <10-300 for extra Max sites vs. 4.5 × 10-7 for Myc sites). The top 1,487 peaks of extra Max (difference values greater than 1.0) were selected for further analysis (Additional file 2: Table S1). These sites were always close to peaks of Myc, Max, and Pol II, but only 417 of these peaks were within 250 bp of an annotated TSS. Gene Ontology (GO) analysis was performed on the associated genes, but no significant enrichment in any type of gene was uncovered. To determine if sites of extra Max might affect gene expression, the mRNA levels of those genes were compared to the mRNA levels of the top 12,000 expressed genes as determined by RNA-Seq. The RNA levels of 351 (of the 417) genes that were identifiable in the RNA-Seq dataset were distributed uniformly across the entire range of top 12,000 expressed genes covering more than three orders of magnitude in RNA levels (Additional file 1: Figure S3B). Thus, the sites with extra Max do not seem to be associated with any particular set of genes and do not correlate with the expression level of the genes they are associated with. Overall, sites with extra Max showed a stronger preference for E-box elements compared to Myc.

Figure 6
figure 6

Examples of sites with more Max than Myc. Genome browser views of normalized Myc, Max, and `Max minus Myc’ occupancy and peaks generated by ChIP-Seq Peak. The distribution of the top 12 6-mers with their relative in vitro occupancies is also displayed. (A) A large region from chromosome 17. (B, C) Close-ups of the two regions with extra Max showing alignment with high scoring 6-mers.

Full size image

Discussion

The results presented here provide evidence supporting a perspective for Myc function where the transcription machinery rather than DNA sequence elements plays a major role in recruiting the Myc-Max heterodimer to genomic sites. Although other studies have found Myc near TSSs [5],[6],[35],[46],[47], previous models for Myc function (Figure 7, top panel) evoked recruitment of the heterodimer to E-box sequences as an initial step (see recent reviews [1],[2]). However, we showed that sites of Myc occupancy were more highly correlated with Pol II rather than specific sequence elements. For sites near annotated TSSs, Myc was found about 100 bp upstream of the promoter proximal paused Pol II. We propose that for a large fraction of genes, the transcription machinery (which includes both initiation and elongation factors) is primarily responsible for recruitment of Myc-Max. The affinity of Myc-Max for DNA (Kd = 10-8 to 10-10 M, depending on the sequence) could then stabilize this interaction with specific sequence elements playing only a minor role (Figure 7, lower panel). The interaction(s) leading to Myc-Max occupancy is likely between the highly unstructured N-terminal transcription activation domain of Myc and factors in the transcription machinery. The model explains how Myc could influence entire transcription programs [5],[6]. One possible mechanism for Myc function that is consistent with what we now know is that Myc could be recruited by Mediator or another factor which is associated with the promoter proximal paused Pol II and then bind with Max relatively non-specifically to the promoter DNA. This could help keep the promoter region free of nucleosomes and primed for preinitiation complex formation if the paused polymerase was released into productive elongation or terminated. This mechanism could explain how Myc leads to universal amplification of gene expression in that it would generally increase accessibility of promoters and, therefore, responsiveness to the signals that regulate transcription of specific genes.

Figure 7
figure 7

Two models of Myc-Max recruitment. The top panel illustrates the prevailing view of E-box recruitment of Myc-Max followed by interactions of Mediator and other factors that affect transcription. The bottom panel is an alternative view supported by the results presented here in which Myc-Max is recruited by the transcription machinery and bound with a relaxed sequence requirement to promoter DNA. We suggest that the Myc-Max occupancy of the promoter region helps keep the promoter free of nucleosomes and the resulting increased accessibility of the promoter is responsible for the amplification of gene expression caused by Myc.

Full size image

Consideration of the size of the mammalian nucleus, the amount of Myc expressed, and the biochemical parameters of Myc-Max binding to DNA determined here leads to the interesting conclusion that almost all the Myc should be bound to DNA regardless of the influence of specific sequences. Quantification of Myc was recently reported for Myc-inducible P493-6 cells. It was found that before induction there were 13,000 molecules of Myc per cell and 77,000 and 362,000 molecules per cell after 1 or 24 h of induction. Taking into account the estimated volume of a mammalian cell nucleus of 4 × 10-13 liters, this means the concentration of Myc would start at 50 nM before induction and reach 1.5 ?M after 24 h induction. If free DNA was available, Myc-Max should be bound before and after induction since the Kd for even low affinity sites is 20 nM. An analysis of DNase I hypersensitivity, FAIRE, and occupancy of DNA binding transcription factors from ENCODE data in five human cell lines led to the estimation that about 3% of the genome is relatively free of histones and available for binding of factors [48]. This would be about 180,000,000 bp considering there are two genomes per cell. This means that even at the highest concentration of Myc, there are 500 bp of available DNA per Myc molecule and, therefore, all Myc should be bound to DNA. Our in vitro binding data demonstrates that the top 16 8-mers cover the top 65% of the normalized occupancy (see Figure 4). One of these 8-mers should occur on average only every 4,000 bp. Therefore, most of the Myc should be bound to lower affinity sites and this is borne out by our bioinformatic analyses. It has been estimated that there are about 180,000 Pol II molecules engaged in transcription in the average HeLa nucleus and it is likely significant that the level of Myc in a high expressing cell is similar, since Myc correlates highly with the position of engaged Pol II. Sites of occupancy of Pol II and Myc rise concomitantly during induction of Myc in P493-6 cells (Additional file 1: Figure S4), further connecting Myc and promoter proximal paused Pol II [5].

The specific interaction(s) that bridges Myc with the transcription machinery is not clear, but Mediator is found in a similar position to Myc upstream of the promoter proximal paused Pol II. Like Mediator, Myc and Max peaked upstream of the TSS. Both the Myc and Mediator peaks extended into the downstream region occupied by Pol II. This could be due to crosslinking between the N-terminal transcription activation domain of Myc with Mediator, which has been shown to interact with both Pol II [49],[50] and Myc [51], leading to indirect crosslinking through the polymerase to downstream DNA. Perhaps in an analogous manner, the N-terminal domain of KLF3 which interacts with other transcription factors but is not involved in DNA binding has been recently demonstrated to play a major role in genome occupancy [52]. Besides Mediator [51] and Pol II [53], Myc has been shown to associate with a number of factors including TRRAP, P-TEFb, chromatin remodeling machinery, DNA repair machinery, and other proteins involved in transcription and other processes [1],[2],[23],[26],[51],[53].

Our results have clarified the DNA binding preferences of the Myc-Max heterodimer, widening the range of core binding sites bound specifically by Myc-Max and suggesting that non-specific interactions could also be significant. In examining binding to all possible 8-mers it became apparent that sequences flanking the core 6-mers were important. Changing the flanking sequences of the canonical E-box 6-mer from GCACGTGC to TCACGTGA resulted in a four- or five-fold reduction in Myc-Max and Max binding. A similar reduction has been previously noted for Pho4p but not Cbf1p, two yeast bHLH proteins [54],[55]. Despite our expanded understanding of the relative interaction potentials of Myc-Max with all possible 8-mers, a strong correlation of Myc occupancy in cells to the underlying sequence could not be made.

How did it become so generally accepted that Myc-Max functioned as a specific DNA binding factor? As shown here and from earlier studies [9],[11]-[13], Myc-Max does have a preference for certain sequences in vitro. Previously, these preferences were determined using EMSAs that routinely included unlabeled carrier DNA. Unfortunately, this leads to an inaccurate determination of the actual range of affinities for different sequences. Many reporter assays have been performed where the presence of a CACGTG E-box sequence enhances the expression of the reporter (for example, [56],[57]). In almost all such experiments, only small (2-fold or less) effects are found and sequence-independent effects of Myc are normalized away. The small effects seen are consistent with our model in which having high affinity sites close to the promoter could fine tune the recruitment of Myc. Another issue arises from the misinterpretation of results obtained by programs, such as Multiple Em for Motif Elicitation (MEME) [58], which discover sequence motifs that are significantly enriched in a population of DNA sequences. However, identification of a motif (with a corresponding low P value) does not mean that sites matching the motif are present in all of the DNA sequences analyzed. In fact, very significant enrichment for a particular motif can be observed even when only a very small fraction of the analyzed sequences contain the motif. A second problem arises because of the length of the sequences analyzed and the resolution of actual binding site. Reasonable quality ChIP-Seq datasets like those analyzed here can determine individual binding sites to within 25 to 50 bps. If the fragments analyzed by a motif discovery algorithm are more than about 100 bp (as they usually are), a particular motif can be present in the sequence analyzed but not be bound (for example, see Figure 1B).

The model we propose for Myc might be applicable to other transcription factors. Significant overlap of occupancy between Myc and other factors such as AP-1 and AP-2 can be observed and two recent ENCODE reports showed co-association of many transcription factors binding to `surprisingly plastic’ DNA sequences [43],[59]. Indeed, our results suggest that it may be time to revisit the significance of sequence-specific binding for many transcription factors in higher eukaryotes. The E. coli lac repressor paradigm [60], in which occupancy of a specific site on a 4 million bp genome is driven by a seven order of magnitude difference between specific and non-specific binding, may not apply to certain human transcription factors. The generally repressive structure of chromatin in eukaryotes can mask most non-specific and specific sites thereby reducing the complexity of DNA that would otherwise be available for binding [61]. In addition, nucleosomes would also obstruct 1-D sliding, one of the important mechanisms used by lac repressor to locate specific sites [62]. The human genome encodes about a thousand DNA binding transcription factors with a wide range of sequence specificity. At one extreme, CTCF is almost always bound to one of several sequence motifs, as determined by differential involvement of its 11 Zinc fingers [63]. TBP on the other hand is found at all promoters regardless of the presence of the TATA sequence it recognizes due to interaction with the transcription machinery [64]. The metazoan genomic landscape may be too complex for all factors to rely solely on sequence specificity for occupancy and function.

Conclusions

We combined an extensive determination of the DNA binding properties of Myc-Max with genome-wide occupancies of Myc, Max, and Pol II and conclude that the affinity of Myc-Max for specific DNA sequences cannot be the main determinant of Myc genomic occupancy in cells. Our results indicate that the range in affinities in vitro for different sequences covers only a little over two orders of magnitude (Kd = 10-8 to 10-10 M). Even taking into account the fact that much of the human genome is not accessible due to chromatin structure, we found that Myc occupancy was not well correlated with affinity for underlying DNA in accessible regions. We found instead a strong correlation of Myc (and Max) occupancy with that of Pol II and on average Myc was located about 100 bp upstream of the promoter proximal paused Pol II. We propose that Myc is a general factor brought to promoters predominately by protein-protein interactions and like TBP, its recruitment to promoters does not strictly require sequence specific binding.

Materials and methods

ChIP-Seq datasets

The alignment files of HeLa, GM12878, and K562 cell lines were downloaded from human ENCODE Project at UCSC [65]. The raw sequence files of H128, H2171, MM1S, P493, and U87 cell lines were obtained from GSE36354. These raw sequences were aligned using ELAND to NCBI Build 36.1 (UCSC hg18) of the human genome. Only sequences that mapped uniquely to the genome with zero or one mismatch were used for further analysis. When multiple sequences mapped to the same genomic position, a maximum of two reads mapping to the same position were used. The sequenced reads were extended 200 bp to account for the size of sequenced fragments and then allocated into 25 base pair bins. The data from each bin were combined to generate the wiggle (WIG) files, which can be uploaded to UCSC genome browser. Max and Myc datasets from HeLa cells were further normalized for the total number of reads and the normalized Myc dataset was subtracted from normalized Max dataset to identify the genomic regions having extra Max compared to Myc.

RefSeq gene list

The complete set of human RefSeq genes was downloaded from the UCSC table browser [66] on 1 January 2012. A custom annotated set of 20,886 RefSeq genes was generated by merging the TSSs from the same gene within 500 bases of each other, and removing the all TSSs within 1,000 bases of each other. This custom list was used for all analyses.

Metagene analyses

The number of reads within 10,000 bases of the TSS of each RefSeq gene was tabulated without binning. The average value of the lowest 2,000 of 20,000 data points was subtracted from each position, and then data were normalized so that the area under each curve was equal. For Figure 2D, the data generated from the previous step were summed at each position using eight cell types from ENCODE (HeLa, GM12878, and K562) and GSE36354 (H128, H2171, MM1S, P493, and U87). Figure 2E was generated from GSE36354 (H2171, MM1S, P493, and U87).

Generation of heatmaps

Heat maps were generated using the program R [67]. Genes were rank-ordered based on the sequence density for Pol II from -2 kb to +2 kb from the TSS. Using this order, base pair resolution sequence density for Pol II, Myc, Max, and CTCF for 20,886 genes was displayed without binning. The raw images were 21,000 × 4,000 pixels each and were adjusted identically using the gamma adjustment in Corel PhotoPaint (version X3) to allow visualization of the wide range of data (Figure 2B).

Correlation of datasets

For Figure 2C, a total of 4 WIG tracks were generated using ENCODE data including Myc, Max, Pol II, and CTCF. The reads at each position were sorted in ascending order, and the cumulative percentage was calculated using R. The lowest 95% of the data points were subtracted as background since they represent mainly isolated single reads and the remaining reads were normalized to reads per million. The difference between two tracks was quantified as the sum of the absolute values of two tracks at each point. The correlation between two transcription factors was calculated within the range of 0 to 1, with a correlation of 1 indicating that the two tracks have a perfect match and 0 denoting total independence.

Peak finding

A peak finding algorithm (ChIP-Seq Peak) [36] was applied to determine precise position and height of each significant peak of Myc, Max, and Pol II. The height is equal to the number of reads contributing to the peak not the highest value at the summit of the peak. The peak positions with heights at least 30 (about 10% of highest peak height) were used to retrieve the sequences within 50 bases of peak locations. To generate Max minus Myc peaks, the Myc and Max datasets were first normalized so that the total reads were 1 million and then subtracted from each other before the ChIP-Seq Peak algorithm was applied. The peaks identified were further annotated with the RefSeq genes and RNA expression levels (GSE23316) across the genome. A total of 351 genes that have a peak of extra Max within 250 bp of their TSSs were classified into three Gene Ontology (GO) categories [68]. The normalized HeLa Myc and Max ChIP-Seq datasets, the Max minus Myc dataset, and its associated ChIP-Seq Peak dataset are available on GEO (GSE43227).

Expression and purification of proteins

Coding sequences of Myc (gi:29839758), MaxL (gi:21704261) and MaxS (gi:21704263) were amplified from HeLa cDNA and individually cloned into pET21a (C-terminal His-tag). All proteins were expressed in BL21 star E. coli after overnight induction of MaxL and MaxS at 18°C or 4 h induction of Myc at 37°C. Soluble Max proteins were purified over Ni-NTA and Mono Q for MaxS or Mono S for MaxL as described for Gdown1 [36]. Yields were approximately 20 mg per liter. Myc was purified over Ni-NTA and Mono Q in the presence of 7 M urea yielding 7 mg of Myc per liter. Both Myc-Max heterodimers were prepared by mixing equal moles of Myc and Max in 6 M urea followed by step dialysis against 4 M, 2 M, 1 M, and 0.5 M urea with 500 mM HGKEDP (25 mM HEPES, pH 7.6, 15% glycerol, indicated KCl, 0.1 mM EDTA, 1 mM DTT, and 0.1% of a solution of saturated PMSF in isopropanol) for 1.5 h each, then against 100 mM HGKEDP overnight. The renatured Myc-Max proteins were finally purified over Mono Q in HGKEDP. All protein samples were aliquotted and kept at -80°C. In the absence of urea, Myc was prone to aggregation, but Max and Myc-Max chromatographed cleanly and behaved well. Myc-Max heterodimers partially dissociated over time generating a small amount (<10%) of Max homodimers.

EMSA

Binding reactions were 25 mM HEPES (pH 7.6), 100 mM KCl, 1 mM DTT, 0.1 mM EDTA, and 0.01% Triton. For reactions analyzed on silver stained gels, 200 ng of protein and 0, 0.1, 0.3, 1, or 3-fold molar excess of unlabeled dsDNA were used in each reaction. For the reactions using radiolabeled probes, each strand of the indicated probe was end-labeled with ?-32P-ATP and then annealed to generate dsDNA. Each reaction contained 0.01 nM labeled dsDNA and 0.01 pM, 0.1 pM, 1 pM, 0.01 nM, 0.1 nM, 1 nM, 10 nM, or 100 nM protein. For competition assays 0.01 nM labeled CACGTG probe, 1 nM of the indicated protein and the indicated amounts of non-labeled competitor DNA were used. After 30 min at room temperature, Ficoll (2% final) was added to each reaction and the samples were run on a 6% polyacrylamide gel in 0.5X Tris/glycine buffer at 10 mA for 2 h. The gels were then analyzed by silver staining or phosphorimaging.

Protein-binding microarrays

Protein-binding microarray (PBM) experiments were performed as described previously [41]. Briefly, 4 x 44 K arrays (Agilent Technologies; AmadID 015681) containing the `all 10-mer’ universal PBM design were used. Arrays were incubated with a PBS buffer based protein mixture of 10 nM His-tagged Myc/Max heterodimer or Max homodimer, 2% milk, 200 ng/?L BSA, 50 ng/?L Salmon Testes DNA, and 0.02% TX-100. Bound protein was tagged with 10 ng/?L anti-His antibody conjugated to Alexa 488 (Qiagen; 35310) in PBS with 2% milk. The microarrays were scanned using a high-resolution GenePix 4400A scanner (Molecular Devices). Data were analyzed to obtain fluorescence intensities for all 8-mers, as described previously [41]. The raw median intensities and normalized in vitro occupancies for Myc-Max and Max-Max across all 8-mers are provided in Additional file 3: Table S2. Genome tracks were created by identifying all 8-mers in the human genome (NCBI build 36, hg18) that contain a core 6-mer that appeared in the top 200 8-mers (these 6-mers were CACGTG, CGCGTG, CACATG, CACGAG, CATGCG, AACGTG, CATGAG, CATATG, CATGCA, CATGCT, GTGCAC, and CTCGAG). Although only the central 6-mer was shown in the genome track, the in vitro occupancy scores were obtained by normalizing the observed intensity for each 8-mer to the intensity of the highest ranked 8-mer (in this case, CCACGTGG). The raw PBM data as well as the 8-mer intensities have been deposited in GEO (GSE58570) and the genome track for the top 12 6-mers has been deposited in GEO (GSE43227).

Statistical analyses

Statistical analyses were performed in R. Fisher’s exact test (dhyper) was used to compare the frequency of CACGTG sites in Myc-occupied regions and in all accessible DNA (DNase I hypersensitive regions [43]). A similar analysis was performed using Pol II ChIP-Seq peaks instead of CACGTG sequences. In Figure 5A, in vitro Myc 8-mer occupancies were compared using the Wilcoxon rank-sum test (wilcox.test with the paired parameter set to false). A receiver operating characteristic (ROC) curve was used to assess the enrichment of sites with high in vitro 8-mer occupancies in the Myc ChIP-Seq peaks. A ROC curve is a plot of false positive rate (1-specificity) versus true positive rate (sensitivity) computed for different 8-mer occupancy cutoffs. Each point on the ROC curve corresponds to one cutoff. Sequences with 8-mer occupancy above or below the cutoff are predicted positives (that is, predicted to be bound by Myc in vivo) or predicted negatives (that is, predicted not bound by Myc in vivo), respectively. In Figure 5B, the correlation between Myc ChIP-Seq and in vitro occupancies was tested using Spearman’s rank correlation coefficient (cor.test with the method parameter set to Spearman).

Authors’ contributions

JG and JJC cloned and expressed the Myc and Max proteins used. JG and FKF performed EMSAs. JS and RG performed the PBM assays and analyzed results with suggestions from DHP. TL, KAN, RG, and JS performed bioinformatic analyses with contributions from JG and DHP. DHP made the original observations that led to the study. JG, KAN, and DHP did most of the writing of the manuscript in consultation with all others. All authors approved the final manuscript.

Additional files

References

  1. Dang CV: MYC on the path to cancer. Cell. 2012, 149: 22-35. 10.1016/j.cell.2012.03.003.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Luscher B, Vervoorts J: Regulation of gene transcription by the oncoprotein MYC. Gene. 2012, 494: 145-160. 10.1016/j.gene.2011.12.027.

    Article  PubMed  Google Scholar 

  3. Eilers M, Eisenman RN: Myc’s broad reach. Genes Dev. 2008, 22: 2755-2766. 10.1101/gad.1712408.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Meyer N, Penn LZ: Reflecting on 25 years with MYC. Nat Rev Cancer. 2008, 8: 976-990. 10.1038/nrc2231.

    Article  PubMed  Google Scholar 

  5. Lin CY, Loven J, Rahl PB, Paranal RM, Burge CB, Bradner JE, Lee TI, Young RA: Transcriptional amplification in tumor cells with elevated c-Myc. Cell. 2012, 151: 56-67. 10.1016/j.cell.2012.08.026.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Nie Z, Hu G, Wei G, Cui K, Yamane A, Resch W, Wang R, Green DR, Tessarollo L, Casellas R, Zhao K, Levens D: c-Myc is a universal amplifier of expressed genes in lymphocytes and embryonic stem cells. Cell. 2012, 151: 68-79. 10.1016/j.cell.2012.08.033.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Blackwell TK, Kretzner L, Blackwood EM, Eisenman RN, Weintraub H: Sequence-specific DNA binding by the c-Myc protein. Science. 1990, 250: 1149-1151. 10.1126/science.2251503.

    Article  PubMed  Google Scholar 

  8. Halazonetis TD, Kandil AN: Determination of the c-MYC DNA-binding site. Proc Natl Acad Sci U S A. 1991, 88: 6162-6166. 10.1073/pnas.88.14.6162.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Blackwood EM, Eisenman RN: Max: a helix-loop-helix zipper protein that forms a sequence-specific DNA-binding complex with Myc. Science. 1991, 251: 1211-1217. 10.1126/science.2006410.

    Article  PubMed  Google Scholar 

  10. Blackwood EM, Luscher B, Eisenman RN: Myc and Max associate in vivo. Genes Dev. 1992, 6: 71-80. 10.1101/gad.6.1.71.

    Article  PubMed  Google Scholar 

  11. Jung KC, Rhee HS, Park CH, Yang CH: Determination of the dissociation constants for recombinant c-Myc, Max, and DNA complexes: the inhibitory effect of linoleic acid on the DNA-binding step. Biochem Biophys Res Commun. 2005, 334: 269-275. 10.1016/j.bbrc.2005.06.088.

    Article  PubMed  Google Scholar 

  12. Nair SK, Burley SK: X-ray structures of Myc-Max and Mad-Max recognizing DNA. Molecular bases of regulation by proto-oncogenic transcription factors. Cell. 2003, 112: 193-205. 10.1016/S0092-8674(02)01284-9.

    Article  PubMed  Google Scholar 

  13. Blackwell TK, Huang J, Ma A, Kretzner L, Alt FW, Eisenman RN, Weintraub H: Binding of myc proteins to canonical and noncanonical DNA sequences. Mol Cell Biol. 1993, 13: 5216-5224.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Jolma A, Yan J, Whitington T, Toivonen J, Nitta KR, Rastas P, Morgunova E, Enge M, Taipale M, Wei G, Palin K, Vaquerizas JM, Vincentelli R, Luscombe NM, Hughes TR, Lemaire P, Ukkonen E, Kivioja T, Taipale J: DNA-binding specificities of human transcription factors. Cell. 2013, 152: 327-339. 10.1016/j.cell.2012.12.009.

    Article  PubMed  Google Scholar 

  15. Littlewood TD, Kreuzaler P, Evan GI: All things to all people. Cell. 2012, 151: 11-13. 10.1016/j.cell.2012.09.006.

    Article  PubMed  Google Scholar 

  16. Zhou Q, Li T, Price DH: RNA polymerase II elongation control. Annu Rev Biochem. 2012, 81: 119-143. 10.1146/annurev-biochem-052610-095910.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Nechaev S, Adelman K: Pol II waiting in the starting gates: Regulating the transition from transcription initiation into productive elongation. Biochim Biophys Acta. 1809, 2011: 34-45.

    Google Scholar 

  18. Li J, Gilmour DS: Promoter proximal pausing and the control of gene expression. Curr Opin Genet Dev. 2011, 21: 231-235. 10.1016/j.gde.2011.01.010.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Guo J, Price DH: RNA polymerase II transcription elongation control. Chem Rev. 2013, 113: 8583-8603. 10.1021/cr400105n.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Peterlin BM, Price DH: Controlling the elongation phase of transcription with P-TEFb. Mol Cell. 2006, 23: 297-305. 10.1016/j.molcel.2006.06.014.

    Article  PubMed  Google Scholar 

  21. Bentley DL, Groudine M: A block to elongation is largely responsible for decreased transcription of c-myc in differentiated HL60 cells. Nature. 1986, 321: 702-706. 10.1038/321702a0.

    Article  PubMed  Google Scholar 

  22. Kanazawa S, Soucek L, Evan G, Okamoto T, Peterlin BM: c-Myc recruits P-TEFb for transcription, cellular proliferation and apoptosis. Oncogene. 2003, 22: 5707-5711. 10.1038/sj.onc.1206800.

    Article  PubMed  Google Scholar 

  23. Eberhardy SR, Farnham PJ: Myc recruits P-TEFb to mediate the final step in the transcriptional activation of the cad promoter. J Biol Chem. 2002, 277: 40156-40162. 10.1074/jbc.M207441200.

    Article  PubMed  Google Scholar 

  24. Eberhardy SR, Farnham PJ: c-Myc mediates activation of the cad promoter via a post-RNA polymerase II recruitment mechanism. J Biol Chem. 2001, 276: 48562-48571.

    PubMed  Google Scholar 

  25. Gargano B, Amente S, Majello B, Lania L: P-TEFb is a crucial co-factor for Myc transactivation. Cell Cycle. 2007, 6: 2031-2037. 10.4161/cc.6.16.4554.

    Article  PubMed  Google Scholar 

  26. Rahl PB, Lin CY, Seila AC, Flynn RA, McCuine S, Burge CB, Sharp PA, Young RA: c-Myc regulates transcriptional pause release. Cell. 2010, 141: 432-445. 10.1016/j.cell.2010.03.030.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Rahl PB, Young RA: MYC and transcription elongation. Cold Spring Harb Perspect Med. 2014, 4: a020990-10.1101/cshperspect.a020990.

    Article  PubMed  Google Scholar 

  28. Guccione E, Martinato F, Finocchiaro G, Luzi L, Tizzoni L, Dall’ Olio V, Zardo G, Nervi C, Bernard L, Amati B: Myc-binding-site recognition in the human genome is determined by chromatin context. Nat Cell Biol. 2006, 8: 764-770. 10.1038/ncb1434.

    Article  PubMed  Google Scholar 

  29. Sabo A, Amati B: Genome recognition by MYC. Cold Spring Harb Perspect Med. 2014, 4: a014191-10.1101/cshperspect.a014191.

    Article  PubMed  Google Scholar 

  30. Guertin MJ, Lis JT: Mechanisms by which transcription factors gain access to target sequence elements in chromatin. Curr Opin Genet Dev. 2013, 23: 116-123. 10.1016/j.gde.2012.11.008.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Kent WJ, Sugnet CW, Furey TS, Roskin KM, Pringle TH, Zahler AM, Haussler D: The human genome browser at UCSC. Genome Res. 2002, 12: 996-1006. 10.1101/gr.229102. Article published online before print in May 2002.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Consortium EP: A user’s guide to the encyclopedia of DNA elements (ENCODE). PLoS Biol. 2011, 9: e1001046-10.1371/journal.pbio.1001046.

    Article  Google Scholar 

  33. Landt SG, Marinov GK, Kundaje A, Kheradpour P, Pauli F, Batzoglou S, Bernstein BE, Bickel P, Brown JB, Cayting P, Chen Y, DeSalvo G, Epstein C, Fisher-Aylor KI, Euskirchen G, Gerstein M, Gertz J, Hartemink AJ, Hoffman MM, Iyer VR, Jung YL, Karmakar S, Kellis M, Kharchenko PV, Li Q, Liu T, Liu XS, Ma L, Milosavljevic A, Myers RM, et al: ChIP-seq guidelines and practices of the ENCODE and modENCODE consortia. Genome Res. 2012, 22: 1813-1831. 10.1101/gr.136184.111.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Core LJ, Waterfall JJ, Lis JT: Nascent RNA sequencing reveals widespread pausing and divergent initiation at human promoters. Science. 2008, 322: 1845-1848. 10.1126/science.1162228.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Lee BK, Bhinge AA, Battenhouse A, McDaniell RM, Liu Z, Song L, Ni Y, Birney E, Lieb JD, Furey TS, Crawford GE, Iyer VR: Cell-type specific and combinatorial usage of diverse transcription factors revealed by genome-wide binding studies in multiple human cells. Genome Res. 2012, 22: 9-24. 10.1101/gr.127597.111.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Cheng B, Li T, Rahl PB, Adamson TE, Loudas NB, Guo J, Varzavand K, Cooper JJ, Hu X, Gnatt A, Young RA, Price DH: Functional association of Gdown1 with RNA polymerase II poised on human genes. Mol Cell. 2012, 45: 38-50. 10.1016/j.molcel.2011.10.022.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Li T, Price D: Gdown1: making a link between mediator and RNA polymerase II elongation control. Transcription. 2012, 3: 177-180. 10.4161/trns.20600.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Farina A, Faiola F, Martinez E: Reconstitution of an E box-binding Myc:Max complex with recombinant full-length proteins expressed in Escherichia coli. Protein Expr Purif. 2004, 34: 215-222. 10.1016/j.pep.2003.11.021.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Sedore SC, Byers SA, Biglione S, Price JP, Maury WJ, Price DH: Manipulation of P-TEFb control machinery by HIV: recruitment of P-TEFb from the large form by Tat and binding of HEXIM1 to TAR. Nucleic Acids Res. 2007, 35: 4347-4358. 10.1093/nar/gkm443.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Ecevit O, Khan MA, Goss DJ: Kinetic analysis of the interaction of b/HLH/Z transcription factors Myc, Max, and Mad with cognate DNA. Biochemistry. 2010, 49: 2627-2635. 10.1021/bi901913a.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Berger MF, Bulyk ML: Universal protein-binding microarrays for the comprehensive characterization of the DNA-binding specificities of transcription factors. Nat Protoc. 2009, 4: 393-411. 10.1038/nprot.2008.195.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Berger MF, Philippakis AA, Qureshi AM, He FS, Estep PW, Bulyk ML: Compact, universal DNA microarrays to comprehensively determine transcription-factor binding site specificities. Nat Biotechnol. 2006, 24: 1429-1435. 10.1038/nbt1246.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Thurman RE, Rynes E, Humbert R, Vierstra J, Maurano MT, Haugen E, Sheffield NC, Stergachis AB, Wang H, Vernot B, Garg K, John S, Sandstrom R, Bates D, Boatman L, Canfield TK, Diegel M, Dunn D, Ebersol AK, Frum T, Giste E, Johnson AK, Johnson EM, Kutyavin T, Lajoie B, Lee BK, Lee K, London D, Lotakis D, Neph S, et al: The accessible chromatin landscape of the human genome. Nature. 2012, 489: 75-82. 10.1038/nature11232.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Grandori C, Cowley SM, James LP, Eisenman RN: The Myc/Max/Mad network and the transcriptional control of cell behavior. Annu Rev Cell Dev Biol. 2000, 16: 653-699. 10.1146/annurev.cellbio.16.1.653.

    Article  PubMed  Google Scholar 

  45. James L, Eisenman RN: Myc and Mad bHLHZ domains possess identical DNA-binding specificities but only partially overlapping functions in vivo. Proc Natl Acad Sci U S A. 2002, 99: 10429-10434. 10.1073/pnas.162369299.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Seitz V, Butzhammer P, Hirsch B, Hecht J, Gutgemann I, Ehlers A, Lenze D, Oker E, Sommerfeld A, von der Wall E, Konig C, Zinser C, Spang R, Hummel M: Deep sequencing of MYC DNA-binding sites in Burkitt lymphoma. PLoS One. 2011, 6: e26837-10.1371/journal.pone.0026837.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Kidder BL, Yang J, Palmer S: Stat3 and c-Myc genome-wide promoter occupancy in embryonic stem cells. PLoS One. 2008, 3: e3932-10.1371/journal.pone.0003932.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Yip KY, Cheng C, Bhardwaj N, Brown JB, Leng J, Kundaje A, Rozowsky J, Birney E, Bickel P, Snyder M, Gerstein M: Classification of human genomic regions based on experimentally determined binding sites of more than 100 transcription-related factors. Genome Biol. 2012, 13: R48-10.1186/gb-2012-13-9-r48.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Poss ZC, Ebmeier CC, Taatjes DJ: The Mediator complex and transcription regulation. Crit Rev Biochem Mol Biol. 2013, 48: 575-608. 10.3109/10409238.2013.840259.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Bernecky C, Grob P, Ebmeier CC, Nogales E, Taatjes DJ: Molecular architecture of the human Mediator-RNA polymerase II-TFIIF assembly. PLoS Biol. 2011, 9: e1000603-10.1371/journal.pbio.1000603.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Liu X, Vorontchikhina M, Wang YL, Faiola F, Martinez E: STAGA recruits Mediator to the MYC oncoprotein to stimulate transcription and cell proliferation. Mol Cell Biol. 2008, 28: 108-121. 10.1128/MCB.01402-07.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Burdach J, Funnell AP, Mak KS, Artuz CM, Wienert B, Lim WF, Tan LY, Pearson RC, Crossley M: Regions outside the DNA-binding domain are critical for proper in vivo specificity of an archetypal zinc finger transcription factor. Nucleic Acids Res. 2014, 42: 276-289. 10.1093/nar/gkt895.

    Article  PubMed  PubMed Central  Google Scholar 

  53. Koch HB, Zhang R, Verdoodt B, Bailey A, Zhang CD, Yates JR, Menssen A, Hermeking H: Large-scale identification of c-MYC-associated proteins using a combined TAP/MudPIT approach. Cell Cycle. 2007, 6: 205-217. 10.4161/cc.6.2.3742.

    Article  PubMed  Google Scholar 

  54. Fisher F, Goding CR: Single amino acid substitutions alter helix-loop-helix protein specificity for bases flanking the core CANNTG motif. EMBO J. 1992, 11: 4103-4109.

    PubMed  PubMed Central  Google Scholar 

  55. Maerkl SJ, Quake SR: A systems approach to measuring the binding energy landscapes of transcription factors. Science. 2007, 315: 233-237. 10.1126/science.1131007.

    Article  PubMed  Google Scholar 

  56. Kretzner L, Blackwood EM, Eisenman RN: Myc and Max proteins possess distinct transcriptional activities. Nature. 1992, 359: 426-429. 10.1038/359426a0.

    Article  PubMed  Google Scholar 

  57. Benvenisty N, Leder A, Kuo A, Leder P: An embryonically expressed gene is a target for c-Myc regulation via the c-Myc-binding sequence. Genes Dev. 1992, 6: 2513-2523. 10.1101/gad.6.12b.2513.

    Article  PubMed  Google Scholar 

  58. Bailey TL, Williams N, Misleh C, Li WW: MEME: discovering and analyzing DNA and protein sequence motifs. Nucleic Acids Res. 2006, 34: W369-W373. 10.1093/nar/gkl198.

    Article  PubMed  PubMed Central  Google Scholar 

  59. Consortium EP: An integrated encyclopedia of DNA elements in the human genome. Nature. 2012, 489: 57-74. 10.1038/nature11247.

    Article  Google Scholar 

  60. von Hippel PH, Revzin A, Gross CA, Wang AC: Non-specific DNA binding of genome regulating proteins as a biological control mechanism: I. The lac operon: equilibrium aspects. Proc Natl Acad Sci U S A. 1974, 71: 4808-4812. 10.1073/pnas.71.12.4808.

    Article  PubMed  Google Scholar 

  61. Kaplan T, Biggin MD: Quantitative models of the mechanisms that control genome-wide patterns of animal transcription factor binding. Methods Cell Biol. 2012, 110: 263-283. 10.1016/B978-0-12-388403-9.00011-4.

    Article  PubMed  Google Scholar 

  62. Halford SE, Marko JF: How do site-specific DNA-binding proteins find their targets?. Nucleic Acids Res. 2004, 32: 3040-3052. 10.1093/nar/gkh624.

    Article  PubMed  PubMed Central  Google Scholar 

  63. Nakahashi H, Kwon KR, Resch W, Vian L, Dose M, Stavreva D, Hakim O, Pruett N, Nelson S, Yamane A, Qian J, Dubois W, Welsh S, Phair RD, Pugh BF, Lobanenkov V, Hager GL, Casellas R: A genome-wide map of CTCF multivalency redefines the CTCF code. Cell Rep. 2013, 3: 1678-1689. 10.1016/j.celrep.2013.04.024.

    Article  PubMed  PubMed Central  Google Scholar 

  64. Burley SK, Roeder RG: Biochemistry and structural biology of transcription factor IID (TFIID). Annu Rev Biochem. 1996, 65: 769-799. 10.1146/annurev.bi.65.070196.004005.

    Article  PubMed  Google Scholar 

  65. Encyclopedia of DNA Elements at UCSC. [], [http://genome.ucsc.edu/ENCODE/]

  66. UCSC Table Browser. [], [http://genome.ucsc.edu/cgi-bin/hgTables]

  67. The R Project for Statistical Computing. [], [www.r-project.org]

  68. Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000, 25: 25-29. 10.1038/75556.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgments

We wish to thank Dan Weeks and Miles Pufall for thoughtful, spirited discussions. This research was supported by National Institutes of Health grant GM35500 to DHP, a PhRMA Foundation Research Starter Grant to RG, an American Heart Association postdoctoral fellowship 12POST12040106 to JG, a Duke Translational Medicine Quality Framework postdoctoral fellowship to JS, a University of Iowa Presidential Graduate Fellowship to KAN, and in part through computational resources provided by the University of Iowa.

Author information

Authors and Affiliations

  1. Department of Biochemistry, University of Iowa, Iowa City, 52242, IA, USA

    Jiannan Guo, Tiandao Li, Francis K Fordjour, Jeffrey J Cooper & David H Price

  2. Molecular and Cellular Biology Program, University of Iowa, Iowa City, 52242, IA, USA

    Kyle A Nilson & David H Price

  3. Institute for Genome Sciences and Policy, Duke University, Durham, 27708, NC, USA

    Joshua Schipper & Raluca Gordân

  4. The Genome Institute, Washington University in St. Louis, St. Louis, 63108, MO, USA

    Tiandao Li

Authors
  1. Jiannan Guo
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Tiandao Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Joshua Schipper
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Kyle A Nilson
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Francis K Fordjour
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Jeffrey J Cooper
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Raluca Gordân
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. David H Price
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to David H Price.

Additional information

Competing interests

The authors declare that they have no competing interests.

Electronic supplementary material

Additional file 1: Supplementary Figures S1-S4. (PDF 1 MB)

Additional file 2: Table S1.: Top 1,487 peaks of extra Max. (XLSX 151 KB)

13059_2014_482_MOESM3_ESM.xlsx

Additional file 3: Table S2.: Raw median intensities and normalized in vitro occupancies for Myc-Max and Max-Max across all 8-mers. (XLSX 2 MB)

Authors’ original submitted files for images

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Guo, J., Li, T., Schipper, J. et al. Sequence specificity incompletely defines the genome-wide occupancy of Myc. Genome Biol 15, 482 (2014). https://doi.org/10.1186/s13059-014-0482-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13059-014-0482-3

Keywords

Genome Biology

ISSN: 1474-760X

Contact us