- Open Access
RISC-mediated control of selected chromatin regulators stabilizes ground state pluripotency of mouse embryonic stem cells
© Pandolfini et al. 2016
- Received: 9 February 2016
- Accepted: 13 April 2016
- Published: 6 May 2016
Embryonic stem cells are intrinsically unstable and differentiate spontaneously if they are not shielded from external stimuli. Although the nature of such instability is still controversial, growing evidence suggests that protein translation control may play a crucial role.
We performed an integrated analysis of RNA and proteins at the transition between naïve embryonic stem cells and cells primed to differentiate. During this transition, mRNAs coding for chromatin regulators are specifically released from translational inhibition mediated by RNA-induced silencing complex (RISC). This suggests that, prior to differentiation, the propensity of embryonic stem cells to change their epigenetic status is hampered by RNA interference. The expression of these chromatin regulators is reinstated following acute inactivation of RISC and it correlates with loss of stemness markers and activation of early cell differentiation markers in treated embryonic stem cells.
We propose that RISC-mediated inhibition of specific sets of chromatin regulators is a primary mechanism for preserving embryonic stem cell pluripotency while inhibiting the onset of embryonic developmental programs.
- Green Fluorescent Protein
- Embryonic Stem
- Embryonic Stem Cell
- Embryonic Stem Cell Line
- Chromatin Regulator
Embryonic stem (ES) cells tend to spontaneously differentiate in the absence of external inductive signals . The first step of ES cell differentiation, commonly reported as “priming”, is mostly associated with changes in the dynamics of chromatin, post-translational modifications of histones, and a general remodeling of nuclear architecture . Priming is considered necessary for lineage specification in the early embryo but the exact mechanisms mediating its action on the transition from pluripotency state to the differentiation of embryonic tissues are not understood. Inhibition of protein translational noise  and transcriptional “leakage” [4, 5] characterize mouse ES cells. This indicates that lineage specification during early embryonic development could be driven by reduction of the transcribed portion of the genome but it also poses the question of how pluripotency can accommodate the transcription of tissue-specific genes.
We speculated that a tight inhibitory control of translation is crucial to maintain pluripotency and that inhibition of protein translation through microRNA (miRNA) and the RNA-induced silencing complex (RISC)  might represent one strategy to avoid a “transcriptional paradox”. There is, indeed, an established body of evidence indicating that release from RISC-mediated translational inhibition, produced through the disruption of components of the miRNA maturation pathway such as Dicer  or DGCR8 , severely impairs pluripotency in ES cells. This observation implies that inhibition of protein translation is necessary for pluripotency. However, while the general involvement of RISC is established, little is known about the families of genes subject to this control. In our investigation, we found that a set of mRNAs encoding chromatin regulators is selectively released from miRNA-mediated protein translation inhibition during priming and we conclude that their inhibition is crucial for the maintenance of ground state pluripotency.
Epiblast-like aggregate cells are equivalent to primed pluripotent cells
Mouse EpiSC can be maintained in a pluripotent state similar to that of human ES cells (hESCs) [17, 18] by FGF2/Activin A treatment [12, 13]. Ground-state mESC cultured for several passages in FGF2/Activin-containing medium (Fig. 2f) acquire a flattened morphology (Fig. 2h) and retain a gene expression profile similar to ELA cells and consistent with epiblast identity [15, 19]. Accordingly, ELA-derived cells could be dissociated and maintained in an EpiSC-like state by adding FGF2/Activin A to culture medium (Fig. 2g). The resulting cells (ELA-EpiSC; Fig. 2i) displayed a high level of correlation with the gene expression profile of mESC-derived EpiSC (Fig. 2j, k).
As well as ground-state mESC, pluripotent stem cells derived from the epiblast spontaneously neuralize when deprived of FGF2 and Activin A. Similarly, more than 90 % of ELA-EpiSC induced the early marker of neuralization Sox1 72 h after FGF2/Activin A withdrawal (Fig. 2l). These cells became nestin-positive neural progenitor cells (Fig. 2m) and differentiating neurons positive to N-tubulin and Pax6 (Fig. 2n) by down-regulating Oct4 and Nanog and up-regulating neuronal differentiation markers (Nestin, β-III-tubulin and NFL; Fig. 2o).
These results thus indicate that ELA cells can be considered to be equivalent to epiblast stem cells and, therefore, the ES–ELA transition could be used as an in vitro model of the transition from inner mass cells to primed post-implantation epiblast cells. Using ELA instead of EpiSC allows the requirement for exogenous FGF and Activin to be avoided, whose action could not parallel the autocrine signaling of the developing embryo. Moreover, since ELA cells are not steady-state cultures as EpiSC are, they retain the transient nature of epiblast cells in vitro and provide a better model than EpiSC for studying spontaneous, rapidly occurring molecular changes during this developmental transition.
Profiling of Argonaute occupancy at early steps of ES cell differentiation
We performed a global survey of the mRNAs that interact with Argonaute (Ago) during the first stages of in vitro differentiation. Towards this goal, we isolated Ago-interacting RNAs by immunoprecipitating Ago proteins (Ago* in Additional file 1: Figure S1a), which are the main RISC components . Ago–RNA and total RNA from cells at 0, 2, 6, and 10 DIV were hybridized to gene expression microarrays (Fig 2a). The release of mRNA from Ago at the transition between sequential steps was evaluated as the variation of the ratio between Ago–RNA and total RNA levels after normalization (Ago Enrichment in Additional file 1: Figure S1b, c; Additional file 2: Table S1; see “Methods”). While most of the transcriptional regulation occurred at the ELA–NPC transition (Fig. 2b), the highest release of mRNA by Ago was observed at the ES–ELA transition (Fig. 2c), with more than 100 mRNA species significantly released from Ago binding. We verified that this could not be explained by either different efficiency in Ago immunoprecipitation (Additional file 1: Figure S1a) or differential expression of Ago- or miRNA-related genes in the two conditions (Additional file 1: Figure S1d). The variation in Ago enrichment was coherently correlated with the variation in total RNA level at the ES–ELA transition for 82.4 % of genes (Additional file 1: Figure S1e), confirming a predominant effect of RISC on mRNA destabilization .
mRNAs of selected chromatin regulator families are preferentially released by Ago during the ES–ELA transition
Chromatin modifications are postulated to play a crucial role in priming [22–24] but miRNA-mediated translational control of chromatin regulators has never been extensively described. Therefore, we focused our attention on the release by Ago of distinct families of chromatin modifiers  (Additional file 3: Table S2).
While the families of mRNAs coding for de novo DNA methyltransferases (DNMT), histone lysine demethylases (KDM), and SWItch/Sucrose NonFermentable nucleosome remodeling complex (SWI/SNF) showed a more than twofold increase of Ago enrichment in ES cells compared with ELA cells, other families of chromatin regulators showed little or no increase, similar to what was observed for families of mRNAs with unrelated functions such as cytochrome P450 (CYP450; Fig. 3c).
Ago enrichment is predictive of protein translation activity
To detect changes in protein translation at the ES–ELA transition, we evaluated ribosome occupancy by means of polysome profiling [26, 27] (see “Methods”; Additional file 1: Figure S2a, b), adapting the translating ribosome affinity purification (TRAP) protocol for our purposes. We measured the amount of mRNAs loaded on translating ribosomes and calculated ribosome enrichment as an estimate of ribosome occupancy. Our analysis showed that the GO categories of genes showing increased ribosome loading during the ES–ELA transition were almost superimposable on those found for genes released from Ago (Fig. 3a, b). Strikingly, DNMT, KDM, and SWI/SNF families showed a ribosome enrichment essentially reciprocal to their Ago enrichment (Fig. 3d).
DNMT, KDM, and SWI/SNF activities are necessary to switch between ground and primed pluripotent cell states
We investigated the role of DNMT in silencing Nanog expression more directly by monitoring Nanog promoter methylation, which was consistently inhibited by AZA treatment at 2 DIV compared with control (Fig. 5g). Interestingly, CHQ treatment and shSA4 transduction blocked the methylation of the Nanog promoter as well (Fig. 5g), indicating the requirement of KDM and SWI/SNF activities in this modification, in agreement with previous reports . At 2 DIV, cells treated with AZA and CHQ showed a global mRNA expression profile almost identical to that of ES cells, as evaluated by hierarchical clustering and fold change analysis (Additional file 1: Figure S3h–n). Cells transduced with shSA4 showed a modulation of their profile consistent with that of AZA and CHQ (Additional file 1: Figure S3l) but, in this case, gene expression clustered with that of ELA cells (Additional file 1: Figure S3h), suggesting that the combined action of several SWI/SNF members might be necessary to control the ES–ELA cell transition.
RISC inactivation de-represses the translation of selected chromatin regulators and destabilizes pluripotency
Ground-state miRNAs affect the levels of chromatin regulators
To get a better insight into the mechanisms of RISC targeting, we performed small RNA-seq of cells during the ES–ELA transition. This experiment had the additional advantage of validating the robustness of Ago-RIP as a method to measure RISC-mediated inhibition of mRNAs. Indeed, we observed that miRNAs predicted to target the most enriched or depleted RNAs in the Ago-RIP were accordingly up- or down-regulated during the ES–ELA transition (Fig. 7c, d). This result suggests that our method of analysis of mRNA/Ago enrichment is indeed predictive of miRNA–mRNA interactions.
More interestingly, the analysis of global miRNA expression at the ES–ELA cell transition highlighted a number of up- and down-regulated miRNAs (Additional file 1: Figure S4d; Additional file 5: Table S5), consistent with the differential pattern of expression observed in naïve and primed pluripotent cells [34, 35]. Specifically, the level of expression of these miRNAs is higher in naïve pluripotent cells and lower in ELA generated from ES cells (Fig. 7e), making them part of a signature of ground-state pluripotency. We found that the 3′ untranslated regions (UTRs) of genes coding for DNMT, KDM, and SWI/SNF are significantly enriched in predicted binding sites for the miRNAs down-regulated at the ES–ELA cell transition compared with control 3′ UTRs (Fig. 7f). High-affinity interactions predicted between 43 miRNAs most changing at the ES–ELA cell transition and distinct members of DNMT, KDM, and SWI/SNF are reported in Additional file 5: Table S5. Altogether, our observations suggest a direct role of miRNAs in inhibiting the translation of DNMT, KDM, and SWI/SNF in ES cells.
Since the high number of predicted miRNA–mRNA interactions impeded their comprehensive experimental validation, we focused on the effects of the miRNAs showing the strongest predicted interaction with Kdm2b (mmu-miR-150-5p and mmu-miR-27b-5p), SmarcA4 (mmu-miR-181c-5p and mmu-miR-425-5p), and Dnmt3b (mmu-miR-295-5p) (Additional file 5: Table S5). To evaluate miRNA–mRNA interactions, we co-transfected mature miRNAs together with reporters carrying enhanced GFP (EGFP) followed by the 3′ UTR of Kdm2b, SmarcA4 or Dnmt3b in ES cells. We then indirectly evaluated the amount of EGFP produced after 48 h in ELA cells or in ES cells maintained in ES cell medium, measuring the ratio between EGFP fluorescence and the fluorescence produced by a DsRed-carrying plasmid (Fig. 7g–i; see “Methods”). EGFP of 3′ UTR-carrying vectors was always expressed at lower levels than EGFP of empty vector. In the presence of Kdm2b, SmarcA4, or Dnmt3b 3′ UTRs, however, EGFP levels were significantly lower in ES cells than in ELA cells, indicating that a translational switch driven by these 3′ UTRs occurs at the ES–ELA transition. Moreover, all the miRNAs expected to bind to the 3′ UTRs with the best efficiency were able to inhibit EGFP expression from a 3′ UTR-carrying vector in ELA cells. According to the prediction of multiple miRNAs interacting with one 3′ UTR, the co-transfection of one single miRNA was not able to inhibit EGFP expression to the levels observed in ES cells. Finally, EGFP expression was less inhibited by miRNA co-transfection in cells maintained in ES medium (ES cells) than in cells differentiated in chemically defined minimal medium (CDMM; ELA cells), in agreement with the evidence that higher amounts of these miRNAs are already present in ES cells and, therefore, their transfection is expected to show a smaller additional effect. These results support the specificity of some of the mRNA–miRNA interactions at the ES–ELA cell transition predicted by in silico tools.
Taking advantage of the spontaneous differentiation program that ES cells undergo when deprived of pluripotency-maintaining signals, we isolated and studied cells primed to differentiation (ELA) just after their escape from the ground-state pluripotency state.
We performed Ago-RIP as a way to estimate RISC-mediated inhibition of gene expression. Two new nuclear roles for Ago2 have been recently found in addition to its main function as a component of RISC: one in pre-mRNA splicing and one in transcriptional repression . We tend to exclude an influence of nuclear RNA in our analysis as Ago–RNA immunoprecipitation was performed on cytosolic protein fractions and strong depletion of nuclear RNA species in Ago–RNA was confirmed by the analysis of gene expression arrays (not shown). Moreover, in light of the strong anti-correlation observed when comparing Ago enrichment with ribosome enrichment (Fig. 4a) and the high consistency between Ago-enriched mRNAs and their predicted targeting miRNAs at the ES–ELA cell transition, we think that the ratio of RNA splice junctions possibly immunoprecipitated by Ago was negligible.
The strong correlation between the decrease of Ago–RNA and the increase of ribosome occupancy and MS detection at the ES–ELA transition strongly supports a protein translation control. However, the protein changes as measured by MS were bigger than ribosome occupancy. We should consider that the normalization method applied for polysome profiling did not take into account the general up-regulation of the translational machinery occurring during cell priming but only measured the modulation of ribosome occupancy above (or below) the overall change (see “Methods”). This, together with the fold change in ribosome occupancy exceeding the level of total RNA up-regulation, would account for the discrepancy.
We observed that the strongest variation, in terms of the number of mRNAs changing their association with Ago, happened at the ES–ELA transition, suggesting an important role for this mechanism of protein translation control in the maintenance of ground-state pluripotency.
A number of mRNAs released by Ago are related to DNA repair and replication. As ES cells and ELA correspond to inner cell mass (ICM) and post-implantation cells, respectively, the requirement of increased proliferation during priming is consistent with the higher cell division rate occurring after embryo implantation .
Chromatin regulators, whose role in pluripotency has been extensively discussed , are the other main class of mRNAs released by Ago at the ES–ELA transition. We found that DNMT, KDM, and SWI/SNF mRNAs are already present in ground-state pluripotent cells but their translation is either inhibited or kept at lower levels.
During early embryonic development, DNMT activity contributes to establish different genomic methylation patterns, which characterize somatic differentiated cells . The observations that a number of pluripotency-related genes are hypomethylated in stem cells  and that global DNA hypermethylation occurs in epiblast cells  suggest that the extensive reprogramming of genomic methylation pattern is a key event in ES cell priming. Global DNA hypomethylation is associated with ground-state pluripotency [40–43] but, surprisingly, appreciable levels of DNMTs are transcribed in naïve ES cells . Our finding that Dnmt1/3 mRNAs are significantly released by RISC and that Dnmt1/3 protein translation is enhanced at the ES–ELA transition might contribute to explain the presence of DNMT mRNAs in ESCs.
Members of the KDM and SWI/SNF families that we found to be released from Ago at the ES–ELA transition can also act as transcriptional repressors. Kdm2b (Jhdm1b) is responsible for the demethylation of H3K36 , which marks the transcribed region of active genes , while SmarcA4 is part of a complex that is necessary to silence Nanog upon differentiation . Indeed, inhibiting the increase of SmarcA4 impedes the ES–ELA transition when LIF and serum are removed from the culture medium. While this fact is not in contrast with the requirement of SmarcA4 for the maintenance of pluripotency , it confirms and accounts for the evidence that SmarcA4 activity is necessary for ES neuralization . Overall, despite the fact that these chromatin regulators are important for the maintenance of pluripotency, we found that their up-regulation is crucial for the control of priming.
Although we do not exclude that transcriptional regulation of DNMT, KDM, and SWI/SNF might contribute to cell priming, our results suggest that the release of selected chromatin regulators from RISC inhibition during priming may play a central role in reducing the pervasive transcription observed in ES cells  by reshaping the nuclear landscape of the cell . Moreover, the release from RISC inhibition of DNMT, KDM, and SWI/SNF members could explain the intrinsic tendency of ES cells to begin differentiating in the absence of external signals .
Naïve and primed mESC have distinct miRNA signatures, with only a few genomic clusters accounting for the majority of the difference in miRNA expression profiles [34, 35]. Global miRNA expression analysis of our in vitro model revealed a consistent modulation during the transition from ES to ELA cells. In agreement with the literature [34, 35], members of the miR-302/367 cluster, which are more expressed in EpiSC, and the miR-290/295 cluster, which are more abundant in ES cells, were up- and down-regulated, respectively, at the ES–ELA transition. The role of miRNAs in ES cell cycle regulation, self-renewal, and their ability to differentiate has been extensively assayed using cell lines in which key components of miRNA processing, Dicer or Dgcr8, were inactivated [49–51]. Moreover, recent studies of global gene expression profiling performed at the single-cell level highlighted that Dgcr8-deficient ES cells are more similar to ES cells cultured in 2i and LIF rather than in serum and LIF  and thus suggested that miRNAs act as key mediators of the transition from ground state pluripotency to primed states, with their absence mimicking the inhibition of the Erk and GSK3 signaling pathways observed in 2i culture. This observation, together with many in vitro studies on the impact of miRNAs in pluripotency, suggest that miRNAs are dispensable for maintaining a ground state of pluripotency in vitro . However, Dicer inactivation dramatically down-regulated Oct4 in the epiblast, causing the loss of stem cells, possibly due to their premature differentiation . In our work, we tried to solve this apparent contradiction between the in vitro and in vivo results by using transient Dicer inactivation assays in vitro. To the best of our knowledge, a key issue of all the in vitro studies of Dicer and Dgcr8 inactivation published so far is that they were performed in cell lines obtained through extensive cell culture selection and stabilization. For this reason, their use neither permitted the role of gene expression buffering exerted by miRNAs in naïve pluripotent cells to be studied nor allowed analysis of the changes that are “naturally” occurring during cell priming. Indeed, by acutely inactivating Dicer or RISC, we found that global miRNA activity is required to maintain pluripotency and that its perturbation triggers the beginning of a differentiation program in ES cell cultures which corresponds to what is observed in vivo.
Our observations do not allow us to evince the exact set of miRNAs that are necessary to maintain pluripotency through the protein translation control of chromatin regulators. By in silico analysis, we identified quite a large number of miRNAs whose modulation at the ES–ELA transition and predicted binding to DNMT, KDM, and SWI/SNF mRNAs might support a direct role in cell priming. For five of these, we showed the ability to inhibit protein translation through the predicted 3′ UTR target. Nevertheless, the high number of DNMT, KDM, and SWI/SNF members possibly involved in the control of cell priming and the consequent abundance of miRNAs predicted to be interactors suggest that a large set of miRNAs might be involved in orchestrating the inhibitory control of priming. The precise characterization of each single miRNA–mRNA interaction is beyond the scope of this work.
ES culture and neuralization
Murine ES cell lines E14Tg2A, 46C, and TNG-A (transgenic Sox1-GFP and Nanog-GFP ES cells, kindly gifted by Dr. A. Smith, University of Cambridge, UK) [16, 48] and the Dicer flox/null ES cell line (generously provided by Dr. G. Hannon, Cancer Research UK Cambridge Institute)  were cultured on gelatin-coated tissue culture dishes at a density of 40,000 cells/cm2. ES cell medium, which was changed daily, contained GMEM (Sigma), 10 % fetal calf serum, 2 mM glutamine, 1 mM sodium pyruvate, 1 mM non-essential amino acids, 0.05 mM β-mercaptoethanol, 100 U/mL penicillin/streptomycin, and 1000 U/mL recombinant mouse LIF (Invitrogen). For ground-state reprogramming of ES cells, medium was switched for three to five passages to 2iL medium: GMEM (Sigma) supplemented with N2/B27 (no vitamin A; Invitrogen), 2 mM glutamine, 1 mM sodium pyruvate, 1 mM non-essential amino acids, 0.05 mM β-mercaptoethanol, 100 U/mL penicillin/streptomycin, 1 μM MEK inhibitor PD0325901, 3 μM GSK3 inhibitor CHIR99021, and 1000 U/mL recombinant mouse LIF .
Chemically defined minimal medium (CDMM) for neural induction consisted of DMEM/F12 (Invitrogen), 2 mM glutamine, 1 mM sodium pyruvate, 0.1 mM non-essential amino acids, 0.05 mM β-mercaptoethanol, and 100 U/mL penicillin/streptomycin supplemented with N2/B27 (no vitamin A; Invitrogen). The protocol of ES neuralization consisted of three steps (Fig. 1a). During step 1, dissociated ES cells were washed with DMEM/F12, aggregated in agar-coated culture dishes (65,000 cells per cm2), and cultured as floating aggregates in CDMM for 2 days. The second day, 75 % of the CDMM was renewed.
In step 2, ES cell aggregates were dissociated and cultured in adhesion (65,000 cells per cm2) on poly-ornithine (Sigma; 20 μg/ml in sterile water, 24 h coating at 37 °C) and natural mouse laminin (Invitrogen; 2.5 μg/ml in PBS, 24-h coating at 37 °C) for 4 days, changing CDMM daily. In step 3, after a second dissociation, ES cells were cultured for an additional 4 days in CDMM devoid of B27 supplement to drive terminal differentiation, using the same type of seeding density and coated surface described for step 2. Serum employed for trypsin inactivation was carefully removed by washing with DMEM/F12. LIF serum withdrawal coupled to culturing as cell aggregates for the first 2 days results in a very low amount of cell death . Cell viability, which was monitored by a trypan blue exclusion test and cell counting, was more than 90 % and did not vary significantly between different steps of the differentiation protocol. Analysis of global gene expression on GO terms of apoptosis between ES and ELA cells showed no significant change (GO 0006917, induction of apoptosis, p = 0.668; GO 1900117, regulation of execution phase of apoptosis, p = 0.373; GO 0097194, execution phase of apoptosis, p = 0.559). The following factors were tested by addition during step 1: 5-azacytidine (Sigma-Aldrich; 50 nM) and 5-carboxy-8-hydroxyquinoline (Sigma-Aldrich; 100 μM).
Ago-bound RNA immunoprecipitation
In order to isolate cellular mRNAs that are bound by Ago proteins, we used a RIP-Assay Kit for microRNA (MBL Japan) according to the manufacturer’s instructions. The protocol includes a step for the isolation of the cytosolic fraction, which allowed us to avoid the co-precipitation of nuclear Ago and RNAs. Anti-pan Ago antibody, clone 2A8 (Ago*; Millipore MABE56) , mouse total IgG (Millipore; 12-371), or Anti-Ago2 antibody, clone 9E8.2 (Millipore; 04-642) were used for immunoprecipitation (IP). For each single experiment, we started from 107 cells collected by pooling three (ES) to ten (ELA, NPC, NPC/NEU Neu) different samples. The binding of RNP to the beads was confirmed after IP by western blotting with anti-Ago* antibody (Additional file 1: Figure S1a). The global screening of Ago–RNA by microarray analysis was performed following IP with Anti-pan Ago antibody (see the next section). Three independent experiments were performed. IP with Anti-Ago2 antibody was used to validate association of selected mRNAs with Ago by RT-PCR using a spike-in as a reference (Additional file 1: Figure S2f, g). Tol2 transposase mRNA (8 ng) synthesized in vitro was added to each pellet as a spike-in before IP. Mouse total IgGs were always used as an IP negative control. Two independent experiments were performed.
Microarray hybridization and data analysis
Total RNA was extracted with NucleoSpin RNA II columns (Macherey-Nagel). RNA quality was assessed with an Agilent Bioanalyzer RNA 6000 Nano kit; 200 ng of RNA was labeled with Low Input Quick Amp Labeling Kit, One-Color (Agilent Technologies, purified and hybridized overnight onto an Agilent SurePrint G3 Mouse Gene Expression Array (8x60K) before detection according to the manufacturer’s instructions. An Agilent DNA microarray scanner (model G2505C) was used for slide acquisition and spot analysis was performed with Feature Extraction software (Agilent Technologies).
For a gene probe to be considered significantly bound to Ago, we assumed that the distribution of signal is bimodal, with a fraction of genes not bound (first peak, background) and another displaying a heterogeneous degree of binding. In fact this is what we observed in the distribution of fluorescence intensity (Additional file 1: Figure S1c). After fitting the distribution to a Gaussian mixture model , we set a threshold T1 equal to (μ1 + 3σ1), which corresponds to p = 0.001 of false positive bound mRNA. Only mRNAs above this threshold were taken into account.
For each step, Ago enrichment (En) was evaluated as log2 ratio between Ago–RNA and total RNA levels after quantile normalization and an arbitrary threshold of T2 ≥ 2 was applied in order to capture only the most enriched genes. Conversely, genes with enrichment less than 0.5 were considered not significantly enriched (T3; this value was relaxed to 1.2 for subsequent analyses; e.g., GO analysis and Fig. 4b–d).
For each transition, genes released from Ago have En ≥ T2 and En+1 ≤ T3 (as depicted in Additional file 1: Figure S1b, right panel) while, conversely, genes loaded by Ago have En ≤ T3 and En+1 ≥ T2.
Finally, a threshold on the differential enrichment between steps (ΔE = En+1 − En) was used as a filter to sort out top regulated genes (|ΔE| ≥ 1.5).
Gene ontologies over-represented in the subset of Ago-released genes during the ES–ELA transition were evaluated using the web service Database for Annotation, Visualization and Integrated Discovery (DAVID) .
Genes that were differentially expressed between ES and ELA cells with |log2FC| ≥1 (where FC is fold change), corrected p value <0.05, and mean log2 expression ≥7.5 (n = 1207; Additional file 6: Table S4) were used for hierarchical clustering in order to assess the effect of manipulations. RNA from three different sets of treatments was pooled.
We changed the cutoffs for different analyses of global mRNA expression because the dynamic ranges of mRNA change were different in the different experimental settings: we used log2 > 2.5 when comparing ES, ELA, NPC, and neuron cells, which are very different from each other, while we used log2 > 1 when comparing ES cells with ELA cells, which are more similar.
Euclidean distance and a complete linkage algorithm were employed for clustering.
To assay the concordance of gene modulation between different conditions and obtain the statistical significance, we used a sign test. The sign test was performed comparing the fold change of ES–ELA cells to the fold change of treated/control ELA cells for each gene. Randomized sets of data used as control were obtained by substituting the fold changes of treated/control datasets with fold changes of random genes (randomization 1) or by scrambling the order of treated/control ELA cell values (randomization 2).
Profiling of ribosome-associated RNA (polysome profiling)
For the pull-down of ribosome-associated RNA we were inspired by the TRAP (translating ribosome affinity purification) protocol developed in the laboratory of Dr. Nathaniel Heintz [26, 27]. In short, this protocol is based on the pull-down of ribosomes using the large subunit protein Rpl10a, which is fused with GFP. We modified this protocol in order to pull down the endogenous Rpl10a protein as this would allow us to sample translation without altering the endogenous translational machinery. Immediately before the analysis, ES and ELA cells were trypsinized and resuspended in the appropriate growth medium supplemented with 100 μg/mL cycloheximide in order to stall ribosomes. After a 15 min incubation at 37 °C, cells were recovered and washed with ice-cold PBS plus 100 μg/mL cycloheximide to remove traces of serum and medium. Each sample (for both ELA and ES cells) included approximately 20 million cells and we conducted all experiments in duplicate. The samples were resuspended in 1 mL ice-cold lysis buffer (20 mM HEPES/KOH pH 7.3, 150 mM KCl, 10 mM MgCl2, 1 % IGEPAL CA 630 (Sigma), 0.5 mM DTT, 100 μg/mL cycloheximide, 1 pill/10 mL complete MINI EDTA free protease inhibitor (Roche), 3 μL/mL RNAsin RNAse inhibitor (Promega), 3 μL/mL SUPERASEin RNAse inhibitor (Life Techologies)) and incubated on ice for 10 min before being triturated in a chilled dounce homogenizer with approximately 20 strokes. After lysis was completed, cellular debris (mostly formed by nuclei and fragments of cell membrane) was removed by centrifugation at 2000 RCF for 10 min. After centrifugation, 0.1 volumes of 300 mM 1,2-diheptanoyl-sn-glycero-3-phosphocholine (DHPC, Generon Ltd) were added to each sample to further solubilize ribosomes, followed by a 10 min incubation on ice and by a second centrifugation at 20,000 RCF for 10 min. The final clarified lysate was transferred to a new low-retention tube and 0.1 volumes were removed to be used as “input” samples for the measurements of total RNA levels. To immunoprecipitate ribosomes, 16 μg of anti-Rpl10a antibody (Sigma, catalogue number WH0004736M1) or an equivalent amount of anti-GFP antibody (control, Sigma catalogue number G1544) were added to each sample and incubated overnight with end-to-end agitation at 4 °C. The following day, 200 μL of 30 mg/mL Protein-G Dynabeads (Life Technologies) were washed in PBS plus 0.1 % Tween 20 (PBST), blocked in PBST plus 1 % bovine serum albumin (BSA) for 1 h, and added to the lysate/antibody mix. Beads were incubated overnight at 4 °C with end-to-end agitation. The following day the precipitated ribosomes were washed four times in 1 mL high-salt buffer (20 mM HEPES/KOH pH 7.3, 350 mM KCl, 10 mM MgCl2, 1 % IGEPAL CA 630 (Sigma), 0.5 mM DTT, 100 μg/mL cycloheximide) and RNA was extracted using the RNEasy Micro kit (Qiagen) as per the supplier’s protocol. The results were quantified by spectrophotometry (Nanodrop) and capillary electrophoresis (Bioanalyzer RNA Nano kit; Additional file 1: Figure S2b). RNA (2 μg) was used to prepare each RNA-seq library. Purified RNA was ribosome-depleted using the RIBO-Zero Gold kit (Illumina), which uses magnetic beads conjugated to rRNA-specific oligonucleotides to remove ribosomal RNA. The final purified RNA was analyzed by capillary electrophoresis (Bioanalyzer Pico Kit) and retro-transcribed and cloned into a sequencing library using the SCRIPTSeq v2 library preparation kit (Illumina) according to the supplier’s protocol. Sequencing libraries were evaluated by Bioanalyzer (to verify the size distribution of fragments) and quantified by quantitative PCR using primers specific for Illumina adaptors (KAPA library quantification kit). All eight libraries (two replicates for ELA and ES cells, both IP and input) were mixed at a similar final concentration before cluster generation and sequencing on a HiSeq 2500 instrument.
Polysome profiling data analysis
All the libraries were sequenced in a multiplexed lane using v4 chemistry on a Illumina HiSeq instrument, producing approximately 250 million reads of raw data. After demultiplexing, each library comprised between 20 and 50 million reads, with the input library overall producing fewer reads than the ribosome IP libraries. This effect was most likely due to an unequal loading and to the different nature of the input RNA samples (including a variety of RNA and small RNA types, while the ribosome IP libraries were mostly composed of messenger RNAs) and raised no concern for the subsequent analysis (Additional file 1: Figure S2a).
All the datasets were subject to quality control using the fastx toolkit software through the following analyses: (1) quality score per cycle, (2) nucleotide distribution per cycle, (3) PCR duplicates (estimated by collapsing the reads). All libraries yielded good quality data, with the exception of one of the replicates for the ELA input sample, which showed a lower complexity and higher amount of duplication (resulting in lower mapping to the transcriptome during the later steps of the analysis; data not shown).
Reads contained in each of the remaining libraries were mapped to the mouse transcriptome (USCS mm10 release) using the software TopHat2 , which accounts for spliced transcripts. Repeated regions (i.e., transposable elements) were excluded and a maximum of two hits on the genome and two mismatches were allowed during mapping. Mapped reads were assigned to genes using the htseq-count script , which counts features overlapping each gene in the released transcriptome. Counts were calculated for each gene_id, collapsing eventual splicing isoforms. The resulting raw counts were normalized using the DESeq2 R package . The DESeq scaling factor for a given lane is computed as the median of the ratio, for each gene, of its read count over its geometric mean across all samples. It is important to notice that this normalization (and, effectively, any other normalization applicable to our data) is based on the assumption that the majority of genes are not differentially expressed among samples. This is not necessarily true during differentiation as there are reports indicating a massive up-regulation of the translational machinery , which would imply a higher read number for most genes during the ES–ELA transition. Unfortunately, it is impossible for any “internal” normalization method to account for global changes. As a consequence, our analysis of differential translation only indicates as significant genes that vary above or below the overall change level. We believe that this further increases the significance of our findings since our results include genes subject to a specific regulation rather than a global effect.
After normalization, we analyzed our results to assess correlation among datasets and among conditions (through hierarchical clustering and principal component analysis; Additional file 1: Figure S2c, d). The results indicated that one of the replicates for the ELA input sample showed poor correlation with both the other input replicate and the other datasets (ELA.input.1; Additional file 1: Figure S2c, d). This, together with the lower mappability, lower read quality, and increased amount of PCR duplicates (as described above), led us to exclude this replicate from further processing. While this reduces the statistical power of our analysis, the fact that we are filtering our results to include only genes presenting high expression and a high degree of change in their ribosome enrichment makes us confident that our conclusions are still significant.
The remaining libraries were filtered to include only genes with an average of 50 counts in at least one of the two conditions (ES or ELA). This threshold was chosen to exclude less-expressed genes, which would be the most affected by the inherent noise of our method. For all the genes passing this first filter, ribosome enrichment was calculated by dividing the average ribosome IP count by the input count. A second filter was then applied, excluding genes for which the standard deviation of the counts accounted for more than 20 % of the average value. This filtering produced a list of approximately 5000 genes, for which we calculated a log2 fold change in ribosome enrichment (ELA average enrichment/ES average enrichment; log2RiboΔE). For our GO enrichment analysis, we considered only genes with log2RiboΔE >1 during differentiation, producing a final list of ~400 differentially regulated genes.
While the GO category enrichment results we obtained were all highly significant according to multiple statistics tests, we couldn’t properly assess the significance of changes in the ribosome enrichment of individual genes due to our decision to remove the faulty replicate in the ELA input dataset. We nevertheless filtered our data to exclude genes with high replicate-to-replicate variation and low expression and reported those with the highest changes, in both directions, in our raw data (Additional file 2: Table S1). Some of the genes individuated by our other analyses (Ago IP and protein MS) fell short of matching our filter in the ribosome enrichment data; nevertheless, they showed a clear increase in enrichment, which was consistent with the modulation highlighted by the other techniques. We decided to show these results since we believe that the fact that we obtained similar results through multiple and complementary analysis types constitute a strong proof of significance.
Semiquantitative real-time PCR
Total RNA was extracted from ES cells or tissue samples with NucleoSpin RNA II columns (Macherey-Nagel). ES cells from at least two to three different wells of 24-well plates were always pooled together to compensate for variability in cell seeding. RNA quantity and RNA quality were assessed by gel electrophoresis. For each sample, 200 ng of total RNA was reverse-transcribed (Eurogentech, RT Core kit). Amplified cDNA was quantified using GoTaq qPCR Master Mix (Promega) on a Rotor-Gene 6000 (Corbett). Primers used for amplification which were not previously reported  were taken from PrimerBank . Amplification take-off values were extracted using the built-in Rotor-Gene 6000 relative quantitation analysis function and relative expression was calculated with the 2-ΔΔCt method, normalizing to the housekeeping gene β-actin. Standard errors shown as error bars in all histograms were obtained from the error propagation formula as described in ; the statistical significance of three independent experiments was probed with a randomization test using the REST software .
Cells prepared for immunocytodetection experiments were cultured on poly-ornithine/laminin-coated round glass coverslips. Cells were fixed using 2 % paraformaldehyde for 10–15 min, washed twice with PBS, permeabilized using 0.1 % Triton X100 in PBS, and blocked using 0.5 % BSA in PBS. Primary antibodies used for microscopy included Oct3/4 (1:200; Santa Cruz Biotechnology C-10; sc-5279), Nanog (1:300; Novus Biologicals; NB100-58842), acetylated N-tubulin (clone 6-11B-1; 1:500; Sigma; T7451), neuronal class III β-tubulin (1:500; Covance; MRB-435P), Nestin (1:200; Millipore; MAB353), Pax6 (1:400; Covance; PRB-278P), and GFP (1:1000; Life Technologies; A-6455). Primary antibodies were incubated for 2 h at room temperature; cells were then washed three times with PBS (10 min each). Alexa Fluor 488 and Alexa Fluor 568 anti-mouse or anti-rabbit IgG conjugates (1:500; Molecular Probes; A-11001 and A-11011, respectively) were incubated for 1 h at room temperature in PBS containing 0.5 % BSA for primary antibody detection, followed by three PBS washes (10 min each). Nuclear staining was done with DAPI (2 μg/mL; Sigma).
The same protocol was applied on sections of embryoid bodies after fixation (2 h at 4 °C in 2 % PFA), dehydration (overnight at 4 °C in 30 % sucrose/PBS), inclusion in Tissue-Tek OCT compound (Sakura), and cryostat sectioning (12–16 μm) on SuperFrost glass slides (Thermo).
Nuclear extract preparation and proteomics sample pre-processing
About 107 cells per sample were dissociated, washed with ice-cold PBS, and resuspended in 400 μL of nuclear extraction buffer (50 mM Tris-Hcl pH 7.9, 10 mM KCl, 0.2 % NP40, 10 % glycerol, and 1 mM PMSF) for 3 min in ice, then nuclei were pelleted by centrifugation at 6000 RPM for 3 min at 4 °C. Nuclei were washed twice with 1 mL of nuclear extraction buffer devoid of NP40; before proceeding to mass spectrometry analysis, a small aliquot was lysed in loading buffer and used for western blotting to confirm that the enrichment in nuclear fraction was homogeneous between samples.
Nuclear extracts were lysed using lysis buffer containing TRIS HCl 50 mM pH 8.1, 0.5 % Triton X-100, and 0.25 % deoxycholic Na and sonicated. Samples were centrifuged at 10,000 g for 10 min at 4 °C to discard cell debris and buffer detergent was removed using a detergent removal spin column (Pierce, Thermo Scientific, USA). Protein concentration was determined by bicinchoninic acid assay (Pierce, Thermo Scientific, USA). We diluted 100 μg of protein in 25 mM of ammonium hydrogen carbonate (pH 8) and reduction was obtained by adding 5 mM dithiothreitol with an incubation of 20 min at 80 °C. Iodoacetamide was added to the samples to a final concentration of 10 mM and incubated in the dark for 30 min at 37 °C. Digestion was performed incubating overnight with 100:1 substrate:trypsin (Roche, Germany) at 37 °C. Nucleic acid contaminants were removed with Amicon Ultra-3 K centrifugal devices (Merck Millipore, Germany), then flow-through containing peptides was loaded on a C18 cartridge in order to eliminate debris and filtered with 0.22 μm filter.
Liquid chromatography-tandem MS analysis
Chromatographic separation of peptides was performed using a nano-HPLC system (Eksigent, ABSciex, USA). The loading pump pre-concentrated the sample in a pre-column cartridge (C18 PepMap-100, 5 μm 100 A, 0.1 x 20 mm; Thermo Scientific, USA).
Chromatographic separation of peptides was performed using a C18 PepMap-100 column (3 μm, 75 μm × 150 mm; Thermo Scientific, USA) at a flow rate of 300 nL/min. Runs were performed with eluent A (Ultrapure water, 0.1 % formic acid) under a 60 min linear gradient from 25 to 40 % of eluent B (acetonitrile (Romil, UK), 0.1 % formic acid) followed by 10 min of a purge step and a 20-min re-equilibration step. The column was directly coupled to a TripleTOF 5600 System (ABSciex, USA), equipped with a DuoSpray ion source (ABSciex, USA). Peptides eluted from chromatography were directly processed using a TripleTOF 5600 mass spectrometer (ABSciex, USA). The mass spectrometer was controlled by Analyst TF 1.6.1 software (ABSciex, USA). For information-dependent acquisition (IDA) analysis, survey scans were acquired in 250 ms and 25 product ion scans were collected if exceeding a threshold of 125 counts per second (counts/s). Dynamic exclusion was set for one-half of peak width (∼8 s), then the precursor was refreshed off the exclusion list. MS/MS data were processed with ProteinPilot Software (ABSciex, USA) using the Paragon and Pro Group Algorithms and UniProt 2013 as protein database for Mus musculus. The false discovery rate (FDR) analysis was done using the integrated tool in the ProteinPilot software and a confidence level of 95 % was set. Label-free statistical comparative analysis  was performed using MarkerView (ABSciex, USA). The ion chromatograms of high confidence peptides identified by ProteinPilot were extracted using PeakView Software (ABSciex, USA), then MS peak areas and identifications were imported into MarkerView. Normalization of the total plaque area (plaque size) was done using a global normalization of profiles (total protein content). Principal component analysis (data not shown) was performed in order to determine groupings among the data set. The two groups (ES and ELA, n = 3) were compared with t-test using a threshold p value ≤0.05 and fold change >2.
Cells were lysed with RIPA buffer (50 mM Tris-HCl pH 7.6, 1 % NP40, 0.5 % deoxycholic acid, 150 mM NaCl, 1 mM EDTA, 1 mM PMSF, 1 % SDS) supplemented with Complete Protease Inhibitor Cocktail (Roche); lysate was incubated for 30 min on ice and sonicated three times for 10 sec each on medium power in order to reduce viscosity. Supernatant was harvested by centrifugation (10 min at 13,000 RPM, 4 °C) and quantified with a Micro BCA Protein Assay Kit (Thermo Scientific). Samples were denatured by adding sample buffer (LDS Sample Buffer, Thermo Scientific) and boiled for 10 min on a thermal block at 99 °C. The total protein extract (10 to 50 μg) was resolved on 8–10 % acrylamide gels, transferred on a nitrocellulose membrane (Hybond-c Extra, GE Healthcare), blocked with 5 % milk proteins in TBST (50 mM Tris pH 7.6, 150 mM NaCl, 0.05 % Tween-20), and probed with primary antibodies, including Dnmt3b (1:2000; Imgenex/Novus Biologicals; IMG-184A), SmarcA4 (1:1000; Abcam; ab4081), Kdm2b (1:2000; Merck Millipore; 09-864), PolII (clone CTD4H8; 1:1000; Millipore), α-tubulin (clone B-5-1-2; 1:5000; Sigma-Aldrich; T6074), GAPDH (Sigma-Aldrich; G9545), panAgo clone 2A8(1:250; Millipore; MABE56) and Ago2 clone 9E8.2 (1:2000; Millipore; 04-642). After 1 h of incubation, the membrane was washed three times with TBST (15 min each) and probed with a HRP-conjugated secondary anti-mouse or anti-rabbit antibody for 1 h (Santa Cruz Biotechnology; sc-2005 and sc-2030, respectively). After three more washes, signal was revealed by means of an enhanced chemiluminescence kit (G&E Healthcare) on a BioMax XAR Film (Kodak).
Lentiviral vector construction and use
Knock-down vectors against SmarcA4 or luciferase (shSA4 and shCtrl, respectively) were constructed in order to sustain high expression of shRNAs in both ES and differentiating conditions: an IRES-PuroR-shRNA cassette from pGIPZ vectors (Open Biosystems) was excised by double restriction with NotI and MluI (NEB), filled in with T4 DNA polymerase (NEB), and inserted into a modified pWPXLd vector (Addgene number 12258; original GFP was swapped with DsRed with restriction sites MluI/XmaI) which was cut with NdeI and blunted as above. This resulted into a self-inactivating transfer plasmid carrying a DsRed-IRES-Puromycin resistance cassette with the shRNA of interest, driven by the strong ubiquitous promoter EF1α. As a reporter of human Nanog promoter activity we employed a PL-SIN-Nanog-EGFP vector  (Addgene, number 21321).
CRISPR/Cas9 vectors were based on LentiCRISPR v2  (Addgene, number 52961) and constructed following the manufacturer’s instructions. Briefly, vector was linearized with BsmBI (NEB), dephosphorylated with calf intestinal phosphatase (NEB), and ligated with an oligo duplex obtained by annealing two single-stranded DNA sequences synthesized in vitro and 5′-phosphorylated by T4 Polynucleotide Kinase (NEB). The following oligos were designed and synthesized for Ago1–2 and control (CTRL) guides, respectively: Ago1–2_fw, CAC CGC GCA TCA TCT TCT ACC GCG A; Ago1–2_rev, AAA CTC GCG GTA GAA GAT GAT GCG C; CTRL_fw, CAC CGG CGA GGT ATT CGG CTC CGC G; CTRL_rev, AAA CCG CGG AGC CGA ATA CCT CGC C. For CRE-GFP vector construction, a CRE-NLS-IRES cassette from CRE-PuroR vector (Addgene, number 30205) was excised by double restriction with NotI and NdeI (NEB), filled in with T4 DNA polymerase (NEB), and inserted into a pWPXLd vector (Addgene, number 12258) which was cut with MluI and blunted as above. This resulted in a self-inactivating transfer plasmid carrying a CRE-IRES-EGFP cassette driven by the strong ubiquitous promoter EF1α.
Lentiviral vectors were produced by transient transfection of 293 T cells using 150 nM polyethylenimine (PEI) reagent (Sigma) with either PL-SIN-Nanog-EGFP, shCtrl, or shSA4 plasmids, together with the Δ8.91 packaging and a VSV-G envelope expressing plasmids  in a ratio of 20 μg:15 μg:5 μg, respectively, per single 100-mm dish. Transfection medium was discarded 24 h after transfection, then viral particles were collected at 48 h and 72 h, pooled and frozen at −80 °C.
For shRNA-mediated knock-down experiments, ES cells underwent two rounds of spinoculation (1 h, 1100 RCF at room temperature) , then transduced cells were selected in puromycin (0.75 μg/mL) for 48 h before ELA formation (Additional file 1: Figure S3a). By 0 DIV more than 90 % of cells were DsRed+, achieving a knock-down of SmarcA4 mRNA to 21.2 % of the expression level in shCtrl-infected cells (Additional file 1: Figure S3b).
For CRISPR/Cas9 experiments, cells were transduced and selected as above and maintained in ES medium (either serum and LIF or 2i plus LIF), achieving genomic editing (Additional file 1: Figure S4b) and strong down-regulation of Ago* protein levels (Fig. 6b) by 4 days after transduction. Rapamycin (100 nM; Abcam) or MG-115 (1 μM; Sigma-Aldrich) activity was tested by adding the drugs to the culture medium 48 h after transduction.
For Dicer conditional knock-out experiments, a Dicer fl/- ES cell line was transduced with CRE-GFP lentiviral vector, then, after 48 h, GFP-positive cells were sorted by cytofluorimetry and reseeded for culturing in ES medium. Control (not excised) cells were transduced with pWPXLd vector devoid of CRE recombinase.
Cells were detached by trypsinization, washed, and resuspended in PBS at room temperature, then analyzed with a FACSCalibur cytometer (BD Bioxciences). At least 25,000 events per sample were collected. Data were processed with Bioconductor package flowCore .
HNP cell line generation
For the generation of a mouse ES cell line carrying human Nanog proximal promoter (Additional file 1: Figure S3d), E14Tg2A cells in FCS plus LIF were infected with PL-SIN-Nanog-EGFP lentiviral vector. After six passages from transduction, cells were detached by trypsinization, washed, and resuspended in PBS, 2 % FBS, and 2 mM EDTA and kept on ice. The Nanog::GFP+ population (Additional file 1: Figure S3e) was sorted on a FACS Aria II Flow Cytometer (BD Biosciences) and reseeded as a polyclonal line for amplification (HNP) to average positional effects of promoter insertion. The resulting cell line remains GFP+ after multiple passages in 2iL (data not shown) but promptly switches off fluorescence upon differentiation (Additional file 1: Figure S3f), with kinetics consistent with developmental Nanog regulation. EGFP mRNA levels correlate well with that of endogenous Nanog, as seen by real-time PCR (Additional file 1: Figure S3g).
Nanog promoter methylation assay
For the characterization of Nanog promoter methylation status, we extracted genomic DNA with the Wizard SV Genomic DNA Purification System (Promega). DNA (500 ng) was bisulfite-treated using an EZ DNA Methylation-Lightning Kit (Zymo Research) and the region of proximal Nanog promoter (Fig. 5g) was amplified by PCR with a primer set from  and cloned into a pGEM-T Easy vector (Promega). For each experiment, DNA obtained from 10 to 15 different colonies transformed with plasmid was sequenced on an ABI3730 DNA analyzer (Applied Biosystems) and data were analyzed with BiQ Analyzer software .
Small RNA-seq and analysis
The sequencing procedure was done using Illumina Sequencing By Synthesis technology. Total RNA (1 μg) extracted with a miRNeasy Mini kit (Qiagen) was used for library preparation (Illumina, TruSeq Small RNA Sample Prep Kit) following the manufacturer’s description. Libraries were sequenced on a miSeq (Illumina) running in 50-bp single-read mode using sequencing chemistry v3 and demultiplexed in FASTQ format using CASAVA v.1.8 (Illumina).
Library adaptors were trimmed with Trimmomatic  and reads were mapped to the mouse genome (NCBI37/mm9) with Bowtie . miRNAs reads were annotated according to miRbase20  and summarized using featureCounts . Normalization was performed as reads per million (RPM) and differential expression was evaluated with the R package edgeR .
In silico analysis of 3′ UTR binding
We predicted the affinity of the specific miRNA subset which is expressed in ES cells and is down-regulated during the ES–ELA transition (log2 mean count per million ≥5; log2 fold change ≤ −1; Additional file 5: Table S5) towards different mRNA families (as in Additional file 3: Table S2). For this purpose, we took advantage of the UCSF Table Browser data retrieval tool , which allows batch downloading of the 3′ UTR sequences of the genes of interest; when genes had several 3′ UTR variants, the longest one was taken into account. We employed the miRanda algorithm, which searches for complementarity matches between miRNAs and 3′ UTRs using dynamic programming alignment and thermodynamic calculation  (Additional file 5: Table S5). We set a minimum score of 140 and a maximum energy of −10 kcal/mol as thresholds. For each gene, miRNA scores were summed and the distributions of sums in different families were tested for statistically significant difference with a Wilcoxon test. 3′ UTR randomization was obtained by scrambling the 3′ UTR sequences . Global miRNA–mRNA binding prediction was performed using the miRVestigator framework  submitting as inputs the lists of either the top Ago-released (ΔE ≤ −2) or Ago-loaded (ΔE ≥ 2) mRNAs during the ES–ELA transition. Only miRNAs with a p value ≤0.002 were retained in the z-score analysis.
3′ UTR transfection assay
The 3′ UTR of mouse Kdm2b, SmarcA4, and Dnmt3b genes were PCR-amplified from E14 ES cell genomic DNA with GoTaq (Promega) using the following primers: Kdm2b_forward, CCA GGA CAA GTA TGT AAA TAT GGA GGG; Kdm2b_reverse, GTA CAA TTG TTT ATA TAA ATC CAA CAA AGG TC; SmarcA4_forward, ACC AGA CAT TCC TGA GTC CTG; SmarcA4_reverse, CCA AGG CAA GTC CTA CTT ATT TAT TTC; Dnmt3b_forward, TTC TAC CCA GGA CTG GGG AGC; Dnmt3b_reverse, GAG AAA TAC AAC TTT AAT CAA CCA GAA AGG. The products of amplification (1035 bp for Kdm2b, 1262 bp for SmarcA4, and 1326 bp for Dnmt3b) were cloned downstream of EGFP in a pEGFPC1 expression vector (Addgene) and used in co-transfection assays together with pDsRed-Express-C1 (Addgene) and mature miRNAs (miRIDIAN microRNA Mimics, Dharmacon). Plasmid DNA and small RNA were co-transfected using Xfect polymer and Xfect RNA transfection reagents (Clontech), respectively, according to the manufacturer’s instructions. E14 cells were plated in 12 multi-well plates (4 × 105/well) in ES cell medium 8 h before transfection. Cells were co-transfected in triplicate with 50 pmol of mature miRNA, 0.5 μg of 3′ UTR/pEGFPC1 vector or pEGFPC1 empty vector, and 0.5 μg of pDsRed-Express-C1, which was used as an internal reference of transfection. Both pEGFPC1 and pDsRed-Express-C1 vectors carried a CMV promoter and SV40_PA_terminator to drive constitutive expression. Four hours after transfection, ES cell medium was replaced with CDMM (ELA cells) or new ES cell medium (ES cells). EGFP and DsRed fluorescence was detected by a GloMax®-Multi + Microplate Multimode Reader (Promega) 48 h after transfection. The EGFP/DsRed fluoresence ratio was compared among different samples of transfected cells to evaluate the effect of the 3′ UTR on EGFP translation.
Availability of data and material
Raw microarray data have been deposited in the Gene Expression Omnibus database with SuperSeries accession GSE79655 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE79655) and SubSeries accessions GSE79649 (Ago-RIP), GSE79650 (ELA-EpiSC profiling), and GSE79652 (chromatin remodeling manipulation). Raw polysome profiling and small RNA sequencing reads are available as SubSeries GSE79653 and GSE79654, respectively.
No ethics approval was required for this study.
We thank Austin Smith for providing 46C and TNG-A cell lines and for critical reading and discussion; Gregory Hannon for providing the Dicer flox/null ES cell line, pGIPZ lentiviral expression system and WPXLd vector; Roberto Ripa for providing the Tol2 transposase mRNA; Paolo Malatesta and Cristina Di Primio for technical support in cell sorting and lentiviral preparation, respectively; Maria Antonietta Calvello and Vania Liverani for technical assistance in the preparation of DNA constructs. We are grateful to Andrew Bannister for comments that greatly improved the manuscript. We thank the Cambridge Institute Genomics Core for Illumina sequencing.
This work was funded by: FIRB RBAP10L8TY (MIUR), Fondazione Roma and PAINCAGE FP7 Collaborative Project number 603191 (RB,MD); Flagship Project InterOmics PB.05 and MIUR-PRIN-2012 (FC); Wellcome Trust Core Grant reference 092096 and Cancer Research UK Grant Reference C6946/A14492 (LP); CRUK-Cambridge Institute Core Grant reference C14303/A17197 (DB).
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Ying QL, Nichols J, Chambers I, Smith A. BMP induction of Id proteins suppresses differentiation and sustains embryonic stem cell self-renewal in collaboration with STAT3. Cell. 2003;115:281–92.View ArticlePubMedGoogle Scholar
- Chen T, Dent SYR. Chromatin modifiers and remodellers: regulators of cellular differentiation. Nat Rev Genet. 2014;15:93–106.View ArticlePubMedPubMed CentralGoogle Scholar
- Schmiedel JM, Klemm SL, Zheng Y, Sahay A, Bluthgen N, Marks DS, et al. MicroRNA control of protein expression noise. Science. 2015;348:128–32.View ArticlePubMedGoogle Scholar
- Efroni S, Duttagupta R, Cheng J, Dehghani H, Hoeppner DJ, Dash C, et al. Global transcription in pluripotent embryonic stem cells. Cell Stem Cell. 2008;2:437–47.View ArticlePubMedPubMed CentralGoogle Scholar
- Marks H, Kalkan T, Menafra R, Denissov S, Jones K, Hofemeister H, et al. The transcriptional and epigenomic foundations of ground state pluripotency. Cell. 2012;149:590–604.View ArticlePubMedPubMed CentralGoogle Scholar
- Bartel DP. MicroRNAs: target recognition and regulatory functions. Cell. 2009;136:215–33.View ArticlePubMedPubMed CentralGoogle Scholar
- Bernstein E, Kim SY, Carmell MA, Murchison EP, Alcorn H, Li MZ, et al. Dicer is essential for mouse development. Nat Genet. 2003;35:215–7.View ArticlePubMedGoogle Scholar
- Melton C, Judson RL, Blelloch R. Opposing microRNA families regulate self-renewal in mouse embryonic stem cells. Nature. 2010;463:621–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Bertacchi M, Pandolfini L, Murenu E, Viegi A, Capsoni S, Cellerino A, et al. The positional identity of mouse ES cell-generated neurons is affected by BMP signaling. Cell Mol Life Sci. 2013;70:1095–111.View ArticlePubMedPubMed CentralGoogle Scholar
- Pelton TA, Sharma S, Schulz TC, Rathjen J, Rathjen PD. Transient pluripotent cell populations during primitive ectoderm formation: correlation of in vivo and in vitro pluripotent cell development. J Cell Sci. 2002;115(Pt 2):329–39.PubMedGoogle Scholar
- Bernemann C, Greber B, Ko K, Sterneckert J, Han DW, Araúzo-Bravo MJ, et al. Distinct developmental ground states of epiblast stem cell lines determine different pluripotency features. Stem Cells. 2011;29:1496–503.View ArticlePubMedGoogle Scholar
- Brons IGM, Smithers LE, Trotter MWB, Rugg-gunn P, Sun B, Chuva SM, et al. Derivation of pluripotent epiblast stem cells from mammalian embryos. Nature. 2007;448:191–6.View ArticlePubMedGoogle Scholar
- Tesar PJ, Chenoweth JG, Brook FA, Davies TJ, Evans EP, Mack DL, et al. New cell lines from mouse epiblast share defining features with human embryonic stem cells. Nature. 2007;448:2–8.View ArticleGoogle Scholar
- Chambers I, Colby D, Robertson M, Nichols J, Lee S, Tweedie S, et al. Functional expression cloning of Nanog, a pluripotency sustaining factor in embryonic stem cells. Cell. 2003;113:643–55.View ArticlePubMedGoogle Scholar
- Ye S, Li P, Tong C, Ying Q-L. Embryonic stem cell self-renewal pathways converge on the transcription factor Tfcp2l1. EMBO J. 2013;32:2548–60.View ArticlePubMedPubMed CentralGoogle Scholar
- Chambers I, Silva J, Colby D, Nichols J, Nijmeijer B, Robertson M, et al. Nanog safeguards pluripotency and mediates germline development. Nature. 2007;450:1230–4.View ArticlePubMedGoogle Scholar
- Thomson JA. Embryonic stem cell lines derived from human blastocysts. Science. 1998;282:1145–7.View ArticlePubMedGoogle Scholar
- Vallier L, Alexander M, Pedersen RA. Activin/Nodal and FGF pathways cooperate to maintain pluripotency of human embryonic stem cells. J Cell Sci. 2005;118(Pt 19):4495–509.View ArticlePubMedGoogle Scholar
- Guo G, Yang J, Nichols J, Hall JS, Eyres I, Mansfield W, et al. Klf4 reverts developmentally programmed restriction of ground state pluripotency. Nature. 2009;1069:1063–9.Google Scholar
- Eichhorn SW, Guo H, McGeary SE, Rodriguez-Mias RA, Shin C, Baek D, et al. mRNA destabilization is the dominant effect of mammalian microRNAs by the time substantial repression ensues. Mol Cell. 2014;56:104–15.View ArticlePubMedPubMed CentralGoogle Scholar
- Snow MHL. Gastrulation in the mouse: growth and regionalization of the epiblast. J Embryol Exp Morphol. 1977;42:293–303.Google Scholar
- Ahmed K, Dehghani H, Rugg-Gunn P, Fussner E, Rossant J, Bazett-Jones DP. Global chromatin architecture reflects pluripotency and lineage commitment in the early mouse embryo. PLoS One. 2010;5:e10531.View ArticlePubMedPubMed CentralGoogle Scholar
- Navarro P, Oldfield A, Legoupi J, Festuccia N, Dubois A, Attia M, et al. Molecular coupling of Tsix regulation and pluripotency. Nature. 2010;468:457–60.View ArticlePubMedGoogle Scholar
- Song J, Saha S, Gokulrangan G, Tesar PJ, Ewing RM. DNA and chromatin modification networks distinguish stem cell pluripotent ground states. Mol Cell Proteomics. 2012;11:1036–47.View ArticlePubMedPubMed CentralGoogle Scholar
- Khare SP, Habib F, Sharma R, Gadewal N, Gupta S, Galande S. HIstome--a relational knowledgebase of human histone proteins and histone modifying enzymes. Nucleic Acids Res. 2012;40:D337–42.View ArticlePubMedPubMed CentralGoogle Scholar
- Heiman M, Schaefer A, Gong S, Peterson JD, Day M, Ramsey KE, et al. A translational profiling approach for the molecular characterization of CNS cell types. Cell. 2008;135:738–48.View ArticlePubMedPubMed CentralGoogle Scholar
- Heiman M, Kulicke R, Fenster RJ, Greengard P, Heintz N. Cell type-specific mRNA purification by translating ribosome affinity purification (TRAP). Nat Protoc. 2014;9:1282–91.View ArticlePubMedPubMed CentralGoogle Scholar
- Silva J, Nichols J, Theunissen TW, Guo G, van Oosten AL, Barrandon O, et al. Nanog is the gateway to the pluripotent ground state. Cell. 2009;138:722–37.View ArticlePubMedPubMed CentralGoogle Scholar
- Santi DV, Norment A, Garrett CE. Covalent bond formation between a DNA-cytosine methyltransferase and DNA containing 5-azacytosine. Proc Natl Acad Sci U S A. 1984;81:6993–7.View ArticlePubMedPubMed CentralGoogle Scholar
- King ONF, Li XS, Sakurai M, Kawamura A, Rose NR, Ng SS, et al. Quantitative high-throughput screening identifies 8- hydroxyquinolines as cell-active histone demethylase inhibitors. PLoS One. 2010;5:e15535.View ArticlePubMedPubMed CentralGoogle Scholar
- Schaniel C, Ang Y-S, Ratnakumar K, Cormier C, James T, Bernstein E, et al. Smarcc1/Baf155 couples self-renewal gene repression with changes in chromatin structure in mouse embryonic stem cells. Stem Cells. 2009;27:2979–91.PubMedGoogle Scholar
- Kamiya D, Banno S, Sasai N, Ohgushi M, Inomata H, Watanabe K, et al. Intrinsic transition of embryonic stem-cell differentiation into neural progenitors. Nature. 2011;470:503–9.View ArticlePubMedGoogle Scholar
- Murchison EP, Partridge JF, Tam OH, Cheloufi S, Hannon GJ. Characterization of Dicer-deficient murine embryonic stem cells. Proc Natl Acad Sci U S A. 2005;102:12135–40.View ArticlePubMedPubMed CentralGoogle Scholar
- Jouneau A, Ciaudo C, Sismeiro O, Brochard V, Jouneau L, Vandormael-Pournin S, et al. Naive and primed murine pluripotent stem cells have distinct miRNA expression profiles. RNA. 2012;18:253–64.View ArticlePubMedPubMed CentralGoogle Scholar
- Greve TS, Judson RL, Blelloch R. microRNA control of mouse and human pluripotent stem cell behavior. Annu Rev Cell Dev Biol. 2013;29:213–39.View ArticlePubMedPubMed CentralGoogle Scholar
- Taliaferro JM, Aspden JL, Bradley T, Marwha D, Blanchette M, Rio DC. Two new and distinct roles for Drosophila Argonaute-2 in the nucleus: alternative pre-mRNA splicing and transcriptional repression. Genes Dev. 2013;27:378–89.View ArticlePubMedPubMed CentralGoogle Scholar
- Reik W, Dean W, Walter J. Epigenetic reprogramming in mammalian development. Science. 2001;293:1089–93.View ArticlePubMedGoogle Scholar
- Farthing CR, Ficz G, Ng RK, Chan CF, Andrews S, Dean W, et al. Global mapping of DNA methylation in mouse promoters reveals epigenetic reprogramming of pluripotency genes. PLoS Genet. 2008;4:e1000116.View ArticlePubMedPubMed CentralGoogle Scholar
- Smith ZD, Chan MM, Humm KC, Karnik R, Mekhoubad S, Regev A, et al. DNA methylation dynamics of the human preimplantation embryo. Nature. 2014;511:611–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Ficz G, Branco MR, Seisenberger S, Santos F, Krueger F, Hore TA, et al. Dynamic regulation of 5-hydroxymethylcytosine in mouse ES cells and during differentiation. Nature. 2011;473:398–402.View ArticlePubMedGoogle Scholar
- Ficz G, Hore TA, Santos F, Lee HJ, Dean W, Arand J, et al. FGF signaling inhibition in ESCs drives rapid genome-wide demethylation to the epigenetic ground state of pluripotency. Cell Stem Cell. 2013;13:351–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Habibi E, Brinkman AB, Arand J, Kroeze LI, Kerstens HHD, Matarese F, et al. Whole-genome bisulfite sequencing of two distinct interconvertible DNA methylomes of mouse embryonic stem cells. Cell Stem Cell. 2013;13:360–9.View ArticlePubMedGoogle Scholar
- Leitch HG, McEwen KR, Turp A, Encheva V, Carroll T, Grabole N, et al. Naive pluripotency is associated with global DNA hypomethylation. Nat Struct Mol Biol. 2013;20:311–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Faddah DA, Wang H, Cheng AW, Katz Y, Buganim Y, Jaenisch R. Single-cell analysis reveals that expression of nanog is biallelic and equally variable as that of other pluripotency factors in mouse escs. Cell Stem Cell. 2013;13:23–9.View ArticlePubMedPubMed CentralGoogle Scholar
- He J, Kallin EM, Tsukada Y-I, Zhang Y. The H3K36 demethylase Jhdm1b/Kdm2b regulates cell proliferation and senescence through p15(Ink4b). Nat Struct Mol Biol. 2008;15:1169–75.View ArticlePubMedPubMed CentralGoogle Scholar
- Barski A, Barski A, Cuddapah S, Cuddapah S, Cui K, Cui K, et al. High-resolution profiling of histone methylations in the human genome. Cell. 2007;129:823–37.View ArticlePubMedGoogle Scholar
- Ho L, Ronan JL, Wu J, Staahl BT, Chen L, Kuo A, et al. An embryonic stem cell chromatin remodeling complex, esBAF, is essential for embryonic stem cell self-renewal and pluripotency. Proc Natl Acad Sci U S A. 2009;106:5181–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Ying Q-L, Stavridis M, Griffiths D, Li M, Smith A. Conversion of embryonic stem cells into neuroectodermal precursors in adherent monoculture. Nat Biotechnol. 2003;21:183–6.View ArticlePubMedGoogle Scholar
- Kanellopoulou C, Muljo SA, Kung AL, Ganesan S, Drapkin R, Jenuwein T, et al. Dicer-deficient mouse embryonic stem cells are defective in differentiation and centromeric silencing. Genes Dev. 2005;19:489–501.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang Y, Medvid R, Melton C, Jaenisch R, Blelloch R. DGCR8 is essential for microRNA biogenesis and silencing of embryonic stem cell self-renewal. Nat Genet. 2007;39:380–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang Y, Baskerville S, Shenoy A, Babiarz JE, Baehner L, Blelloch R. Embryonic stem cell-specific microRNAs regulate the G1-S transition and promote rapid proliferation. Nat Genet. 2008;40:1478–83.View ArticlePubMedPubMed CentralGoogle Scholar
- Kumar RM, Cahan P, Shalek AK, Satija R, DaleyKeyser AJ, Li H, et al. Deconstructing transcriptional heterogeneity in pluripotent stem cells. Nature. 2014;516:56–61.View ArticlePubMedPubMed CentralGoogle Scholar
- Dunn S-J, Martello G, Yordanov B, Emmott S, Smith AG. Defining an essential transcription factor program for naïve pluripotency. Science. 2014;344:1156–60.View ArticlePubMedPubMed CentralGoogle Scholar
- Sampath P, Pritchard DK, Pabon L, Reinecke H, Schwartz SM, Morris DR, et al. A hierarchical network controls protein translation during murine embryonic stem cell self-renewal and differentiation. Cell Stem Cell. 2008;2:448–60.View ArticlePubMedGoogle Scholar
- Ying Q-L, Wray J, Nichols J, Batlle-Morera L, Doble B, Woodgett J, et al. The ground state of embryonic stem cell self-renewal. Nature. 2008;453:519–23.View ArticlePubMedGoogle Scholar
- Nelson PT, De Planell-Saguer M, Lamprinaki S, Kiriakidou M, Zhang P, O’Doherty U, et al. A novel monoclonal antibody against human Argonaute proteins reveals unexpected characteristics of miRNAs in human blood cells. RNA. 2007;13:1787–92.View ArticlePubMedPubMed CentralGoogle Scholar
- Smyth GK, Speed T. Normalization of cDNA microarray data. Methods. 2003;31:265–73.View ArticlePubMedGoogle Scholar
- Benaglia T, Chauveau D, Hunter DR, Young DS. mixtools: An R package for analyzing finite mixture models. J Stat Softw. 2009;32:1–29.View ArticleGoogle Scholar
- Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4:44–57.View ArticleGoogle Scholar
- Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013;14:R36.View ArticlePubMedPubMed CentralGoogle Scholar
- Anders S, Pyl PT, Huber W. HTSeq A Python framework to work with high-throughput sequencing data. bioRxiv. 2014;31:002824.Google Scholar
- Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550.View ArticlePubMedPubMed CentralGoogle Scholar
- Spandidos A, Wang X, Wang H, Seed B. PrimerBank: a resource of human and mouse PCR primer pairs for gene expression detection and quantification. Nucleic Acids Res. 2009;38:D792–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Nordgård O, Kvaløy JT, Farmen RK, Heikkilä R. Error propagation in relative real-time reverse transcription polymerase chain reaction quantification models: The balance between accuracy and precision. Anal Biochem. 2006;356:182–93.View ArticlePubMedGoogle Scholar
- Pfaffl MW, Horgan GW, Dempfle L. Relative expression software tool (REST) for group-wise comparison and statistical analysis of relative expression results in real-time PCR. Nucleic Acids Res. 2002;30:e36.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhu W, Smith JW, Huang CM. Mass spectrometry-based label-free quantitative proteomics. J Biomed Biotechnol. 2010;2010:1–6.Google Scholar
- Hotta A, Cheung AYL, Farra N, Vijayaragavan K, Séguin CA, Draper JS, et al. Isolation of human iPS cells using EOS lentiviral vectors to select for pluripotency. Nat Methods. 2009;6:370–6.View ArticlePubMedGoogle Scholar
- Sanjana NE, Shalem O, Zhang F. Improved vectors and genome-wide libraries for CRISPR screening. bioRxiv. 2014;11:006726.Google Scholar
- Zufferey R, Nagy D, Mandel RJ, Naldini L, Trono D. Multiply attenuated lentiviral vector achieves efficient gene delivery in vivo. Nat Biotechnol. 1997;15:871–5.View ArticlePubMedGoogle Scholar
- O’Doherty U, Swiggard WJ, Malim MH. Human immunodeficiency virus type 1 spinoculation enhances infection through virus binding. J Virol. 2000;74:10074–80.View ArticlePubMedPubMed CentralGoogle Scholar
- Hahne F, LeMeur N, Brinkman RR, Ellis B, Haaland P, Sarkar D, et al. flowCore: a Bioconductor package for high throughput flow cytometry. BMC Bioinformatics. 2009;10:106.View ArticlePubMedPubMed CentralGoogle Scholar
- Lienert F, Wirbelauer C, Som I, Dean A, Mohn F, Schübeler D. Identification of genetic elements that autonomously determine DNA methylation states. Nat Genet. 2011;43:1091–7.View ArticlePubMedGoogle Scholar
- Bock C, Reither S, Mikeska T, Paulsen M, Walter J, Lengauer T. BiQ Analyzer: Visualization and quality control for DNA methylation data from bisulfite sequencing. Bioinformatics. 2005;21:4067–8.View ArticlePubMedGoogle Scholar
- Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible read trimming tool for Illumina NGS data. Bioinformatics. 2014;30:2114–20.View ArticlePubMedPubMed CentralGoogle Scholar
- Landmead B, Trapnell C, Pot M, Salzberg S. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10:R25.View ArticleGoogle Scholar
- Kozomara A, Griffiths-Jones S. MiRBase: Annotating high confidence microRNAs using deep sequencing data. Nucleic Acids Res. 2014;42:D68–73.View ArticlePubMedPubMed CentralGoogle Scholar
- Liao Y, Smyth GK, Shi W. FeatureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014;30:923–30.View ArticlePubMedGoogle Scholar
- Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2009;26:139–40.View ArticlePubMedPubMed CentralGoogle Scholar
- Karolchik D, Hinrichs AS, Furey TS, Roskin KM, Sugnet CW, Haussler D, et al. The UCSC Table Browser data retrieval tool. Nucleic Acids Res. 2004;32(Database issue):D493–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Enright AJ, John B, Gaul U, Tuschl T, Sander C, Marks DS. MicroRNA targets in Drosophila. Genome Biol. 2003;5:R1.View ArticlePubMedPubMed CentralGoogle Scholar
- Plaisier CL, Bare JC, Baliga NS. MiRvestigator: web application to identify miRNAs responsible for co-regulated gene expression patterns discovered through transcriptome profiling. Nucleic Acids Res. 2011;39 Suppl 2:W125–31.View ArticlePubMedPubMed CentralGoogle Scholar