ChIP on Chip: surprising results are often artifacts
© Waldminghaus and Skarstad; licensee BioMed Central Ltd. 2010
Received: 9 February 2010
Accepted: 5 July 2010
Published: 5 July 2010
The method of chromatin immunoprecipitation combined with microarrays (ChIP-Chip) is a powerful tool for genome-wide analysis of protein binding. However, a high background signal is a common phenomenon.
Reinvestigation of the chromatin immunoprecipitation procedure led us to discover four causes of high background: i) non-unique sequences, ii) incomplete reversion of crosslinks, iii) retention of protein in spin-columns and iv) insufficient RNase treatment. The chromatin immunoprecipitation method was modified and applied to analyze genome-wide binding of SeqA and σ32 in Escherichia coli.
False positive findings originating from these shortcomings of the method could explain surprising and contradictory findings in published ChIP-Chip studies. We present a modified chromatin immunoprecipitation method greatly reducing the background signal.
We set out to investigate the genome-wide binding of the sequestration protein SeqA in E. coli. This task can be considered especially challenging because SeqA has been shown to bind selectively to hemimethylated GATC sites . Although there are about 20.000 GATCs around the Escherichia coli chromosome only about 2% will be hemimethylated in unsynchronized cells . Such cell-to-cell variation increases the amount of cell material needed and therefore potentially the level of background signals. In fact, we found that application of a published ChIP-Chip method produced a background signal exceeding the specific signal. However, we were able to reduce the background significantly by modifying the protocol. The new protocol allowed us to uncover the genome-wide binding of SeqA and to reinvestigate σ32 binding to the E. coli chromosome.
High background signal in ChIP-Chip experiments
The background signal we observed seemed to correspond to highly transcribed regions, i.e. DNA with many RNA polymerase molecules bound (Fig. 2). Protein-rich DNA is segregated into the organic phase during phenol-chloroform extraction of crosslinked DNA . However, this phenomenon should not have affected a ChIP-Chip experiment, because the crosslinking is reversed before extraction is performed. The appearance of protein-rich gene regions as background might indicate an incomplete reversion of crosslinking at these sites. To clarify this question we compared DNA that was crosslinked and reversed with DNA that was not crosslinked. (Fig. 5B; see materials and methods for details). If the reversion of the crosslinking in this protocol is complete one would expect the two signals to be the same. This was indeed the case for the uvrD region. However, the rpsD DNA was more than seven-fold reduced in the crosslinked-reversed sample compared to the non-crosslinked DNA. To analyze the effect of crosslinking and reversion on a global scale we differentially labeled the DNA and applied it to a microarray. Ratios of the crosslinked-reversed versus the non-crosslinked DNA are shown in Fig. 5C (blue signal). The results show that the same regions that gave a high background signal in the SeqA ChIP-Chip yielded a reduced signal if the DNA is crosslinked and reversed (Fig. 5C; compare blue and red signal, Fig. 3D-E).
qPCR ratios of crosslinked-reversed versus non crosslinked DNAa
37°C proteinase K
0.36 ± 0.07
0.03 ± 0.01
42/65°C proteinase K
0.59 ± 0.25
0.03 ± 0.01
65°C o. n.
0.47 ± 0.12
0.03 ± 0.01
100°C 10 min
1.03 ± 0.44
0.04 ± 0.02
Modification of the ChIP-Chip procedure allows genome-wide analysis of SeqA binding
Considering the identified weaknesses of the ChIP-Chip protocol it was possible to make appropriate modifications (see material and methods for details). The first change was the omission of spin-columns in the washing of agarose beads. Second, the input DNA was taken from the supernatant resulting from centrifugation of the immunoprecipitated chromatin beads. In addition, we included RNase digestion of immunoprecipitated DNA and excluded signals originating from microarray-probes to non-unique sequences during data analysis. The reasoning behind the latter two will be described in detail below.
To put the results in a biological context we calculated the SeqA binding signal for a 60.000 bp moving window (Fig. 7, inner ring). The reasoning behind this is that SeqA has been shown to bind specifically to hemimethylated DNA "trailing" the replication fork. We estimated the stretch of hemimethylated DNA following the replication fork to be 60.000 bp (based on a replication speed of 1000 bp/sec and an average hemimethylation time of 1 min). The result shows that SeqA binding is not evenly distributed over the chromosome. Instead there are regions with strong binding, such as the origin of replication (oriC) and areas with low binding, such as to the left and right of oriC (Fig. 7). The most extended area with low SeqA binding is about one-fourth of the chromosome around the replication terminus with distinct borders rather than smooth transitions to the neighboring high SeqA binding regions. A clear correlation was observed between the number of GATC sites in the probe region and the corresponding ChIP signal (Fig. 7B). In summary, we have shown that the revised ChIP-Chip protocol can be successfully used to gain insight into the challenging question of chromosome-wide SeqA binding in E. coli.
Reinvestigation of σ32 binding to the E. coli genome
Target detection for previously reported σ32-sitesa
Target detection for σ32-sites found only by Wade et al., 2006a
σ32 target candidates not detected in Wade et al.
Limited RNase treatment is an additional source of false positives in ChIP-Chip studies
Non-unique sequences can cause false positives in ChIP-Chip analysis
One important source of false positive findings in ChIP-Chip studies is the inclusion of non-unique sequences. For the 40.000 probes on the microarray used in this study we examined the number of complementary sequences on the E. coli chromosome. 889 probes were found to match multiple loci on the chromosome, the numbers ranging from 2 to 11 (data not shown). Note that signals obtained with these probes and the surrounding probes were routinely excluded from all results shown above as mentioned. However, to investigate the effect of these non-unique probes we reanalyzed the σ32-ChIP-Chip experiment of 30°C cells described above including the non-unique sequences (Fig. 8C). Some of these probes gave an elevated ChIP signal. Consequently, six new peaks were detected by our search algorithm in addition to the 15 peaks detected before (Fig. 8C). Also the published σ32-study includes two target sites in non-unique sequence regions. These are the yibA promoter close to the rhsA gene and the yrdA promoter downstream of the ribosomal RNA gene rrsD. In summary, our data demonstrate the potential of non-unique sequences to cause false-positive findings in ChIP-Chip studies.
Multiple sources of false positives in ChIP-Chip studies
Here we present four sources of high background signals that caused false positive target site detection in our experiments as well as in many published studies. In the following we discuss how this unspecific background might occur. The first two problems, namely the selective enrichment of some DNA fragments during spin-column washing and the variability in reversion of crosslinking, might actually be due to the same circumstance. Both affected chromosomal regions with high transcription activity, such as the ribosomal protein gene rpsD (Fig. 5). In such regions crosslinking of RNA polymerase, DNA and transcribed mRNA will form large complexes. Concerning the washing of immunoprecipitated DNA with spin-columns it is easy to imagine that such highly crosslinked fragments could be trapped in the column matrix. A release of these bound complexes in the elution step would explain the enrichment of protein-rich DNA through washing with spin-columns. This would be limited to the IP DNA in a ChIP experiment because usually no beads are used to purify the input DNA. The logical improvement of the protocol in this case was to wash the immunoprecipitated DNA without spin-columns. Another possibility would be to use systems which separate beads by magnetism instead of centrifugation.
In contrast, the difference in crosslinking/reversion efficiencies at genomic loci could not be reduced by leaving out the crosslinking because it is an essential part of the protocol. The incomplete reversion of crosslinking led to depletion of protein-rich chromosomal regions during DNA preparation (Fig. 5). If this depletion were similar in the IP and input DNA it would not appear as ChIP signal because the corresponding ratio would be one. However, different rates of depletion in IP and input DNA would let this ratio go up or down. If for example 60% of a crosslinked site is reversed in the IP DNA but only 30% in the input DNA this would appear as two-fold enrichment and potentially as false positive target. Thus, transcriptionally active regions of the chromosome are more likely to show a high background signal. This problem could not be solved by variation of reversion conditions (Table 1). However, as one way to better separate the real targets from such background we increased the specific signal by using the supernatant of the immuno-precipitation as input DNA. This should amplify the specific signal because it will be enriched in the immuno-precipitated DNA and at the same time reduced in the reference DNA.
A high background signal originating from non-digested RNA may also occur in ChIP-Chip experiments. This will for example be high if the Klenow fragment is used for labelling of immunoprecipitated DNA, since it can use RNA as primer to incorporate labelled nucleotides. If a linker-mediated PCR is used to amplify the immunoprecipitated DNA the amount of RNA relative to DNA will be reduced, potentially reducing the RNA-caused background. Here we show that a thorough RNase digestion is a suitable way to eliminate the RNA background, allowing a free choice of subsequent labelling and amplification techniques.
An additional origin of high background signals in microarray analysis is caused by the occurrence of non-unique sequences on the chromosomes. A systematic evaluation of labeling and microarray hybridization of predefined DNA targets revealed such genome redundancy as one major cause of false positives . A probe to a non-unique sequence will bind a mix of DNA fragments originating from different chromosomal loci. The chromosomal position can influence the protein binding to the different copies of a non-unique sequence and may therefore lead to erroneous ChIP-Chip results. If for example one copy is located downstream an active promoter and the other copy not, a RNA polymerase ChIP would enrich the first locus but not the second. On the microarray this would appear as a medium enrichment at both chromosomal positions. Additional errors might occur at non-unique sequences with multiple copies and some sequence variation. In this case one probe might be complementary to for example two copies and the neighboring one to seven copies. Genes that are typically non-unique are the ribosomal and transfer RNA genes or transposons but also for example the rhsABCD gene family or gadAB in E. coli.
To estimate the degree of false positives caused by non-unique sequences we screened the literature for occurrence of the mentioned genes as target sites in microarray studies. Appearance of non-unique sequence false positives turned out to be quite frequent. For example, 36 out of 269 'extended protein occupancy domains' in a recent study from Vora et al. are in regions with non-unique sequences . Some studies even draw major conclusions from the appearance of non-unique sequence false positives. For example, the heat shock regulator HspR was suggested to be involved in regulation of tRNA and rRNA genes in Streptomyces coelicolor, the B. subtilis condensin SMC was proposed to be recruited to rRNA and tRNA genes  and tRNA genes were described to be cohesin loading sites both in budding and fission yeast [16, 17]. All of the mentioned gene loci are non-unique in the respective genomes. Note that in principle the described conclusions could be right; it is just that the results of microarray experiments can say nothing about it and might actually be misleading instead. Fortunately, non-unique sequences can be easily detected and corresponding probes be excluded from data sets. Even better would be elimination during array design.
Beside the causes of high background described in this study other factors have been shown to affect the background level. For example Lee and colleagues point out that ChIP-Chip experiments are highly dependent on the antibody used for the immunoprecipitation . The background signal will be high if the antibody performs poorly or if it binds other proteins unspecifically. In this context the salt concentration of the IP and wash buffer is critical and can be adjusted to optimize immunoprecipitation . In addition to the experimental procedure improper data processing can lead to false positive findings. How the data are analyzed will depend on different factors such as probe density and the relative number of binding sites . Correct normalization regarding the dye bias in two color microarrays has been shown to be essential for ChIP-Chip experiments .
How frequent are false positives in published ChIP-Chip data?
The presence of non-unique sequence false positives might indicate that a high number of false positives are the rule, rather than exception in published ChIP-Chip studies. A false positive rate about 50% was found by our reinvestigation of a published σ32-study . The conclusion that the targets found in the published ChIP-Chip experiment but not in our study are false positives is supported by findings from others [11, 19]. While almost all of the targets we detected have been found with other methods then ChIP-Chip, the only evidence for the supposed new targets by Wade et al. is their ChIP-Chip analysis [10, 11, 19]. It is noteworthy that this analysis was done with the protocol used in the first experiment of our study producing a high background . In addition the supposed new targets lacked a typical σ32-recognition site .
Further evidence for a frequent false-positive rate in ChIP-Chip studies comes from large differences of binding site detections in parallel studies. For example, FIS was found to bind all regions on the E. coli genome that are bound by RNA polymerase despite the absence of consensus binding sites . A later study showed very different results with data that nicely fit the distribution of FIS binding motifs . In two independent studies the binding of the estrogen receptor to the human chromosome 17 of MCF-7 breast cancer cells was analyzed [22, 23]. We compared the 389 binding sites described in the Gevry study to the 390 sites detected in the Carroll study and found only about 50% overlap (binding sites were considered the same when not more than 2000 bp apart, data not shown). Interestingly, others have also suggested an extended degree of false positives as explanation for contradictory results in parallel ChIP-Chip studies. Highly dissimilar binding patterns of the Mediator complex in yeast were reported [24–26]. Fan and Struhl reinvestigated the contradictory results and suggested that the differences were caused by a high degree of false positives due to the experimental set-up of Andrau and colleagues . These supposed false positives are mainly located in transcriptionally active coding regions as is also the case in our study.
A high number of false positives would make systematic approaches to analyze ChIP-Chip-derived binding sites especially difficult. Indeed, a recent analysis of yeast ChIP-Chip data revealed that only 48% of detected transcription factor binding sites could be explained by direct binding and an additional 16% by indirect binding . The remaining 36% of the data set could not be explained by either direct or indirect transcription factor binding and were suggested to be noise. Taken together, high false positive rates seem to be common in ChIP-Chip studies. In some cases it actually seems to be an accepted fact. For example, Partridge and colleagues removed over one third of ChIP-Chip detected NsrR target sites just because they did not fit their expectations of lying in promoter regions . However, this high false-positive rate was not investigated any further.
How to deal with the background
Beside the need for technical improvements, the high level of ChIP-Chip false positives emphasizes the great importance of suitable control experiments. Good controls are ChIP-Chip experiments with cells lacking the IP epitope (for instance ΔseqA; Fig. 6), mock IPs without antibody (Fig. 5A) or IPs with preimmune serum or IPs from cells growing under conditions that are expected to give no or reduced binding of the respective protein (such as 30° for the heat shock sigma factor σ32; Fig. 8C). A suitable control experiment has two important functions. First, it allows estimation of the experimental quality. In this study the ΔseqA control was the key to understanding that the ChIP-Chip method gave high background (Fig. 4 and 6). Second, a control experiment can help to detect targets in the actual experiment. We used the σ32 control ChIP-Chip at 30°C to find significant targets in the corresponding data set of heat shocked cells (see materials and methods).
It has been suggested that DNA from control experiments should be used as a hybridization reference, meaning that for example the IP DNA from a wt strain and a deletion strain are differentially labeled and hybridized to the same array . However, others point out that a control should never be used as hybridization reference . We agree with the latter opinion because use of control DNA as hybridization reference would not allow assessment of the experimental quality as outlined above. For instance, bad quality DNA from experiments with limited digestion of RNA (Fig. 8A-B) might not be detected if used as hybridization reference. Taken together, appropriate control experiments should be included in every ChIP-Chip study. Submission of the raw and processed control data to the public should be self-evident but is an exception in published studies so far.
Recently, chromatin immunoprecipitation has been combined with high throughput sequencing methods (ChIP-Seq). Interestingly, an analysis of different types of control DNA resulted in a variable pattern of background distributed over the chromosomes [31, 32]. The pattern of background peaks varied between input DNA, non-crosslinked DNA and mock-IP DNA and lead to the conclusion that the type of reference DNA directly influence the number of sites deemed significant when scoring ChIP-Seq data. This underlines that the described problems apply to chromatin immunoprecipitation based methods in general.
Revised ChIP-Chip method reveals new biological insights
The revised ChIP-Chip method we developed enabled us to analyze binding of the sequestration protein SeqA to the E. coli chromosome. SeqA is involved in regulation of replication initiation and also proposed to play a role in chromosome organization and segregation . It was found to exhibit prolonged binding to hemimethylated GATC sites at oriC and thereby hindering reinitiation [7, 33]. Enhanced binding of SeqA at oriC was also found in our ChIP-Chip analysis, in fact it was the highest peak detected (Fig. 7). The second-highest peak was in the dnaA promoter region which has been shown to have an exceptionally long hemimethylation period . While our data support SeqA binding as proposed for oriC and the dnaA promoter it contradicts published suggestions on chromosome-wide binding. Brendler and colleagues found an even distribution of potential SeqA binding sites over the chromosome . Our data suggest that SeqA structures retain specific DNA tracts for varying amounts of time. Most striking is the relatively short duration of SeqA binding to the left and right of oriC and to the DNA at about one-quarter of the chromosome surrounding the replication terminus. The latter finding is in contrast to results from ChIP-PCR experiments with synchronized cells which suggested a prolonged SeqA binding in the terminus region . Clearly, further analysis and additional experiments are needed to understand the biological meaning of the SeqA binding pattern.
We describe here a revised ChIP-Chip method and show its potential to greatly reduce false positive target site detection, which seems to be a widespread problem. Although we present many examples of high false positive rates in published studies, it has to be pointed out that this will vary greatly with the exact experimental details as outlined above. Since method details such as the duration of the RNase treatment or the use of spin columns have a major impact on the background signal, it is of high importance t give an accurate description of the procedure used. The results reported here should allow critical reviewing of published ChIP-Chip studies as well as assessment and potential modification of other variants of the ChIP-Chip method and related methods.
Cell growth, crosslinking and preparation of cell extracts
Strains and oligonucleotides used in this study
Strain or oligonucleotide
Relevant characteristic(s) or sequence
Source or reference
E. coli MG1655
F-λ rph-1 (wild type)
Guyer et al., 1981; Jensen, 1993
Torheim et al., 2000
AGACGCCCGCCTTCATCCAG (5' FAM - 3' TAMRA)
AACTCGCGTCCTTCGGCAGC (5' FAM - 3' TAMRA)
TCAGAACGCTCCGGCTTACGC (5' FAM - 3' TAMRA)
The ChIP protocol initially used in this study was as described in Grainger et al., 2004 except that DNA was purified with phenol/chloroform instead of a PCR clean up kit. 800 μl of sonicated cell extract (see above) was incubated with 20 μl protein A/G agarose beads (Ultralink) and 5 μl of SeqA antiserum or antibody against RNA polymerase subunit β (Neoclone) at 4°C over night. Samples were transferred to a Spin-X centrifuge column (Costar), centrifuged for 2 min at 4.000 rpm to collect the beads. The flow through was removed. Washing was done by adding 500 μl buffer to the beads on the spin column and rotation at room temperature for three minutes with subsequent collection of the beads by centrifugation as above. Washing was performed with the following buffers (IP buffer two times all others one time): IP buffer (50 mM HEPES-KOH pH 7.5, 150 mM NaCl, 1 mM EDTA, 1% Triton × 100, 0.1% Sodium deoxycholate, 0.1% SDS), IP buffer with 500 mM NaCl, wash buffer (10 mM Tris pH 8.0, 250 mM LiCl, 1 mM EDTA, 0.5% Nonidet-P40, 0.5% Sodium deoxycholate) and TE. For elution, 100 μl elution buffer (50 mM Tris (pH 7.5), 10 mM EDTA, 1% SDS) was added to the column with the beads, incubated in a 65°C water bath for 10 min and centrifuged as above. To reverse the cross link 80 μl TE and 20 μl proteinase K (20 mg/ml) were added and samples incubated for 2 h at 42 and 6 h at 65°C. DNA was purified with phenol/chloroform. To prepare the control DNA, 800 μl of sonicated cell extract was incubated at 65°C over night. 1 μl RNase A (20 mg/ml) were added and samples incubated 30 min at 65°C before extraction with phenol/chloroform. The ChIP protocol as described above resulted in the high background signal (Fig. 2 and 4).
The following modifications were applied for the other ChIP-Chip experiments. First, agarose beads were not collected on a spin column but instead at the bottom of a usual 1.5 ml eppendorf tube. The supernatant was then removed by pipetting. Second, the control DNA was taken from the supernatant resulting from centrifugation of the precipitated chromatin beads processed further as the immuno precipitated DNA after elution. Third, before addition of proteinase K, sample and control DNA were incubated with RNase A (50 μg/ml) for at least 90 min at 42°C (except in the σ32-analysis shown in Fig. 8A where incubation was 30 min as indicated). Incubation of 800 μl cell extract with 15 μl σ32- or 5 μl SeqA antiserum was for 1 h at 4°C.
Labeling and array hybridisation
Usually DNA from six parallel immuno-precipitations (each with 800 μl extract as described) were joined and labeled with Cy3-dCTP using the Klenow fragment and random primers of the BioPrime kit from Invitrogen. An equal amount of hybridization control DNA was labeled with Cy5-dCTP. Hybridization was for about 36 h at 55°C to E. coli whole genome microarrays from Oxford Gene Technology. The arrays have a probe length of 60 bases and a start to start spacing of about 150 bases. ChIP-Chip analysis were made in duplicates, except the crosslink-reversion array (Fig. 5), the ΔseqA arrays (Fig. 6) and the shorter RNase incubation array (Fig. 8A). Please note that the array hybridized with the SeqA ChIP of the ΔseqA strain with the unmodified method was of poor quality but regarded sufficient for its purpose described above.
Microarray data processing
Arrays were scanned on an Agilent SureScan High-Resolution Scanner. Spot intensities were extracted using the Feature Extraction software 10.5.1.1 from Applied Biosystems with a linear dye normalization correction method. The data were further analyzed with the statistics software R, in particular the Bioconductor package and the limma library [36, 37]. The background was subtracted and data points with a value below 0 after background subtraction were removed. Ratios of g (sample) to r (control) were calculated and normalized to the array wide average. For arrays performed in duplicates the mean of the two normalized values was calculated. Probes in gene regions with non unique sequences were deleted (a list is available on request). For σ32-target detection data obtained from heat-shocked cells were searched for two or more neighboring probes with a log2 signal > 0.5 in both replicates. This resulted in 74 potential targets (34 previously described, 9 described exclusively by Wade et al., 2006, 31 not found by Wade et al.). After subtraction of log2 signals of the corresponding replicates from non-heat-shocked cells, 47 potential targets remained (Tables 2, 3, 4; 34 previously described, 7 described exclusively by Wade et al., 2006, 6 not found by Wade et al.). For peak detection in σ32-data of non-heat-shocked cells (Fig. 8C) we searched for probes with a log2 ratio > 1 and the one to the left and right > 0.5.
GenomeViz was used for visualization of ChiP-Chip data . Data points with log2 ratios > 0.5 were extracted and the corresponding genome locus assigned as 1000 bp up- and down-stream of the respective probe middle. For the moving window calculation of SeqA binding the sum of positive log2 ratios of 60.000 bp windows were calculated with a step size of 1000 bps. Raw as well as processed data are available at the Genome Omnibus Database, accession number GSE19053. To analyze the overlap of ChIP-Chip experiments a cut-off was chosen for each data set to select ~1000 probes with the highest ChIP signal (or the lowest signal for the crosslinking experiment). The overlap is the number of probes were the signal is beyond this cut-off at similar positions in the two compared data sets.
ChIP washing comparison
For the comparison of washing methods (Fig. 5A) 2 × 800 μl of crosslinked, sonicated MG1655 cell extract were incubated with 20 μl protein A/G agarose beads (Ultralink) without antibody for 1 h at 4°C. One of these mock IP samples was then processed with the use of spin-columns and one without as described above. Eluted DNA was purified with phenol/chloroform and analysed by quantitative PCR as described below. Note that purification of the DNA with a Qiagen PCR cleanup kit gave the same results as the phenol extraction (data not shown).
To compare crosslinked-reversed with non crosslinked DNA 100 ml E. coli MG1655 LB culture was grown at 37°C to an OD600 of 0.15. After collecting 50 ml as 'non crosslinked' sample, crosslinking was done as described above. Crosslinked and non crosslinked cells were washed and sonicated corresponding to the ChIP-Chip protocol above. For experiments presented in Fig. 5B and 5D, 400 μl of the sonicated extracts were mixed with 400 μl TE and incubated with 2 μl RNase A (20 mg/ml) at 42°C for 1 h. Next, 200 μl proteinase K (20 mg/ml) were added and samples incubated for 2 h at 42 and 6 h at 65°C. For experiments without proteinase K shown in table 1, 200 μl of crosslinked extract was mixed with 200 μl TE and incubated at 65°C over night or 10 min at 100°C. For the other experiments 200 μl were mixed with 160 μl TE plus 40 μl proteinase K (20 mg/ml) and incubated at 37°C over night or for 2 h at 42°C followed by 65°C for 6 h. DNA was extracted with phenol and chlorophorm and analyzed by microarray hybridization (as above) or qPCR as described below.
Quantitative RT PCR
Reactions were carried out in triplicates of 25 μl volume each. About 10 ng DNA was used as template in 10 μl ddH2O and added to a mix of 12.5 TaqMan Gene Expression mix (Applied Biosystems) and 2.5 μl primer mix (9 μM each forward and reverse primer and 2.5 μM probe) in 96 well PCR plates. For a primer list see Table 5. Reactions were carried out with a 7500 Real Time PCR System (Applied Biosystems). The system software was used to calculate Ct values which were transformed to relative values of template DNA. qPCR values for the yahEF gene region were used for normalization.
chromatin immunoprecipitation combined with microarrays
chromatin immunoprecipitation combined with next generation sequencing
We thank Franz Narberhaus (Bochum) for the σ32 antiserum and Douglas Hurd (Oxford Gene Technology) for instruction in DNA labeling and microarray hybridization. We are grateful for the support from the Helse Sør-Øst/University of Oslo Microarray Core Facility, supported by the functional genomics programme (FUGE) of the Research Council of Norway. We thank Erik Boye, Frank Führer and Leonardo A. Meza-Zepeda for critical reading of the manuscript and the Skarstad-group for helpful discussions. Irene Kim is acknowledged for her help with submission of the microarray data to the Genome Omnibus Database. Supported by the Norwegian Research Council FUGE program and the German Research Foundation (WA 2713/1-1).
- Ren B, Robert F, Wyrick JJ, Aparicio O, Jennings EG, Simon I, Zeitlinger J, Schreiber J, Hannett N, Kanin E: Genome-wide location and function of DNA binding proteins. Science . 2000, 290 (5500): 2306-2309.PubMedView ArticleGoogle Scholar
- Buck MJ, Lieb JD: ChIP-chip: considerations for the design, analysis, and application of genome-wide chromatin immunoprecipitation experiments. Genomics. 2004, 83 (3): 349-360. 10.1016/j.ygeno.2003.11.004.PubMedView ArticleGoogle Scholar
- Grainger DC, Overton TW, Reppas N, Wade JT, Tamai E, Hobman JL, Constantinidou C, Struhl K, Church G, Busby SJ: Genomic studies with Escherichia coli MelR protein: applications of chromatin immunoprecipitation and microarrays. J Bacteriol. 2004, 186 (20): 6938-6943. 10.1128/JB.186.20.6938-6943.2004.PubMed CentralPubMedView ArticleGoogle Scholar
- Lee TI, Johnstone SE, Young RA: Chromatin immunoprecipitation and microarray-based analysis of protein location. Nat Protoc. 2006, 1 (2): 729-748. 10.1038/nprot.2006.98.PubMed CentralPubMedView ArticleGoogle Scholar
- Sandmann T, Jakobsen JS, Furlong EE: ChIP-on-chip protocol for genome-wide analysis of transcription factor binding in Drosophila melanogaster embryos. Nat Protoc. 2006, 1 (6): 2839-2855. 10.1038/nprot.2006.383.PubMedView ArticleGoogle Scholar
- Waldminghaus T, Skarstad K: The Escherichia coli SeqA protein. Plasmid. 2009, 61 (3): 141-150. 10.1016/j.plasmid.2009.02.004.PubMedView ArticleGoogle Scholar
- Slater S, Wold S, Lu M, Boye E, Skarstad K, Kleckner N: E. coli SeqA protein binds oriC in two different methyl-modulated reactions appropriate to its roles in DNA replication initiation and origin sequestration. Cell. 1995, 82 (6): 927-936. 10.1016/0092-8674(95)90272-4.PubMedView ArticleGoogle Scholar
- Campbell JL, Kleckner N: E. coli oriC and the dnaA gene promoter are sequestered from dam methyltransferase following the passage of the chromosomal replication fork. Cell. 1990, 62 (5): 967-979. 10.1016/0092-8674(90)90271-F.PubMedView ArticleGoogle Scholar
- Nagy PL, Cleary ML, Brown PO, Lieb JD: Genomewide demarcation of RNA polymerase II transcription units revealed by physical fractionation of chromatin. P Natl Acad Sci USA . 2003, 100 (11): 6364-6369. 10.1073/pnas.1131966100.View ArticleGoogle Scholar
- Wade JT, Roa DC, Grainger DC, Hurd D, Busby SJ, Struhl K, Nudler E: Extensive functional overlap between sigma factors in Escherichia coli. Nat Struct Biol. 2006, 13 (9): 806-814. 10.1038/nsmb1130.View ArticleGoogle Scholar
- Nonaka G, Blankschien M, Herman C, Gross CA, Rhodius VA: Regulon and promoter analysis of the E. coli heat-shock factor, σ32, reveals a multifaceted cellular response to heat stress. Genes Dev. 2006, 20 (13): 1776-1789. 10.1101/gad.1428206.PubMed CentralPubMedView ArticleGoogle Scholar
- Johnson DS, Li W, Gordon DB, Bhattacharjee A, Curry B, Ghosh J, Brizuela L, Carroll JS, Brown M, Flicek P: Systematic evaluation of variability in ChIP-chip experiments using predefined DNA targets. Genome Res. 2008, 18 (3): 393-403. 10.1101/gr.7080508.PubMed CentralPubMedView ArticleGoogle Scholar
- Vora T, Hottes AK, Tavazoie S: Protein occupancy landscape of a bacterial genome. Mol Cell. 2009, 35 (2): 247-253. 10.1016/j.molcel.2009.06.035.PubMed CentralPubMedView ArticleGoogle Scholar
- Bucca G, Laing E, Mersinias V, Allenby N, Hurd D, Holdstock J, Brenner V, Harrison M, Smith CP: Development and application of versatile high density microarrays for genome-wide analysis of Streptomyces coelicolor: characterization of the HspR regulon. Genome Biol. 2009, 10 (1): R5-10.1186/gb-2009-10-1-r5.PubMed CentralPubMedView ArticleGoogle Scholar
- Gruber S, Errington J: Recruitment of condensin to replication origin regions by ParB/SpoOJ promotes chromosome segregation in B. subtilis. Cell. 2009, 137 (4): 685-696. 10.1016/j.cell.2009.02.035.PubMedView ArticleGoogle Scholar
- D'Ambrosio C, Schmidt CK, Katou Y, Kelly G, Itoh T, Shirahige K, Uhlmann F: Identification of cis-acting sites for condensin loading onto budding yeast chromosomes. Genes Dev. 2008, 22 (16): 2215-2227. 10.1101/gad.1675708.PubMed CentralPubMedView ArticleGoogle Scholar
- Schmidt CK, Brookes N, Uhlmann F: Conserved features of cohesin binding along fission yeast chromosomes. Genome Biol. 2009, 10 (5): R52-10.1186/gb-2009-10-5-r52.PubMed CentralPubMedView ArticleGoogle Scholar
- Peng S, Alekseyenko AA, Larschan E, Kuroda MI, Park PJ: Normalization and experimental design for ChIP-chip data. BMC Bioinformatics. 2007, 8: 219-10.1186/1471-2105-8-219.PubMed CentralPubMedView ArticleGoogle Scholar
- Zhao K, Liu M, Burgess RR: The global transcriptional response of Escherichia coli to induced σ32 protein involves σ32 regulon activation followed by inactivation and degradation of σ32 in vivo. J Biol Chem. 2005, 280 (18): 17758-17768. 10.1074/jbc.M500393200.PubMedView ArticleGoogle Scholar
- Grainger DC, Hurd D, Goldberg MD, Busby SJ: Association of nucleoid proteins with coding and non-coding segments of the Escherichia coli genome. Nucleic Acids Res. 2006, 34 (16): 4642-4652. 10.1093/nar/gkl542.PubMed CentralPubMedView ArticleGoogle Scholar
- Cho BK, Knight EM, Barrett CL, Palsson BO: Genome-wide analysis of Fis binding in Escherichia coli indicates a causative role for A-/AT-tracts. Genome Res. 2008, 18 (6): 900-910. 10.1101/gr.070276.107.PubMed CentralPubMedView ArticleGoogle Scholar
- Carroll JS, Meyer CA, Song J, Li W, Geistlinger TR, Eeckhoute J, Brodsky AS, Keeton EK, Fertuck KC, Hall GF: Genome-wide analysis of estrogen receptor binding sites. Nat Gen. 2006, 38 (11): 1289-1297. 10.1038/ng1901.View ArticleGoogle Scholar
- Gevry N, Hardy S, Jacques PE, Laflamme L, Svotelis A, Robert F, Gaudreau L: Histone H2A.Z is essential for estrogen receptor signaling. Genes Dev. 2009, 23 (13): 1522-1533. 10.1101/gad.1787109.PubMed CentralPubMedView ArticleGoogle Scholar
- Andrau JC, van de Pasch L, Lijnzaad P, Bijma T, Koerkamp MG, van de Peppel J, Werner M, Holstege FC: Genome-wide location of the coactivator mediator: Binding without activation and transient Cdk8 interaction on DNA. Mol Cell. 2006, 22 (2): 179-192. 10.1016/j.molcel.2006.03.023.PubMedView ArticleGoogle Scholar
- Fan X, Chou DM, Struhl K: Activator-specific recruitment of Mediator in vivo. Nat Struct Mol Biol. 2006, 13 (2): 117-120. 10.1038/nsmb1049.PubMedView ArticleGoogle Scholar
- Zhu X, Wiren M, Sinha I, Rasmussen NN, Linder T, Holmberg S, Ekwall K, Gustafsson CM: Genome-wide occupancy profile of mediator and the Srb8-11 module reveals interactions with coding regions. Mol Cell. 2006, 22 (2): 169-178. 10.1016/j.molcel.2006.03.032.PubMedView ArticleGoogle Scholar
- Fan X, Struhl K: Where does mediator bind in vivo?. PLoS ONE. 2009, 4 (4): e5029-10.1371/journal.pone.0005029.PubMed CentralPubMedView ArticleGoogle Scholar
- Gordan R, Hartemink AJ, Bulyk ML: Distinguishing direct versus indirect transcription factor-DNA interactions. Genome Res. 2009, 19 (11): 2090-100. 10.1101/gr.094144.109.PubMed CentralPubMedView ArticleGoogle Scholar
- Partridge JD, Bodenmiller DM, Humphrys MS, Spiro S: NsrR targets in the Escherichia coli genome: new insights into DNA sequence requirements for binding and a role for NsrR in the regulation of motility. Mol Microbiol. 2009, 73 (4): 680-94. 10.1111/j.1365-2958.2009.06799.x.PubMedView ArticleGoogle Scholar
- Efromovich S, Grainger D, Bodenmiller D, Spiro S: Genome-wide identification of binding sites for the nitric oxide-sensitive transcriptional regulator NsrR. Methods Enzymol. 2008, 437: 211-233. full_text.PubMedView ArticleGoogle Scholar
- Auerbach RK, Euskirchen G, Rozowsky J, Lamarre-Vincent N, Moqtaderi Z, Lefrancois P, Struhl K, Gerstein M, Snyder M: Mapping accessible chromatin regions using Sono-Seq. Proc Natl Acad Sci USA. 2009, 106 (35): 14926-14931. 10.1073/pnas.0905443106.PubMed CentralPubMedView ArticleGoogle Scholar
- Lefrancois P, Euskirchen GM, Auerbach RK, Rozowsky J, Gibson T, Yellman CM, Gerstein M, Snyder M: Efficient yeast ChIP-Seq using multiplex short-read DNA sequencing. BMC Genomics. 2009, 10: 37-10.1186/1471-2164-10-37.PubMed CentralPubMedView ArticleGoogle Scholar
- Lu M, Campbell JL, Boye E, Kleckner N: SeqA: a negative modulator of replication initiation in E. coli. Cell. 1994, 77 (3): 413-426. 10.1016/0092-8674(94)90156-2.PubMedView ArticleGoogle Scholar
- Brendler T, Sawitzke J, Sergueev K, Austin S: A case for sliding SeqA tracts at anchored replication forks during Escherichia coli chromosome replication and segregation. EMBO J. 2000, 19 (22): 6249-6258. 10.1093/emboj/19.22.6249.PubMed CentralPubMedView ArticleGoogle Scholar
- Yamazoe M, Adachi S, Kanaya S, Ohsumi K, Hiraga S: Sequential binding of SeqA protein to nascent DNA segments at replication forks in synchronized cultures of Escherichia coli. Mol Microbiol. 2005, 55 (1): 289-298. 10.1111/j.1365-2958.2004.04389.x.PubMedView ArticleGoogle Scholar
- Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004, 5 (10): R80-10.1186/gb-2004-5-10-r80.PubMed CentralPubMedView ArticleGoogle Scholar
- Smyth GK, Speed T: Normalization of cDNA microarray data. Methods. 2003, 31 (4): 265-273.PubMedView ArticleGoogle Scholar
- Ghai R, Hain T, Chakraborty T: GenomeViz: visualizing microbial genomes. BMC Bioinformatics. 2004, 5: 198-10.1186/1471-2105-5-198.PubMed CentralPubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.