Description

This track shows probable sites of STAT1 binding in HeLa cells as determined by chromatin immunoprecipitation followed by microarray analysis. STAT1 (Signal Transducer and Activator of Transcription) is a transcription factor that moves to the nucleus and binds DNA only in response to a cytokine signal such as interferon-gamma. HeLa cells are a common cell line derived from a cervical cancer. Each of the four subtracks represents a different microarray platform. The track as a whole can be used to compare results across microarray platforms.

The first three platforms are custom maskless photolithographic arrays with oligonucleotides tiling most of the non-repetitive DNA sequence of the ENCODE regions:

• Maskless design #1: 50-mer oligonucleotides tiled every 38 bps (overlapping by 12 nts)
• Maskless design #2: 36-mer oligonucleotides tiled end to end
• Maskless design #3: 50-mer oligonucleotides tiled end to end

The subtracks show the ratio of immunoprecipitated DNA from cytokine-stimulated cells vs. unstimulated cells in each of the four platforms. The ratio is calculated as -log₁₀(p-value) in a 501-base window. The data shown is the combined result of multiple biological replicates: five for the first maskless array (50-mer every 38 bp), two for the second maskless array (36-mer every 36 bp), three for the third maskless array (50-mer every 50 bp) and six for the PCR Amplicon array.

These data are available at NCBI GEO as GSE2714, which also provides additional information about the experimental protocols.

Display Conventions and Configuration

This annotation follows the display conventions for composite "wiggle" tracks. The subtracks within this annotation may be configured in a variety of ways to highlight different aspects of the displayed data. The graphical configuration options are shown at the top of the track description page, followed by a list of subtracks. To display only selected subtracks, uncheck the boxes next to the tracks you wish to hide. For more information about the graphical configuration options, click the Graph configuration help link.

Methods

For all arrays, the STAT1 ChIP DNA was labeled with Cy5 and the control DNA was labeled with Cy3.

Maskless photolithographic arrays

The data from replicates were median-scaled and quantile-normalized to each other. After normalization, replicates were condensed to a single value. Using a 501 bp sliding window centered on each oligonucleotide probe, a signal map (estimating the fold enrichment [log₂ scale] of ChIP DNA) is generated by computing the pseudomedian signal of all log₂(Cy5/Cy3) ratios (median of pairwise averages) within the window (including replicates). Using the same procedure, a -log₁₀(p-value) map (measuring significance of enrichment of oligonucleotide probes in the window) for all sliding windows can be made by computing P-values using the Wilcoxon paired signed rank test comparing fluorensent intensity between Cy5 and Cy3 for each oligonucleotide probe (Cy5 and Cy3 signals from the same array). A binding site is determined by thresholding both on fold enrichment and -log₁₀(p-value) and requiring a maximum gap and a minimum run between oligonucleotide positions.

For the first maskless array (50-mer every 38 bp):
   log₂(Cy5/Cy3) >= 1.25, -log₁₀(p-value) >=8.0, MaxGap <= 100 bp, MinRun >= 180 bp

For the second maskless array (36-mer every 36 bp):
   log₂(Cy5/Cy3) >= 0.25, -log₁₀(p-value) >=4.0, MaxGap <= 250 bp, MinRun >= 0 bp

For the third maskless array (50-mer every 50 bp):
   log₂(Cy5/Cy3) >= 0.25, -log₁₀(p-value) >=4.0, MaxGap <= 250 bp, MinRun >= 0 bp

PCR Amplicon Arrays

The Cy5 and Cy3 array data were loess-normalized between channels on the same slide and then between slides. A z-score was then determined for each PCR amplicon from the distribution of log(Cy5/Cy3) in a local log(Cy5*Cy3) intensity window (see Quackenbush, 2002 and the Express Yourself website for more details). From the z-score, a P-value was then associated with each PCR amplicon. Hits were determined using a 3 sigma threshold and requiring a spot to be present on three out of six arrays.

Verification

ChIP-chip binding sites were verified by comparing "hit lists" generated from combinations of different biological replicates. Only experiments that yielded a significant overlap (greater than 50 percent) were accepted. As an independent check (for maskless arrays), data on the microarray were randomized with respect to position and re-scored; significantly fewer hits (consistent with random noise) were generated this way.

Credits

This data was generated and analyzed by the labs of Michael Snyder, Mark Gerstein and Sherman Weissman at Yale University. The PCR Amplicon arrays were manufactured by Bing Ren's lab at UCSD.

References

Cawley, S., Bekiranov, S., Ng, H.H., Kapranov, P., Sekinger, E.A., Kampa, D., Piccolboni, A., Sementchenko, V., Cheng, J. et al. Unbiased mapping of transcription factor binding sites along human chromosomes 21 and 22 points to widespread regulation of noncoding RNAs. Cell 116(4), 499-509 (2004).

Euskirchen, G., Royce, T.E., Bertone, P., Martone, R., Rinn, J.L., Nelson, F.K., Sayward, F., Luscombe, N.M., Miller, P. et al. CREB binds to multiple loci on human chromosome 22, Mol Cell Biol. 24(9), 3804-14 (2004).

Luscombe, N.M., Royce, T.E., Bertone, P., Echols, N., Horak, C.E., Chang, J.T., Snyder, M. and Gerstein, M. ExpressYourself: A modular platform for processing and visualizing microarray data. Nucleic Acids Res. 31(13), 3477-82 (2003).

Martone, R., Euskirchen, G., Bertone, P., Hartman, S., Royce, T.E., Luscombe, N.M., Rinn, J.L., Nelson, F.K., Miller, P. et al. Distribution of NF-kappaB-binding sites across human chromosome 22. Proc Natl Acad Sci U S A. 100(21), 12247-52 (2003).

Quackenbush, J.. Microarray data normalization and transformation, Nat Genet. 32(Suppl), 496-501 (2002).