Core Transcriptional Regulatory Circuitry in Human Embryonic Stem Cells

Mapping Transcription Factor Occupancy in Human Embryonic Stem Cells

• Growth Conditions
• Human ES Cell Quality Control
• Technology and Protocols

• Download Data
• Data Analysis
• Array Scanning and Data Extraction
• Normalization and Analysis
• Identification of Bound Regions
• Data Quality Control

We used the Genepix software to obtain background-subtracted intensity values for each fluorophore for every feature on the array. To obtain set-normalized intensities, we first calculated, for each slide, the median intensities in each channel for the set of 2,043 control probes described above and included on each array. We then calculated the average of these median intensities for the set of 10 slides. Intensities were then normalized such that the median intensity of each channel for an individual slide equaled the average of the median intensities of that channel across all slides.

Each slide contains a set of negative control spots that contain 60-mer sequences that do not cross-hybridize to human genomic DNA. We calculated the median intensity of these negative control spots in each channel and then subtracted this number from the set-normalized intensities of all other features.

To correct for different amounts of genomic and immunoprecipitated DNA hybridized to the chip, the set-normalized, negative control-subtracted median intensity value of the IP-enriched DNA channel was then divided by the median of the genomic DNA channel. This yielded a normalization factor that was applied to each intensity in the genomic DNA channel.

Next, we calculated the log of the ratio of intensity in the IP-enriched channel to intensity in the genomic DNA channel for each probe and used a whole chip error model (Hughes et al., 2000) to calculate confidence values for each spot on each array (single probe p-value). This error model functions by converting the intensity information in both channels to an X score which is dependent on both the absolute value of intensities and background noise in each channel. The X scores for an array are assumed to be normally distributed which allows for calculation of a p-value for the enrichment ratio seen at each feature.