The MOV10 RNA helicase is a dosage-dependent host restriction factor for LINE1 retrotransposition in mice.

Transposable elements constitute nearly half of the mammalian genome and play important roles in genome evolution. While a multitude of both transcriptional and post-transcriptional mechanisms exist to silence transposable elements, control of transposition in vivo remains poorly understood. MOV10, an RNA helicase, is an inhibitor of mobilization of retrotransposons and retroviruses in cell culture assays. Here we report that MOV10 restricts LINE1 retrotransposition in mice. Although MOV10 is broadly expressed, its loss causes only incomplete penetrance of embryonic lethality, and the surviving MOV10-deficient mice are healthy and fertile. Biochemically, MOV10 forms a complex with UPF1, a key component of the nonsense-mediated mRNA decay pathway, and primarily binds to the 3' UTR of somatically expressed transcripts in testis. Consequently, loss of MOV10 results in an altered transcriptome in testis. Analyses using a LINE1 reporter transgene reveal that loss of MOV10 leads to increased LINE1 retrotransposition in somatic and reproductive tissues from both embryos and adult mice. Moreover, the degree of LINE1 retrotransposition inhibition is dependent on the Mov10 gene dosage. Furthermore, MOV10 deficiency reduces reproductive fitness over successive generations. Our findings demonstrate that MOV10 attenuates LINE1 retrotransposition in a dosage-dependent manner in mice.

TERA-Seq: true end-to-end sequencing of native RNA molecules for transcriptome characterization.

Direct sequencing of single, native RNA molecules through nanopores has a strong potential to transform research in all aspects of RNA biology and clinical diagnostics. The existing platform from Oxford Nanopore Technologies is unable to sequence the very 5' ends of RNAs and is limited to polyadenylated molecules. Here, we develop True End-to-end RNA Sequencing (TERA-Seq), a platform that addresses these limitations, permitting more thorough transcriptome characterization. TERA-Seq describes both poly- and non-polyadenylated RNA molecules and accurately identifies their native 5' and 3' ends by ligating uniquely designed adapters that are sequenced along with the transcript. We find that capped, full-length mRNAs in human cells show marked variation of poly(A) tail lengths at the single molecule level. We report prevalent capping downstream of canonical transcriptional start sites in otherwise fully spliced and polyadenylated molecules. We reveal RNA processing and decay at single molecule level and find that mRNAs decay cotranslationally, often from their 5' ends, while frequently retaining poly(A) tails. TERA-Seq will prove useful in many applications where true end-to-end direct sequencing of single, native RNA molecules and their isoforms is desirable.

The RNA helicase MOV10L1 binds piRNA precursors to initiate piRNA processing.

Piwi-piRNA (Piwi-interacting RNA) ribonucleoproteins (piRNPs) enforce retrotransposon silencing, a function critical for preserving the genome integrity of germ cells. The molecular functions of most of the factors that have been genetically implicated in primary piRNA biogenesis are still elusive. Here we show that MOV10L1 exhibits 5'-to-3' directional RNA-unwinding activity in vitro and that a point mutation that abolishes this activity causes a failure in primary piRNA biogenesis in vivo. We demonstrate that MOV10L1 selectively binds piRNA precursor transcripts and is essential for the generation of intermediate piRNA processing fragments that are subsequently loaded to Piwi proteins. Multiple analyses suggest an intimate coupling of piRNA precursor processing with elements of local secondary structures such as G quadruplexes. Our results support a model in which MOV10L1 RNA helicase activity promotes unwinding and funneling of the single-stranded piRNA precursor transcripts to the endonuclease that catalyzes the first cleavage step of piRNA processing.

Sequence-dependent but not sequence-specific piRNA adhesion traps mRNAs to the germ plasm.

The conserved Piwi family of proteins and piwi-interacting RNAs (piRNAs) have a central role in genomic stability, which is inextricably linked to germ-cell formation, by forming Piwi ribonucleoproteins (piRNPs) that silence transposable elements. In Drosophila melanogaster and other animals, primordial germ-cell specification in the developing embryo is driven by maternal messenger RNAs and proteins that assemble into specialized messenger ribonucleoproteins (mRNPs) localized in the germ (pole) plasm at the posterior of the oocyte. Maternal piRNPs, especially those loaded on the Piwi protein Aubergine (Aub), are transmitted to the germ plasm to initiate transposon silencing in the offspring germ line. The transport of mRNAs to the oocyte by midoogenesis is an active, microtubule-dependent process; mRNAs necessary for primordial germ-cell formation are enriched in the germ plasm at late oogenesis via a diffusion and entrapment mechanism, the molecular identity of which remains unknown. Aub is a central component of germ granule RNPs, which house mRNAs in the germ plasm, and interactions between Aub and Tudor are essential for the formation of germ granules. Here we show that Aub-loaded piRNAs use partial base-pairing characteristics of Argonaute RNPs to bind mRNAs randomly in Drosophila, acting as an adhesive trap that captures mRNAs in the germ plasm, in a Tudor-dependent manner. Notably, germ plasm mRNAs in drosophilids are generally longer and more abundant than other mRNAs, suggesting that they provide more target sites for piRNAs to promote their preferential tethering in germ granules. Thus, complexes containing Tudor, Aub piRNPs and mRNAs couple piRNA inheritance with germline specification. Our findings reveal an unexpected function for piRNP complexes in mRNA trapping that may be generally relevant to the function of animal germ granules.

A MicroRNA precursor surveillance system in quality control of MicroRNA synthesis.

MicroRNAs (miRNAs) are essential for regulation of gene expression. Though numerous miRNAs have been identified by high-throughput sequencing, few precursor miRNAs (pre-miRNAs) are experimentally validated. Here we report a strategy for constructing high-throughput sequencing libraries enriched for full-length pre-miRNAs. We find widespread and extensive uridylation of Argonaute (Ago)-bound pre-miRNAs, which is primarily catalyzed by two terminal uridylyltransferases: TUT7 and TUT4. Uridylation by TUT7/4 not only polishes pre-miRNA 3' ends, but also facilitates their degradation by the exosome, preventing clogging of Ago with defective species. We show that the exosome exploits distinct substrate preferences of DIS3 and RRP6, its two catalytic subunits, to distinguish productive from defective pre-miRNAs. Furthermore, we identify a positive feedback loop formed by the exosome and TUT7/4 in triggering uridylation and degradation of Ago-bound pre-miRNAs. Our study reveals a pre-miRNA surveillance system that comprises TUT7, TUT4, and the exosome in quality control of miRNA synthesis.

SamQL: a structured query language and filtering tool for the SAM/BAM file format.

BACKGROUND: The Sequence Alignment/Map Format Specification (SAM) is one of the most widely adopted file formats in bioinformatics and many researchers use it daily. Several tools, including most high-throughput sequencing read aligners, use it as their primary output and many more tools have been developed to process it. However, despite its flexibility, SAM encoded files can often be difficult to query and understand even for experienced bioinformaticians. As genomic data are rapidly growing, structured, and efficient queries on data that are encoded in SAM/BAM files are becoming increasingly important. Existing tools are very limited in their query capabilities or are not efficient. Critically, new tools that address these shortcomings, should not be able to support existing large datasets but should also do so without requiring massive data transformations and file infrastructure reorganizations. RESULTS: Here we introduce SamQL, an SQL-like query language for the SAM format with intuitive syntax that supports complex and efficient queries on top of SAM/BAM files and that can replace commonly used Bash one-liners employed by many bioinformaticians. SamQL has high expressive power with no upper limit on query size and when parallelized, outperforms other substantially less expressive software. CONCLUSIONS: SamQL is a complete query language that we envision as a step to a structured database engine for genomics. SamQL is written in Go, and is freely available as standalone program and as an open-source library under an MIT license, https://github.com/maragkakislab/samql/ .

Kc167, a widely used Drosophila cell line, contains an active primary piRNA pathway.

PIWI family proteins bind to small RNAs known as PIWI-interacting RNAs (piRNAs) and play essential roles in the germline by silencing transposons and by promoting germ cell specification and function. Here we report that the widely used Kc167 cell line, derived from Drosophila melanogaster embryos, expresses piRNAs that are loaded to Aub and Piwi. Kc167 piRNAs are produced by a canonical, primary piRNA biogenesis pathway, from phased processing of precursor transcripts by the Zuc endonuclease, Armi helicase, and dGasz mitochondrial scaffold protein. Kc167 piRNAs derive from cytoplasmic transcripts, notably tRNAs and mRNAs, and their abundance correlates with that of parent transcripts. The expression of Aub is robust in Kc167, that of Piwi is modest, while Ago3 is undetectable, explaining the lack of transposon-related piRNA amplification by the Aub-Ago3, ping-pong mechanism. We propose that the default state of the primary piRNA biogenesis machinery is random transcript sampling to allow generation of piRNAs from any transcript, including newly acquired retrotransposons. This state is unmasked in Kc167, likely because they do not express piRNA cluster transcripts in sufficient amounts and do not amplify transposon piRNAs. We use Kc167 to characterize an inactive isoform of Aub protein. Since most Kc167 piRNAs are genic, they can be mapped uniquely to the genome, facilitating computational analyses. Furthermore, because Kc167 is a widely used and well-characterized cell line that is easily amenable to experimental manipulations, we expect that it will serve as an excellent system to study piRNA biogenesis and piRNA-related factors.

CLIPSeqTools--a novel bioinformatics CLIP-seq analysis suite.

Immunoprecipitation of RNA binding proteins (RBPs) after in vivo crosslinking, coupled with sequencing of associated RNA footprints (HITS-CLIP, CLIP-seq), is a method of choice for the identification of RNA targets and binding sites for RBPs. Compared with RNA-seq, CLIP-seq analysis is widely diverse and depending on the RBPs that are analyzed, the approaches vary significantly, necessitating the development of flexible and efficient informatics tools. In this study, we present CLIPSeqTools, a novel, highly flexible computational suite that can perform analysis from raw sequencing data with minimal user input. It contains a wide array of tools to provide an in-depth view of CLIP-seq data sets. It supports extensive customization and promotes improvization, a critical virtue, since CLIP-seq analysis is rarely well defined a priori. To highlight CLIPSeqTools capabilities, we used the suite to analyze Ago-miRNA HITS-CLIP data sets that we prepared from human brains.

FUS regulates genes coding for RNA-binding proteins in neurons by binding to their highly conserved introns.

Dominant mutations and mislocalization or aggregation of Fused in Sarcoma (FUS), an RNA-binding protein (RBP), cause neuronal degeneration in Amyotrophic Lateral Sclerosis (ALS) and Frontotemporal Lobar Degeneration (FTLD), two incurable neurological diseases. However, the function of FUS in neurons is not well understood. To uncover the impact of FUS in the neuronal transcriptome, we used high-throughput sequencing of immunoprecipitated and cross-linked RNA (HITS-CLIP) of FUS in human brains and mouse neurons differentiated from embryonic stem cells, coupled with RNA-seq and FUS knockdowns. We report conserved neuronal RNA targets and networks that are regulated by FUS. We find that FUS regulates splicing of genes coding for RBPs by binding to their highly conserved introns. Our findings have important implications for understanding the impact of FUS in neurodegenerative diseases and suggest that perturbations of FUS can impact the neuronal transcriptome via perturbations of RBP transcripts.

Mitochondrial protein BmPAPI modulates the length of mature piRNAs.

PIWI proteins and their associated PIWI-interacting RNAs (piRNAs) protect genome integrity by silencing transposons in animal germlines. The molecular mechanisms and components responsible for piRNA biogenesis remain elusive. PIWI proteins contain conserved symmetrical dimethylarginines (sDMAs) that are specifically targeted by TUDOR domain-containing proteins. Here we report that the sDMAs of PIWI proteins play crucial roles in PIWI localization and piRNA biogenesis in Bombyx mori-derived BmN4 cells, which harbor fully functional piRNA biogenesis machinery. Moreover, RNAi screenings for Bombyx genes encoding TUDOR domain-containing proteins identified BmPAPI, a Bombyx homolog of Drosophila PAPI, as a factor modulating the length of mature piRNAs. BmPAPI specifically recognized sDMAs and interacted with PIWI proteins at the surface of the mitochondrial outer membrane. BmPAPI depletion resulted in 3'-terminal extensions of mature piRNAs without affecting the piRNA quantity. These results reveal the BmPAPI-involved piRNA precursor processing mechanism on mitochondrial outer membrane scaffolds.

DIANA-LncBase: experimentally verified and computationally predicted microRNA targets on long non-coding RNAs.

Recently, the attention of the research community has been focused on long non-coding RNAs (lncRNAs) and their physiological/pathological implications. As the number of experiments increase in a rapid rate and transcriptional units are better annotated, databases indexing lncRNA properties and function gradually become essential tools to this process. Aim of DIANA-LncBase (www.microrna.gr/LncBase) is to reinforce researchers' attempts and unravel microRNA (miRNA)-lncRNA putative functional interactions. This study provides, for the first time, a comprehensive annotation of miRNA targets on lncRNAs. DIANA-LncBase hosts transcriptome-wide experimentally verified and computationally predicted miRNA recognition elements (MREs) on human and mouse lncRNAs. The analysis performed includes an integration of most of the available lncRNA resources, relevant high-throughput HITS-CLIP and PAR-CLIP experimental data as well as state-of-the-art in silico target predictions. The experimentally supported entries available in DIANA-LncBase correspond to >5000 interactions, while the computationally predicted interactions exceed 10 million. DIANA-LncBase hosts detailed information for each miRNA-lncRNA pair, such as external links, graphic plots of transcripts' genomic location, representation of the binding sites, lncRNA tissue expression as well as MREs conservation and prediction scores.

DIANA miRPath v.2.0: investigating the combinatorial effect of microRNAs in pathways.

MicroRNAs (miRNAs) are key regulators of diverse biological processes and their functional analysis has been deemed central in many research pipelines. The new version of DIANA-miRPath web server was redesigned from the ground-up. The user of DNA Intelligent Analysis (DIANA) DIANA-miRPath v2.0 can now utilize miRNA targets predicted with high accuracy based on DIANA-microT-CDS and/or experimentally verified targets from TarBase v6; combine results with merging and meta-analysis algorithms; perform hierarchical clustering of miRNAs and pathways based on their interaction levels; as well as elaborate sophisticated visualizations, such as dendrograms or miRNA versus pathway heat maps, from an intuitive and easy to use web interface. New modules enable DIANA-miRPath server to provide information regarding pathogenic single nucleotide polymorphisms (SNPs) in miRNA target sites (SNPs module) or to annotate all the predicted and experimentally validated miRNA targets in a selected molecular pathway (Reverse Search module). DIANA-miRPath v2.0 is an efficient and yet easy to use tool that can be incorporated successfully into miRNA-related analysis pipelines. It provides for the first time a series of highly specific tools for miRNA-targeted pathway analysis via a web interface and can be accessed at http://www.microrna.gr/miRPathv2.

TarBase 6.0: capturing the exponential growth of miRNA targets with experimental support.

As the relevant literature and the number of experiments increase at a super linear rate, databases that curate and collect experimentally verified microRNA (miRNA) targets have gradually emerged. These databases attempt to provide efficient access to this wealth of experimental data, which is scattered in thousands of manuscripts. Aim of TarBase 6.0 (http://www.microrna.gr/tarbase) is to face this challenge by providing a significant increase of available miRNA targets derived from all contemporary experimental techniques (gene specific and high-throughput), while incorporating a powerful set of tools in a user-friendly interface. TarBase 6.0 hosts detailed information for each miRNA-gene interaction, ranging from miRNA- and gene-related facts to information specific to their interaction, the experimental validation methodologies and their outcomes. All database entries are enriched with function-related data, as well as general information derived from external databases such as UniProt, Ensembl and RefSeq. DIANA microT miRNA target prediction scores and the relevant prediction details are available for each interaction. TarBase 6.0 hosts the largest collection of manually curated experimentally validated miRNA-gene interactions (more than 65,000 targets), presenting a 16.5-175-fold increase over other available manually curated databases.

Full-length direct RNA sequencing uncovers stress-granule dependent RNA decay upon cellular stress.

Cells react to stress by triggering response pathways, leading to extensive alterations in the transcriptome to restore cellular homeostasis. The role of RNA metabolism in shaping the cellular response to stress is vital, yet the global changes in RNA stability under these conditions remain unclear. In this work, we employ direct RNA sequencing with nanopores, enhanced by 5' end adaptor ligation, to comprehensively interrogate the human transcriptome at single-molecule and nucleotide resolution. By developing a statistical framework to identify robust RNA length variations in nanopore data, we find that cellular stress induces prevalent 5' end RNA decay that is coupled to translation and ribosome occupancy. Unlike typical RNA decay models in normal conditions, we show that stress-induced RNA decay is dependent on XRN1 but does not depend on removal of the poly(A) tail. We observed that RNAs undergoing decay are predominantly enriched in the stress granule transcriptome. Inhibition of stress granule formation via genetic ablation of G3BP1 and G3BP2 fully rescues RNA length and suppresses stress-induced decay. Our findings reveal RNA decay as a key determinant of RNA metabolism upon cellular stress and dependent on stress-granule formation.

Gene body DNA hydroxymethylation restricts the magnitude of transcriptional changes during aging.

DNA hydroxymethylation (5hmC), the most abundant oxidative derivative of DNA methylation, is typically enriched at enhancers and gene bodies of transcriptionally active and tissue-specific genes. Although aberrant genomic 5hmC has been implicated in age-related diseases, its functional role in aging remains unknown. Here, using mouse liver and cerebellum as model organs, we show that 5hmC accumulates in gene bodies associated with tissue-specific function and restricts the magnitude of gene expression changes with age. Mechanistically, 5hmC decreases the binding of splicing associated factors and correlates with age-related alternative splicing events. We found that various age-related contexts, such as prolonged quiescence and senescence, drive the accumulation of 5hmC with age. We provide evidence that this age-related transcriptionally restrictive function is conserved in mouse and human tissues. Our findings reveal that 5hmC regulates tissue-specific function and may play a role in longevity.

Mis-spliced transcripts generate de novo proteins in TDP-43-related ALS/FTD.

Functional loss of TDP-43, an RNA binding protein genetically and pathologically linked to amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD), leads to the inclusion of cryptic exons in hundreds of transcripts during disease. Cryptic exons can promote the degradation of affected transcripts, deleteriously altering cellular function through loss-of-function mechanisms. Here, we show that mRNA transcripts harboring cryptic exons generated de novo proteins in TDP-43-depleted human iPSC-derived neurons in vitro, and de novo peptides were found in cerebrospinal fluid (CSF) samples from patients with ALS or FTD. Using coordinated transcriptomic and proteomic studies of TDP-43-depleted human iPSC-derived neurons, we identified 65 peptides that mapped to 12 cryptic exons. Cryptic exons identified in TDP-43-depleted human iPSC-derived neurons were predictive of cryptic exons expressed in postmortem brain tissue from patients with TDP-43 proteinopathy. These cryptic exons produced transcript variants that generated de novo proteins. We found that the inclusion of cryptic peptide sequences in proteins altered their interactions with other proteins, thereby likely altering their function. Last, we showed that 18 de novo peptides across 13 genes were present in CSF samples from patients with ALS/FTD spectrum disorders. The demonstration of cryptic exon translation suggests new mechanisms for ALS/FTD pathophysiology downstream of TDP-43 dysfunction and may provide a potential strategy to assay TDP-43 function in patient CSF.

Functional microRNA targets in protein coding sequences.

MOTIVATION: Experimental evidence has accumulated showing that microRNA (miRNA) binding sites within protein coding sequences (CDSs) are functional in controlling gene expression. RESULTS: Here we report a computational analysis of such miRNA target sites, based on features extracted from existing mammalian high-throughput immunoprecipitation and sequencing data. The analysis is performed independently for the CDS and the 3(')-untranslated regions (3(')-UTRs) and reveals different sets of features and models for the two regions. The two models are combined into a novel computational model for miRNA target genes, DIANA-microT-CDS, which achieves higher sensitivity compared with other popular programs and the model that uses only the 3(')-UTR target sites. Further analysis indicates that genes with shorter 3(')-UTRs are preferentially targeted in the CDS, suggesting that evolutionary selection might favor additional sites on the CDS in cases where there is restricted space on the 3(')-UTR.

DIANA-microT Web server upgrade supports Fly and Worm miRNA target prediction and bibliographic miRNA to disease association.

microRNAs (miRNAs) are small endogenous RNA molecules that are implicated in many biological processes through post-transcriptional regulation of gene expression. The DIANA-microT Web server provides a user-friendly interface for comprehensive computational analysis of miRNA targets in human and mouse. The server has now been extended to support predictions for two widely studied species: Drosophila melanogaster and Caenorhabditis elegans. In the updated version, the Web server enables the association of miRNAs to diseases through bibliographic analysis and provides insights for the potential involvement of miRNAs in biological processes. The nomenclature used to describe mature miRNAs along different miRBase versions has been extensively analyzed, and the naming history of each miRNA has been extracted. This enables the identification of miRNA publications regardless of possible nomenclature changes. User interaction has been further refined allowing users to save results that they wish to analyze further. A connection to the UCSC genome browser is now provided, enabling users to easily preview predicted binding sites in comparison to a wide array of genomic tracks, such as single nucleotide polymorphisms. The Web server is publicly accessible in www.microrna.gr/microT-v4.

Ribonucleic Acid-Mediated Control of Protein Translation Under Stress.

Significance: The need of cells to constantly respond to endogenous and exogenous stress has necessitated the evolution of pathways to counter the deleterious effects of stress and to restore cellular homeostasis. The inability to activate a timely and adequate response can lead to disease and is a hallmark of aging. Besides protein-coding genes, cells contain a plethora of noncoding regulatory elements that allow cells to respond rapidly and efficiently to external stimuli by activating highly specific and tightly controlled mechanisms. Many of these programs converge on the regulation of translation, one of the most energy-consuming processes in cells. Recent Advances: The noncoding dimension of translational regulation includes short and long noncoding ribonucleic acids (ncRNAs), as well as messenger RNA features, such as the sequence and modification status of the 5' and 3' untranslated regions (UTRs), that do not change the amino acid sequence of the produced protein. Critical Issues: In this review, we discuss the regulatory role of the nonprotein-coding components of translation under stress, particularly oxidative stress. We conclude that the regulation of translation through ncRNAs, UTRs, and nucleotide modifications is emerging as a critical component of the stress response. Future Directions: Further areas of study using long-read sequencing technologies will be discussed. Antioxid. Redox Signal. 39, 374-389.

Deep learning and direct sequencing of labeled RNA captures transcriptome dynamics.

Quantification of the dynamics of RNA metabolism is essential for understanding gene regulation in health and disease. Existing methods rely on metabolic labeling of nascent RNAs and physical separation or inference of labeling through PCR-generated mutations, followed by short-read sequencing. However, these methods are limited in their ability to identify transient decay intermediates or co-analyze RNA decay with cis-regulatory elements of RNA stability such as poly(A) tail length and modification status, at single molecule resolution. Here we use 5-ethynyl uridine (5EU) to label nascent RNA followed by direct RNA sequencing with nanopores. We developed RNAkinet, a deep convolutional and recurrent neural network that processes the electrical signal produced by nanopore sequencing to identify 5EU-labeled nascent RNA molecules. RNAkinet demonstrates generalizability to distinct cell types and organisms and reproducibly quantifies RNA kinetic parameters allowing the combined interrogation of RNA metabolism and cis-acting RNA regulatory elements.