Streamlining differential exon and 3' UTR usage with diffUTR.
Ontology highlight
ABSTRACT: Background
Despite the importance of alternative poly-adenylation and 3' UTR length for a variety of biological phenomena, there are limited means of detecting UTR changes from standard transcriptomic data.Results
We present the diffUTR Bioconductor package which streamlines and improves upon differential exon usage (DEU) analyses, and leverages existing DEU tools and alternative poly-adenylation site databases to enable differential 3' UTR usage analysis. We demonstrate the diffUTR features and show that it is more flexible and more accurate than state-of-the-art alternatives, both in simulations and in real data.Conclusions
diffUTR enables differential 3' UTR analysis and more generally facilitates DEU and the exploration of their results.
SUBMITTER: Gerber S
PROVIDER: S-EPMC8045333 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Alternative usage of exons provides genomes with plasticity to produce different transcripts from the same gene, modulating the function, localization, and life cycle of gene products. It affects most human genes. For a limited number of cases, alternative functions and tissue-specific roles are known. However, recent high-throughput sequencing studies have suggested that much alternative isoform usage across tissues is nonconserved, raising the question of the extent of its functional importance. We address this question in a genome-wide manner by analyzing the transcriptomes of five tissues for six primate species, focusing on exons that are 1:1 orthologous in all six species. Our results support a model in which differential usage of exons has two major modes: First, most of the exons show only weak differences, which are dominated by interspecies variability and may reflect neutral drift and noisy splicing. These cases dominate the genome-wide view and explain why conservation appears to be so limited. Second, however, a sizeable minority of exons show strong differences between tissues, which are mostly conserved. We identified a core set of 3,800 exons from 1,643 genes that show conservation of strongly tissue-dependent usage patterns from human at least to macaque. This set is enriched for exons encoding protein-disordered regions and untranslated regions. Our findings support the theory that isoform regulation is an important target of evolution in primates, and our method provides a powerful tool for discovering potentially functional tissue-dependent isoforms.
| S-EPMC3780897 | biostudies-other
Project description:Large alternatively spliced internal exons are uncommon in vertebrate genes, and the mechanisms governing their usage are unknown. In this report, we examined alternative splicing of a 1-kb internal exon from the human caldesmon gene containing two regulated 5' splice sites that are 687 nucleotides apart. In cell lines normally splicing caldesmon RNA via utilization of the exon-internal 5' splice site, inclusion of the differential exon required a long purine-rich sequence located between the two competing 5' splice sites. This element consisted of four identical 32-nucleotide purine-rich repeats that resemble exon-splicing enhancers (ESE) identified in other genes. One 32-nucleotide repeat supported exon inclusion, repressed usage of the terminal 5' splice site, and functioned in a heterologous exon dependent on exon enhancers for inclusion, indicating that the caldesmon purine-rich sequence can be classified as an ESE. The ESE was required for utilization of the internal 5' splice site only in the presence of the competing 5' splice site and had no effect when placed downstream of the terminal 5' splice site. In the absence of the internal 5' splice site, the ESE activated a normally silent cryptic 5' splice site near the natural internal 5' splice site, indicating that the ESE stimulates upstream 5' splice site selection. We propose that the caldesmon ESE functions to regulate competition between two 5' splice sites within a differential internal exon.
| S-EPMC230637 | biostudies-other
Project description:BACKGROUND:With the advent of the GeneChip Exon Arrays, it is now possible to extract "exon-level" expression estimates, allowing for detection of alternative splicing events, one of the primary mechanisms of transcript diversity. In the context of (1) a complex trait use case and (2) a human cerebellum vs. heart comparison on previously validated data, we present a transcript-based statistical model and validation framework to allow detection of alternative exon usage (AEU) between different groups. To illustrate the approach, we detect and confirm differences in exon usage in the two of the most widely studied mouse genetic models (the C57BL/6J and DBA/2J inbred strains) and in a human dataset. RESULTS:We developed a computational framework that consists of probe level annotation mapping and statistical modeling to detect putative AEU events, as well as visualization and alignment with known splice events. We show a dramatic improvement (?25 fold) in the ability to detect these events using the appropriate annotation and statistical model which is actually specified at the transcript level, as compared with the transcript cluster/gene-level annotation used on the array. An additional component of this workflow is a probe index that allows ranking AEU candidates for validation and can aid in identification of false positives due to single nucleotide polymorphisms. DISCUSSION:Our work highlights the importance of concordance between the functional unit interrogated (e.g., gene, transcripts) and the entity (e.g., exon, probeset) within the statistical model. The framework we present is broadly applicable to other platforms (including RNAseq).
| S-EPMC3097375 | biostudies-literature
Project description:Neuronal activity-regulated gene transcription underlies plasticity-dependent changes in the molecular composition and structure of neurons. A large number of genes regulated by different neuronal plasticity inducing pathways have been identified, but altered gene expression levels represent only part of the complexity of the activity-regulated transcriptional program. Alternative splicing, the differential inclusion and exclusion of exonic sequence in mRNA, is an additional mechanism that is thought to define the activity-dependent transcriptome. Here, we present a genome wide microarray-based survey to identify exons with increased expression levels at 1, 4 or 8 h following neuronal activity in the murine hippocampus provoked by generalized seizures. We used two different bioinformatics approaches to identify alternative activity-induced exon usage and to predict alternative splicing, ANOSVA (ANalysis Of Splicing VAriation) which we here adjusted to accommodate data from different time points and FIRMA (Finding Isoforms using Robust Multichip Analysis). RNA sequencing, in situ hybridization and reverse transcription PCR validate selected activity-dependent splicing events of previously described and so far undescribed activity-regulated transcripts, including Homer1a, Homer1d, Ania3, Errfi1, Inhba, Dclk1, Rcan1, Cda, Tpm1 and Krt75. Taken together, our survey significantly adds to the comprehensive understanding of the complex activity-dependent neuronal transcriptomic signature. In addition, we provide data sets that will serve as rich resources for future comparative expression analyses.
| S-EPMC7656758 | biostudies-literature
Project description:This article contains data on RDM16 and STA1 regulate differential usage of exon/intron in RNA directed DNA Methylation pathway (RdDM) (Sharma et al., 2016) [5]. This data include expression profiles of top 100 genes that has at least one exon or intron differentially expressed in three different contrast, i.e., WT (Wild type) vs RDM16, WT vs STA1, and RDM16 vs STA1. Also we included the alignment of MORC6 protein to the ATPase-C family members that have conserved three ATP binding sites and conserved Mg2+ binding sites in the spliced exon.
| S-EPMC5403763 | biostudies-literature
Project description:Shrew-1, also called AJAP1, is a transmembrane protein associated with E-cadherin-mediated adherence junctions and a putative tumor suppressor. Apart from its interaction with ?-catenin and involvement in E-cadherin internalization, little structure or function information exists. Here we explored shrew-1 expression during postnatal differentiation of mammary gland as a model system. Immunohistological analyses with antibodies against either the extracellular or the cytoplasmic domains of shrew-1 consistently revealed the expression of full-length shrew-1 in myoepithelial cells, but only part of it in luminal cells. While shrew-1 localization remained unaltered in myoepithelial cells, nuclear localization occurred in luminal cells during lactation. Based on these observations, we identified two unknown shrew-1 transcript variants encoding N-terminally truncated proteins. The smallest shrew-1 protein lacks the extracellular domain and is most likely the only variant present in luminal cells. RNA analyses of human tissues confirmed that the novel transcript variants of shrew-1 exist in vivo and exhibit a differential tissue expression profile. We conclude that our findings are essential for the understanding and interpretation of future functional and interactome analyses of shrew-1 variants.
| S-EPMC5155531 | biostudies-literature
Project description:Epigenomics and transcriptomics data from high-throughput sequencing techniques such as RNA-seq and ChIP-seq have been successfully applied in predicting gene transcript expression. However, the locations of chromatin loops in the genome identified by techniques such as Chromatin Interaction Analysis with Paired End Tag sequencing (ChIA-PET) have never been used for prediction tasks. Here, we developed machine learning models to investigate if ChIA-PET could contribute to transcript and exon usage prediction. In doing so, we used a large set of transcription factors as well as ChIA-PET data. We developed different Gradient Boosting Trees models according to the different tasks with the integrated datasets from three cell lines, including GM12878, HeLaS3 and K562. We validated the models via 10-fold cross validation, chromosome-split validation and cross-cell validation. Our results show that both transcript and splicing-derived exon usage can be effectively predicted with at least 0.7512 and 0.7459 of accuracy, respectively, on all cell lines from all kinds of validations. Examining the predictive features, we found that RNA Polymerase II ChIA-PET was one of the most important features in both transcript and exon usage prediction, suggesting that chromatin loop anchors are predictive of both transcript and exon usage.
| S-EPMC8575016 | biostudies-literature
Project description:Kanate was used to induce seizures in mice, and gene expression in the hippocampus was measured
2020-11-18 | GSE148028 | GEO
Project description:The methyl-CpG-binding protein 2 (MECP2) is a critical global regulator of gene expression. Mutations in MECP2 cause neurodevelopmental disorders including Rett syndrome (RTT). MECP2 exon 2 is spliced into two alternative messenger ribonucleic acid (mRNA) isoforms encoding MECP2-E1 or MECP2-E2 protein isoforms that differ in their N-termini. MECP2-E2, isolated first, was used to define the general roles of MECP2 in methyl-deoxyribonucleic acid (DNA) binding, targeting of transcriptional regulatory complexes, and its disease-causing impact in RTT. It was later found that MECP2-E1 is the most abundant isoform in the brain and its exon 1 is also mutated in RTT. MECP2 transcripts undergo alternative polyadenylation generating mRNAs with four possible 3'untranslated region (UTR) lengths ranging from 130 to 8600 nt. Together, the exon and 3'UTR isoforms display remarkable abundance disparity across cell types and tissues during development. These findings indicate discrete means of regulation and suggest that protein isoforms perform non-overlapping roles. Multiple regulatory programs have been explored to explain these disparities. DNA methylation patterns of the MECP2 promoter and first intron impact MECP2-E1 and E2 isoform levels. Networks of microRNAs and RNA-binding proteins also post-transcriptionally regulate the stability and translation efficiency of MECP2 3'UTR isoforms. Finally, distinctions in biophysical properties in the N-termini between MECP2-E1 and E2 lead to variable protein stabilities and DNA binding dynamics. This review describes the steps taken from the discovery of MECP2, the description of its key functions, and its association with RTT, to the emergence of evidence revealing how MECP2 isoforms are differentially regulated at the transcriptional, post-transcriptional and post-translational levels.
| S-EPMC7530521 | biostudies-literature
Project description:Alternative gene splicing is a common phenomenon in which a single gene gives rise to multiple transcript isoforms. The process is strictly guided and involves a multitude of proteins and regulatory complexes. Unfortunately, aberrant splicing events have been linked to genetic disorders. Therefore, understanding mechanisms of alternative splicing regulation and differences in splicing events between diseased and healthy tissues is crucial in advancing personalized medicine and drug developments. We propose a linear mixed model, Random Effects for the Identification of Differential Splicing (REIDS), for the identification of alternative splicing events using Human Transcriptome Arrays (HTA). For each exon, a splicing score is calculated based on two scores, an exon score and an array score. The junction information is used to rank the identified exons from strongly confident to less confident candidates for alternative splicing. The design of junctions was also discussed to highlight the complexity of exon-exon and exon-junction interactions. Based on a list of Rt-PCR validated probe sets, REIDS outperforms AltAnalyze and iGems in the % recall rate.
| S-EPMC5974242 | biostudies-literature