The transcription factor TCFL5 responds to A-MYB to elaborate the male meiotic program in mice.

In brief: The testis-specific transcription factor, TCFL5, expressed in pachytene spermatocytes regulates the meiotic gene expression program in collaboration with the transcription factor A-MYB. Abstract: In male mice, the transcription factors STRA8 and MEISON initiate meiosis I. We report that STRA8/MEISON activates the transcription factors A-MYB and TCFL5, which together reprogram gene expression after spermatogonia enter into meiosis. TCFL5 promotes the transcription of genes required for meiosis, mRNA turnover, miR-34/449 production, meiotic exit, and spermiogenesis. This transcriptional architecture is conserved in rhesus macaque, suggesting TCFL5 plays a central role in meiosis and spermiogenesis in placental mammals. Tcfl5em1/em1 mutants are sterile, and spermatogenesis arrests at the mid- or late-pachytene stage of meiosis. Moreover, Tcfl5+/em1 mutants produce fewer motile sperm.

A-MYB/TCFL5 regulatory architecture ensures the production of pachytene piRNAs in placental mammals.

In male mice, the transcription factor A MYB initiates the transcription of pachytene piRNA genes during meiosis. Here, we report that A MYB activates the transcription factor Tcfl5 produced in pachytene spermatocytes. Subsequently, A MYB and TCFL5 reciprocally reinforce their own transcription to establish a positive feedback circuit that triggers pachytene piRNA production. TCFL5 regulates the expression of genes required for piRNA maturation and promotes transcription of evolutionarily young pachytene piRNA genes, whereas A-MYB activates the transcription of older pachytene piRNA genes. Intriguingly, pachytene piRNAs from TCFL5-dependent young loci initiates the production of piRNAs from A-MYB-dependent older loci ensuring the self-propagation of pachytene piRNAs. A MYB and TCFL5 act via a set of incoherent feedforward loops that drive regulation of gene expression by pachytene piRNAs during spermatogenesis. This regulatory architecture is conserved in rhesus macaque, suggesting that it was present in the last common ancestor of placental mammals.

Estimating RNA dynamics using one time point for one sample in a single-pulse metabolic labeling experiment.

BACKGROUND: Over the past decade, experimental procedures such as metabolic labeling for determining RNA turnover rates at the transcriptome-wide scale have been widely adopted and are now turning to single cell measurements. Several computational methods to estimate RNA synthesis, processing and degradation rates from such experiments have been suggested, but they all require several RNA sequencing samples. Here we present a method that can estimate those three rates from a single sample. METHODS: Our method relies on the analytical solution to the Zeisel model of RNA dynamics. It was validated on metabolic labeling experiments performed on mouse embryonic stem cells. Resulting degradation rates were compared both to previously published rates on the same system and to a state-of-the-art method applied to the same data. RESULTS: Our method is computationally efficient and outputs rates that correlate well with previously published data sets. Using it on a single sample, we were able to reproduce the observation that dynamic biological processes tend to involve genes with higher metabolic rates, while stable processes involve genes with lower rates. This supports the hypothesis that cells control not only the mRNA steady-state abundance, but also its responsiveness, i.e., how fast steady state is reached. Moreover, degradation rates obtained with our method compare favourably with the other tested method. CONCLUSIONS: In addition to saving experimental work and computational time, estimating rates for a single sample has several advantages. It does not require an error-prone normalization across samples and enables the use of replicates to estimate uncertainty and assess sample quality. Finally the method and theoretical results described here are general enough to be useful in other contexts such as nucleotide conversion methods and single cell metabolic labeling experiments.

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.

The Contribution of lincRNAs at the Interface between Cell Cycle Regulation and Cell State Maintenance.

Cell cycle progression is controlled by the interplay of established cell cycle regulators. Changes in these regulators' activity underpin differences in cell cycle kinetics between cell types. We investigated whether long intergenic noncoding RNAs (lincRNAs) contribute to embryonic stem cell cycle adaptations. Using single-cell RNA sequencing data for mouse embryonic stem cells (mESCs) staged as G1, S, or G2/M we found differentially expressed lincRNAs are enriched among cell cycle-regulated genes. These lincRNAs (CC-lincRNAs) are co-expressed with genes involved in cell cycle regulation. We tested the impact of two CC-lincRNA candidates and show using CRISPR activation that increasing their expression is associated with deregulated cell cycle progression. Interestingly, CC-lincRNAs are often differentially expressed between G1 and S, their promoters are enriched in pluripotency transcription factor (TF) binding sites, and their transcripts are frequently co-regulated with genes involved in the maintenance of pluripotency, suggesting a contribution of CC-lincRNAs to mESC cell cycle adaptations.

A Protocol for Transcriptome-Wide Inference of RNA Metabolic Rates in Mouse Embryonic Stem Cells.

The relative ease of mouse Embryonic Stem Cells (mESCs) culture and the potential of these cells to differentiate into any of the three primary germ layers: ectoderm, endoderm and mesoderm (pluripotency), makes them an ideal and frequently used ex vivo system to dissect how gene expression changes impact cell state and differentiation. These efforts are further supported by the large number of constitutive and inducible mESC mutants established with the aim of assessing the contributions of different pathways and genes to cell homeostasis and gene regulation. Gene product abundance is controlled by the modulation of the rates of RNA synthesis, processing, and degradation. The ability to determine the relative contribution of these different RNA metabolic rates to gene expression control using standard RNA-sequencing approaches, which only capture steady state abundance of transcripts, is limited. In contrast, metabolic labeling of RNA with 4-thiouridine (4sU) coupled with RNA-sequencing, allows simultaneous and reproducible inference of transcriptome wide synthesis, processing, and degradation rates. Here we describe, a detailed protocol for 4sU metabolic labeling in mESCs that requires short 4sU labeling times at low concentration and minimally impacts cellular homeostasis. This approach presents a versatile method for in-depth characterization of the gene regulatory strategies governing gene steady state abundance in mESC.

Splicing of enhancer-associated lincRNAs contributes to enhancer activity.

Transcription is common at active mammalian enhancers sometimes giving rise to stable enhancer-associated long intergenic noncoding RNAs (elincRNAs). Expression of elincRNA is associated with changes in neighboring gene product abundance and local chromosomal topology, suggesting that transcription at these loci contributes to gene expression regulation in cis Despite the lack of evidence supporting sequence-dependent functions for most elincRNAs, splicing of these transcripts is unexpectedly common. Whether elincRNA splicing is a mere consequence of cognate enhancer activity or if it directly impacts enhancer function remains unresolved. Here, we investigate the association between elincRNA splicing and enhancer activity in mouse embryonic stem cells. We show that multi-exonic elincRNAs are enriched at conserved enhancers, and the efficient processing of elincRNAs is strongly associated with their cognate enhancer activity. This association is supported by their enrichment in enhancer-specific chromatin signatures; elevated binding of co-transcriptional regulators; increased local intra-chromosomal DNA contacts; and strengthened cis-regulation on target gene expression. Our results support the role of efficient RNA processing of enhancer-associated transcripts to cognate enhancer activity.

Dicer and Hsp104 function in a negative feedback loop to confer robustness to environmental stress.

Epigenetic mechanisms can be influenced by environmental cues and thus evoke phenotypic variation. This plasticity can be advantageous for adaptation but also detrimental if not tightly controlled. Although having attracted considerable interest, it remains largely unknown if and how environmental cues such as temperature trigger epigenetic alterations. Using fission yeast, we demonstrate that environmentally induced discontinuous phenotypic variation is buffered by a negative feedback loop that involves the RNase Dicer and the protein disaggregase Hsp104. In the absence of Hsp104, Dicer accumulates in cytoplasmic inclusions and heterochromatin becomes unstable at elevated temperatures, an epigenetic state inherited for many cell divisions after the heat stress. Loss of Dicer leads to toxic aggregation of an exogenous prionogenic protein. Our results highlight the importance of feedback regulation in building epigenetic memory and uncover Hsp104 and Dicer as homeostatic controllers that buffer environmentally induced stochastic epigenetic variation and toxic aggregation of prionogenic proteins.