il-4(+)
DESCRIPTION
A. B. (fold). (fold). IL-4(+). IL-4(+). IL-4(-). IL-4(-). Enrichment against background. Enrichment against background. C. D. Supplementary Figure S1. (fold). (fold). E. F. Enrichment against background. Enrichment against background. H3K4me3. Pol II. H3Ac. H3Ac. H3K4me3. - PowerPoint PPT PresentationTRANSCRIPT
A B
IL-4(+)IL-4(-)
IL-4(+)IL-4(-)
ChIP-Seq (STAT6)Ramos IL-4 (+)P-value
Ramos IL-4 (-)P-value
BEAS2B IL-4 (+)P-value
BEASB IL-4 (-)P-value
fold (IP/WCE in 121 bp bin)= 1 0.24 0.24 0.46 0.572 0.07 0.07 0.20 0.33 0.02 0.02 0.07 0.13 4 0.003 0.003 0.02 0.046 5 4E-04 4E-04 0.005 0.0146 5.7E-05 5.7E-05 1.0E-03 3.6E-037 6.6E-06 6.6E-06 2.0E-04 8.5E-04 8 6.8E-07 6.8E-07 3.4E-05 1.8E-049 6.3E-08 6.3E-08 5.1E-06 3.3E-0510 5.4E-09 5.4E-09 7.1E-07 5.7E-06
Ramos IL-4 (+)#total peak
Ramos IL-4 (-)#total peak
BEAS2B IL-4 (+)#total peak
BEAS2B IL-4 (-)#total peak
window size= 60 bp 1319 171 1324 491
121 bp 600 73 773 314
180 bp 259 20 504 168
C
D
(fold) (fold)
Enric
hmen
t aga
inst
bac
kgro
und
Enric
hmen
t aga
inst
bac
kgro
und
Supplementary Figure S1
False Detection Rate (fold=10) #Total Peak (121bp)
Ramos BEAS2B Ramos BEAS2BIL-4(+) IL-4(-) IL-4(+) IL-4(-) IL-4(+) IL-4(+)
H3K4me3 5.5E-07 1.7E-07 7.3E-08 3.1E-08 27364 24281H3Ac 3.1E-08 1.9E-07 1.9E-07 2.9E-08 46294 37480Pol II 8.1E-05 8.9E-05 1.6E-05 1.4E-05 17431 13457
H3K4me3 H3Ac Pol II
(fold)
H3K4me3 H3Ac Pol II
(fold)E F
G
Enric
hmen
t aga
inst
bac
kgro
und
Enric
hmen
t aga
inst
bac
kgro
und
Supplementary Figure S1.Validation of the detected binding sites by independent real time PCR in Ramos cells (A) and BEAS2B (B) cells. Fold enrichment (y-axis) was calculated based on delta Ct method against the non-binding intergenic genomic region for detected STAT6 binding site in each gene. Red and blue bars represent the results from IL-4 (-) and IL-4 (+) experiments, respectively. PCR primers used for each region are shown in (Supplementary Table S2). Estimated false detection rate (C) and number of detected binding sites of STAT6 (D), when different cutoffs were used. (E-G) Similar validation analysis for ChIP Seq analyses of H3K4me3, H3Ac and pol II. Genomic coordinate used for the analyses are shown in Supplementary Table S2. The parameters used in the text are indicated by a red box. WCE: whole cell extract genome used for the background control. Error bars are standard deviations of the triplicate experiments.
Supplementary Figure S1
- - - + - - - - + - Competitor- + + + + - + + + + Nuclear protein- - + + + - - + + + IL4 stimulation- - - - + - - - - + anti-STAT6
NM_003791 NM_022913
1 2 3 4 5 6 7 8 9 10
Super Shift
Complex
Supplementary Figure S2.STAT6 bind to the TCTCGCG sequence in BEAS2B cells. Nuclear protein was extracted from BEAS2B cells with or without IL-4.Left side was used NM_003791 (MBTPS1) probe and right side was used NM_0022913 (GPBP1) probe.
Supplementary Figure S2
A
TSSseq Fold
4.9 7.9 9.1 6.1 11.9 13.7TSSseq Fold
2.8 6.1 28.4 13.2 15.0 26.0
B
cell typeRamos IL-4 (+)
Ramos IL-4 (-)
BEAS2B IL-4 (+)
BEAS2B IL-4 (-)
total
outside of the Refseq region 2,453,479 2,764,432 501,068 546,381 16,814,181
within the Refseq region 12,815,014 13,875,727 14,317,748 11,082,366 68,573,421
upstream 3,373,426 4,121,651 3,655,518 3,215,510 22,462,370
first exon 7,129,630 7,602,531 7,831,972 6,261,221 33,317,005
second or later exon
762,470 738,366 1,680,299 962,612 6,463,399
intron 1,549,488 1,413,179 1,149,959 643,023 6,330,647
%full 94 95 88 91 81
cell typeRamos IL-4 (+)
Ramos IL-4 (-)
BEAS2B IL-4 (+)
BEAS2B IL-4 (-)
#TSC (total) 247640 213901 141051 84198
#TSC (>1ppm) 33196 32128 41359 37335
C
D
Enric
hmen
t aga
inst
bac
kgro
und
Enric
hmen
t aga
inst
bac
kgro
und
Supplementary Figure S3
si-Statsi-GAPDH
cell typesi ControlIL-4 (+)
si ControlIL-4 (-)
siSTAT6IL-4 (+)
siSTAT6IL-4 (-)
outside of the Refseq region 3,501,064 3,189,838 2,317,653 2,421,010
within the Refseq region 4,742,036 4,667,671 3,562,124 3,510,735
upstream 1,960,000 2,122,649 2,040,194 1,973,422
first exon 1,070,921 1,422,152 937,398 1,061,180
second or later exon
992,511 690,723 341,180 295,238
intron 718,604 432,147 243,352 180,895
%full 79 85 90 91
cell typesi ControlIL-4 (+)
si ControlIL-4 (-)
siSTAT6IL-4 (+)
siSTAT6IL-4 (-)
#TSC (total) 228766 125664 155016 103794
#TSC (>1ppm) 42511 45753 41509 44407
E
F
G
Enric
hmen
t aga
inst
bac
kgro
und
Supplementary Figure S3
Supplementary Figure S3.Validation Analyses of TSS Seq. Real time RT-PCR validation of the transcripts which were identified as being induced by IL-4 stimulation in Ramos cells (A) and BEAS2B (B) cells. Fold induction of the transcript levels were calculated against GAPDH by delta Ct method. Fold difference evaluated by digital TSS Seq tag counts are shown in the margin for each gene. (C) Mapped positions of the TSS Seq tags relative to the RefSeq genes are shown. %full was calculated as (#“upstream”+#”first exon”+#“intron”)/(#“total tags”) within the RefSeq region. (D) Number of TSCs having the indicated expression levels. (E) Knockdown experiments of STAT6. Knockdown efficiency (y-axis) of each gene was evaluated by real time RT-PCR using delta Ct method against non-targeting siRNA. siRNA used for each experiment is shown in the inset. Sequences of the siRNAs are shown in Supplementary Table S2. (F, G) Evaluation of the TSS Seq libraries constructed using STAT6-knowckdown BEAS2B cells. Error bars are standard deviations of the triplicate experiments.
Supplementary Figure S3
ChIP (STAT6)
Ramos IgM # total reads 80,673,092# IL-4(+) Peaks 7887
# IL-4(-) Peaks 9969
110
cell typeRamos Ramos
IgM(+) IL-4(+) IgM(+) IL-4(-)outside of the Refseq region 639,267 399,820within the Refseq region 10,719,287 6,301,263
upstream 3,013,506 1,797,419first exon 6,770,078 3,971,466second or later exon 266,685 165,318intron 669,018 367,060
%full 92 92
A B
D
cell typeRamos Ramos
IgM(+) IL-4 (+) IgM (+) IL-4 (-)#TSC (total) 84,662 74,417 #TSC (>1ppm) 16,090 17,722
C
2
ChIP positive
ChIP & TSS positive
IL4
IL4
IgMIL4
IgMIL4
Supplementary Figure S4
# TSS tag (ppm)
Ramos IgM(+) IL-4(+) 59 (5.2)
Ramos IgM(+) IL-4(-) 17 (2.5)
IL4(+) STAT6 (IP)
IL4(+)STAT6 (WCE)
ChIP Seq
TSS Seq
IgM(+)IL-4 (+)
IgM(+)IL-4 (-)
NM_016483 ( PHF7)E
Supplementary Figure S4
Supplementary Figure S4.Identification of STAT6 target in antigen-stimulated Ramos cells. (A) Summary of the ChIP Seq analyses of STAT6. The numbers of peaks detected using the described parameters are shown. (B) Mapped positions of the TSS Seq tags relative to the RefSeq genes are shown. %full was calculated as (#“upstream”+#”first exon”+#“intron”)/(#“total tags”) within the RefSeq region. (C) Number of TSCs having the indicated expression levels. (D) Comparison Ramos IL4(+) and Ramos IgM(+) IL-4(+) peaks. Upper panel was ChIP positive peaks and lower panel was ChIP & TSS positive peaks. (E) Examples of identified STAT6 targets in authentic RefSeq promoter regions in Ramos cells. STAT6 binding sites are indicated by red lines. Expression changes evaluated by digital TSS tag counts are also shown in the bottom margin. Discussion: Interestingly, we compared the targets identified from antigen-stimulated and non-stimulated conditions and found that most of them did not overlap. Particularly, we observed that 20% of ChIP Seq positive sites overlapped, though with rare exceptions, they did not induce the downstream transcriptional initiations. These results also suggested that transcriptional regulations controlled by STAT6 are complex, depending on cellular conditions. Further details will be published elsewhere.
Supplementary Figure S4
NM_016077(PTRH2)
# TSS tag (ppm)
BEAS2B IL-4(+) 1574 (103.7)
BEAS2B IL-4(-) 331 (28.4)
NM_207123(GAB1)
# TSS tag (ppm)
Ramos IL-4(+) 388 (25.4)
Ramos IL-4(-) 36 (2.2)
STAT6 (IP)
H3Ac (IP)
H3K4me3 (IP)
Pol II (IP)
STAT6 (WCE)
H3Ac (WCE)
H3K4me3 (WCE)
Pol II (WCE)
ChIP Seq
TSS Seq
IL-4 (+)
IL-4 (-)
Supplementary Figure S5
A
Supplementary Figure S5.Examples of identified STAT6 targets. (A) Other examples of the identified STAT6 targets in authentic promoter regions (left panel) and putative alternative regions (right panel) in Ramos cells. Also see the main text.
Supplementary Figure S3
Active target Silent targetRamos
IL-4 (-)
IL-4 (+)
BEAS2B
IL-4 (-)
IL-4 (+)
STAT6
pol IIactive chromatin
additional factor 1
transcriptionalinduction
silent chromatin
additional factor 2
Intermediate
STAT6
STAT6
A
Supplementary Figure S6
# TSS tag (ppm)
Ramos IL-4(+) 0 (0.0)
Ramos IL-4(-) 3 (0.1)
NM_078483 ( SLC36A1)
# TSS tag (ppm)
BEAS2B IL-4(+) 10 (0.6)
BEAS2B IL-4(-) 0 (0.0)
NM_017677 ( MTMR8)
B
STAT6 (IP)
H3Ac (IP)
H3K4me3 (IP)
Pol II (IP)
STAT6 (WCE)
H3Ac (WCE)
H3K4me3 (WCE)
Pol II (WCE)
ChIP Seq
TSS Seq
IL-4 (+)
IL-4 (-)
Supplementary Figure S6
NM_078483(Ramos)
NM_017677(BEAS2B)
IL-4(+)IL-4(-)
IL-4(+)IL-4(-)
WCEIL-4(+)
WCEIL-4(-)
TSSIL-4(+)
TSSIL-4(-)
H3K4me3IL-4(+)
H3K4me3IL-4(-)
H3AcIL4(+)
H3AcIL4(-)
Pol IIIL-4(+)
Pol IIIL-4(-)
Fold 1 1 2.9E-5 2.6E-5 497.4 148.1 17.2 20.1 10.94 44.55
WCEIL-4(+)
WCEIL-4(-)
TSSIL-4(+)
TSSIL-4(-)
H3K4me3IL-4(+)
H3K4me3IL-4(-)
H3AcIL-4(+)
H3AcIL-4(-)
Pol IIIL-4(+)
Pol IIIL-4(-)
Fold 1 1 6.4E-5 7.9E-5 3.1 43.3 3 1.2 6.3 6.1
(1) (1)(2.9E-05)
(2.6E-05)
(497.4)(148.1)
(17.2)(20.1) (10.9)
(1) (1) (6.4E-05)(7.9E-05) (3.1)
(43.3)
(3) (1.2) (6.3) (6.1)
(44.5)
C
Supplementary Figure S6
Supplementary Figure S6.Chromatin, pol II binding and transcriptional statuses around non-active STAT6 binding sites. (A) Schematic representation of the cases in which the chromosome status and the binding status of pol II are intermediate between active and silent target genes. Right and left panels in (B) exemplify “intermediate” binding sites represented by yellow and red shade boxes in (A), respectively. (C) Upper and lower panels represent results of validation analyses for histone modification patterns, pol II binding statuses and TSS inductions in right and left panels in (B), respectively. Numbers in parentheses represent fold induction compared to background (WCE). Note that we used the same primer set for the real time RT-PCR validations of TSS Seq and observed no positive signals from 1 ng of first strand cDNA used in Supplementary Figure S3.
Supplementary Figure S6
Active Target in Th2 (Ramos) Active Target in Th2 (BEAS2B)
H3Ac
H3K4me3
Pol II
H3Ac
H3K4me3
Pol II
Distance from TSS(bp)
Ave
rage
tag
coun
ts (
ppm
)
Distance from TSS(bp)
Distance from TSS(bp)
Distance from TSS(bp)
Distance from TSS(bp)
Distance from TSS(bp)
WCE IL-4(+)ChIP IL-4(-)ChIP IL-4(+)
WCE IL-4(-)WCE IL-4(+)
ChIP IL-4(-)ChIP IL-4(+)
WCE IL-4(-)
Ave
rage
tag
coun
ts (
ppm
)A
vera
ge ta
g co
unts
(pp
m)
Ave
rage
tag
coun
ts (
ppm
)A
vera
ge ta
g co
unts
(pp
m)
Ave
rage
tag
coun
ts (
ppm
)
A
Supplementary Figure S7
Target in BEAS2B Target in Th2#total 358 394Ramos H3K4me3 (+) 322 (90%) 241 (61%)Ramos H3Ac (+) 318 (89%) 236 (60%)
Ramos Pol II (+) 299 (84%) 192 (49%)
Target in Ramos Target in Th2#total 150 394BEAS2B H3K4me3 (+) 101(67%) 229 (58%)BEAS2B H3Ac (+) 90 (60%) 210 (53%)
BEAS2B Pol II (+) 93 (62%) 178 (45%)
B CChromatin/pol II status in Ramos Chromatin/pol II status in BEAS2B
Supplementary Figure S7.Results of the similar analysis as shown in Figure 3 (A) and Table 2B (B) in Th2 cells. (A) Chromatin and binding statuses of pol II in Ramos cells (left panels) and BEAS2B cells (right panels) around the STAT6 targets which are active in Th2 cells. (B) nd; not determined. Note transcriptional consequences of the STAT6 binding were not specified in the previous paper in Th2 cells.
Supplementary Figure S7
Transfac ID P-value Definition Refseq IDTSS Seq (ppm)(IL-4 (+)/IL-4 (-))
V$FXR_Q3 0.00077 FXR1 NM_005087 61.3/58.6
V$CEBPGAMMA_Q6 0.0016CCAAT/enhancer binding
protein NM_001806 3.2/0.1
V$CETS1P54_02 0.038v-ets erythroblastosis virus E26 oncogene homolog 1
NM_005238 2.0/0.5
Transfac ID P-value Definition Refseq IDTSS Seq (ppm)(IL-4 (+)/IL-4 (-))
V$FAC1_01 2.9E-05bromodomain PHD finger
transcription factor NM_004459 4.0/0.2
V$FXR_Q3 0.0025 FXR1 NM_005087 122.6/66.2
V$CEBPGAMMA_Q6 0.0025CCAAT/enhancer binding
protein NM_001806 2.0/1.9
A
B
Putative cooperative transcription factors in Ramos
Putative cooperative transcription factors in BEAS2B
Supplementary Figure S8.List of the transcription factor candidates which are likely to play roles in cell-type specific transcriptional induction in Ramos cells (A) and BEAS2B cells (B). Statistical significance of the enrichment of the consensus binding sequences in active targets is shown in the second column. Expression level evaluated by digital TSS tag counts is shown in the 5th column.
Supplementary Figure S8
Supplementary Table S1.Sequence primers used for the validation analyses.
Supplementary Table S1
Fw Primer Rv Primer
ChIP for STAT6 in Ramos NM_031479 GCAGATGGCTGGGGTATTTA CCCTACCTCCCCACTGACAC
NM_021950 GCCCTAAAAGTGAAGCCAGA AAGCAGGTGTGGATGGTAGG
NM_004015 AGTTGCGCACCTGTGTGAT TGTTTTCAGCAGCACTGGAC
NM_002002 GCCAAGACTCAGCTCAAACA CCACATCCGCACAATTTTCT
NM_021150 TCAAACAGACTGGAAACAACTCA CACAACAGGTGCCTGTCACT
NM_015226 GCTGGTTTGGCTGTTCTTTC CCTCGGAAGGAAGTGCATTA
NM_001006605 GTGTTTTTGCTGTGGCATTG TCCCATATGGCAAGGTGAAG
NM_198153 GCAGGGATTCATTGAGCCTA GCCCTCTGACTGTGAGTCCT
ChIP for STAT6 in BEAS2B NM_021102 AGCATGAGCCCAGTGAACTT TCTTGGAAGCCTGTGGTTTT
NM_001039775 AAGAGACCTGACTCTGGGCA ATTCCACCAGCTTCTCAGGA
NM_198951 ATTCCAGTCCCAGTGACCAG TGGGAAAACTGAAACCCAGA
NM_145806 CGTCATTTTGTGGCTTCAGA CTTCACTTCCTGCCACCACT
TSS in Ramos NM_021950 CCAACCAAAGAATGAAGAAG TCTGGAAAGTTCGTCTCTGT
NM_002002 TTCTGCGACCGTAAGCTG GGGTCTGGTCTTGAATCAG
NM_003745 CAGCGGAACTGCTTTTT GAAGGAGCTCAGGTAGTCG
NM_001006605 TTGGGATAATCTACCAGGTG TTCCTCTAGTTGGCTTATCAA
NM_207435 GTTGATGGGAACCATTTTC GAAGAACTTGTCTGGCTGTC
NM_004233 CTGCACAGCGTAAAGAAGAG CAAGTGAAAATGATGAGTGTT
TSS in BEAS2B NM_021102 AACAGCTACCGCTCTGAGGA CCAGGAAGAGGATCAACACC
NM_032177 AAGAGGAAAATGGGCAAGGT TCGGGCTATCAGGTCTTTCT
NM_014739 AGGAGTTTTTGCTGGGACAA CCACACCATCATCTCTGTCG
NM_016530 AGACAAGCGCAAAATCCAGT CTCCTGCTGAATTGCTGTCA
NM_198951 CTACCAGGGATCCAGCTTTG ATACTCCTGCCGCTCCTCTT
NM_017977 CCCCACGCAGACTGTGAGTA GTCCTTGTCCCACTGTGGTC
TSS in si-BEAS2B STAT6 CCTTTTGGCAGTGGTTTGAT ATCTGTGGAGAGCCATCCTG
GAPDH TCGGAGTCAACGGATTTGGT TGACGGTGCCATGGAATTTG
Ramos NM_078483 CTCAAGGGAGCACGTGACCT CTCGCCATCCTCCTTTGCT
BEAS2B NM_017677 GTTGAACCTGAACCAGAATAA CTCTCGTCACCTGATCCTC
BEAS2B EMSA NM_003791 TACTCCCAACTCTCGCGAGACTGGGCGACC GGTCGCCCAGTCTCGCGAGAGTTGGGAGTA
NM_022913 CTGTGCACGATCTCGCGAGACTTGGCCCAG CTGGGCCAAGTCTCGCGAGATCGTGCACAG