microRNAs are important regulators of gene expression that guide translational repression and degradation of target mRNAs. Only relatively few miRNA targets have been characterized, and computational prediction is hampered by the relatively small number of nucleotides that seem to be involved in target recognition. Argonaute (Ago) crosslinking and immunoprecipitation (CLIP) in combination with next-generation sequencing proved to be a successful method for identifying targets of endogenous cellular miRNAs on a transcriptome-wide scale. Here we review various approaches to Ago CLIP, describe in detail the PAR-CLIP method and provide an outline of the necessary computational analysis for identification of in vivo miRNA binding sites.

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.

Most of what is presently known about how miRNAs regulate gene expression comes from studies that characterized the regulatory effect of miRNA binding sites located in the 3' untranslated regions (UTR) of mRNAs. In recent years, there has been increasing evidence that miRNAs also bind in the coding region (CDS), but the implication of these interactions remains obscure because they have a smaller impact on mRNA stability compared with miRNA-target interactions that involve 3' UTRs. Here we show that miRNA-complementary sites that are located in both CDS and 3'-UTRs are under selection pressure and share the same sequence and structure properties. Analyzing recently published data of ribosome-protected fragment profiles upon miRNA transfection from the perspective of the location of miRNA-complementary sites, we find that sites located in the CDS are most potent in inhibiting translation, while sites located in the 3' UTR are more efficient at triggering mRNA degradation. Our study suggests that miRNAs may combine targeting of CDS and 3' UTR to flexibly tune the time scale and magnitude of their post-transcriptional regulatory effects.

The stability, localization and translation rate of mRNAs are regulated by a multitude of RNA-binding proteins (RBPs) that find their targets directly or with the help of guide RNAs. Among the experimental methods for mapping RBP binding sites, cross-linking and immunoprecipitation (CLIP) coupled with deep sequencing provides transcriptome-wide coverage as well as high resolution. However, partly due to their vast volume, the data that were so far generated in CLIP experiments have not been put in a form that enables fast and interactive exploration of binding sites. To address this need, we have developed the CLIPZ database and analysis environment. Binding site data for RBPs such as Argonaute 1-4, Insulin-like growth factor II mRNA-binding protein 1-3, TNRC6 proteins A-C, Pumilio 2, Quaking and Polypyrimidine tract binding protein can be visualized at the level of the genome and of individual transcripts. Individual users can upload their own sequence data sets while being able to limit the access to these data to specific users, and analyses of the public and private data sets can be performed interactively. CLIPZ, available at http://www.clipz.unibas.ch, aims to provide an open access repository of information for post-transcriptional regulatory elements.

Small nucleolar RNAs (snoRNAs) are localized within the nucleolus, a sub-nuclear compartment, in which they guide ribosomal or spliceosomal RNA modifications, respectively. Up until now, snoRNAs have only been identified in eukaryal and archaeal genomes, but are notably absent in bacteria. By screening B lymphocytes for expression of non-coding RNAs (ncRNAs) induced by the Epstein-Barr virus (EBV), we here report, for the first time, the identification of a snoRNA gene within a viral genome, designated as v-snoRNA1. This genetic element displays all hallmark sequence motifs of a canonical C/D box snoRNA, namely C/C'- as well as D/D'-boxes. The nucleolar localization of v-snoRNA1 was verified by in situ hybridisation of EBV-infected cells. We also confirmed binding of the three canonical snoRNA proteins, fibrillarin, Nop56 and Nop58, to v-snoRNA1. The C-box motif of v-snoRNA1 was shown to be crucial for the stability of the viral snoRNA; its selective deletion in the viral genome led to a complete down-regulation of v-snoRNA1 expression levels within EBV-infected B cells. We further provide evidence that v-snoRNA1 might serve as a miRNA-like precursor, which is processed into 24 nt sized RNA species, designated as v-snoRNA1(24pp). A potential target site of v-snoRNA1(24pp) was identified within the 3'-UTR of BALF5 mRNA which encodes the viral DNA polymerase. V-snoRNA1 was found to be expressed in all investigated EBV-positive cell lines, including lymphoblastoid cell lines (LCL). Interestingly, induction of the lytic cycle markedly up-regulated expression levels of v-snoRNA1 up to 30-fold. By a computational approach, we identified a v-snoRNA1 homolog in the rhesus lymphocryptovirus genome. This evolutionary conservation suggests an important role of v-snoRNA1 during gamma-herpesvirus infection.

MicroRNAs (miRNAs) are short RNAs that act as guides for the degradation and translational repression of protein-coding mRNAs. A large body of work showed that miRNAs are involved in the regulation of a broad range of biological functions, from development to cardiac and immune system function, to metabolism, to cancer. For most of the over 500 miRNAs that are encoded in the human genome the functions still remain to be uncovered. Identifying miRNAs whose expression changes between cell types or between normal and pathological conditions is an important step towards characterizing their function as is the prediction of mRNAs that could be targeted by these miRNAs. To provide the community the possibility of exploring interactively miRNA expression patterns and the candidate targets of miRNAs in an integrated environment, we developed the MirZ web server, which is accessible at www.mirz.unibas.ch. The server provides experimental and computational biologists with statistical analysis and data mining tools operating on up-to-date databases of sequencing-based miRNA expression profiles and of predicted miRNA target sites in species ranging from Caenorhabditis elegans to Homo sapiens.

In animals, microRNAs frequently form families with related sequences. The functional relevance of miRNA families and the relative contribution of family members to target repression have remained, however, largely unexplored. Here, we used the Caenorhabditis elegans miR-58 miRNA family, composed primarily of the four highly abundant members miR-58.1, miR-80, miR-81, and miR-82, as a model to investigate the redundancy of miRNA family members and their impact on target expression in an in vivo setting. We found that miR-58 family members repress largely overlapping sets of targets in a predominantly additive fashion. Progressive deletions of miR-58 family members lead to cumulative up-regulation of target protein and RNA levels. Phenotypic defects could only be observed in the family quadruple mutant, which also showed the strongest change in target protein levels. Interestingly, although the seed sequences of miR-80 and miR-58.1 differ in a single nucleotide, predicted canonical miR-80 targets were efficiently up-regulated in the mir-58.1 single mutant, indicating functional redundancy of distinct members of this miRNA family. At the aggregate level, target binding leads mainly to mRNA degradation, although we also observed some degree of translational inhibition, particularly in the single miR-58 family mutants. These results provide a framework for understanding how miRNA family members interact to regulate target mRNAs.

Preserving skeletal muscle function is essential to maintain life quality at high age. Calorie restriction (CR) potently extends health and lifespan, but is largely unachievable in humans, making "CR mimetics" of great interest. CR targets nutrient-sensing pathways centering on mTORC1. The mTORC1 inhibitor, rapamycin, is considered a potential CR mimetic and is proven to counteract age-related muscle loss. Therefore, we tested whether rapamycin acts via similar mechanisms as CR to slow muscle aging. Here we show that long-term CR and rapamycin unexpectedly display distinct gene expression profiles in geriatric mouse skeletal muscle, despite both benefiting aging muscles. Furthermore, CR improves muscle integrity in mice with nutrient-insensitive, sustained muscle mTORC1 activity and rapamycin provides additive benefits to CR in naturally aging mouse muscles. We conclude that rapamycin and CR exert distinct, compounding effects in aging skeletal muscle, thus opening the possibility of parallel interventions to counteract muscle aging.

The speed of translation elongation is primarily determined by the abundance of tRNAs. Thus, the codon usage influences the rate with which individual mRNAs are translated. As the nature of tRNA pools and modifications can vary across biological conditions, codon elongation rates may also vary, leading to fluctuations in the protein production from individual mRNAs. Although it has been observed that functionally related mRNAs exhibit similar codon usage, presumably to provide an effective way to coordinate expression of multiple proteins, experimental evidence for codon-mediated translation efficiency modulation of functionally related mRNAs in specific conditions is scarce and the associated mechanisms are still debated.

Understanding the molecular signatures of colorectal cancer progression under chemotherapeutic treatment will be crucial for the success of future therapy improvements. Here, we used a xenograft-based mouse model to investigate, how whole transcriptome signatures change during metastatic colorectal cancer progression and how such signatures are affected by LDM chemotherapy using RNA sequencing. We characterized mRNAs as well as non-coding RNAs such as microRNAs, long non-coding RNAs and circular RNAs in colorectal-cancer bearing mice with or without LDM chemotherapy. Furthermore, we found that circZNF609 functions as oncogene, since over-expression studies lead to an increased tumor growth while specific knock down results in smaller tumors. Our data represent novel insights into the relevance of non-coding and circRNAs in colorectal cancer and provide a comprehensive resource of gene expression changes in primary tumors and metastases. In addition, we present candidate genes that could be important modulators for successful LDM chemotherapy.

Mutations in the RNA-binding protein Fused in Sarcoma (FUS) cause early-onset amyotrophic lateral sclerosis (ALS). However, a detailed understanding of central RNA targets of FUS and their implications for disease remain elusive. Here, we use a unique blend of crosslinking and immunoprecipitation (CLIP) and NMR spectroscopy to identify and characterise physiological and pathological RNA targets of FUS. We find that U1 snRNA is the primary RNA target of FUS via its interaction with stem-loop 3 and provide atomic details of this RNA-mediated mode of interaction with the U1 snRNP. Furthermore, we show that ALS-associated FUS aberrantly contacts U1 snRNA at the Sm site with its zinc finger and traps snRNP biogenesis intermediates in human and murine motor neurons. Altogether, we present molecular insights into a FUS toxic gain-of-function involving direct and aberrant RNA-binding and strengthen the link between two motor neuron diseases, ALS and spinal muscular atrophy (SMA).

The behavior of cells in vivo is complex and highly dynamic, as it results from an interplay between intercellular matrix proteins with surface receptors and other microenvironmental cues. Although the effects of the cellular niche have been investigated for a number of cell types using different molecular approaches, comprehensive assessments of how the global transcriptome responds to 3D scaffolds composed of various extracellular matrix (ECM) constituents at different concentrations are still lacking.

The tremendous rate with which data is generated and analysis methods emerge makes it increasingly difficult to keep track of their domain of applicability, assumptions, limitations, and consequently, of the efficacy and precision with which they solve specific tasks. Therefore, there is an increasing need for benchmarks, and for the provision of infrastructure for continuous method evaluation. APAeval is an international community effort, organized by the RNA Society in 2021, to benchmark tools for the identification and quantification of the usage of alternative polyadenylation (APA) sites from short-read, bulk RNA-sequencing (RNA-seq) data. Here, we reviewed 17 tools and benchmarked eight on their ability to perform APA identification and quantification, using a comprehensive set of RNA-seq experiments comprising real, synthetic, and matched 3'-end sequencing data. To support continuous benchmarking, we have incorporated the results into the OpenEBench online platform, which allows for continuous extension of the set of methods, metrics, and challenges. We envisage that our analyses will assist researchers in selecting the appropriate tools for their studies, while the containers and reproducible workflows could easily be deployed and extended to evaluate new methods or data sets.

Sarcopenia, the age-related loss of skeletal muscle mass and function, can dramatically impinge on quality of life and mortality. While mitochondrial dysfunction and imbalanced proteostasis are recognized as hallmarks of sarcopenia, the regulatory and functional link between these processes is underappreciated and unresolved. We therefore investigated how mitochondrial proteostasis, a crucial process that coordinates the expression of nuclear- and mitochondrial-encoded mitochondrial proteins with supercomplex formation and respiratory activity, is affected in skeletal muscle aging. Intriguingly, a robust mitochondrial translation impairment was observed in sarcopenic muscle, which is regulated by the peroxisome proliferator-activated receptor γ coactivator 1 α (PGC-1α) with the estrogen-related receptor α (ERRα). Exercise, a potent inducer of PGC-1α activity, rectifies age-related reduction in mitochondrial translation, in conjunction with quality control pathways. These results highlight the importance of mitochondrial proteostasis in muscle aging, and elucidate regulatory interactions that underlie the powerful benefits of physical activity in this context.

Alternative polyadenylation is a main driver of transcriptome diversity in mammals, generating transcript isoforms with different 3' ends via cleavage and polyadenylation at distinct polyadenylation (poly(A)) sites. The regulation of cell type-specific poly(A) site choice is not completely resolved, and requires quantitative poly(A) site usage data across cell types. 3' end-based single-cell RNA-seq can now be broadly used to obtain such data, enabling the identification and quantification of poly(A) sites with direct experimental support. We propose SCINPAS, a computational method to identify poly(A) sites from scRNA-seq datasets. SCINPAS modifies the read deduplication step to favor the selection of distal reads and extract those with non-templated poly(A) tails. This approach improves the resolution of poly(A) site recovery relative to standard software. SCINPAS identifies poly(A) sites in genic and non-genic regions, providing complementary information relative to other tools. The workflow is modular, and the key read deduplication step is general, enabling the use of SCINPAS in other typical analyses of single cell gene expression. Taken together, we show that SCINPAS is able to identify experimentally-supported, known and novel poly(A) sites from 3' end-based single-cell RNA sequencing data.

The RNA-binding protein Roquin is required to prevent autoimmunity. Roquin controls T-helper cell activation and differentiation by limiting the induced expression of costimulatory receptors such as tumor necrosis factor receptor superfamily 4 (Tnfrs4 or Ox40). A constitutive decay element (CDE) with a characteristic triloop hairpin was previously shown to be recognized by Roquin. Here we use SELEX assays to identify a novel U-rich hexaloop motif, representing an alternative decay element (ADE). Crystal structures and NMR data show that the Roquin-1 ROQ domain recognizes hexaloops in the SELEX-derived ADE and in an ADE-like variant present in the Ox40 3'-UTR with identical binding modes. In cells, ADE-like and CDE-like motifs cooperate in the repression of Ox40 by Roquin. Our data reveal an unexpected recognition of hexaloop cis elements for the posttranscriptional regulation of target messenger RNAs by Roquin.

Expression of mature messenger RNAs (mRNAs) requires appropriate transcription initiation and termination, as well as pre-mRNA processing by capping, splicing, cleavage, and polyadenylation. A core 3'-end processing complex carries out the cleavage and polyadenylation reactions, but many proteins have been implicated in the selection of polyadenylation sites among the multiple alternatives that eukaryotic genes typically have. In recent years, high-throughput approaches to map both the 3'-end processing sites as well as the binding sites of proteins that are involved in the selection of cleavage sites and in the processing reactions have been developed. Here, we review these approaches as well as the insights into the mechanisms of polyadenylation that emerged from genome-wide studies of polyadenylation across a range of cell types and states.

The conserved pre-mRNA splicing factor SF1 is implicated in 3' splice site recognition by binding directly to the intron branch site. However, because SF1 is not essential for constitutive splicing, its role in pre-mRNA processing has remained mysterious. Here, we used crosslinking and immunoprecipitation (CLIP) to analyze short RNAs directly bound by human SF1 in vivo. SF1 bound mainly pre-mRNAs, with 77% of target sites in introns. Binding to target RNAs in vitro was dependent on the newly defined SF1 binding motif ACUNAC, strongly resembling human branch sites. Surprisingly, the majority of SF1 binding sites did not map to the expected position near 3' splice sites. Instead, target sites were distributed throughout introns, and a smaller but significant fraction occurred in exons within coding and untranslated regions. These data suggest a more complex role for SF1 in splicing regulation. Indeed, SF1 silencing affected alternative splicing of endogenous transcripts, establishing a previously unexpected role for SF1 and branch site-like sequences in splice site selection.

Accurate reconstruction of the regulatory networks that control gene expression is one of the key current challenges in molecular biology. Although gene expression and chromatin state dynamics are ultimately encoded by constellations of binding sites recognized by regulators such as transcriptions factors (TFs) and microRNAs (miRNAs), our understanding of this regulatory code and its context-dependent read-out remains very limited. Given that there are thousands of potential regulators in mammals, it is not practical to use direct experimentation to identify which of these play a key role for a particular system of interest. We developed a methodology that models gene expression or chromatin modifications in terms of genome-wide predictions of regulatory sites and completely automated it into a web-based tool called ISMARA (Integrated System for Motif Activity Response Analysis). Given only gene expression or chromatin state data across a set of samples as input, ISMARA identifies the key TFs and miRNAs driving expression/chromatin changes and makes detailed predictions regarding their regulatory roles. These include predicted activities of the regulators across the samples, their genome-wide targets, enriched gene categories among the targets, and direct interactions between the regulators. Applying ISMARA to data sets from well-studied systems, we show that it consistently identifies known key regulators ab initio. We also present a number of novel predictions including regulatory interactions in innate immunity, a master regulator of mucociliary differentiation, TFs consistently disregulated in cancer, and TFs that mediate specific chromatin modifications.

Genome-wide identification of mRNAs regulated by RNA-binding proteins is crucial to uncover post-transcriptional gene regulatory systems. The conserved PUF family RNA-binding proteins repress gene expression post-transcriptionally by binding to sequence elements in 3'-UTRs of mRNAs. Despite their well-studied implications for development and neurogenesis in metazoa, the mammalian PUF family members are only poorly characterized and mRNA targets are largely unknown. We have systematically identified the mRNAs associated with the two human PUF proteins, PUM1 and PUM2, by the recovery of endogenously formed ribonucleoprotein complexes and the analysis of associated RNAs with DNA microarrays. A largely overlapping set comprised of hundreds of mRNAs were reproducibly associated with the paralogous PUM proteins, many of them encoding functionally related proteins. A characteristic PUF-binding motif was highly enriched among PUM bound messages and validated with RNA pull-down experiments. Moreover, PUF motifs as well as surrounding sequences exhibit higher conservation in PUM bound messages as opposed to transcripts that were not found to be associated, suggesting that PUM function may be modulated by other factors that bind conserved elements. Strikingly, we found that PUF motifs are enriched around predicted miRNA binding sites and that high-confidence miRNA binding sites are significantly enriched in the 3'-UTRs of experimentally determined PUM1 and PUM2 targets, strongly suggesting an interaction of human PUM proteins with the miRNA regulatory system. Our work suggests extensive connections between the RBP and miRNA post-transcriptional regulatory systems and provides a framework for deciphering the molecular mechanism by which PUF proteins regulate their target mRNAs.