1 Cell Stem Cell, Volume 3 Supplemental Data Sall4 Regulates Distinct Transcription Circuitries in Different Blastocyst-Derived Stem Cell Lineages Chin Yan Lim, Wai-Leong Tam, Jinqiu Zhang, Haw Siang Ang, Hui Jia, Leonard Lipovich, Huck-Hui Ng, Chia-Lin Wei, Wing Kin Sung, Paul Robson, Henry Yang, and Bing Lim Experimental Procedures RNAi design and construction of plasmids for shRNA synthesis. For RNAi validation of gene function, two short hairpin RNA constructs were made for each targeted gene. These oligonucleotides were annealed and cloned into pSuper.puro (Oligoengine). Gata6 shRNA 1 GATCCCCGAGGACCTGTTGCTCTTCAttcaagagaTGAAGAGCAACAGGTCCTCTTTTTA AGCTTAAAAAGAGGACCTGTTGCTCTTCAtctcttgaaTGAAGAGCAACAGGTCCTCGGG Gata 6 shRNA 2 gatccccAtGCAGACAtAACAttCCTttcaagagaAGGAAtGttAtGtCtGCATttttta agcttaaaaaAtGCAGACAtAACAttCCTtctcttgaaAGGAAtGttAtGtCtGCATggg Sox17 shRNA 1 GATCCCCGGACCCGGCTTTCTTTGCAttcaagagaTGCAAAGAAAGCCGGGTCCTTTTTA AGCTTAAAAAGGACCCGGCTTTCTTTGCAtctcttgaaTGCAAAGAAAGCCGGGTCCGGG Sox17 shRNA 2 GATCCCCGCACGGAATTCGAACAGTAttcaagagaTACTGTTCGAATTCCGTGCTTTTTA AGCTTAAAAAGCACGGAATTCGAACAGTAtctcttgaaTACTGTTCGAATTCCGTGCGGG Sox 7 shRNA 1 GATCCCCGAACACGCTGCCTGAGAAAttcaagagaTTTCTCAGGCAGCGTGTTCTTTTTA AGCTTAAAAAGAACACGCTGCCTGAGAAAtctcttgaaTTTCTCAGGCAGCGTGTTCGGG Sox 7 shRNA 2 GATCCCCGTCCTCTGTTTCTGGATGTttcaagagaACATCCAGAAACAGAGGACTTTTTA AGCTTAAAAAGTCCTCTGTTTCTGGATGTtctcttgaaACATCCAGAAACAGAGGACGGG Ezh2 shRNA 1 GATCCCCGGATGGCACTTTCATTGAAttcaagagaTTCAATGAAAGTGCCATCCTTTTTA AGCTTAAAAAGGATGGCACTTTCATTGAAtctcttgaaTTCAATGAAAGTGCCATCCGGG Ezh2 shRNA 3 GATCCCCGCAAATTCTCGGTGTCAAAttcaagagaTTTGACACCGAGAATTTGCTTTTTA AGCTTAAAAAGCAAATTCTCGGTGTCAAAtctcttgaaTTTGACACCGAGAATTTGCGGG Slc25a36 shRNA 1
20
Embed
Supplemental Data Sall4 Regulates Distinct Transcription Circuitries … · Sall4 Regulates Distinct Transcription Circuitries in Different Blastocyst-Derived Stem ... Chia-Lin Wei,
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
in Different Blastocyst-Derived Stem Cell Lineages Chin Yan Lim, Wai-Leong Tam, Jinqiu Zhang, Haw Siang Ang, Hui Jia, Leonard Lipovich, Huck-Hui Ng, Chia-Lin Wei, Wing Kin Sung, Paul Robson, Henry Yang, and Bing Lim Experimental Procedures RNAi design and construction of plasmids for shRNA synthesis. For RNAi validation of gene function, two short hairpin RNA constructs were made for each targeted gene. These oligonucleotides were annealed and cloned into pSuper.puro (Oligoengine). Gata6 shRNA 1 GATCCCCGAGGACCTGTTGCTCTTCAttcaagagaTGAAGAGCAACAGGTCCTCTTTTTA AGCTTAAAAAGAGGACCTGTTGCTCTTCAtctcttgaaTGAAGAGCAACAGGTCCTCGGG Gata 6 shRNA 2 gatccccAtGCAGACAtAACAttCCTttcaagagaAGGAAtGttAtGtCtGCATttttta agcttaaaaaAtGCAGACAtAACAttCCTtctcttgaaAGGAAtGttAtGtCtGCATggg Sox17 shRNA 1 GATCCCCGGACCCGGCTTTCTTTGCAttcaagagaTGCAAAGAAAGCCGGGTCCTTTTTA AGCTTAAAAAGGACCCGGCTTTCTTTGCAtctcttgaaTGCAAAGAAAGCCGGGTCCGGG Sox17 shRNA 2 GATCCCCGCACGGAATTCGAACAGTAttcaagagaTACTGTTCGAATTCCGTGCTTTTTA AGCTTAAAAAGCACGGAATTCGAACAGTAtctcttgaaTACTGTTCGAATTCCGTGCGGG Sox 7 shRNA 1 GATCCCCGAACACGCTGCCTGAGAAAttcaagagaTTTCTCAGGCAGCGTGTTCTTTTTA AGCTTAAAAAGAACACGCTGCCTGAGAAAtctcttgaaTTTCTCAGGCAGCGTGTTCGGG Sox 7 shRNA 2 GATCCCCGTCCTCTGTTTCTGGATGTttcaagagaACATCCAGAAACAGAGGACTTTTTA AGCTTAAAAAGTCCTCTGTTTCTGGATGTtctcttgaaACATCCAGAAACAGAGGACGGG Ezh2 shRNA 1 GATCCCCGGATGGCACTTTCATTGAAttcaagagaTTCAATGAAAGTGCCATCCTTTTTA AGCTTAAAAAGGATGGCACTTTCATTGAAtctcttgaaTTCAATGAAAGTGCCATCCGGG Ezh2 shRNA 3 GATCCCCGCAAATTCTCGGTGTCAAAttcaagagaTTTGACACCGAGAATTTGCTTTTTA AGCTTAAAAAGCAAATTCTCGGTGTCAAAtctcttgaaTTTGACACCGAGAATTTGCGGG Slc25a36 shRNA 1
buckleyc
Cross-Out
2
GATCCCCGGCATATCAGAGACTGTTAttcaagagaTAACAGTCTCTGATATGCCTTTTTA AGCTTAAAAAGGCATATCAGAGACTGTTAtctcttgaaTAACAGTCTCTGATATGCCGGG Slc25a36 shRNA 1 GATCCCCGCAGTCTTCTTCTGTGACAttcaagagaTGTCACAGAAGAAGACTGCTTTTTA AGCTTAAAAAGCAGTCTTCTTCTGTGACAtctcttgaaTGTCACAGAAGAAGACTGCGGG Zfp27 shRNA 1 GATCCCCGGACACGTCTTATGTAATTttcaagagaAATTACATAAGACGTGTCCTTTTTA AGCTTAAAAAGGACACGTCTTATGTAATTtctcttgaaAATTACATAAGACGTGTCCGGG Zfp27 shRNA 2 gatccccttACtGtCtCtAtAGGACAttcaagagatGtCCtAtAGAGACAGtAAttttta agcttaaaaattACtGtCtCtAtAGGACAtctcttgaatGtCCtAtAGAGACAGtAAggg Control shRNA GATCCCCGAACGGCATCAAGGTGAACttcaagagaGTTCACCTTGATGCCGTTCTTTTTA AGCTTAAAAAGAACGGCATCAAGGTGAACtctcttgaaGTTCACCTTGATGCCGTTCGGG Chromatin Immunoprecipitation and DNA Microarray Analysis. ChIP assays with feeder-free E14 mouse ESCs were carried out as described previously (Zhang et al., 2006). Briefly, cells were cross-linked with 1 % formaldehyde for 10 min at room temperature and formaldehyde was then inactivated by the addition of 125 mM glycine. Chromatin extracts containing DNA fragments with an average size of 500 bp were immunoprecipitated using antibodies. Approximately 5x107 cells were used for each ChIP reaction. Amounts of antibodies used for immunoprecipitation were 10 ug for anti-Sall4, 10 ug for anti-H3K4me3 (Abcam) and 15 ug for anti-H3K27me3 (Upstate biotechnology). For all ChIP experiments, quantitative PCR analyses were performed in real-time using the ABI PRISM 7900 Sequence Detection System and SYBR Green Master Mix as described previously (Zhang et al., 2006). Relative occupancy values were calculated by determining the apparent immunoprecipitation efficiency (ratios of the amount of immunoprecipitated DNA to that of the input sample) and normalized to the level observed. ChIP combined with DNA microarray was carried out according to Agilent Mammalian ChIP-on-chip protocol (version 3) (Boyer et al., 2005). Briefly, immunoprecipitated (IP) DNA was blunted with T4 polymerase and ligated to linkers with T4 DNA ligase. For reference, 200 ng of input DNA was used. Ligated DNA was amplified with two rounds of PCR cycles. Two µg of each amplified sample were labeled with Cy5 for IP DNA, and Cy3 for input DNA using CGH labeling kit (Invitrogen). Four µg of each labeled sample was blocked with mouse Cot-1 DNA (Invitogen) and blocking agent (Agilent), and hybridized to DNA microarrays (G4490A, Agilent) for 40 h at 65 oC. The microarrays were washed and stabilized with acetonitrile (Sigma), followed by stabilization and drying solution (Agilent). The arrays were scanned with Agilent DNA microarray scanner (G2565BA), and probe features were extracted using Feature Extraction Software (G2567AA, version 9.1). The data was then analysed with ChIP Analytics 1.3 (G4477AA), according to the Whitehead Neighbourhood Model, P(Xbar)<0.001, with intra-array median normalization (Boyer et al., 2005). GEO: GSE12482 (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE12482). Selected
3
bound targets, determined by the algorithm, were validated by qPCR to determine the false discovery rate. Quantitative RT-PCR - List of ABI Taqman probes Mapped ID Assay ID Actb Mm00607939_s1 Afp Mm00431715_m1 Bmp4 Mm00432087_m1 Cdx2 Mm00432449_m1 Eomes Mm01351984_m1 Esx1 Mm01297208_m1 Fgf5 Mm00438919_m1 Foxa2 Mm00839704_mH Gata4 Mm00484689_m1 Gata6 Mm00802636_m1 Hand1 Mm00433931_m1 Hnf4a Mm00433964_m1 Ihh Mm00439613_m1 MyoD1 Mm00440387_m1 Nanog Mm02019550_s1 Otx2 Mm00446859_m1 Pax6 Mm00443072_m1 Pdgfra Mm01211694_m1 Pou5f1 Mm00658129_gH Psx1/Rhox6 Mm00655990_g1 Pthr1 Mm00441046_m1 Sall4 Mm00614351_m1 Sox1 Mm00486299_s1 Sox17 Mm00488363_m1 Sox2 Mm00488369_s1 Sox7 Mm00776876_m1 T Mm00436877_m1 Tcf2/Hnf1b Mm00447452_m1 Utf1 Mm00447703_g1 Promoter reporter assays. The promoters of selected genes were PCR amplified from mouse genomic DNA and cloned into pGL4.10 (Promega). Primers used for cloning are listed in table. The Oct4 promoter was previous described in Zhang et al., 2006. To perform the promoter-luciferase assays, HEK293T cells were seeded 24 h before transfection at a density of 2.5 × 104 cells per well in 96-well culture plates. For testing the effect of Sall4 activity on these promoters, Sall4 over-expressing plasmid, pCAG-Sall4-IRES-GFP, or pCAG-IRES-GFP vector (200 ng) was co-transfected with individual promoter-luciferase (50ng) and pRL-SV40 (1ng) (Promega). Firefly and Renilla luciferase activities were measured with the Dual-Luciferase Reporter system (Promega) 48h post-transfection on Centro LB960 Luminometer (Berthold Technologies). The data were expressed as relative to corresponding control vector transfections, after normalization to Renilla luciferase readings.
4
Primers for promoter reporter constructs Gene Promoter
Intersection of Sall4 ChIP-Seq data with MTL in ESCs. To determine the common target (co-regulated) genes of Oct4, Sox2, Nanog and Sall4, co-localization of the binding loci in the neighborhood of a target gene was analyzed. Firstly, the co-localized binding loci of all these four factors were extracted and mapped to genes using CEAS (Ji et al., 2006) with the mouse genome assembly version Build 36 (mm8). Sall4 binding loci generated from Sall4 ChIP-seq data with enrichment ≥ 9 intensity reads were used to intersect with the 3582 MTL obtained from Chen et al., 2008. Amongst these 3582 MTL, clustering of 14 transcription factors (13 TFs + Sall4) was performed based on co-occurrences of the binding loci to evaluate the similarity of the TF targeting. The co-occurrences between 14 TFs within the MTL were first computed and the correlation coefficients between each pair of co-occurrence vector were then calculated and modeled to a 14×14 correlation matrix. With the matrix, a heatmap reflecting the hierarchical clustering of the correlation coefficients was generated. De novo motif analysis. To determine the binding motif of Sall4, we extracted sequences (+/- 200) for all (132) ChIP-seq peaks in the Sall4 library containing at least 25 overlapping extended tags. An enriched motif (TCGCCATA) was identified in the 132 sequences using the motif discovery algorithm Weeder (Pavesi et al., 2001). Supplemental References Boyer, L. A., Lee, T. I., Cole, M. F., Johnstone, S. E., Levine, S. S., Zucker, J. P., Guenther, M. G., Kumar, R. M., Murray, H. L., Jenner, R. G., et al. (2005). Core transcriptional regulatory circuitry in human embryonic stem cells. Cell 122, 947-956. Chen, X., Xu, H., Yuan, P., Fang, F., Huss, M., Vega, V. B., Wong, E., Orlov, Y. L., Zhang, W., Jiang, J., et al. (2008). Integration of external signaling pathways with the core transcriptional network in embryonic stem cells. Cell 133, 1106-1117. Ji, X., Li, W., Song, J., Wei, L., and Liu, X. S. (2006). CEAS: cis-regulatory element annotation system. Nucleic Acids Res 34, W551-554. Pavesi, G., Mauri, G., and Pesole, G. (2001). An algorithm for finding signals of unknown length in DNA sequences. Bioinformatics 17 Suppl 1, S207-214.
6
Zhang, J., Tam, W. L., Tong, G. Q., Wu, Q., Chan, H. Y., Soh, B. S., Lou, Y., Yang, J., Ma, Y., Chai, L., et al. (2006). Sall4 modulates embryonic stem cell pluripotency and early embryonic development by the transcriptional regulation of Pou5f1. Nat Cell Biol 8, 1114-1123.
Figure S1. Expression of XEN marker genes in blastocyst-derived XEN cells. Immunofluorescence analysis of Gata6, Gata4, Sox7 and Sox17 expression and nuclear localization in XEN cells. Scale bars = 50 µm.
7
Figure S2. Blastocyst-derived XEN-GIS cell line is highly similar to previously characterized XEN cells. (A) Dendrogram showing hierarchical clustering of the transcriptome profile of the XEN-GIS line with that of three characterized XEN lines XEN1-3, IM8A1-I, IM8A-II, four ES cell lines and four mouse embryonic fibroblasts (MEF) lines. (B) Heatmap representation of gene expression profiles in the indicated XEN and ES cell lines. (C) Marker gene expression changes in differentiated XEN-GIS cells analyzed by qRT-PCR. Values are normalized to β-actin and plotted relative to expression levels in undifferentiated XEN-GIS cells, except Afp (*) which was not expressed in the undifferentiated cells. For Afp, an arbitrary CT value at the threshold of real-time PCR detection was used in the calculations. Error bars indicate s.e.m of four replicates.
8
Figure S3A. Validation of ES-specific Sall4 ChIP-chip binding sites in ES cells by qPCR. (A) Specific enrichment of Sall4 at 109 ES-specific loci in ES cells was determined by ChIP using anti-Sall4 or anti-GST control antibodies. Fold enrichment was calculated as the ratio of signals at Sall4-bound locus to the mean of three negative control genomic regions. (B) ChIP-qPCR analysis of the ES-specific Sall4 binding sites in XEN cells. Low or no enrichment of Sall4 was observed. Error bars indicate standard deviation of two technical replicates.
9
Figure S3B. Validation of XEN-specific Sall4 ChIP-chip binding sites in XEN cells by qPCR. (A) Specific enrichment of Sall4 at 118 XEN-specific loci in ES cells was determined by ChIP using anti-Sall4 or anti-GST control antibodies. Fold enrichment was calculated as the ratio of signals at Sall4-bound locus to the mean of three negative control genomic regions. (B) ChIP-qPCR analysis of the XEN-specific Sall4 binding sites in ES cells. Low or no enrichment of Sall4 was observed. Error bars indicate standard deviation of two technical replicates.
10
Figure S3C. Validation of Sall4 ChIP-chip binding sites common to both ES and XEN cells by qPCR. (A) Enrichment of Sall4 at 28 loci in ES cells was determined by ChIP using anti-Sall4 or anti-GST control antibodies. Fold enrichment was calculated as the ratio of signals at Sall4-bound locus to the mean of three negative control genomic regions. (B) ChIP-qPCR analysis of the same 28 Sall4 binding sites in XEN cells. Error bars indicate standard deviation of two technical replicates.
11
Figure S4. Sall4 transcriptionally activates the promoters of its target genes. Promoter reporter constructs of Sall4 target genes were co-transfected with Sall4 over-expression or control vector into HEK293T cells. Luciferase reporter gene activities were measured 48hrs post-transfection. Reporter activity with Sall4 over-expression are plotted relative to control vector. The error bars represent standard error of three independent experiments, each performed in triplicates. Primers used in cloning these promoters are available in Supplemental Information.
12
Figure S5. Changes in lineage marker gene expression upon knockdown of Sall4 targets in ES cells by RNAi. Quantitative RT-PCR analysis of marker gene expression was performed four days after the cells were transfected with gene-specific shRNA or control shRNA plasmids and cultured under puromycin selection. All values were normalized to levels of β-actin, and plotted relative to expression levels in control shRNA-transfected cells. Error bars indicate standard error of three replicates.
13
Figure S6. RNAi knockdown of Slc25a36, Zfp27 and Ezh2 in ES cells leads to the emergence of Cdx2-positive cells. Gene-specific or control shRNA plasmids were transfected into ES cells cultured in trophoblast stem cell medium to promote differentiation into trophoblast stem cells. Cells were fixed and immunostained for Cdx2 three days after transfection. Scale bars = 50 µm.
14
Figure S7. Analysis of Sall4 ChIP-seq data. (A) Heatmap showing the cluster analysis of the co-occurrence of Sall4 with other transcription factors within multiple transcription factor-binding loci (MTL) in ES cells. The factors were clustered along both axes based on the degree of colocalization with other factors. Colors reflect the colocalization frequency of each pair of factors in MTL; yellow indicates more frequent colocalization, red indicate less colocalization. (B) Putative Sall4 binding motif identified by de novo sequence analysis of Sall4 ChIP-seq data in ES cells. The TCGCCATA sequence was found to be enriched in 132 high-intensity Sall4 binding peak sequences. Methodology for the analysis is available in Supplemental Information.
15
Table S1. Expression of known genes in different XEN cell lines Expression levels (Affymetrix) (Illumina) Gene Symbol Unigene ID XEN1-3 IM8A1-I IM8A1-II XEN-GIS
Correlation between XEN
lines
PE and VE genes Sox7 Sox17 Tcf2 Stra6 Gata6 Gata4 Krt2-8 Dab2 Vegfa Krt1-18 PE genes Pthr1 Sparc Pdgfra Snai1 Plat Thbd Fst VE genes Ihh Acvr1 Furin Cited1 Hnf4a Foxa2 Afp Cyp26a1 Plau ES and TS genes Pou5f1 Nanog Zfp42 Eomes Cdx2
2344 4528 235 346 3903 1455 ND 1400 1696 4650 538 2134 2739 ND 3873 178 1029 183 519 178 497 1110 ND A 90 A A A A 97 A
Good Good Good Good Good Good N/A
Good Good Good
Good Good Good N/A
Good Good Good
Good No
Good Good No* N/A
Good No
Good
Good Good Good Good Good
Correlation in the expression of a subset of known genes in XEN-GIS and three characterized XEN lines, XEN1-3, IM8A-I, IM8A-II (Kunath et al. 2005). Microarray expression profiling of XEN-GIS was performed on Illumina Beadarrays, while the Affymetrix microarray data of XEN1-3, IM8A1-I and IM8A1-II were obtained from Kunath et al. 2005. “A” indicates absence of gene expression and “ND” indicates genes not detected on the microarray. * Expression of Hnf4a was detected by RT-PCR in XEN1-3, IM8A1-I and IM8A1-II (Kunath et al. 2005).
16
Table S2. Expression of genes in XEN cells upon Sall4 knockdown Log2 Fold difference
(compared to control KD)
Gene Symbol Unigene ID TargetID Sall4 KD 3 Days
Sall4 KD 4 Days
Sall4 KD 5 Days
PE and VE genes Sox7 Sox17 Tcf2 Stra6 Gata6 Gata4 Krt2-8 Dab2 Vegfa Krt1-18 PE genes Pthr1 Sparc Pdgfra Snai1 Plat Thbd Fst VE genes Ihh Acvr1 Furin Cited1 Hnf4a Foxa2 Afp Cyp26a1 Plau Top 40 XEN expressed genes (compared to expression in ES) Glipr1 Foxq1 Clic6 Ctsh Hs3st1 Plat Gata6 Hprt Pdgfra Ankrd1 Sox17
Changes in expression of a subset of genes after RNAi-mediated Sall4 knockdown in XEN cells. The log2 fold difference in gene expression between Sall4 knockdown and control knockdown XEN cells after RNAi treatment for 3, 4, and 5 days are shown for genes known to be expressed in the extra-embroynic endoderm lineages, and for genes with the highest differential expression in XEN compared to ES cells. The log2 fold difference in gene expression is represented as ND (not determined) for genes that are not expressed. The log2 fold difference for each time point was obtained from the average of triplicate experiments. Tables S3 and S4 are available as Excel files.
18
Table S5. Sall4 target genes with highly differential expression in Sall4 knockdown ES cells.
The list of Sall4-bound genes that exhibited >1.75 fold change in expression upon RNAi-mediated Sall4 knockdown in ES cells. Expression microarray profiling was performed on triplicate sets of ES cells transfected with Sall4- or control- shRNA and cultured in puromycin drug selection for four days. Average fold changes in gene expression between Sall4 knockdown and control knockdown cells were determined for each Sall4-bound gene.
Gene Symbol Fold Change (Downregulated) Gene Symbol Fold Change
Table S6. Sall4 target genes with highly differential expression in Sall4 knockdown XEN cells.
The list of Sall4-bound genes that exhibited >1.75 fold change in expression upon RNAi-mediated Sall4 knockdown in XEN cells. Expression microarray profiling was performed on triplicate sets of XEN cells transfected with Sall4- or control- shRNA and cultured under puromycin drug selection for five days. Average fold changes in gene expression between Sall4 knockdown and control knockdown cells were determined for each Sall4-bound gene.