Alternative polyadenylation (APA) is a general mechanism of transcript diversification in mammals, which has been recently linked to proliferative states and cancer. Different 3' untranslated region ...(3' UTR) isoforms interact with different RNA-binding proteins (RBPs), which modify the stability, translation, and subcellular localization of the corresponding transcripts. Although the heterogeneity of pre-mRNA 3' end processing has been established with high-throughput approaches, the mechanisms that underlie systematic changes in 3' UTR lengths remain to be characterized. Through a uniform analysis of a large number of 3' end sequencing data sets, we have uncovered 18 signals, six of which are novel, whose positioning with respect to pre-mRNA cleavage sites indicates a role in pre-mRNA 3' end processing in both mouse and human. With 3' end sequencing we have demonstrated that the heterogeneous ribonucleoprotein C (HNRNPC), which binds the poly(U) motif whose frequency also peaks in the vicinity of polyadenylation (poly(A)) sites, has a genome-wide effect on poly(A) site usage. HNRNPC-regulated 3' UTRs are enriched in ELAV-like RBP 1 (ELAVL1) binding sites and include those of the CD47 gene, which participate in the recently discovered mechanism of 3' UTR-dependent protein localization (UDPL). Our study thus establishes an up-to-date, high-confidence catalog of 3' end processing sites and poly(A) signals, and it uncovers an important role of HNRNPC in regulating 3' end processing. It further suggests that U-rich elements mediate interactions with multiple RBPs that regulate different stages in a transcript's life cycle.
Understanding the regulation of gene expression, including transcription start site usage, alternative splicing, and polyadenylation, requires accurate quantification of expression levels down to the ...level of individual transcript isoforms. To comparatively evaluate the accuracy of the many methods that have been proposed for estimating transcript isoform abundance from RNA sequencing data, we have used both synthetic data as well as an independent experimental method for quantifying the abundance of transcript ends at the genome-wide level.
We found that many tools have good accuracy and yield better estimates of gene-level expression compared to commonly used count-based approaches, but they vary widely in memory and runtime requirements. Nucleotide composition and intron/exon structure have comparatively little influence on the accuracy of expression estimates, which correlates most strongly with transcript/gene expression levels. To facilitate the reproduction and further extension of our study, we provide datasets, source code, and an online analysis tool on a companion website, where developers can upload expression estimates obtained with their own tool to compare them to those inferred by the methods assessed here.
As many methods for quantifying isoform abundance with comparable accuracy are available, a user's choice will likely be determined by factors such as the memory and runtime requirements, as well as the availability of methods for downstream analyses. Sequencing-based methods to quantify the abundance of specific transcript regions could complement validation schemes based on synthetic data and quantitative PCR in future or ongoing assessments of RNA-seq analysis methods.
Through alternative polyadenylation, human mRNAs acquire longer or shorter 3′ untranslated regions, the latter typically associated with higher transcript stability and increased protein production. ...To understand the dynamics of polyadenylation site usage, we performed transcriptome-wide mapping of both binding sites of 3′ end processing factors CPSF-160, CPSF-100, CPSF-73, CPSF-30, Fip1, CstF-64, CstF-64τ, CF Im25, CF Im59, and CF Im68 and 3′ end processing sites in HEK293 cells. We found that although binding sites of these factors generally cluster around the poly(A) sites most frequently used in cleavage, CstF-64/CstF-64τ and CFIm proteins have much higher positional specificity compared to CPSF components. Knockdown of CF Im68 induced a systematic use of proximal polyadenylation sites, indicating that changes in relative abundance of a single 3′ end processing factor can modulate the length of 3′ untranslated regions across the transcriptome and suggesting a mechanism behind the previously observed increase in tumor cell invasiveness upon CF Im68 knockdown.
Display omitted
► Thousands of poly(A) sites are identified by the A-seq method in a human cell line ► Binding sites of pre-mRNA 3′ end processing factors are mapped by PAR-CLIP ► CstF-64 and CF Im68 exhibit the highest positional binding specificity ► CF Im 68 siRNA treatment causes a global shift toward proximal poly(A) sites
Alternative cleavage and polyadenylation generate mRNAs with 3′ untranslated regions of different lengths. Keller, Zavolan, and colleagues mapped the binding sites of ten 3′ end processing proteins by PAR-CLIP and identified the 3′ end processing sites by A-seq. CstF-64 and CF Im68 proteins showed the highest positional specificity. Knockdown of CF Im68 induced a systematic shift toward proximal polyadenylation sites, indicating that changes in the relative abundance of a single 3′ end processing factor can modulate the length of 3′ untranslated regions transcriptome-wide.
3' Untranslated regions (3' UTRs) length is regulated in relation to cellular state. To uncover key regulators of poly(A) site use in specific conditions, we have developed PAQR, a method for ...quantifying poly(A) site use from RNA sequencing data and KAPAC, an approach that infers activities of oligomeric sequence motifs on poly(A) site choice. Application of PAQR and KAPAC to RNA sequencing data from normal and tumor tissue samples uncovers motifs that can explain changes in cleavage and polyadenylation in specific cancers. In particular, our analysis points to polypyrimidine tract binding protein 1 as a regulator of poly(A) site choice in glioblastoma.
Alternative polyadenylation is a cellular mechanism that generates mRNA isoforms differing in their 3' untranslated regions (3' UTRs). Changes in polyadenylation site usage have been described upon ...induction of proliferation in resting cells, but the underlying mechanism and functional significance of this phenomenon remain largely unknown. To understand the functional consequences of shortened 3' UTR isoforms in a physiological setting, we used 3' end sequencing and quantitative mass spectrometry to determine polyadenylation site usage, mRNA and protein levels in murine and human naive and activated T cells. Although 3' UTR shortening in proliferating cells is conserved between human and mouse, orthologous genes do not exhibit similar expression of alternative 3' UTR isoforms. We generally find that 3' UTR shortening is not accompanied by a corresponding change in mRNA and protein levels. This suggests that although 3' UTR shortening may lead to changes in the RNA-binding protein interactome, it has limited effects on protein output.
Small nucleolar RNAs (snoRNAs) are a class of non-coding RNAs that guide the post-transcriptional processing of other non-coding RNAs (mostly ribosomal RNAs), but have also been implicated in ...processes ranging from microRNA-dependent gene silencing to alternative splicing. In order to construct an up-to-date catalog of human snoRNAs we have combined data from various databases, de novo prediction and extensive literature review. In total, we list more than 750 curated genomic loci that give rise to snoRNA and snoRNA-like genes. Utilizing small RNA-seq data from the ENCODE project, our study characterizes the plasticity of snoRNA expression identifying both constitutively as well as cell type specific expressed snoRNAs. Especially, the comparison of malignant to non-malignant tissues and cell types shows a dramatic perturbation of the snoRNA expression profile. Finally, we developed a high-throughput variant of the reverse-transcriptase-based method for identifying 2'-O-methyl modifications in RNAs termed RimSeq. Using the data from this and other high-throughput protocols together with previously reported modification sites and state-of-the-art target prediction methods we re-estimate the snoRNA target RNA interaction network. Our current results assign a reliable modification site to 83% of the canonical snoRNAs, leaving only 76 snoRNA sequences as orphan.
The prediction of a consensus structure for a set of related RNAs is an important first step for subsequent analyses. RNAalifold, which computes the minimum energy structure that is simultaneously ...formed by a set of aligned sequences, is one of the oldest and most widely used tools for this task. In recent years, several alternative approaches have been advocated, pointing to several shortcomings of the original RNAalifold approach.
We show that the accuracy of RNAalifold predictions can be improved substantially by introducing a different, more rational handling of alignment gaps, and by replacing the rather simplistic model of covariance scoring with more sophisticated RIBOSUM-like scoring matrices. These improvements are achieved without compromising the computational efficiency of the algorithm. We show here that the new version of RNAalifold not only outperforms the old one, but also several other tools recently developed, on different datasets.
The new version of RNAalifold not only can replace the old one for almost any application but it is also competitive with other approaches including those based on SCFGs, maximum expected accuracy, or hierarchical nearest neighbor classifiers.
Celotno besedilo
Dostopno za:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Cleavage and polyadenylation specificity factor (CPSF) is the central component of the 3' processing machinery for polyadenylated mRNAs in metazoans: CPSF recognizes the polyadenylation signal ...AAUAAA, providing sequence specificity in both pre-mRNA cleavage and polyadenylation, and catalyzes pre-mRNA cleavage. Here we show that of the seven polypeptides that have been proposed to constitute CPSF, only four (CPSF160, CPSF30, hFip1, and WDR33) are necessary and sufficient to reconstitute a CPSF subcomplex active in AAUAAA-dependent polyadenylation, whereas CPSF100, CPSF73, and symplekin are dispensable. WDR33 is required for binding of reconstituted CPSF to AAUAAA-containing RNA and can be specifically UV cross-linked to such RNAs, as can CPSF30. Transcriptome-wide identification of WDR33 targets by photoactivatable ribonucleoside-enhanced cross-linking and immunoprecipitation (PAR-CLIP) showed that WDR33 binds in and very close to the AAUAAA signal in vivo with high specificity. Thus, our data indicate that the large CPSF subunit participating in recognition of the polyadenylation signal is WDR33 and not CPSF160, as suggested by previous studies.
Vienna RNA Websuite Gruber, Andreas R; Lorenz, Ronny; Bernhart, Stephan H ...
Nucleic acids research,
07/2008, Letnik:
36, Številka:
suppl-2
Journal Article
Recenzirano
Odprti dostop
The Vienna RNA Websuite is a comprehensive collection of tools for folding, design and analysis of RNA sequences. It provides a web interface to the most commonly used programs of the Vienna RNA ...package. Among them, we find folding of single and aligned sequences, prediction of RNA-RNA interactions, and design of sequences with a given structure. Additionally, we provide analysis of folding landscapes using the barriers program and structural RNA alignments using LocARNA. The web server together with software packages for download is freely accessible at http://rna.tbi.univie.ac.at/.
Recruitment of the UPF1 nonsense-mediated mRNA decay (NMD) factor to target mRNAs was initially proposed to occur through interaction with release factors at terminating ribosomes. However, recently ...emerging evidence points toward translation-independent interaction with the 3' untranslated region (UTR) of mRNAs. We mapped transcriptome-wide UPF1-binding sites by individual-nucleotide-resolution UV cross-linking and immunoprecipitation in human cells and found that UPF1 preferentially associated with 3' UTRs in translationally active cells but underwent significant redistribution toward coding regions (CDS) upon translation inhibition, thus indicating that UPF1 binds RNA before translation and gets displaced from the CDS by translating ribosomes. Corroborated by RNA immunoprecipitation and by UPF1 cross-linking to long noncoding RNAs, our evidence for translation-independent UPF1-RNA interaction suggests that the triggering of NMD occurs after UPF1 binding to mRNA, presumably through activation of RNA-bound UPF1 by aberrant translation termination.
Celotno besedilo
Dostopno za:
DOBA, IJS, IZUM, KILJ, NUK, PILJ, PNG, SAZU, UILJ, UKNU, UL, UM, UPUK