Connecting microRNA genes to the core transcriptional regulatory circuitry of embryonic stem cells - PubMed (original) (raw)
. 2008 Aug 8;134(3):521-33.
doi: 10.1016/j.cell.2008.07.020.
Stuart S Levine, Megan F Cole, Garrett M Frampton, Tobias Brambrink, Sarah Johnstone, Matthew G Guenther, Wendy K Johnston, Marius Wernig, Jamie Newman, J Mauro Calabrese, Lucas M Dennis, Thomas L Volkert, Sumeet Gupta, Jennifer Love, Nancy Hannett, Phillip A Sharp, David P Bartel, Rudolf Jaenisch, Richard A Young
Affiliations
- PMID: 18692474
- PMCID: PMC2586071
- DOI: 10.1016/j.cell.2008.07.020
Connecting microRNA genes to the core transcriptional regulatory circuitry of embryonic stem cells
Alexander Marson et al. Cell. 2008.
Abstract
MicroRNAs (miRNAs) are crucial for normal embryonic stem (ES) cell self-renewal and cellular differentiation, but how miRNA gene expression is controlled by the key transcriptional regulators of ES cells has not been established. We describe here the transcriptional regulatory circuitry of ES cells that incorporates protein-coding and miRNA genes based on high-resolution ChIP-seq data, systematic identification of miRNA promoters, and quantitative sequencing of short transcripts in multiple cell types. We find that the key ES cell transcription factors are associated with promoters for miRNAs that are preferentially expressed in ES cells and with promoters for a set of silent miRNA genes. This silent set of miRNA genes is co-occupied by Polycomb group proteins in ES cells and shows tissue-specific expression in differentiated cells. These data reveal how key ES cell transcription factors promote the ES cell miRNA expression program and integrate miRNAs into the regulatory circuitry controlling ES cell identity.
Figures
Figure 1. High-resolution genome-wide mapping of core ES cell transcription factors with ChIP-seq
A. Summary of binding data for Oct4, Sox2, Nanog and Tcf3. 14,230 sites are co-bound genome wide and mapped to either promoter proximal (TSS +/− 8kb, dark green, 27% of binding sites), genic (>8kb from TSS, middle green, 30% of binding sites), or intergenic (light green, 43% of binding sites). The promoter proximal binding sites are associated with 3,289 genes. B. (upper) Binding of Oct4 (blue), Sox2 (purple), Nanog (orange) and Tcf3 (red) across 37.5kb of mouse chromosome 3 surrounding the Sox2 gene (black box below the graph, arrow indicates transcription start site). Short sequences uniquely and perfectly mapping to the genome were extended to 200bp (maximum fragment length) and scored in 25bp bins. The score of the bins were then normalized to the total number of reads mapped. Oct4/Sox2 DNA binding motifs (Loh et al., 2006) were mapped across the genome and are shown as grey boxes below the graph. Height of the box reflects the quality of the motif. (lower) Detailed analysis of three enriched regions (Chromosome 3: 4,837,600-34,838,300, 34,845,300-34,846,000, and 34,859,900–34,860,500) at the Sox2 gene indicated with dotted boxes above. The 5’ most base from ChIP-seq were separated by strand and binned into 25bp regions. Sense (darker tone) and anti-sense (light tone) of each of the four factors tested are directed towards the binding site, which in each case occurs at a high-confidence Oct/Sox2 DNA binding motif indicated below.
Figure 2. Identification of miRNA promoters
A. Description of algorithm for miRNA promoter identification. A library of candidate transcriptional start sites was generated with histone H3 lysine 4 tri-methyl (H3K4me3) location analysis data from multiple tissues (Barski et al., 2007; Guenther et al., 2007; Mikkelsen et al., 2007). Candidates were scored to assess likelihood that they represent true miRNA promoters. Based on scores, a list of mouse and human miRNA promoters was assembled. Additional details can be found in Supplemental Text. B. Examples of identified miRNA promoter regions. A map of H3K4me3 enrichment is displayed in regions neighbouring selected human and mouse miRNAs for multiple cell types: human ES cells (hES), REH human pro-B cell line (B cell), primary human hepatocytes (Liver), primary human T cells (T cell), mouse ES cells (mES), neural precursor cells (NPCs) and mouse embryonic fibroblasts (MEFs). miRNA promoter coordinates were confirmed by distance to mature miRNA genomic sequence, conservation and EST data (shown as solid line where available). Predicted transcriptional start site and direction of transcription are noted by an arrow, with mature miRNA sequences indicated (red). CpG islands, commonly found at promoters, are indicated (green). Dotted lines denote presumed transcripts. C. Confirmation of predicted transcription start sites for active miRNAs using chromatin modifications. Normalized ChIP-seq counts for H3K4me3 (red), H3K79me2 (blue) and H3K36me3 (green) are shown for two miRNA genes where EST data was unavailable. Predicted start site (arrow), CpG islands (green bar), presumed transcript (dotted lines) and miRNA positions (red bar) are shown. D. Most human and mouse miRNA promoters show evidence of H3K4me3 enrichment in multiple tissues.
Figure 3. Oct4, Sox2, Nanog and Tcf3 occupancy and regulation of miRNA promoters
A. Oct4 (blue), Sox2 (purple), Nanog (orange) and Tcf3 (red) binding is shown at four murine miRNA genes as in Figure 1A. H3K4me3 enrichment in ES cells is indicated by shading across genomic region. Presumed transcripts are shown as dotted lines. Coordinates for the mmu-mir-290-295 cluster are derived from NCBI build 37. B. Oct4 ChIP enrichment ratios (ChIP-enriched versus total genomic DNA) are shown across human miRNA promoter region for the hsa-mir-302 cluster. H3K4me3 enrichment in ES cells is indicated by shading across genomic region. C. Schematic of miRNAs with conserved binding by the core transcription factors in ES cells. Transcription factors are represented by dark blue circles and miRNAs are represented by purple hexagons. D. Quantitative RT-PCR analysis of RNA extracted from ZHBTc4 cells in the presence or absence of doxycycline treatment. Fold change was calculated for each pri-miRNA for samples from 12 hours and 24 hours of doxycyline treatment relative to those from untreated cells. Transcript levels were normalized to Gapdh levels. Error bars indicate standard deviation derived from triplicate PCR reactions.
Figure 4. Regulation of Oct4/Sox2/Nanog/TCF3-bound miRNAs during differentiation
A. Pie charts showing relative contributions of miRNAs to the complete population of miRNAs in mES cells (red), MEFs (blue) and neural precursors (NPCs, green) based on quantification of miRNAs from by small RNA sequencing. A full list of the miRNAs identified can be found in Table S9. B. Normalized frequency of detection of individual mature miRNAs whose promoters are occupied by Oct4/Sox2/Nanog/ Tcf3 in mouse. Red line in center and right panel show the level of detection in ES cells. C. Histogram of changes in frequency of detection. Changes for miRNAs whose promoters are occupied by Oct4, Sox2, Nanog and Tcf3 in mouse are shown as bars (red for ES enriched, blue for MEF enriched and green for NPC enriched). The background frequency for non-occupied miRNAs is shown as a grey line.
Figure 5. Polycomb represses lineage-specific miRNAs in ES cells
A. Suz12 (light green) and H3K27me3 (dark green, Mikkelsen et al., 2007) binding are shown for two miRNA genes in murine ES cells. Predicted start sites (arrow), CpG islands (green bar), presumed miRNA primary transcript (dotted line) and mature miRNA (red bar) are shown. B. Expression analysis of miRNAs from mES cells based on quantitative small RNA sequencing. Cumulative distributions for Polycomb-bound miRNAs (green line) and all miRNAs (grey line) are shown. C. Expression analysis of miRNAs occupied by Suz12 in mES cells. Relative counts are shown for mES (red), NPCs (orange) and MEFs (yellow). H3K27me3 (green line) and H3K4me3 (blue line) mapped reads are shown for mES cells, MEFs and NPCs (Mikkelsen et al., 2007). D. Schematic of a subset of miRNA genes occupied by Suz12 in both mES and hES cells as in Figure 3C. miRNA genes where Oct4/Sox2/Nanog/Tcf3 are also present are indicated. Cells known to selectively express these miRNAs based on computation predictions (Farh et al., 2005) or experimental confirmation (Landgraf et al., 2007) are indicated. The Polycomb group (PcG) protein Suz12 is represented by a green circle.
Figure 6. miRNA modulation of the gene regulatory network in ES cells
A. An incoherent feed-forward motif (Alon, 2007) involving a miRNA repression of a transcription factor target gene is illustrated (left). Transcription factors are represented by dark blue circles, miRNAs in purple hexagons, protein-coding gene in pink rectangles and proteins in orange ovals. Selected instances of this network motif identified in ES cells based on data from Sinkkonen et al., 2008 or data in Figure S11 are shown (right). B. Second model of incoherent feed-forward motif (Alon 2007) involving protein repression of a miRNA is illustrated (left). In ES cells, Lin28 blocks the maturation of primary Let-7g (Visiwanthan et al., 2008). Lin28 and the Let-7g gene are occupied by Oct4/Sox2/Nanog/Tcf3. Targetscan prediction (Grimson et al., 2007), of Lin28 by mature Let-7g is noted (purple dashed line, right). C. A coherent feed-forward motif (Alon 2007) involving miRNA repression of a transcriptional repressor that regulates a transcription factor target gene is illustrated (left). This motif is found in ES cells, where mir-290-295 miRNAs repress Rbl2 indirectly maintaining the expression of Dnmt3a and Dnmt3a, which are also occupied at their promoters by Oct4/Sox2/Nanog/Tcf3 (right).
Figure 7. Multi-level regulatory network controlling ES cell identity
Updated map of ES cell regulatory circuitry is shown. Interconnected auto-regulatory loop is shown to the left. Active genes are shown at the top right, and inactive genes are shown at the bottom right. Transcription factors are represented by dark blue circles, and Suz12 by a green circle. Gene promoters are represented by red rectangles, gene products by orange circles, and miRNA promoters are represented by purple hexagons.
Similar articles
- Embryonic stem cell-specific miR302-367 cluster: human gene structure and functional characterization of its core promoter.
Barroso-delJesus A, Romero-López C, Lucena-Aguilar G, Melen GJ, Sanchez L, Ligero G, Berzal-Herranz A, Menendez P. Barroso-delJesus A, et al. Mol Cell Biol. 2008 Nov;28(21):6609-19. doi: 10.1128/MCB.00398-08. Epub 2008 Aug 25. Mol Cell Biol. 2008. PMID: 18725401 Free PMC article. - Stat3 and c-Myc genome-wide promoter occupancy in embryonic stem cells.
Kidder BL, Yang J, Palmer S. Kidder BL, et al. PLoS One. 2008;3(12):e3932. doi: 10.1371/journal.pone.0003932. Epub 2008 Dec 11. PLoS One. 2008. PMID: 19079543 Free PMC article. - Embryonic stem cell microRNAs: defining factors in induced pluripotent (iPS) and cancer (CSC) stem cells?
Gunaratne PH. Gunaratne PH. Curr Stem Cell Res Ther. 2009 Sep;4(3):168-77. doi: 10.2174/157488809789057400. Curr Stem Cell Res Ther. 2009. PMID: 19492978 Review. - Signed weighted gene co-expression network analysis of transcriptional regulation in murine embryonic stem cells.
Mason MJ, Fan G, Plath K, Zhou Q, Horvath S. Mason MJ, et al. BMC Genomics. 2009 Jul 20;10:327. doi: 10.1186/1471-2164-10-327. BMC Genomics. 2009. PMID: 19619308 Free PMC article. - The Expression and Functional Roles of miRNAs in Embryonic and Lineage-Specific Stem Cells.
Farzaneh M, Alishahi M, Derakhshan Z, Sarani NH, Attari F, Khoshnam SE. Farzaneh M, et al. Curr Stem Cell Res Ther. 2019;14(3):278-289. doi: 10.2174/1574888X14666190123162402. Curr Stem Cell Res Ther. 2019. PMID: 30674265 Review.
Cited by
- Overexpression of Let-7a mitigates diploidization in mouse androgenetic haploid embryonic stem cells.
He W, Tang H, Li Y, Wang M, Li Y, Chen J, Gao S, Han Z. He W, et al. iScience. 2024 Apr 17;27(5):109769. doi: 10.1016/j.isci.2024.109769. eCollection 2024 May 17. iScience. 2024. PMID: 38711447 Free PMC article. - NANOG controls testicular germ cell tumour stemness through regulation of MIR9-2.
Cardenas RP, Zyoud A, McIntyre A, Alberio R, Mongan NP, Allegrucci C. Cardenas RP, et al. Stem Cell Res Ther. 2024 May 1;15(1):128. doi: 10.1186/s13287-024-03724-1. Stem Cell Res Ther. 2024. PMID: 38693576 Free PMC article. - MiR-290 Family Maintains Pluripotency and Self-Renewal by Regulating MAPK Signaling Pathway in Intermediate Pluripotent Stem Cells.
Liu Y, Li X, Ma X, Du Q, Wang J, Yu H. Liu Y, et al. Int J Mol Sci. 2024 Feb 26;25(5):2681. doi: 10.3390/ijms25052681. Int J Mol Sci. 2024. PMID: 38473927 Free PMC article. - Efficient derivation of transgene-free porcine induced pluripotent stem cells enables in vitro modeling of species-specific developmental timing.
Conrad JV, Meyer S, Ramesh PS, Neira JA, Rusteika M, Mamott D, Duffin B, Bautista M, Zhang J, Hiles E, Higgins EM, Steill J, Freeman J, Ni Z, Liu S, Ungrin M, Rancourt D, Clegg DO, Stewart R, Thomson JA, Chu LF. Conrad JV, et al. Stem Cell Reports. 2023 Dec 12;18(12):2328-2343. doi: 10.1016/j.stemcr.2023.10.009. Epub 2023 Nov 9. Stem Cell Reports. 2023. PMID: 37949072 Free PMC article. - Pioneer activity distinguishes activating from non-activating SOX2 binding sites.
Maresca M, van den Brand T, Li H, Teunissen H, Davies J, de Wit E. Maresca M, et al. EMBO J. 2023 Oct 16;42(20):e113150. doi: 10.15252/embj.2022113150. Epub 2023 Sep 11. EMBO J. 2023. PMID: 37691488 Free PMC article.
References
- Alon U. Network motifs: theory and experimental approaches. Nat Rev Genet. 2007;8:450–461. - PubMed
- Barski A, Cuddapah S, Cui K, Roh TY, Schones DE, Wang Z, Wei G, Chepelev I, Zhao K. High-resolution profiling of histone methylations in the human genome. Cell. 2007;129:823–837. - PubMed
- Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004;116:281–297. - PubMed
Publication types
MeSH terms
Substances
Grants and funding
- R01 HG002668-05/HG/NHGRI NIH HHS/United States
- 5-R01-HDO45022/PHS HHS/United States
- HHMI/Howard Hughes Medical Institute/United States
- P30 CA014051/CA/NCI NIH HHS/United States
- P30-CA14051/CA/NCI NIH HHS/United States
- 5-R37-CA084198/CA/NCI NIH HHS/United States
- 5-R01-CA087869/CA/NCI NIH HHS/United States
- R37 CA084198/CA/NCI NIH HHS/United States
- R01 GM067031/GM/NIGMS NIH HHS/United States
- HG002668/HG/NHGRI NIH HHS/United States
- R01 HG002668/HG/NHGRI NIH HHS/United States
- R01 CA087869/CA/NCI NIH HHS/United States
LinkOut - more resources
Full Text Sources
Other Literature Sources