The limits of human microRNA annotation have been met.

Over the last few years, the number of microRNAs in the human genome has become a controversially debated issue. Several publications reported thousands of putative novel microRNAs not included in the curated microRNA gene database MirGeneDB and the repository miRBase. Recently, by using sequencing of ∼300 human tissues and cell lines, the human RNA atlas, an expanded inventory of human RNA annotations, was published, reporting thousands of putative microRNAs. We, the developers of established microRNA prediction tools and hosts of MirGeneDB, raise concerns about the frequently applied prediction and functional validation strategies, briefly discussing the drawbacks of false positive detections. By means of quantifying well-established biogenesis-derived features, we show that the reported novel microRNAs essentially represent false-positives and argue that the human microRNA complement, at about 550 microRNA genes, is already near complete. Output of available tools must be curated as false predictions will misguide scientists looking for biomarkers or therapeutic targets.

MirGeneDB 2.1: toward a complete sampling of all major animal phyla.

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. MirGeneDB 2.1 is part of RNAcentral and Elixir Norway, publicly and freely available at http://www.mirgenedb.org/.

The SWI/SNF subunit BRG1 affects alternative splicing by changing RNA binding factor interactions with nascent RNA.

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.

Ancient microRNA profiles of a 14,300-year-old canid samples confirm taxonomic origin and give glimpses into tissue-specific gene regulation from the Pleistocene.

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 years, these microRNA sequences are by far the oldest ever reported. The authenticity of the sequences is further supported by a) the presence of canid / Caniformia-specific sequences that never evolved outside of this clade, b) tissue-specific expression patterns (cartilage, liver and muscle) that resemble those of modern dogs and c) 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. .

Next-generation sequencing reveals two populations of damage-induced small RNAs at endogenous DNA double-strand breaks.

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.

Malat1 affects transcription and splicing through distinct pathways in mouse embryonic stem cells.

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.

A combined approach for single-cell mRNA and intracellular protein expression analysis.

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.

Nuclear gene proximity and protein interactions shape transcript covariations in mammalian single cells.

Single-cell RNA sequencing studies on gene co-expression patterns could yield important regulatory and functional insights, but have so far been limited by the confounding effects of differentiation and cell cycle. We apply a tailored experimental design that eliminates these confounders, and report thousands of intrinsically covarying gene pairs in mouse embryonic stem cells. These covariations form a network with biological properties, outlining known and novel gene interactions. We provide the first evidence that miRNAs naturally induce transcriptome-wide covariations and compare the relative importance of nuclear organization, transcriptional and post-transcriptional regulation in defining covariations. We find that nuclear organization has the greatest impact, and that genes encoding for physically interacting proteins specifically tend to covary, suggesting importance for protein complex formation. Our results lend support to the concept of post-transcriptional RNA operons, but we further present evidence that nuclear proximity of genes may provide substantial functional regulation in mammalian single cells.

Circulating miRNA Spaceflight Signature Reveals Targets for Countermeasure Development.

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.

microRNA-205-5p is a modulator of insulin sensitivity that inhibits FOXO function.

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.

Control of hepatic gluconeogenesis by Argonaute2.

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 α1 (AMPKα1) by de-repressing its targeting by miR-148a, an abundant microRNA in the liver. Deletion of Ago2 from hyperglycemic, obese, and insulin-resistant Lepob/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 AMPKα1 by miR-148a, thereby identifying a regulatory link between non-coding RNAs and a key stress regulator in the hepatocyte.

Differential Impact of Glucose Administered Intravenously and Orally on Circulating miR-375 Levels in Human Subjects.

Background: To date, numerous nucleic acid species have been detected in the systemic circulation including microRNAs (miRNAs); however, their functional role in this compartment remains unclear. Objective: The aim of this study was to determine whether systemic levels of miRNAs abundant in blood, including the neuroendocrine tissue-enriched miR-375, are altered in response to a glucose challenge. Design: Twelve healthy males were recruited for an acute crossover study that consisted of two tests each following an 8-hour fasting period. An oral glucose tolerance test (OGTT) was performed, and blood samples were collected over a 3-hour period. Following a period of at least 1 week, the same participants were administered an isoglycemic intravenous glucose infusion (IIGI) with the same blood-collection protocol. Results: The glucose response curve following the IIGI mimicked that obtained after the OGTT, but as expected, systemic insulin levels were lower during the IIGI compared with the OGTT (P < 0.05). miR-375 levels in circulation were increased only in response to an OGTT and not during an IIGI. In addition, the response to the OGTT also coincided with the transient increase of circulating glucagon-like peptide (GLP)-1, GLP-2, and glucose-dependent insulinotropic polypeptide. Conclusions: The present findings show levels of miR-375 increase following administration of an OGTT and, in light of its enrichment in cells of the gut, suggest that the gastrointestinal tract may play an important role in the abundance and function of this miRNA in the blood.