Malat1 is a long-noncoding RNA with critical roles in gene regulation and cancer metastasis, however its functional role in stem cells is largely unexplored. We here perform a nuclear knockdown of Malat1 in mouse embryonic stem cells, causing the de-regulation of 320 genes and aberrant splicing of 90 transcripts, some of which potentially affecting the translated protein sequence. We find evidence that Malat1 directly interacts with gene bodies and aberrantly spliced transcripts, and that it locates upstream of down-regulated genes at their putative enhancer regions, in agreement with functional genomics data. Consistent with this, we find these genes affected at both exon and intron levels, suggesting that they are transcriptionally regulated by Malat1. Besides, the down-regulated genes are regulated by specific transcription factors and bear both activating and repressive chromatin marks, suggesting that some of them might be regulated by bivalent promoters. We propose a model in which Malat1 facilitates the transcription of genes involved in chromatid dynamics and mitosis in one pathway, and affects the splicing of transcripts that are themselves involved in RNA processing in a distinct pathway. Lastly, we compare our findings with Malat1 perturbation studies performed in other cell systems and in vivo.

We describe an update of MirGeneDB, the manually curated microRNA gene database. Adhering to uniform and consistent criteria for microRNA annotation and nomenclature, we substantially expanded MirGeneDB with 30 additional species representing previously missing metazoan phyla such as sponges, jellyfish, rotifers and flatworms. MirGeneDB 2.1 now consists of 75 species spanning over ∼800 million years of animal evolution, and contains a total number of 16 670 microRNAs from 1549 families. Over 6000 microRNAs were added in this update using ∼550 datasets with ∼7.5 billion sequencing reads. By adding new phylogenetically important species, especially those relevant for the study of whole genome duplication events, and through updating evolutionary nodes of origin for many families and genes, we were able to substantially refine our nomenclature system. All changes are traceable in the specifically developed MirGeneDB version tracker. The performance of read-pages is improved and microRNA expression matrices for all tissues and species are now also downloadable. Altogether, this update represents a significant step toward a complete sampling of all major metazoan phyla, and a widely needed foundation for comparative microRNA genomics and transcriptomics studies.

BRG1 and BRM are ATPase core subunits of the human SWI/SNF chromatin remodelling complexes mainly associated with transcriptional initiation. They also have a role in alternative splicing, which has been shown for BRM-containing SWI/SNF complexes at a few genes. Here, we have identified a subset of genes which harbour alternative exons that are affected by SWI/SNF ATPases by expressing the ATPases BRG1 and BRM in C33A cells, a BRG1- and BRM-deficient cell line, and analysed the effect on splicing by RNA sequencing. BRG1- and BRM-affected sub-sets of genes favouring both exon inclusion and exon skipping, with only a minor overlap between the ATPase. Some of the changes in alternative splicing induced by BRG1 and BRM expression did not require the ATPase activity. The BRG1-ATPase independent included exons displayed an exon signature of a high GC content. By investigating three genes with exons affected by the BRG-ATPase-deficient variant, we show that these exons accumulated phosphorylated RNA pol II CTD, both serine 2 and serine 5 phosphorylation, without an enrichment of the RNA polymerase II. The ATPases were recruited to the alternative exons, together with both core and signature subunits of SWI/SNF complexes, and promoted the binding of RNA binding factors to chromatin and RNA at the alternative exons. The interaction with the nascent RNP, however, did not reflect the association to chromatin. The hnRNPL, hnRNPU and SAM68 proteins associated with chromatin in cells expressing BRG1 and BRM wild type, but the binding of hnRNPU to the nascent RNP was excluded. This suggests that SWI/SNF can regulate alternative splicing by interacting with splicing-RNA binding factor and influence their binding to the nascent pre-mRNA particle.

Combined measurements of mRNA and protein expression in single cells enable in-depth analysis of cellular states. We present SPARC, an approach that combines single-cell RNA-sequencing with proximity extension essays to simultaneously measure global mRNA and 89 intracellular proteins in individual cells. We show that mRNA expression fails to accurately reflect protein abundance at the time of measurement, although the direction of changes is in agreement during neuronal differentiation. Moreover, protein levels of transcription factors better predict their downstream effects than do their corresponding transcripts. Finally, we highlight that protein expression variation is overall lower than mRNA variation, but relative protein variation does not reflect the mRNA level. Our results demonstrate that mRNA and protein measurements in single cells provide different and complementary information regarding cell states. SPARC presents a state-of-the-art co-profiling method that overcomes current limitations in throughput and protein localization, including removing the need for cell fixation.

DNA sequencing is the current key technology for historic or ancient biological samples and has led to many exciting discoveries in the field of paleogenomics. However, functional insights into tissue identity, cellular composition, or gene regulation cannot be gained from DNA. Recent analyses have shown that, under favorable conditions, RNA can also be sequenced from ancient samples, enabling studies at the transcriptomic and regulatory level. Analyzing ancient RNA data from a Pleistocene canid, we find hundreds of intact microRNAs that are taxonomically informative, show tissue specificity and have functionally predictive characteristics. With an extraordinary age of 14,300 yr, these microRNA sequences are by far the oldest ever reported. The authenticity of the sequences is further supported by (i) the presence of canid/Caniformia-specific sequences that never evolved outside of this Glade, (ii) tissue-specific expression patterns (cartilage, liver, and muscle) that resemble those of modern dogs, and (iii) RNA damage patterns that are clearly distinct from those of fresh samples. By performing computational microRNA-target enrichment analyses on the ancient sequences, we predict microRNA functions consistent with their tissue pattern of expression. For instance, we find a liver-specific microRNA that regulates carbohydrate metabolism and starvation responses in canids. In summary, we show that straightforward paleotranscriptomic microRNA analyses can give functional glimpses into tissue identity, cellular composition, and gene regulatory activity of ancient samples and biological processes that took place in the Pleistocene, thus holding great promise for deeper insights into gene regulation in extinct animals based on ancient RNA sequencing.

We have identified and validated a spaceflight-associated microRNA (miRNA) signature that is shared by rodents and humans in response to simulated, short-duration and long-duration spaceflight. Previous studies have identified miRNAs that regulate rodent responses to spaceflight in low-Earth orbit, and we have confirmed the expression of these proposed spaceflight-associated miRNAs in rodents reacting to simulated spaceflight conditions. Moreover, astronaut samples from the NASA Twins Study confirmed these expression signatures in miRNA sequencing, single-cell RNA sequencing (scRNA-seq), and single-cell assay for transposase accessible chromatin (scATAC-seq) data. Additionally, a subset of these miRNAs (miR-125, miR-16, and let-7a) was found to regulate vascular damage caused by simulated deep space radiation. To demonstrate the physiological relevance of key spaceflight-associated miRNAs, we utilized antagomirs to inhibit their expression and successfully rescue simulated deep-space-radiation-mediated damage in human 3D vascular constructs.

The ability to profile transcriptomes and proteomes in a high-throughput fashion in single cells has truly revolutionized functional genomics, and countless functional and regulatory insights have been based on these technologies. While major applications include the discovery of new cell types and the a posteriori sorting of cell populations, studies of gene expression noise and gene co-expression have made use of this inter-cellular heterogeneity in a genuine quantitative fashion. Yet, there are still major limitations to overcome.

First, strong dynamic processes, such as cell cycle or differentiation axes, tend to overshadow more subtle underlying regulatory processes. While this has sparked the development of tools that can identify and correct these biases at large, few insights into the subtleties of gene regulation have been published thus far. The majority of studies still focus on drastic changes such as differentiation or disease. We address this issue in paper I and to a limited extend in paper II and paper III through the elimination of major confounders during experimental design. In these papers, we show that variation and covariation of miRNAs, mRNAs and proteins between individual cells of a homogeneous non-dynamic population are informative of gene regulation.

Second, while single-cell technologies are booming, with new technologies being published every day, the co-profiling of RNA and protein in the same single cells still remains a major challenge. All current technologies are limited either by protein location or throughput, or require invasive cell fixation that can compromise mRNA stability. We overcome these limitations in paper II through the combination of quantitative single-cell RNA sequencing with proximity extension assays for protein detection. Using this technology, SPARC, we show that transcription factor protein, but not transcription factor RNA, covaries with the RNA expression of its targets. We also show that translation is a major mediator of the shift in variation from the RNA to the protein level.

Third, some technologies still suffer from limited sensitivity. While, for instance, the first single-cell miRNA detection already succeeded in 2006 and the first single-cell small RNA sequencing technique was published in 2016, few insights into miRNA dynamics or function have been gained from single-cell data since. Using an optimized single-cell small RNA sequencing protocol, we quantify the miRNA transcriptome of close to 200 single cells in paper III. We show that variation and covariation can be linked to miRNA transcription and turnover. Integrating miRNA and miRNA target data from all three papers, we present evidence that the induction of variation on the RNA level and the buffering of protein expression noise are naturally occurring for many miRNAs.

In summary, we present new strategies and new protocols that overcome existing limitations in the field, and we present regulatory insights that were enabled by quantitative measurements of single-cell gene expression variation and covariation.

Objective: The liver performs a central role in regulating energy homeostasis by increasing glucose output during fasting. Recent studies on Argonaute2 (Ago2), a key RNA-binding protein mediating the microRNA pathway, have illustrated its role in adaptive mechanisms according to changes in metabolic demand. Here we sought to characterize the functional role of Ago2 in the liver in the maintenance of systemic glucose homeostasis. Methods: We first analyzed Ago2 expression in mouse primary hepatocyte cultures after modulating extracellular glucose concentrations and in the presence of activators or inhibitors of glucokinase activity. We then characterized a conditional loss-of-function mouse model of Ago2 in liver for alterations in systemic energy metabolism. Results: Here we show that Ago2 expression in liver is directly correlated to extracellular glucose concentrations and that modulating glucokinase activity is adequate to affect hepatic Ago2 levels. Conditional deletion of Ago2 in liver resulted in decreased fasting glucose levels in addition to reducing hepatic glucose production. Moreover, loss of Ago2 promoted hepatic expression of AMP-activated protein kinase alpha 1 (AMPK alpha 1) by de-repressing its targeting by miR-148a, an abundant microRNA in the liver. Deletion of Ago2 from hyperglycemic, obese, and insulin-resistant Lep(ob/ob) mice reduced both random and fasted blood glucose levels and body weight and improved insulin sensitivity. Conclusions: These data illustrate a central role for Ago2 in the adaptive response of the liver to fasting. Ago2 mediates the suppression of AMPKa1 by miR-148a, thereby identifying a regulatory link between non-coding RNAs and a key stress regulator in the hepatocyte.

Objectives: Hepatic insulin resistance is a hallmark of type 2 diabetes and obesity. Insulin receptor signaling through AKT and FOXO has important metabolic effects that have traditionally been ascribed to regulation of gene expression. However, whether all the metabolic effects of FOXO arise from its regulation of protein-encoding mRNAs is unknown. Methods: To address this question, we obtained expression profiles of FOXO-regulated murine hepatic microRNAs (miRNAs) during fasting and refeeding using mice lacking Foxo1, 3a, and 4 in liver (L-Foxo1,3a, 4). Results: Out of 439 miRNA analyzed, 175 were differentially expressed in Foxo knockouts. Their functions were associated with insulin, Wnt, Mapk signaling, and aging. Among them, we report a striking increase of miR-205-5p expression in L-Foxo1,3a,4 knockouts, as well as in obese mice. We show that miR-205-5p gain-of-function increases AKT phosphorylation and decreases SHIP2 in primary hepatocytes, resulting in FOXO inhibition. This results in decreased hepatocyte glucose production. Consistent with these observations, miR-205-5p gain-of-function in mice lowered glucose levels and improved pyruvate tolerance. Conclusions: These findings reveal a homeostatic miRNA loop regulating insulin signaling, with potential implications for in vivo glucose metabolism.

Recent studies suggest that transcription takes place at DNA double-strand breaks (DSBs), that transcripts at DSBs are processed by Drosha and Dicer into damage-induced small RNAs (diRNAs), and that diRNAs are required for DNA repair. However, diRNAs have been mostly detected in reporter constructs or repetitive sequences, and their existence at endogenous loci has been questioned by recent reports. Using the homing endonuclease I-PpoI, we have investigated diRNA production in genetically unperturbed human and mouse cells. I-PpoI is an ideal tool to clarify the requirements for diRNA production because it induces DSBs in different types of loci: the repetitive 28S locus, unique genes and intergenic loci. We show by extensive sequencing that the rDNA locus produces substantial levels of diRNAs, whereas unique genic and intergenic loci do not. Further characterization of diRNAs emerging from the 28S locus reveals the existence of two diRNA subtypes. Surprisingly, Drosha and its partner DGCR8 are dispensable for diRNA production and only one diRNAs subtype depends on Dicer processing. Furthermore, we provide evidence that diRNAs are incorporated into Argonaute. Our findings provide direct evidence for diRNA production at endogenous loci in mammalian cells and give insights into RNA processing at DSBs.