- Research article
- Open Access
Novel sequential ChIP and simplified basic ChIP protocols for promoter co-occupancy and target gene identification in human embryonic stem cells
© Medeiros et al; licensee BioMed Central Ltd. 2009
- Received: 16 April 2009
- Accepted: 29 June 2009
- Published: 29 June 2009
The investigation of molecular mechanisms underlying transcriptional regulation, particularly in embryonic stem cells, has received increasing attention and involves the systematic identification of target genes and the analysis of promoter co-occupancy. High-throughput approaches based on chromatin immunoprecipitation (ChIP) have been widely used for this purpose. However, these approaches remain time-consuming, expensive, labor-intensive, involve multiple steps, and require complex statistical analysis. Advances in this field will greatly benefit from the development and use of simple, fast, sensitive and straightforward ChIP assay and analysis methodologies.
We initially developed a simplified, basic ChIP protocol that combines simplicity, speed and sensitivity. ChIP analysis by real-time PCR was compared to analysis by densitometry with the ImageJ software. This protocol allowed the rapid identification of known target genes for SOX2, NANOG, OCT3/4, SOX17, KLF4, RUNX2, OLIG2, SMAD2/3, BMI-1, and c-MYC in a human embryonic stem cell line. We then developed a novel Sequential ChIP protocol to investigate in vivo promoter co-occupancy, which is basically characterized by the absence of antibody-antigen disruption during the assay. It combines centrifugation of agarose beads and magnetic separation. Using this Sequential ChIP protocol we found that c-MYC associates with the SOX2/NANOG/OCT3/4 complex and identified a novel RUNX2/BMI-1/SMAD2/3 complex in BG01V cells. These two TF complexes associate with two distinct sets of target genes. The RUNX2/BMI-1/SMAD2/3 complex is associated predominantly with genes not expressed in undifferentiated BG01V cells, consistent with the reported role of those TFs as transcriptional repressors.
These simplified basic ChIP and novel Sequential ChIP protocols were successfully tested with a variety of antibodies with human embryonic stem cells, generated a number of novel observations for future studies and might be useful for high-throughput ChIP-based assays.
- Densitometric Analysis
- ChIP Assay
- Fold Enrichment
- Sequential ChIP
- Complex Statistical Analysis
Regulatory transcription factors (TFs) are encoded by approximately 10% of the human genome . The search for an accurate and complete list of target genes for thousands of TFs and the elucidation of their complex interactions at promoter sites, particularly in embryonic stem (ES) cells, has gained increasing interest. However, only a small fraction of the in vivo target genes and relatively few TF-TF interactions have been elucidated [2–4]. Chromatin immunoprecipitation (ChIP) and its derivatives (ChIP-chip, ChIP-seq, ChIP-SAGE, ChIP-PET, Sequential ChIP, etc) have been widely used for the investigation of TF-DNA interactions [4–9]. High-throughput approaches, such as ChIP-chip and ChIP-SAGE, are necessary for genome-wide analysis and the systematic identification of new DNA-binding sequences. Real-time (rt) PCR remains extensively used for validation of genome-wide data and for analysis of ChIP results in general. High-throughput approaches are time-consuming, expensive, labor-intensive, involve multiple steps that facilitate error introduction, and require complex statistical analysis [7, 10]. Therefore, advances in this field will greatly benefit from the development and use of faster and straightforward ChIP assay and analysis methodologies.
Here, we present data obtained with a simplified, basic ChIP assay and analysis protocol that allowed the rapid identification of known target genes for SOX2, NANOG, OCT3/4, SOX17, KLF4, RUNX2, OLIG2, SMAD2/3, BMI-1, and c-MYC in the human ES cell line BG01V. We used rtPCR to initially validate the protocol/antibodies and densitometric analysis of PCR results with the ImageJ software as a more practical, less expansive, less time-consuming readout alternative. In addition, we developed a novel, non-disruptive, highly sensitive Sequential ChIP protocol for the identification of promoter co-occupancy, based on our simplified basic ChIP protocol. The data obtained with this Sequential ChIP protocol are consistent with data previously obtained with more labor-intensive, expensive, time-consuming ChIP-chip platforms. Furthermore, Sequential ChIP analysis led to the identification of two TF complexes in BG01V ES cells: SOX2/NANOG/OCT3/4/c-MYC and RUNX2/BMI-1/SMAD2/3 complexes. These two TF complexes associate with two different sets of target genes. The RUNX2/BMI-1/SMAD2/3 complex is associated predominantly with genes not expressed in undifferentiated BG01V cells, consistent with the reported role of those TFs as transcriptional repressors. These simplified basic ChIP and novel Sequential ChIP protocols were successfully tested with a variety of antibodies with BG01V ES cells, generated a number of novel observations for future studies and might be useful for high-throughput ChIP-based assays.
Development of an improved basic ChIP protocol
Primer sets used in PCR/rtPCR reactions.
Primer sequences (Forward/Reverse)
BMI-1, implicated in cell proliferation and stem cell self-renewal [46, 47], binds to the HOXC13 promoter , which was detected by rtPCR and densitometric analysis (Fig. 3D) of BMI-1 ChIP assays using undifferentiated BG01V cells. SMAD2/3, involved in epithelial-to-mesenchymal transition, neuronal differentiation, and chondrogenesis [29, 48–50], have been shown to regulate c-MYC and GLI1 , which were also detected by rtPCR and densitometric analysis of SMAD2/3 ChIP assays using undifferentiated BG01V cells (Fig. 3E). OLIG2, implicated in neuron, astrocyte and oligodendrocyte differentiation, binds to the P21 promoter , which was detected by rtPCR and densitometric analysis (Fig. 3F) of OLIG2 ChIP assays using undifferentiated BG01V cells.
Taken together, these results confirm that our simplified, basic ChIP assay can be used with a variety of antibodies and that the densitometric quantification of PCR results using the ImageJ software can be used as an alternative ChIP readout.
SOX2/OCT3/4/NANOG/c-MYC complex co-occupies several promoters in BG01V cells
SOX2/OCT3/4NANOG and c-MYC have been found to share a large number of target genes in ChIP-chip studies [10, 39, 53–55]. We then investigated if c-MYC co-occupies promoters with the SOX2/OCT3/4/NANOG complex. Indeed, immunoprecipitation with the SOX2 antibody followed by the c-MYC antibody led to the detection of a number of promoters in undifferentiated BG01V cells (Fig. 6B, top panel, "Sequential ChIP: SOX2 > c-MYC"). The same promoters were detected if the NANOG antibody preceded the c-MYC antibody (Fig. 6B, middle panel, "Sequential ChIP: NANOG > c-MYC") or if the OCT3/4 antibody preceded the c-MYC antibody (data not shown). In contrast, none of these promoters were detected when normal IgG was used as the first antibody (Fig. 6B, last panel).
These results indicate that, at the indicated promoters (Fig. 6B), c-MYC is a subunit of the SOX2/OCT3/4/NANOG complex in undifferentiated BG01V cells. To further test this hypothesis we immunodepleted SOX2 from undifferentiated BG01V cells lysates (Fig. 6C, ERK1/2 as negative control) and then analyzed c-MYC targets by ChIP. Fig. 6D shows which c-MYC-bound promoters were detected in those SOX2-depleted lysates and shows a significant reduction in the number of c-MYC-associated promoters. Most of the promoters detected in Fig. 6B top panel ("Sequential ChIP: SOX2 > c-MYC") were not detected by c-MYC ChIP after SOX2 depletion (Fig. 6D), including P21, FGF4, NANOG, VEGF, PTEN, BAX, MUC4, OCT3/4 and SUZ12 (densitometric analysis confirmed that these signals were above input DNA levels in Fig. 6B but not in Fig. 6D, data not shown).
Collectively these results corroborate the applicability of the novel Sequential ChIP assay. In addition these results indicate that c-MYC associates with the SOX2/OCT3/4/NANOG complex in BG01V cells.
RUNX2/BMI-1/SMAD2/3 complex co-occupies a second set of promoters in BG01V cells
We then sought to investigate the presence of other TF complexes in BG01V cells using the Sequential ChIP assay. In a screening with 25 potential target genes, we found that BMI-1, RUNX2, SMAD2/3, KLF4 and OLIG2 bind to a similar set of promoters in BG01V cells (data not shown). Therefore, we use this group of TFs in Sequential ChIP assays. Seven of the promoters that were detected by most or none of those TFs were selected for this set of experiments. Our results indicate that IFNG, E-CADHERIN, and BAX promoters are co-occupied by a BMI-1/RUNX2/SMAD2/3 complex, and that the SMA promoter is co-occupied by BMI-1/SMAD2/3 (Fig. 6E). The GATA4 promoter was not detected in any of the individual ChIP assays for this group of TFs (data not shown) and was used as an additional negative control. Consistently GATA4 was not detected by Sequential ChIP using any combination of antibodies (Fig. 6E). In addition, MUC4 was not detected when BMI-1 was the first antibody in the Sequential ChIP assay (Fig. 6E), as expected, since MUC4 is not detected in BMI-1 ChIP assays (data not shown). Furthermore, Fig. 6E shows that SNAIL was not detected when the BMI-1 antibody was followed by RUNX2 or SMAD2/3 antibodies, which was also expected, since SNAIL was detected in BMI-1 ChIP assays, but not detected in RUNX2 or SMAD2/3 ChIP (data not shown). None of the above mentioned promoters were detected when the BMI-1 antibody was followed by the SOX2 antibody (Fig. 6E) or by the KLF4 or OLIG2 antibodies (data not shown). Finally, none of the above mentioned promoters were detected when the SMAD2/3 antibody was followed by the KLF4 or OLIG2 antibodies (data not shown).
Overall, these experiments suggest the presence of a RUNX2/BMI-1/SMAD2/3 complex at the above mentioned promoters in undifferentiated BG01V cells.
Results shown here indicate that our simplified, basic ChIP and the novel Sequential ChIP assays are suitable for systematic target gene identification and detection of TF interactions at specific DNA sites. The suitability of the basic ChIP assay was confirmed by our ability to rapidly identify known targets of SOX2/OCT3/4/NANOG (FGF4, LEFTY, NANOG, VEGF, BCL2, GLI1, E-CADHERIN, OCT3/4, c-MYC, HESX1, ZFP206, and SUZ12) [10, 35, 36, 39], SOX17 (mouse Lama1) , RUNX2 (VEGF and BAX) [31, 43], KLF4 (OCT3/4, NANOG, B2R) , BMI-1 (HOXC13) , SMAD2/3 (VEGF, IFNG and FOS) , and OLIG2 (P21) .
The basic ChIP protocol used here is a simplified, combined version of fast ChIP protocols previously described [11, 12]. The main advantages of this basic ChIP assay over previously published conventional protocols [[20, 31, 37], for example] is the combination of speed (only about 2 hours to obtain ChIP samples, via the elimination of some common steps) and sensitivity (known targets could be detected with only about 20,000 cell equivalents were used per PCR reaction). Our basic ChIP protocol combines biotinylated antibodies, streptavidin magnetic beads, chelating resin for elimination of PCR inhibitors , and DNA purification with commercially available silica columns. It differs from previously described fast ChIP protocols [11, 12] by not involving: antibody capture with protein A, cross-linking reversal, proteinase K treatment and phenol-chloroform-ethanol treatment for DNA isolation. In preliminary experiments we noticed that the combination of chelating resin  and silica columns for DNA purification leads to greater PCR products yields than with chelating resin alone (data not shown). The biotin-streptavidin interaction represents one of the strongest non-covalent interactions known  and therefore it represents an advantage over protein A-IgG interaction . Our preliminary observations (data not shown) and those reported by other investigators [11, 12] suggest that cross-linking reversal and proteinase K treatment can be avoided, although it might result in a slightly lower DNA yield. The boiling of the samples with chelating resin and subsequent DNA purification with silica columns seem to replace the need for cross-linking reversal and proteinase K treatment.
The use of densitometric analysis of PCR results qualitatively matched the rtPCR results for all the antibodies used here since all above mentioned known targets were detected above background by both densitometry and rtPCR (Figs. 2 and 3). Densitometric analysis with the ImageJ software, due to its simplicity, practicality and less time and resources used in preliminary calibration experiments, offers a viable alternative to rtPCR when absolute quantification is not required. In addition, it can be applied for pre-screenings prior to rtPCR. It was expected, however, that densitometric analysis with the ImageJ software would show a quantification limitation, probably coming from pixel saturation in DNA bands from PCR reactions and consequent lower measurement capacity. Such limitation probably caused the lower levels of fold enrichment observed in the densitometric analysis as compared to rtPCR for some of the ChIP assays (Fig. 3; KLF4, BMI-1 and OLIG2 ChIP assay). However, our comparisons indicate that when the goal is simply the identification of target genes, densitometric analysis provides a reliable alternative, besides, for some antibodies, such limitation was minimum (Fig. 3; SOX17 and SMAD2/3 ChIP assays).
We then analyzed promoter co-occupancy using the newly developed Sequential ChIP assay (diagram in Fig. 5). This assay provides a significant improvement from previously described Sequential ChIP assays [57–60] and it differs from those in two main aspects. First, during the Sequential ChIP assay the antibody-antigen interaction is never disrupted by epitope-specific peptides or other treatments, and consequently immunoprecipitations can be carried-out in the same, intact samples. This feature might contribute for lower background, higher specificity and the possibility of using smaller samples. Second, the modifications introduced in the basic ChIP assay to make it a straightforward, fast, sensitive assay, were maintained in the Sequential ChIP protocol. In addition, the Sequential ChIP combines two methods of physical separation of the targeted DNA fragments. The first antibody, non-biotinylated, binds to agarose beads, via protein G-IgG interaction, and is brought down by centrifugation. The second antibody, this time a biotinylated one, binds to magnetic beads, via biotin-streptavidin interaction, and is separated from the rest of the sample with a magnet. These two forms of physical separation (gravity versus magnetism) are quite complementary and exert minimum interference over each other. The results are even visual with the reagents used here: when the antibodies used in the assay recognize TFs that co-occupy a significant number of target genes, as it is the case for SOX2/OCT3/4/NANOG/c-MYC, one can see the agarose beads slurry attaching to the magnet in the last steps of the assay, since they are carried by the streptavidin-magnetic beads. This visual effect is not observed with any of the negative controls. In proof-of-principle experiments we investigated promoter co-occupancy in undifferentiated BG01V cells using SOX2, OCT3/4 and NANOG antibodies (Fig. 6A) and the obtained results are consistent with those obtained in ChIP-chip studies [10, 39].
Furthermore, our Sequential ChIP assay identified c-MYC as a potential interacting partner of the SOX2/OCT3/4/NANOG complex at several gene promoters in undifferentiated BG01V cells (Fig. 6). It has been suggested that c-MYC is involved in the regulation of up to 15% of all genes. It has been reported that c-MYC modulates global chromatin structure via its interaction with histone acetyl-transferase (HAT) and SWI/SNF chromatin remodeling complexes [61–63], which might be c-MYC function in the SOX2/OCT3/4/NANOG complex. Supporting this hypothesis, SOX2/OCT3/4/NANOG and c-MYC have been found to share a large number of target genes in ChIP-chip studies [10, 39, 53–55].
Furthermore, data shown in Fig. 6E suggests the presence of a RUNX2/BMI-1/SMAD2/3 complex in undifferentiated BG01V cells. Interestingly, the association of SMAD2/3 with RUNX2 has been reported previously in B-lymphocytes and human breast cancer cells [64, 65]. BMI-1 is part of the Polycomb-group complex [46, 47] and a known transcriptional repressor. SMAD2/3 and RUNX2 have been reported to act as both activators and repressors [50, 66]. Accordingly, many of the genes associated with this RUNX2/BMI-1/SMAD2/3 complex in BG01V cells (for example SMA, IFNG, GATA4 and 6, MUC4, FOS, SNAIL, P21; data not shown) are not expressed in undifferentiated ES cells [38, 54], suggesting that this RUNX2/BMI-1/SMAD2/3 complex plays a transcriptional repressor role in undifferentiated BG01V cells.
In summary, our simplified basic ChIP and novel Sequential ChIP protocols were successfully tested with a variety of antibodies and generated a number of novel observations for future studies. The basic ChIP assay combines simplicity (avoidance of several commonly used steps), speed (whole assay performed in 2 hours) and sensitivity (using 20,000 cells). More importantly, promoter co-occupancy data obtained with our novel Sequential ChIP protocol indicates the existence of two TF complexes in human ES cells: SOX2/OCT3/4/NANOG/c-MYC and RUNX2/BMI-1/SMAD2/3. This novel Sequential ChIP protocol provides a significant improvement from previously described Sequential ChIP assays by combining two methods of physical separation of the targeted DNA (centrifugation and magnetism) and avoiding antibody-antigen interaction disruption during the assay. Both protocols might be useful for high-throughput ChIP-based assays, which are expensive, labor intensive, require complex statistical analysis, and might greatly benefit from faster and straightforward ChIP assay and analysis methodologies.
Cells and cell culture
Mouse ES cell line D3 was purchased from American Type Culture Collection (ATCC) and cultured according to the manufacturer's instructions. Undifferentiated human ES cell line BG01V (Novocell) were cultured according to the manufacturer's instructions, without feeder cells . Endoderm differentiation of D3 cells was performed as described previously .
Primary antibodies were used at 10 μg/ml. BG01V cells were cultured in 24-well plates, fixed with 4% paraformaldehyde/PBS at room temperature (RT) for 20 min, blocked and permeabilized with 0.1% Triton X-100, 1% BSA, and 10% normal donkey serum in PBS, at RT, for 45 min. After blocking, cells were incubated with diluted primary antibody overnight at 4°C, followed by 5 μg/ml secondary antibody (anti-IgG-NL557, R&D Systems) at RT, in the dark, for 1 hour. Cells were washed with PBS/0.1% BSA between each step. No staining was obtained with secondary antibody only, used as a negative control (data not shown).
Reagents and antibodies
Anti-SOX2, NANOG, OCT3/4, c-MYC, PDX1, OLIG2, RUNX2, SMAD2/3, KLF4, ERK1/2 and BMI-1 antibodies, streptavidin magnetic beads (R&D Systems); antibody specificity was determined previously with appropriate controls by Western blotting and immunocytochemistry [ref.  and data not shown]; Odyssey western blot blocker buffer (LI-COR Biosciences); IDRdye 800 anti-goat IgG antibody (Rockland Immunochemicals); phorbol-12-myristate-13-acetate (PMA), Chelex-100, lithium chloride, sodium deoxycholate, leupeptin, phenyl-methanessulfonyl-floride (PMSF), 37% formaldehyde, NP-40, glycine, dimethyl-sulfoxide (DMSO), DAPI, aprotinin (Sigma); AmpliTaq® Gold Polymerase PCR kit, dNTP mix (Applied Biosystems); EXPRESS SYBR® GreenER™ qPCR kit (Invitrogen); QIAQuick® DNA purification kit (Qiagen); 96-well PCR plates, 96-well plate foil (Eppendorf), protein A-agarose beads (Thermo Scientific), and protein G-agarose beads (Pierce).
Western blotting and Immunoprecipitation
Cells were lysed with 2× lysis buffer (2% NP-40, 100 mM Tris-HCl, 300 mM NaCl, 4 mM EDTA, 2 mM sodium vanadate, 10 μg/ml leupeptin, 2 mM PMSF, 10 μg/ml aprotinin) and equal number of cell equivalents were separated by 10% SDS-PAGE. After Western transfer, nitrocellulose membranes were blocked with Odyssey blocker/PBS (1:1) for 1 hour, incubated with primary antibodies for 30 min in Odyssey blocker/PBS/0.2% Tween, washed with PBS/0.2% Tween 3× for 5 min, incubated for 30 min with secondary fluorescent antibody (IDRdye 800 donkey anti-goat IgG) in Odyssey blocker/PBS/0.2%Tween/0.02%SDS, and washed with PBS/0.2%Tween 3× for 5 min. Fluorescence was detected and images captured with an Odyssey infrared imager (LI-COR Biosciences). For immunoprecipitation (IP) lysates were incubated with 5 μg of biotinylated anti-SOX2 antibody for 4 hrs at 4°C, for each round of IP. Samples were incubated with 20 μL of streptavidin magnetic beads (R&D Systems) for 30 min at 4°C, beads were collected with a magnet, and washed with wash buffer (20 mM HEPES pH 7.6, 50 mM NaCl, 2.5 mM MgCl2, 0.1 mM EDTA, and 0.05% Triton X-100).
Basic ChIP Protocol
The basic ChIP assay (diagram shown in Fig. 1) is a combination of sensitive and fast ChIP assays described previously [11, 12] with some modifications. Cells were fixed with 1% formaldehyde for 15 min (RT), quenched with 125 mM glycine for 5 min (RT), centrifuged and resuspended in Lysis Buffer (1% SDS, 10 mM EDTA, 50 mM Tris, pH 8) containing protease inhibitors (10 μg/mL leupeptin, 10 μg/mL aprotinin, and 1 mM PMSF) and incubated on ice (10 min). Initial samples consisted of about 1 × 106 cells in 0.5 mL of Lysis Buffer. Next samples were sonicated (Heat Systems-Ultrasonic device) to shear chromatin to an average length of about 1 kb and transferred to 1.5 mL tubes, microcentrifuged for 10 min (max speed). Supernatants were collected in 1.5 mL tubes containing 1 mL of the Dilution Buffer (0.01% SDS, 1.1% Triton, 1.2 mM EDTA, 167 mM NaCl, 17 mM Tris, pH 8). Five μg of biotinylated primary antibodies were added, samples were incubated (15 min) in an ultrasonic bath (Branson), followed by addition of 50 μL of streptavidin magnetic beads, incubated (30 min, 4°C) on a rotator (Labnet Int.). Beads were collected with a magnet (R&D Systems), washed 4× with 1 mL of each of four Wash Buffers (Wash Buffer 1: 0.1% SDS, 1% Triton, 2 mM EDTA, 150 mM NaCl, 20 mM Tris, pH 8; Wash Buffer 2: 0.1% SDS, 1% Triton, 2 mM EDTA, 500 mM NaCl, 20 mM Tris, pH 8; Wash Buffer 3: 0.25 M LiCl, 1% NP-40, 1% deoxycholate, 1 mM EDTA, 10 mM Tris, pH 8; Wash Buffer 4: 10 mM Tris, pH 8, 1 mM EDTA). After the last wash, 100 μL of a 10% Chelex-100 (Sigma)/PBS resin solution was added to the beads, samples were boiled (10 min) in a heat block, microcentrifuged (1 min, max speed), supernatants transferred to a 1.5 mL tube, followed by the addition of 120 μL of MilliQ water back to the beads, microcentrifuged again (1 min, max speed), and the new supernatant pooled with the previous one. DNA samples were then cleaned up with QIAQuick kit®, resuspended in 50 μL of Wash Buffer 4 and 1 μL (containing about 20,000 cell equivalents) was used per PCR reaction. Experiments were performed 3×, with independent samples.
The Sequential ChIP protocol (diagram shown in Fig. 5) was performed as follows. Samples were treated as above until addition of the first antibody. Five μg of the first antibody was added, this time a non-biotinylated antibody, incubated for 2 hrs to overnight at 4°C, followed by 50 μL of protein G-agarose beads for 30 min at 4°C on a rotator (Labnet Int.). Samples were then washed as described above, using a microcentrifuge instead of a magnet to pull-down the agarose beads. Subsequently the agarose beads were resuspended in 1 mL of Dilution buffer. Five μg of the second antibody was added, this time a biotinylated antibody, incubated for 2 hrs to overnight at 4°C, followed by incubation with magnetic streptavidin beads. Beads were washed as described above, using a magnet, and the final DNA samples were obtained as described above (basic ChIP protocol). Experiments were performed 3×, with independent samples.
PCR and Densitometric Image Analysis
DNA samples, from ChIP assays or input DNA (0.1% of the sample used for the ChIP assay), primers (1 μM, IDT) and DMSO (5%) were added to the additional PCR reagents (AmpliTaq® Gold Polymerase PCR kit and dNTP mix) in a 50 μL reaction volume and subjected to the following cycle: 95°C/9 min, 43× (95°C/1 min, 60°C/1 min), and 60°C/10 min, using 96-well plates and an Eppendorf Mastercycler®. PCR products were subjected to a 1.5% agarose gel electrophoresis, visualized with Ethidium bromide staining and photographed. Images were saved as TIFF files for analysis with ImageJ. All primers have been validated previously (Table 1). Signal intensities from PCR data obtained from ChIP assays or from whole cell lysates (Input DNA) were quantified from the TIFF images with ImageJ software (National Institutes of Health; http://rsbweb.nih.gov/ij/), as previously described  with some modifications. Initially images were transformed to 16-bit-type images. The threshold function was set to black and white type of image to eliminate background interference. Threshold limits were set to compensate possible saturation interference. The rectangle tool was used to draw the area for measurement around the PCR bands and the same rectangle tool was moved to each band in the blot. Area, area fraction and mean gray values were taken. Mean gray values were used as a measure of relative number of pixels. Fold enrichment = ChIP/Input DNA. For each TF three independent ChIP experiments were performed. Average and standard deviations (SD) were calculated based on the signal intensities from each experiment.
Real-time PCR (rtPCR)
DNA samples, from ChIP assays or input DNA (0.1% of the sample used for the ChIP assay), and primers (200 nM, IDT) were added to the additional rtPCR reagents (EXPRESS SYBR® GreenER™ qPCR kit) in a 20 μL reaction volume and subjected to the following cycle: 50°C/2 min, 95°C/2 min, 40× (95°C/15 sec, 60°C/1 min), using ABI 7900 HT or Cepheid Smart rtPCR instruments. No template (primers only) reactions were used for normalization. Each sample was amplified in triplicate. Three independent experiments were performed. Average and SD were calculated based on relative copy numbers from each experiment.
We thank Dr. Y. Shimizu (Univ. Minnesota) for generously providing several reagents, and Drs. M. Donohoe (Cornell Univ.) and D. Mueller (Univ. Minnesota) for critical suggestions.
- Levine M, Tjian R: Transcription regulation and animal diversity. Nature. 2003, 424: 147-151. 10.1038/nature01763.View ArticleGoogle Scholar
- Meng X, Brodsky MH, Wolfe SA: A bacterial one-hybrid system for determining the DNA-binding specificity of transcription factors. Nat Biotechnol. 2005, 23: 988-994. 10.1038/nbt1120.View ArticleGoogle Scholar
- Mukherjee S, Berger MF, Jona G, Wang XS, Muzzey D, Snyder M, Young RA, Bulyk ML: Rapid analysis of the DNA-binding specificities of transcription factors with DNA microarrays. Nat Genet. 2004, 36: 1331-1339. 10.1038/ng1473.View ArticleGoogle Scholar
- Walhout AJ: Unraveling transcription regulatory networks by protein-DNA and protein-protein interaction mapping. Genome Res. 2006, 16: 1445-1454. 10.1101/gr.5321506.View ArticleGoogle Scholar
- Solomon MJ, Larsen PL, Varshavsky A: Mapping protein-DNA interactions in vivo with formaldehyde: evidence that histone H4 is retained on a highly transcribed gene. Cell. 1988, 53: 937-947. 10.1016/S0092-8674(88)90469-2.View ArticleGoogle Scholar
- Ren B, Robert F, Wyrick JJ, Aparicio O, Jennings EG, Simon I, Zeitlinger J, Schreiber J, Hannett N, Kanin E, Volkert TL, Wilson CJ, Bell SP, Young RA: Genome-wide location and function of DNA binding proteins. Science. 2000, 290: 2306-2309. 10.1126/science.290.5500.2306.View ArticleGoogle Scholar
- Kim TH, Ren B: Genome-wide analysis of protein-DNA interactions. Annu Rev Genomics Hum Genet. 2006, 7: 81-102. 10.1146/annurev.genom.7.080505.115634.View ArticleGoogle Scholar
- Robertson G, Hirst M, Bainbridge M, Bilenky M, Zhao Y, Zeng T, Euskirchen G, Bernier B, Varhol R, Delaney A, Thiessen N, Griffith OL, He A, Marra M, Snyder M, Jones S: Genome-wide profiles of STAT1 DNA association using chromatin immunoprecipitation and massively parallel sequencing. Nat Methods. 2007, 4: 651-657. 10.1038/nmeth1068.View ArticleGoogle Scholar
- Collas P, Dahl JA: Chop it, ChIP it, check it: the current status of chromatin immunoprecipitation. Front Biosci. 2008, 13: 929-943. 10.2741/2733.View ArticleGoogle Scholar
- Sharov AA, Masui S, Sharova LV, Piao Y, Aiba K, Matoba R, Xin L, Niwa H, Ko MS: Identification of Pou5f1, Sox2, and Nanog downstream target genes with statistical confidence by applying a novel algorithm to time course microarray and genome-wide chromatin immunoprecipitation data. BMC Genomics. 2008, 9: 269-10.1186/1471-2164-9-269.View ArticleGoogle Scholar
- Nelson JD, Denisenko O, Bomsztyk K: Protocol for the fast chromatin immunoprecipitation (ChIP) method. Nat Protoc. 2006, 1: 179-185. 10.1038/nprot.2006.27.View ArticleGoogle Scholar
- Dahl JA, Collas P: A rapid micro chromatin immunoprecipitation assay (microChIP). Nat Protoc. 2008, 3: 1032-1045. 10.1038/nprot.2008.68.View ArticleGoogle Scholar
- Cai J, Olson JM, Rao MS, Stanley M, Taylor E, Ni HT: Development of antibodies to human embryonic stem cell antigens. BMC Dev Biol. 2005, 5: 26-10.1186/1471-213X-5-26.View ArticleGoogle Scholar
- Masui S, Nakatake Y, Toyooka Y, Shimosato D, Yagi R, Takahashi K, Okochi H, Okuda A, Matoba R, Sharov AA, Ko MS, Niwa H: Pluripotency governed by Sox2 via regulation of Oct3/4 expression in mouse embryonic stem cells. Nat Cell Biol. 2007, 9: 625-635. 10.1038/ncb1589.View ArticleGoogle Scholar
- Su L, David M: Distinct mechanisms of STAT phosphorylation via the interferon-alpha/beta receptor. Selective inhibition of STAT3 and STAT5 by piceatannol. J Biol Chem. 2000, 275: 12661-12666. 10.1074/jbc.275.17.12661.View ArticleGoogle Scholar
- Nelson EA, Walker SR, Alvarez JV, Frank DA: Isolation of unique STAT5 targets by chromatin immunoprecipitation-based gene identification. J Biol Chem. 2004, 279: 54724-54730. 10.1074/jbc.M408464200.View ArticleGoogle Scholar
- Barlev NA, Liu L, Chehab NH, Mansfield K, Harris KG, Halazonetis TD, Berger SL: Acetylation of p53 activates transcription through recruitment of coactivators/histone acetyltransferases. Mol Cell. 2001, 8: 1243-1254. 10.1016/S1097-2765(01)00414-2.View ArticleGoogle Scholar
- Kirmizis A, Bartley SM, Farnham PJ: Identification of the polycomb group protein SU(Z)12 as a potential molecular target for human cancer therapy. Mol Cancer Ther. 2003, 2: 113-121.Google Scholar
- Jonckheere N, Vincent A, Perrais M, Ducourouble MP, Male AK, Aubert JP, Pigny P, Carraway KL, Freund JN, Renes IB, Van Seuningen I: The human mucin MUC4 is transcriptionally regulated by caudal-related homeobox, hepatocyte nuclear factors, forkhead box A, and GATA endodermal transcription factors in epithelial cancer cells. J Biol Chem. 2007, 282: 22638-22650. 10.1074/jbc.M700905200.View ArticleGoogle Scholar
- Nandiwada SL, Li W, Zhang R, Mueller DL: p300/Cyclic AMP-responsive element binding-binding protein mediates transcriptional coactivation by the CD28 T cell costimulatory receptor. J Immunol. 2006, 177: 401-413.View ArticleGoogle Scholar
- Cao R, Tsukada Y, Zhang Y: Role of Bmi-1 and Ring1A in H2A ubiquitylation and Hox gene silencing. Mol Cell. 2005, 20: 845-854. 10.1016/j.molcel.2005.12.002.View ArticleGoogle Scholar
- Hu MC, Mo R, Bhella S, Wilson CW, Chuang PT, Hui CC, Rosenblum ND: GLI3-dependent transcriptional repression of Gli1, Gli2 and kidney patterning genes disrupts renal morphogenesis. Development. 2006, 133: 569-578. 10.1242/dev.02220.View ArticleGoogle Scholar
- Regl G, Kasper M, Schnidar H, Eichberger T, Neill GW, Philpott MP, Esterbauer H, Hauser-Kronberger C, Frischauf AM, Aberger F: Activation of the BCL2 promoter in response to Hedgehog/GLI signal transduction is predominantly mediated by GLI2. Cancer Res. 2004, 64: 7724-7731. 10.1158/0008-5472.CAN-04-1085.View ArticleGoogle Scholar
- Peinado H, Olmeda D, Cano A: Snail, Zeb and bHLH factors in tumour progression: an alliance against the epithelial phenotype?. Nat Rev Cancer. 2007, 7: 415-428. 10.1038/nrc2131.View ArticleGoogle Scholar
- Yu P, Kodadek T: Dynamics of the hypoxia-inducible factor-1-vascular endothelial growth factor promoter complex. J Biol Chem. 2007, 282: 35035-35045. 10.1074/jbc.M707557200.View ArticleGoogle Scholar
- Peiro S, Escriva M, Puig I, Barbera MJ, Dave N, Herranz N, Larriba MJ, Takkunen M, Franci C, Munoz A, Virtanen I, Baulida J, Garcia de Herreros A: Snail1 transcriptional repressor binds to its own promoter and controls its expression. Nucleic Acids Res. 2006, 34: 2077-2084. 10.1093/nar/gkl141.View ArticleGoogle Scholar
- Escriva M, Peiro S, Herranz N, Villagrasa P, Dave N, Montserrat-Sentis B, Murray SA, Franci C, Gridley T, Virtanen I, Garcia de Herreros A: Repression of PTEN phosphatase by Snail1 transcriptional factor during gamma radiation-induced apoptosis. Mol Cell Biol. 2008, 28: 1528-1540. 10.1128/MCB.02061-07.View ArticleGoogle Scholar
- Tang Y, Urs S, Liaw L: Hairy-related transcription factors inhibit Notch-induced smooth muscle alpha-actin expression by interfering with Notch intracellular domain/CBF-1 complex interaction with the CBF-1-binding site. Circ Res. 2008, 102: 661-668. 10.1161/CIRCRESAHA.107.165134.View ArticleGoogle Scholar
- Furumatsu T, Tsuda M, Taniguchi N, Tajima Y, Asahara H: Smad3 induces chondrogenesis through the activation of SOX9 via CREB-binding protein/p300 recruitment. J Biol Chem. 2005, 280: 8343-8350. 10.1074/jbc.M413913200.View ArticleGoogle Scholar
- Caslini C, Capo-chichi CD, Roland IH, Nicolas E, Yeung AT, Xu XX: Histone modifications silence the GATA transcription factor genes in ovarian cancer. Oncogene. 2006, 25: 5446-5461. 10.1038/sj.onc.1209533.View ArticleGoogle Scholar
- Eliseev RA, Dong YF, Sampson E, Zuscik MJ, Schwarz EM, O'Keefe RJ, Rosier RN, Drissi MH: Runx2-mediated activation of the Bax gene increases osteosarcoma cell sensitivity to apoptosis. Oncogene. 2008, 27: 3605-3614. 10.1038/sj.onc.1211020.View ArticleGoogle Scholar
- Liu G, Ding W, Neiman J, Mulder KM: Requirement of Smad3 and CREB-1 in mediating transforming growth factor-beta (TGF beta) induction of TGF beta 3 secretion. J Biol Chem. 2006, 281: 29479-29490. 10.1074/jbc.M600579200.View ArticleGoogle Scholar
- Saifudeen Z, Dipp S, Fan H, El-Dahr SS: Combinatorial control of the bradykinin B2 receptor promoter by p53, CREB, KLF-4, and CBP: implications for terminal nephron differentiation. Am J Physiol Renal Physiol. 2005, 288: F899-909. 10.1152/ajprenal.00370.2004.View ArticleGoogle Scholar
- Chen C, Rowell EA, Thomas RM, Hancock WW, Wells AD: Transcriptional regulation by Foxp3 is associated with direct promoter occupancy and modulation of histone acetylation. J Biol Chem. 2006, 281: 36828-36834. 10.1074/jbc.M608848200.View ArticleGoogle Scholar
- Wang ZX, Teh CH, Kueh JL, Lufkin T, Robson P, Stanton LW: Oct4 and Sox2 directly regulate expression of another pluripotency transcription factor, Zfp206, in embryonic stem cells. J Biol Chem. 2007, 282: 12822-12830. 10.1074/jbc.M611814200.View ArticleGoogle Scholar
- Chakravarthy H, Boer B, Desler M, Mallanna SK, McKeithan TW, Rizzino A: Identification of DPPA4 and other genes as putative Sox2:Oct-3/4 target genes using a combination of in silico analysis and transcription-based assays. J Cell Physiol. 2008, 216: 651-662. 10.1002/jcp.21440.View ArticleGoogle Scholar
- Patterson ES, Addis RC, Shamblott MJ, Gearhart JD: SOX17 directly activates Zfp202 transcription during in vitro endoderm differentiation. Physiol Genomics. 2008, 34: 277-284. 10.1152/physiolgenomics.90236.2008.View ArticleGoogle Scholar
- Zeng X, Miura T, Luo Y, Bhattacharya B, Condie B, Chen J, Ginis I, Lyons I, Mejido J, Puri RK, Rao MS, Freed WJ: Properties of pluripotent human embryonic stem cells BG01 and BG02. Stem Cells. 2004, 22: 292-312. 10.1634/stemcells.22-3-292.View ArticleGoogle Scholar
- Boyer LA, Lee TI, Cole MF, Johnstone SE, Levine SS, Zucker JP, Guenther MG, Kumar RM, Murray HL, Jenner RG, Gifford DK, Melton DA, Jaenisch R, Young RA: Core transcriptional regulatory circuitry in human embryonic stem cells. Cell. 2005, 122: 947-956. 10.1016/j.cell.2005.08.020.View ArticleGoogle Scholar
- Kanai-Azuma M, Kanai Y, Gad JM, Tajima Y, Taya C, Kurohmaru M, Sanai Y, Yonekawa H, Yazaki K, Tam PP, Hayashi Y: Depletion of definitive gut endoderm in Sox17-null mutant mice. Development. 2002, 129: 2367-2379.Google Scholar
- Niimi T, Hayashi Y, Futaki S, Sekiguchi K: SOX7 and SOX17 regulate the parietal endoderm-specific enhancer activity of mouse laminin alpha1 gene. J Biol Chem. 2004, 279: 38055-38061. 10.1074/jbc.M403724200.View ArticleGoogle Scholar
- Komori T: Regulation of bone development and maintenance by Runx2. Front Biosci. 2008, 13: 898-903. 10.2741/2730.View ArticleGoogle Scholar
- Zelzer E, Glotzer DJ, Hartmann C, Thomas D, Fukai N, Soker S, Olsen BR: Tissue specific regulation of VEGF expression during bone development requires Cbfa1/Runx2. Mech Dev. 2001, 106: 97-106. 10.1016/S0925-4773(01)00428-2.View ArticleGoogle Scholar
- Liu Y, Sinha S, McDonald OG, Shang Y, Hoofnagle MH, Owens GK: Kruppel-like factor 4 abrogates myocardin-induced activation of smooth muscle gene expression. J Biol Chem. 2005, 280: 9719-9727. 10.1074/jbc.M412862200.View ArticleGoogle Scholar
- Jiang J, Chan YS, Loh YH, Cai J, Tong GQ, Lim CA, Robson P, Zhong S, Ng HH: A core Klf circuitry regulates self-renewal of embryonic stem cells. Nat Cell Biol. 2008, 10: 353-360. 10.1038/ncb1698.View ArticleGoogle Scholar
- Jacobs JJ, Kieboom K, Marino S, DePinho RA, van Lohuizen M: The oncogene and Polycomb-group gene bmi-1 regulates cell proliferation and senescence through the ink4a locus. Nature. 1999, 397: 164-168. 10.1038/16476.View ArticleGoogle Scholar
- Molofsky AV, He S, Bydon M, Morrison SJ, Pardal R: Bmi-1 promotes neural stem cell self-renewal and neural development but not mouse growth and survival by repressing the p16Ink4a and p19Arf senescence pathways. Genes Dev. 2005, 19: 1432-1437. 10.1101/gad.1299505.View ArticleGoogle Scholar
- Roberts AB, Tian F, Byfield SD, Stuelten C, Ooshima A, Saika S, Flanders KC: Smad3 is key to TGF-beta-mediated epithelial-to-mesenchymal transition, fibrosis, tumor suppression and metastasis. Cytokine Growth Factor Rev. 2006, 17: 19-27. 10.1016/j.cytogfr.2005.09.008.View ArticleGoogle Scholar
- Garcia-Campmany L, Marti E: The TGFbeta intracellular effector Smad3 regulates neuronal differentiation and cell fate specification in the developing spinal cord. Development. 2007, 134: 65-75. 10.1242/dev.02702.View ArticleGoogle Scholar
- Brown KA, Pietenpol JA, Moses HL: A tale of two proteins: differential roles and regulation of Smad2 and Smad3 in TGF-beta signaling. J Cell Biochem. 2007, 101: 9-33. 10.1002/jcb.21255.View ArticleGoogle Scholar
- Ligon KL, Fancy SP, Franklin RJ, Rowitch DH: Olig gene function in CNS development and disease. Glia. 2006, 54: 1-10. 10.1002/glia.20273.View ArticleGoogle Scholar
- Fellous TG, Guppy NJ, Brittan M, Alison MR: Cellular pathways to beta-cell replacement. Diabetes/Metabolism Res and Rev. 2007, 23: 87-99. 10.1002/dmrr.692.View ArticleGoogle Scholar
- Zeller KI, Jegga AG, Aronow BJ, O'Donnell KA, Dang CV: An integrated database of genes responsive to the Myc oncogenic transcription factor: identification of direct genomic targets. Genome Biol. 2003, 4: R69-10.1186/gb-2003-4-10-r69.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Kim J, Chu J, Shen X, Wang J, Orkin SH: An extended transcriptional network for pluripotency of embryonic stem cells. Cell. 2008, 132: 1049-1061. 10.1016/j.cell.2008.02.039.View ArticleGoogle Scholar
- van Werven FJ, Marc Timmers HTh: The use of biotin tagging in Sccharomyces cerevisiae improves the sensitivity of chromatin immunoprecipitation. Nucleic Acids Res. 2006, 34: e33-10.1093/nar/gkl003.View ArticleGoogle Scholar
- Scully KM, Jacobson EM, Jepsen K, Lunyak V, Viadiu H, Carriere C, Rose DW, Hooshmand F, Aggarwal AK, Rosenfeld MG: Allosteric effects of Pit-1 DNA sites on long-term repression in cell type specification. Science. 2000, 290: 1127-1131. 10.1126/science.290.5494.1127.View ArticleGoogle Scholar
- Metivier R, Penot G, Hubner MR, Reid G, Brand H, Kos M, Gannon F: Estrogen receptor-alpha directs ordered, cyclical, and combinatorial recruitment of cofactors on a natural target promoter. Cell. 2003, 115: 751-763. 10.1016/S0092-8674(03)00934-6.View ArticleGoogle Scholar
- Geisberg JV, Struhl K: Quantitative sequential chromatin immunoprecipitation, a method for analyzing co-occupancy of proteins at genomic regions in vivo. Nucleic Acids Res. 2004, 32: e151-10.1093/nar/gnh148.View ArticleGoogle Scholar
- Geisberg JV, Struhl K: Cellular stress alters the transcriptional properties of promoter-bound Mot1-TBP complexes. Mol Cell. 2004, 14: 479-489. 10.1016/j.molcel.2004.05.003.View ArticleGoogle Scholar
- Patel JH, Du Y, Ard PG, Phillips C, Carella B, Chen CJ, Rakowski C, Chatterjee C, Lieberman PM, Lane WS, Blobel GA, McMahon SB: The c-MYC oncoprotein is a substrate of the acetyltransferases hGCN5/PCAF and TIP60. Mol Cell Biol. 2004, 24: 10826-10834. 10.1128/MCB.24.24.10826-10834.2004.View ArticleGoogle Scholar
- Fernandez PC, Frank SR, Wang L, Schroeder M, Liu S, Greene J, Cocito A, Amati B: Genomic targets of the human c-Myc protein. Genes Dev. 2003, 17: 1115-1129. 10.1101/gad.1067003.View ArticleGoogle Scholar
- Amati B, Frank SR, Donjerkovic D, Taubert S: Function of the c-Myc oncoprotein in chromatin remodeling and transcription. Biochim Biophys Acta. 2001, 1471: M135-45.Google Scholar
- Selvamurugan N, Kwok S, Alliston T, Reiss M, Partridge NC: Transforming growth factor-beta 1 regulation of collagenase-3 expression in osteoblastic cells by cross-talk between the Smad and MAPK signaling pathways and their components, Smad2 and Runx2. J Biol Chem. 2004, 279: 19327-19334. 10.1074/jbc.M314048200.View ArticleGoogle Scholar
- Hanai J, Chen LF, Kanno T, Ohtani-Fujita N, Kim WY, Guo WH, Imamura T, Ishidou Y, Fukuchi M, Shi MJ, Stavnezer J, Kawabata M, Miyazono K, Ito Y: Interaction and functional cooperation of PEBP2/CBF with Smads. Synergistic induction of the immunoglobulin germline Calpha promoter. J Biol Chem. 1999, 274: 31577-31582. 10.1074/jbc.274.44.31577.View ArticleGoogle Scholar
- Jensen ED, Nair AK, Westendorf JJ: Histone deacetylase co-repressor complex control of Runx2 and bone formation. Crit Rev Eukaryot Gene Expr. 2007, 17: 187-96.View ArticleGoogle Scholar
- D'Amour KA, Bang AG, Eliazer S, Kelly OG, Agulnick AD, Smart NG, Moorman MA, Kroon E, Carpenter MK, Baetge EE: Production of pancreatic hormone-expressing endocrine cells from human embryonic stem cells. Nat Biotechnol. 2006, 24: 1392-1401. 10.1038/nbt1259.View ArticleGoogle Scholar
- Chapman MD, Keir G, Petzold A, Thompson EJ: Measurement of high affinity antibodies on antigen-immunoblots. J Immunol Methods. 2006, 310: 62-66. 10.1016/j.jim.2005.12.004.View 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.