APAlyzer: a bioinformatics package for analysis of alternative polyadenylation isoforms.
Ontology highlight
ABSTRACT: Summary
Most eukaryotic genes produce alternative polyadenylation (APA) isoforms. APA is dynamically regulated under different growth and differentiation conditions. Here, we present a bioinformatics package, named APAlyzer, for examining 3'UTR APA, intronic APA and gene expression changes using RNA-seq data and annotated polyadenylation sites in the PolyA_DB database. Using APAlyzer and data from the GTEx database, we present APA profiles across human tissues.Availability and implementation
APAlyzer is freely available at https://bioconductor.org/packages/release/bioc/html/APAlyzer.html as an R/Bioconductor package.Supplementary information
Supplementary data are available at Bioinformatics online.
SUBMITTER: Wang R
PROVIDER: S-EPMC7320624 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:The use of alternative promoters, splicing, and cleavage and polyadenylation (APA) generates mRNA isoforms that expand the diversity and complexity of the transcriptome. Here, we uncovered thousands of previously undescribed 5' uncapped and polyadenylated transcripts (5' UPTs). We show that these transcripts resist exonucleases due to a highly structured RNA and N6-methyladenosine modification at their 5' termini. 5' UPTs appear downstream of APA sites within their host genes and are induced upon APA activation. Strong enrichment in polysomal RNA fractions indicates 5' UPT translational potential. Indeed, APA promotes downstream translation initiation, non-canonical protein output, and consistent changes to peptide presentation at the cell surface. Lastly, we demonstrate the biological importance of 5' UPTs using Bcl2, a prominent anti-apoptotic gene whose entire coding sequence is a 5' UPT generated from 5' UTR-embedded APA sites. Thus, APA is not only accountable for terminating transcripts, but also for generating downstream uncapped RNAs with translation potential and biological impact.
| S-EPMC9636002 | biostudies-literature
Project description:Acetylcholinesterase (AChE), encoded by the ACHE gene, hydrolyzes the neurotransmitter acetylcholine to terminate synaptic transmission. Alternative splicing close to the 3΄ end generates three distinct isoforms of AChET, AChEH and AChER. We found that hnRNP H binds to two specific G-runs in exon 5a of human ACHE and activates the distal alternative 3΄ splice site (ss) between exons 5a and 5b to generate AChET. Specific effect of hnRNP H was corroborated by siRNA-mediated knockdown and artificial tethering of hnRNP H. Furthermore, hnRNP H competes for binding of CstF64 to the overlapping binding sites in exon 5a, and suppresses the selection of a cryptic polyadenylation site (PAS), which additionally ensures transcription of the distal 3΄ ss required for the generation of AChET. Expression levels of hnRNP H were positively correlated with the proportions of the AChET isoform in three different cell lines. HnRNP H thus critically generates AChET by enhancing the distal 3΄ ss and by suppressing the cryptic PAS. Global analysis of CLIP-seq and RNA-seq also revealed that hnRNP H competitively regulates alternative 3΄ ss and alternative PAS in other genes. We propose that hnRNP H is an essential factor that competitively regulates alternative splicing and alternative polyadenylation.
| S-EPMC5388418 | biostudies-literature
Project description:The human ether-a-go-go-related gene 1 (hERG1) encodes the pore-forming subunit of the rapidly activating delayed rectifier potassium channel. Several hERG1 isoforms with different N- and C-terminal ends have been identified. The hERG1a, hERG1b, and hERG1-3.1 isoforms contain the full-length C terminus, whereas the hERG1(USO) isoforms, hERG1a(USO) and hERG1b(USO), lack most of the C-terminal domain and contain a unique C-terminal end. The mechanisms underlying the generation of hERG1(USO) isoforms are not understood. We show that hERG1 isoforms with different C-terminal ends are generated by alternative splicing and polyadenylation of hERG1 pre-mRNA. We identified an intrinsically weak, noncanonical poly(A) signal, AGUAAA, within intron 9 of hERG1 that modulates the expression of hERG1a and hERG1a(USO). Replacing AGUAAA with the strong, canonical poly(A) signal AAUAAA resulted in the predominant production of hERG1a(USO) and a marked decrease in hERG1 current. In contrast, eliminating the intron 9 poly(A) signal or increasing the strength of 5' splice site led to the predominant production of hERG1a and a significant increase in hERG1 current. We found significant variation in the relative abundance of hERG1 C-terminal isoforms in different human tissues. Taken together, these findings suggest that post-transcriptional regulation of hERG1 pre-mRNA may represent a novel mechanism to modulate the expression and function of hERG1 channels.
| S-EPMC2952224 | biostudies-literature
Project description:Messenger RNA 3'-end formation is an essential posttranscriptional processing step for most eukaryotic genes. Different from plants and animals where AAUAAA and its variants routinely are found as the main poly(A) signal, Chlamydomonas reinhardtii uses UGUAA as the major poly(A) signal. The advance of sequencing technology provides an enormous amount of sequencing data for us to explore the variations of poly(A) signals, alternative polyadenylation (APA), and its relationship with splicing in this algal species. Through genome-wide analysis of poly(A) sites in C. reinhardtii, we identified a large number of poly(A) sites: 21,041 from Sanger expressed sequence tags, 88,184 from 454, and 195,266 from Illumina sequence reads. In comparison with previous collections, more new poly(A) sites are found in coding sequences and intron and intergenic regions by deep-sequencing. Interestingly, G-rich signals are particularly abundant in intron and intergenic regions. The prevalence of different poly(A) signals between coding sequences and a 3'-untranslated region implies potentially different polyadenylation mechanisms. Our data suggest that the APA occurs in about 68% of C. reinhardtii genes. Using Gene Ontolgy analysis, we found most of the APA genes are involved in RNA regulation and metabolic process, protein synthesis, hydrolase, and ligase activities. Moreover, intronic poly(A) sites are more abundant in constitutively spliced introns than retained introns, suggesting an interplay between polyadenylation and splicing. Our results support that APA, as in higher eukaryotes, may play significant roles in increasing transcriptome diversity and gene expression regulation in this algal species. Our datasets also provide useful information for accurate annotation of transcript ends in C. reinhardtii.
| S-EPMC4025486 | biostudies-literature
Project description:MotivationThe length of the 3' untranslated region (3' UTR) of an mRNA is essential for many biological activities such as mRNA stability, sub-cellular localization, protein translation, protein binding and translation efficiency. Moreover, correlation between diseases and the shortening (or lengthening) of 3' UTRs has been reported in the literature. This length is largely determined by the polyadenylation cleavage site in the mRNA. As alternative polyadenylation (APA) sites are common in mammalian genes, several tools have been published recently for detecting APA sites from RNA-Seq data or performing shortening/lengthening analysis. These tools consider either up to only two APA sites in a gene or only APA sites that occur in the last exon of a gene, although a gene may generally have more than two APA sites and an APA site may sometimes occur before the last exon. Furthermore, the tools are unable to integrate the analysis of shortening/lengthening events with APA site detection.ResultsWe propose a new tool, called TAPAS, for detecting novel APA sites from RNA-Seq data. It can deal with more than two APA sites in a gene as well as APA sites that occur before the last exon. The tool is based on an existing method for finding change points in time series data, but some filtration techniques are also adopted to remove change points that are likely false APA sites. It is then extended to identify APA sites that are expressed differently between two biological samples and genes that contain 3' UTRs with shortening/lengthening events. Our extensive experiments on simulated and real RNA-Seq data demonstrate that TAPAS outperforms the existing tools for APA site detection or shortening/lengthening analysis significantly.Availability and implementationhttps://github.com/arefeen/TAPAS.Supplementary informationSupplementary data are available at Bioinformatics online.
| S-EPMC6454472 | biostudies-literature
Project description:BackgroundAlternative polyadenylation (APA) is emerging as an important mechanism in the post-transcriptional regulation of gene expression across eukaryotic species. Recent studies have shown that APA plays key roles in biological processes, such as cell proliferation and differentiation. Single-cell RNA-seq technologies are widely used in gene expression heterogeneity studies; however, systematic studies of APA at the single-cell level are still lacking.ResultsHere, we described a novel computational framework, SAPAS, that utilizes 3'-tag-based scRNA-seq data to identify novel poly(A) sites and quantify APA at the single-cell level. Applying SAPAS to the scRNA-seq data of phenotype characterized GABAergic interneurons, we identified cell type-specific APA events for different GABAergic neuron types. Genes with cell type-specific APA events are enriched for synaptic architecture and communications. In further, we observed a strong enrichment of heritability for several psychiatric disorders and brain traits in altered 3' UTRs and coding sequences of cell type-specific APA events. Finally, by exploring the modalities of APA, we discovered that the bimodal APA pattern of Pak3 could classify chandelier cells into different subpopulations that are from different laminar positions.ConclusionsWe established a method to characterize APA at the single-cell level. When applied to a scRNA-seq dataset of GABAergic interneurons, the single-cell APA analysis not only identified cell type-specific APA events but also revealed that the modality of APA could classify cell subpopulations. Thus, SAPAS will expand our understanding of cellular heterogeneity.
| S-EPMC8299648 | biostudies-literature
Project description:Alternative polyadenylation (APA) is an important post-transcriptional mechanism that has major implications in biological processes and diseases. Although specialized sequencing methods for polyadenylation exist, availability of these data are limited compared to RNA-sequencing data. We developed REPAC, a framework for the analysis of APA from RNA-sequencing data. Using REPAC, we investigate the landscape of APA caused by activation of B cells. We also show that REPAC is faster than alternative methods by at least 7-fold and that it scales well to hundreds of samples. Overall, the REPAC method offers an accurate, easy, and convenient solution for the exploration of APA.
| S-EPMC9912678 | biostudies-literature
Project description:BackgroundAlternative polyadenylation (APA) has emerged as a pervasive mechanism that contributes to the transcriptome complexity and dynamics of gene regulation. The current tsunami of whole genome poly(A) site data from various conditions generated by 3' end sequencing provides a valuable data source for the study of APA-related gene expression. Cluster analysis is a powerful technique for investigating the association structure among genes, however, conventional gene clustering methods are not suitable for APA-related data as they fail to consider the information of poly(A) sites (e.g., location, abundance, number, etc.) within each gene or measure the association among poly(A) sites between two genes.ResultsHere we proposed a computational framework, named PASCCA, for clustering genes from replicated or unreplicated poly(A) site data using canonical correlation analysis (CCA). PASCCA incorporates multiple layers of gene expression data from both the poly(A) site level and gene level and takes into account the number of replicates and the variability within each experimental group. Moreover, PASCCA characterizes poly(A) sites in various ways including the abundance and relative usage, which can exploit the advantages of 3' end deep sequencing in quantifying APA sites. Using both real and synthetic poly(A) site data sets, the cluster analysis demonstrates that PASCCA outperforms other widely-used distance measures under five performance metrics including connectivity, the Dunn index, average distance, average distance between means, and the biological homogeneity index. We also used PASCCA to infer APA-specific gene modules from recently published poly(A) site data of rice and discovered some distinct functional gene modules. We have made PASCCA an easy-to-use R package for APA-related gene expression analyses, including the characterization of poly(A) sites, quantification of association between genes, and clustering of genes.ConclusionsBy providing a better treatment of the noise inherent in repeated measurements and taking into account multiple layers of poly(A) site data, PASCCA could be a general tool for clustering and analyzing APA-specific gene expression data. PASCCA could be used to elucidate the dynamic interplay of genes and their APA sites among various biological conditions from emerging 3' end sequencing data to address the complex biological phenomenon.
| S-EPMC6343338 | biostudies-literature
Project description:Alternative polyadenylation (APA) is increasingly recognized to regulate gene expression across different cell types, but obtaining APA maps from individual cell types typically requires prior purification, a stressful procedure that can itself alter cellular states. Here, we describe a new platform, cTag-PAPERCLIP, that generates APA profiles from single cell populations in intact tissues; cTag-PAPERCLIP requires no tissue dissociation and preserves transcripts in native states. Applying cTag-PAPERCLIP to profile four major cell types in the mouse brain revealed common APA preferences between excitatory and inhibitory neurons distinct from astrocytes and microglia, regulated in part by neuron-specific RNA-binding proteins NOVA2 and PTBP2. We further identified a role of APA in switching Araf protein isoforms during microglia activation, impacting production of downstream inflammatory cytokines. Our results demonstrate the broad applicability of cTag-PAPERCLIP and a previously undiscovered role of APA in contributing to protein diversity between different cell types and cellular states within the brain.
| S-EPMC5637551 | biostudies-literature
Project description:RNA stability plays an important role in gene expression. Here, using 3' end sequencing of newly made and pre-existing poly(A)+ RNAs, we compare transcript stability in multiple human cell lines, including HEK293T, HepG2, and SH-SY5Y. We show that while mRNA stability is generally conserved across the cell lines, specific transcripts having a high GC content and possibly more stable secondary RNA structures are relatively more stable in SH-SY5Y cells compared to the other 2 cell lines. These features also differentiate stability levels of alternative polyadenylation (APA) 3'UTR isoforms in a cell type-specific manner. Using differentiation of a neural stem cell line as a model, we show that mRNA stability difference could contribute to gene expression changes in neurogenesis and confirm the neuronal identity of SH-SY5Y cells at both gene expression and APA levels. In addition, compared to transcripts using 3'-most exon cleavage/polyadenylation sites (PASs), those using intronic PASs are generally less stable, especially when the PAS-containing intron is large and has a strong 5' splice site, suggesting that intronic polyadenylation mostly plays a negative role in gene expression. Interestingly, the differential mRNA stability among APA isoforms appears to buffer PAS choice in these cell lines. Moreover, we found that several other poly(A)+ RNA species, including promoter-associated long noncoding RNAs and transcripts encoded by the mitochondrial genome, are more stable in SH-SY5Y cells than the other 2 cell lines, further highlighting distinct RNA metabolism in neuronal cells. Together, our results indicate that distinct RNA stability control in neuronal cells may contribute to the gene expression and APA programs that define their cell identity.
| S-EPMC9159357 | biostudies-literature