Translation is required for miRNA-dependent decay of endogenous transcripts.

Post-transcriptional repression of gene expression by miRNAs occurs through transcript destabilization or translation inhibition. mRNA decay is known to account for most miRNA-dependent repression. However, because transcript decay occurs co-translationally, whether target translation is a requirement for miRNA-dependent transcript destabilization remains unknown. To decouple these two molecular processes, we used cytosolic long noncoding RNAs (lncRNAs) as models for endogenous transcripts that are not translated. We show that, despite interacting with the miRNA-loaded RNA-induced silencing complex, the steady-state abundance and decay rates of these transcripts are minimally affected by miRNA loss. To further validate the apparent requirement of translation for miRNA-dependent decay, we fused two lncRNA candidates to the 3'-end of a protein-coding gene reporter and found this results in their miRNA-dependent destabilization. Further analysis revealed that the few natural lncRNAs whose levels are regulated by miRNAs in mESCs tend to associate with translating ribosomes, and possibly represent misannotated micropeptides, further substantiating the necessity of target translation for miRNA-dependent transcript decay. In summary, our analyses suggest that translation is required for miRNA-dependent transcript destabilization, and demonstrate that the levels of coding and noncoding transcripts are differently affected by miRNAs.

LINE-1 activity in facultative heterochromatin formation during X chromosome inactivation.

During X chromosome inactivation (XCI), Xist RNA coats and silences one of the two X chromosomes in female cells. Little is known about how XCI spreads across the chromosome, although LINE-1 elements have been proposed to play a role. Here we show that LINEs participate in creating a silent nuclear compartment into which genes become recruited. A subset of young LINE-1 elements, however, is expressed during XCI, rather than being silenced. We demonstrate that such LINE expression requires the specific heterochromatic state induced by Xist. These LINEs often lie within escape-prone regions of the X chromosome, but close to genes that are subject to XCI, and are associated with putative endo-siRNAs. LINEs may thus facilitate XCI at different levels, with silent LINEs participating in assembly of a heterochromatic nuclear compartment induced by Xist, and active LINEs participating in local propagation of XCI into regions that would otherwise be prone to escape.

Sequestration of LINE-1 in cytosolic aggregates by MOV10 restricts retrotransposition.

LINE-1 (L1) retroelements have retained their ability to mobilize. Mechanisms regulating L1 mobility include DNA methylation in somatic cells and the piRNA pathway in the germline. During preimplantation stages of mouse embryonic development, however, both pathways are inactivated leading to a window necessitating alternate means of L1 regulation. We previously reported an increase in L1 levels in Dicer_KO mouse embryonic stem cells (mESCs), which was accompanied by only a marginal increase in retrotransposition, suggesting additional mechanisms suppressing L1 mobility. Here, we demonstrate that L1 ribonucleoprotein complexes (L1 RNP) accumulate as aggregates in the cytoplasm of Dicer_KO mESCs along with the RNA helicase MOV10. The combined overexpression of L1 ORF1p and MOV10 is sufficient to create L1 RNP aggregates. In Dicer_KO mESCs, MOV10 is upregulated due to the loss of its direct regulation by miRNAs. The newly discovered posttranscriptional regulation of Mov10, and its role in preventing L1 retrotransposition by driving cytosolic aggregation, provides routes to explore for therapy in disease conditions where L1s are upregulated.

Global and precise identification of functional miRNA targets in mESCs by integrative analysis.

MicroRNA (miRNA) loaded Argonaute (AGO) complexes regulate gene expression via direct base pairing with their mRNA targets. Previous works suggest that up to 60% of mammalian transcripts might be subject to miRNA-mediated regulation, but it remains largely unknown which fraction of these interactions are functional in a specific cellular context. Here, we integrate transcriptome data from a set of miRNA-depleted mouse embryonic stem cell (mESC) lines with published miRNA interaction predictions and AGO-binding profiles. Using this integrative approach, combined with molecular validation data, we present evidence that < 10% of expressed genes are functionally and directly regulated by miRNAs in mESCs. In addition, analyses of the stem cell-specific miR-290-295 cluster target genes identify TFAP4 as an important transcription factor for early development. The extensive datasets developed in this study will support the development of improved predictive models for miRNA-mRNA functional interactions.

Comparative analysis of differential gene expression tools for RNA sequencing time course data.

RNA sequencing (RNA-seq) has become a standard procedure to investigate transcriptional changes between conditions and is routinely used in research and clinics. While standard differential expression (DE) analysis between two conditions has been extensively studied, and improved over the past decades, RNA-seq time course (TC) DE analysis algorithms are still in their early stages. In this study, we compare, for the first time, existing TC RNA-seq tools on an extensive simulation data set and validated the best performing tools on published data. Surprisingly, TC tools were outperformed by the classical pairwise comparison approach on short time series (<8 time points) in terms of overall performance and robustness to noise, mostly because of high number of false positives, with the exception of ImpulseDE2. Overlapping of candidate lists between tools improved this shortcoming, as the majority of false-positive, but not true-positive, candidates were unique for each method. On longer time series, pairwise approach was less efficient on the overall performance compared with splineTC and maSigPro, which did not identify any false-positive candidate.

Embryonic stem cell-specific microRNAs contribute to pluripotency by inhibiting regulators of multiple differentiation pathways.

The findings that microRNAs (miRNAs) are essential for early development in many species and that embryonic miRNAs can reprogram somatic cells into induced pluripotent stem cells suggest that these miRNAs act directly on transcriptional and chromatin regulators of pluripotency. To elucidate the transcription regulatory networks immediately downstream of embryonic miRNAs, we extended the motif activity response analysis approach that infers the regulatory impact of both transcription factors (TFs) and miRNAs from genome-wide expression states. Applying this approach to multiple experimental data sets generated from mouse embryonic stem cells (ESCs) that did or did not express miRNAs of the ESC-specific miR-290-295 cluster, we identified multiple TFs that are direct miRNA targets, some of which are known to be active during cell differentiation. Our results provide new insights into the transcription regulatory network downstream of ESC-specific miRNAs, indicating that these miRNAs act on cell cycle and chromatin regulators at several levels and downregulate TFs that are involved in the innate immune response.

RNAi-dependent and independent control of LINE1 accumulation and mobility in mouse embryonic stem cells.

In most mouse tissues, long-interspersed elements-1 (L1s) are silenced via methylation of their 5'-untranslated regions (5'-UTR). A gradual loss-of-methylation in pre-implantation embryos coincides with L1 retrotransposition in blastocysts, generating potentially harmful mutations. Here, we show that Dicer- and Ago2-dependent RNAi restricts L1 accumulation and retrotransposition in undifferentiated mouse embryonic stem cells (mESCs), derived from blastocysts. RNAi correlates with production of Dicer-dependent 22-nt small RNAs mapping to overlapping sense/antisense transcripts produced from the L1 5'-UTR. However, RNA-surveillance pathways simultaneously degrade these transcripts and, consequently, confound the anti-L1 RNAi response. In Dicer(-/-) mESC complementation experiments involving ectopic Dicer expression, L1 silencing was rescued in cells in which microRNAs remained strongly depleted. Furthermore, these cells proliferated and differentiated normally, unlike their non-complemented counterparts. These results shed new light on L1 biology, uncover defensive, in addition to regulatory roles for RNAi, and raise questions on the differentiation defects of Dicer(-/-) mESCs.

Naive and primed murine pluripotent stem cells have distinct miRNA expression profiles.

Over the last years, the microRNA (miRNA) pathway has emerged as a key component of the regulatory network of pluripotency. Although clearly distinct states of pluripotency have been described in vivo and ex vivo, differences in miRNA expression profiles associated with the developmental modulation of pluripotency have not been extensively studied so far. Here, we performed deep sequencing to profile miRNA expression in naive (embryonic stem cell [ESC]) and primed (epiblast stem cell [EpiSC]) pluripotent stem cells derived from mouse embryos of identical genetic background. We developed a graphical representation method allowing the rapid identification of miRNAs with an atypical profile including mirtrons, a small nucleolar RNA (snoRNA)-derived miRNA, and miRNAs whose biogenesis may differ between ESC and EpiSC. Comparison of mature miRNA profiles revealed that ESCs and EpiSCs exhibit very different miRNA signatures with one third of miRNAs being differentially expressed between the two cell types. Notably, differential expression of several clusters, including miR290-295, miR17-92, miR302/367, and a large repetitive cluster on chromosome 2, was observed. Our analysis also showed that differentiation priming of EpiSC compared to ESC is evidenced by changes in miRNA expression. These dynamic changes in miRNAs signature are likely to reflect both redundant and specific roles of miRNAs in the fine-tuning of pluripotency during development.

ncPRO-seq: a tool for annotation and profiling of ncRNAs in sRNA-seq data.

SUMMARY: Non-coding RNA (ncRNA) PROfiling in small RNA (sRNA)-seq (ncPRO-seq) is a stand-alone, comprehensive and flexible ncRNA analysis pipeline. It can interrogate and perform detailed profiling analysis on sRNAs derived from annotated non-coding regions in miRBase, Rfam and RepeatMasker, as well as specific regions defined by users. The ncPRO-seq pipeline performs both gene-based and family-based analyses of sRNAs. It also has a module to identify regions significantly enriched with short reads, which cannot be classified under known ncRNA families, thus enabling the discovery of previously unknown ncRNA- or small interfering RNA (siRNA)-producing regions. The ncPRO-seq pipeline supports input read sequences in fastq, fasta and color space format, as well as alignment results in BAM format, meaning that sRNA raw data from the three current major platforms (Roche-454, Illumina-Solexa and Life technologies-SOLiD) can be analyzed with this pipeline. The ncPRO-seq pipeline can be used to analyze read and alignment data, based on any sequenced genome, including mammals and plants. AVAILABILITY: Source code, annotation files, manual and online version are available at http://ncpro.curie.fr/. CONTACT: bioinfo.ncproseq@curie.fr or cciaudo@ethz.ch SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Highly dynamic and sex-specific expression of microRNAs during early ES cell differentiation.

Embryonic stem (ES) cells are pluripotent cells derived from the inner cell mass of the mammalian blastocyst. Cellular differentiation entails loss of pluripotency and gain of lineage-specific characteristics. However, the molecular controls that govern the differentiation process remain poorly understood. We have characterized small RNA expression profiles in differentiating ES cells as a model for early mammalian development. High-throughput 454 pyro-sequencing was performed on 19-30 nt RNAs isolated from undifferentiated male and female ES cells, as well as day 2 and 5 differentiating derivatives. A discrete subset of microRNAs (miRNAs) largely dominated the small RNA repertoire, and the dynamics of their accumulation could be readily used to discriminate pluripotency from early differentiation events. Unsupervised partitioning around meloids (PAM) analysis revealed that differentiating ES cell miRNAs can be divided into three expression clusters with highly contrasted accumulation patterns. PAM analysis afforded an unprecedented level of definition in the temporal fluctuations of individual members of several miRNA genomic clusters. Notably, this unravelled highly complex post-transcriptional regulations of the key pluripotency miR-290 locus, and helped identify miR-293 as a clear outlier within this cluster. Accordingly, the miR-293 seed sequence and its predicted cellular targets differed drastically from those of the other abundant cluster members, suggesting that previous conclusions drawn from whole miR-290 over-expression need to be reconsidered. Our analysis in ES cells also uncovered a striking male-specific enrichment of the miR-302 family, which share the same seed sequence with most miR-290 family members. Accordingly, a miR-302 representative was strongly enriched in embryonic germ cells derived from primordial germ cells of male but not female mouse embryos. Identifying the chromatin remodelling and E2F-dependent transcription repressors Ari4a and Arid4b as additional targets of miR-302 and miR-290 supports and possibly expands a model integrating possible overlapping functions of the two miRNA families in mouse cell totipotency during early development. This study demonstrates that small RNA sampling throughout early ES cell differentiation enables the definition of statistically significant expression patterns for most cellular miRNAs. We have further shown that the transience of some of these miRNA patterns provides highly discriminative markers of particular ES cell states during their differentiation, an approach that might be broadly applicable to the study of early mammalian development.

Argonaute proteins regulate a specific network of genes through KLF4 in mouse embryonic stem cells.

The Argonaute proteins (AGOs) are well known for their role in post-transcriptional gene silencing in the microRNA (miRNA) pathway. Here we show that in mouse embryonic stem cells, AGO1&2 serve additional functions that go beyond the miRNA pathway. Through the combined deletion of both Agos, we identified a specific set of genes that are uniquely regulated by AGOs but not by the other miRNA biogenesis factors. Deletion of Ago2&1 caused a global reduction of the repressive histone mark H3K27me3 due to downregulation at protein levels of Polycomb repressive complex 2 components. By integrating chromatin accessibility, prediction of transcription factor binding sites, and chromatin immunoprecipitation sequencing data, we identified the pluripotency factor KLF4 as a key modulator of AGO1&2-regulated genes. Our findings revealed a novel axis of gene regulation that is mediated by noncanonical functions of AGO proteins that affect chromatin states and gene expression using mechanisms outside the miRNA pathway.

AGO1 regulates pericentromeric regions in mouse embryonic stem cells.

Argonaute proteins (AGOs), which play an essential role in cytosolic post-transcriptional gene silencing, have been also reported to function in nuclear processes like transcriptional activation or repression, alternative splicing and, chromatin organization. As most of these studies have been conducted in human cancer cell lines, the relevance of AGOs nuclear functions in the context of mouse early embryonic development remains uninvestigated. Here, we examined a possible role of the AGO1 protein on the distribution of constitutive heterochromatin in mouse embryonic stem cells (mESCs). We observed a specific redistribution of the repressive histone mark H3K9me3 and the heterochromatin protein HP1α, away from pericentromeric regions upon Ago1 depletion. Furthermore, we demonstrated that major satellite transcripts are strongly up-regulated in Ago1_KO mESCs and that their levels are partially restored upon AGO1 rescue. We also observed a similar redistribution of H3K9me3 and HP1α in Drosha_KO mESCs, suggesting a role for microRNAs (miRNAs) in the regulation of heterochromatin distribution in mESCs. Finally, we showed that specific miRNAs with complementarity to major satellites can partially regulate the expression of these transcripts.

The solution structure of Dead End bound to AU-rich RNA reveals an unusual mode of tandem RRM-RNA recognition required for mRNA regulation.

Dead End (DND1) is an RNA-binding protein essential for germline development through its role in post-transcriptional gene regulation. The molecular mechanisms behind selection and regulation of its targets are unknown. Here, we present the solution structure of DND1's tandem RNA Recognition Motifs (RRMs) bound to AU-rich RNA. The structure reveals how an NYAYUNN element is specifically recognized, reconciling seemingly contradictory sequence motifs discovered in recent genome-wide studies. RRM1 acts as a main binding platform, including atypical extensions to the canonical RRM fold. RRM2 acts cooperatively with RRM1, capping the RNA using an unusual binding pocket, leading to an unusual mode of tandem RRM-RNA recognition. We show that the consensus motif is sufficient to mediate upregulation of a reporter gene in human cells and that this process depends not only on RNA binding by the RRMs, but also on DND1's double-stranded RNA binding domain (dsRBD), which is dispensable for binding of a subset of targets in cellulo. Our results point to a model where DND1 target selection is mediated by a non-canonical mode of AU-rich RNA recognition by the tandem RRMs and a role for the dsRBD in the recruitment of effector complexes responsible for target regulation.

Girafe--an R/Bioconductor package for functional exploration of aligned next-generation sequencing reads.

UNLABELLED: The R/Bioconductor package girafe facilitates the functional exploration of alignments of sequence reads from next-generation sequencing data to a genome. It allows users to investigate the genomic intervals together with the aligned reads and to work with, visualise and export these intervals. Moreover, the package operates within and extends the ever-growing Bioconductor framework and thus enables users to leverage a multitude of methods for their data in order to answer specific research questions. AVAILABILITY AND IMPLEMENTATION: The R package girafe is available from the Bioconductor web site: http://www.bioconductor.org/packages/release/bioc/html/girafe.html. An extensive vignette and the Bioconductor mailing lists provide additional documentation and help for using the package.

Prediction of the miRNA interactome - Established methods and upcoming perspectives.

MicroRNAs (miRNAs) are well-studied small noncoding RNAs involved in post-transcriptional gene regulation in a wide range of organisms, including mammals. Their function is mediated by base pairing with their target RNAs. Although many features required for miRNA-mediated repression have been described, the identification of functional interactions is still challenging. In the last two decades, numerous Machine Learning (ML) models have been developed to predict their putative targets. In this review, we summarize the biological knowledge and the experimental data used to develop these ML models. Recently, Deep Neural Network-based models have also emerged in miRNA interaction modeling. We thus outline established and emerging models to give a perspective on the future developments needed to improve the identification of genes directly regulated by miRNAs.

Fast In Vitro Procedure to Identify Extraembryonic Differentiation Defect of Mouse Embryonic Stem Cells.

Mouse embryonic stem cells (mESCs) are a powerful model to study early mouse development. These blastocyst-derived cells can maintain pluripotency and differentiate into the three embryonic germ layers and an extraembryonic layer, the extraembryonic endoderm (ExEn), which shares similar molecular markers to the definitive endoderm. Here, we present a fast procedure to identify a differentiation defect of mESCs toward ExEn in vitro through the molecular and cellular characterization of embryoid bodies, followed by direct differentiation of mESCs into ExEn. For complete details on the use and execution of this protocol, please refer to Ngondo et al. (2018).

Inhibition of FGF and TGF-ß Pathways in hESCs Identify STOX2 as a Novel SMAD2/4 Cofactor.

The fibroblast growth factor (FGF) and the transforming growth factor-ß (TGF-ß) pathways are both involved in the maintenance of human embryonic stem cells (hESCs) and regulate the onset of their differentiation. Their converging functions have suggested that these pathways might share a wide range of overlapping targets. Published studies have focused on the long-term effects (24-48 h) of FGF and TGF-ß inhibition in hESCs, identifying direct and indirect target genes. In this study, we focused on the earliest transcriptome changes occurring between 3 and 9 h after FGF and TGF-ß inhibition to identify direct target genes only. Our analysis clearly shows that only a handful of target transcripts are common to both pathways. This is surprising in light of the previous literature, and has implications for models of cell signaling in human pluripotent cells. In addition, we identified STOX2 as a novel primary target of the TGF-ß signaling pathway. We show that STOX2 might act as a novel SMAD2/4 cofactor. Taken together, our results provide insights into the effect of cell signaling on the transcription profile of human pluripotent cells.