scholarly journals Terminitor: Cleavage Site Prediction Using Deep Learning Models

2019 ◽  
Author(s):  
Chen Yang ◽  
Chenkai Li ◽  
Ka Ming Nip ◽  
René L Warren ◽  
Inanc Birol
Keyword(s):  
Experimental Data ◽  
Gene Regulation ◽  
Deep Learning ◽  
Cleavage Site ◽  
Differential Expression Analysis ◽  
Underlying Mechanism ◽  
Library Preparation ◽  
Rna Seq ◽  
Sequencing Data ◽  
Reference Transcriptome

AbstractAs a widespread RNA processing machinery, alternative polyadenylation plays a crucial role in gene regulation. To help decipher its underlying mechanism and understand its impact, it is desirable to comprehensively profile 3’-untranslated region cleavage and associated polyadenylation sites. State-of-the-art polyadenylation site detection tools are known to be influenced by library preparation artefacts or manually selected features. Moreover, recently published machine learning methods have only been tested on pre-constructed datasets, thus lacking validation on experimental data. Here we present Terminitor, the first deep neural network-based profiling pipeline to make predictions from RNA-seq data. We show how Terminitor outperforms competing tools in sensitivity and precision on experimental transcriptome sequencing data, and demonstrate its use with data from short- and long-read sequencing technologies. For species without a good reference transcriptome annotation, Terminitor is still able to pass on the information learnt from a related species and make reasonable predictions. We used Terminitor to showcase how single nucleotide variations can create or destroy polyadenylated cleavage sites in human RNA-seq samples.Author Summary3’ cleavage and polyadenylation of pre-mRNA is part of RNA maturation process. One gene can be cleaved at different positions at its 3’ end, namely alternatively polyadenylation, thus identifying the correct polyadenylated cleavage site (poly(A) CS) is essential to unveil its role in gene regulation under different physiological and pathological conditions. The current poly(A) CS prediction tools are either heavily influenced by RNA-Seq library preparation artefacts or have only been designed and tested on ad hoc datasets, lacking association with real world applications. In this study, we present a deep learning model, Terminitor, that predicts the probability of a nucleotide sequence containing a poly(A) CS, and validated its performance on human and mouse data. Along with the model, we propose a poly(A) CS profiling pipeline for RNA-seq data. We benchmarked our pipeline against competing tools and achieved higher sensitivity and precision in experimental data. The usage of Terminitor is not limited to genome and transcriptome annotation and we expect it to facilitate the identification of novel isoforms, improve the accuracy of transcript quantification and differential expression analysis, and contribute to the repertoire of reference transcriptome annotation.

scholarly journals The long and the short of it: unlocking nanopore long-read RNA sequencing data with short-read differential expression analysis tools

2021 ◽  
Vol 3 (2) ◽  
Author(s):  
Xueyi Dong ◽  
Luyi Tian ◽  
Quentin Gouil ◽  
Hasaru Kariyawasam ◽  
Shian Su ◽  
...  
Keyword(s):  
Differential Expression ◽  
Expression Analysis ◽  
Differential Expression Analysis ◽  
Transcriptomic Analysis ◽  
Statistical Testing ◽  
Rna Seq ◽  
Sequencing Data ◽  
Short Read ◽  
Sequencing Platform ◽  
Long Read

Abstract Application of Oxford Nanopore Technologies’ long-read sequencing platform to transcriptomic analysis is increasing in popularity. However, such analysis can be challenging due to the high sequence error and small library sizes, which decreases quantification accuracy and reduces power for statistical testing. Here, we report the analysis of two nanopore RNA-seq datasets with the goal of obtaining gene- and isoform-level differential expression information. A dataset of synthetic, spliced, spike-in RNAs (‘sequins’) as well as a mouse neural stem cell dataset from samples with a null mutation of the epigenetic regulator Smchd1 was analysed using a mix of long-read specific tools for preprocessing together with established short-read RNA-seq methods for downstream analysis. We used limma-voom to perform differential gene expression analysis, and the novel FLAMES pipeline to perform isoform identification and quantification, followed by DRIMSeq and limma-diffSplice (with stageR) to perform differential transcript usage analysis. We compared results from the sequins dataset to the ground truth, and results of the mouse dataset to a previous short-read study on equivalent samples. Overall, our work shows that transcriptomic analysis of long-read nanopore data using long-read specific preprocessing methods together with short-read differential expression methods and software that are already in wide use can yield meaningful results.

scholarly journals HTSeq - A Python framework to work with high-throughput sequencing data

2014 ◽  
Author(s):  
Simon Anders ◽  
Paul Theodor Pyl ◽  
Wolfgang Huber
Keyword(s):  
High Throughput ◽  
High Throughput Sequencing ◽  
Rapid Development ◽  
Differential Expression Analysis ◽  
Rna Seq ◽  
Sequencing Data ◽  
Standard Work ◽  
Data Formats ◽  
High Throughput Sequencing Data ◽  
Python Package

Motivation: A large choice of tools exists for many standard tasks in the analysis of high-throughput sequencing (HTS) data. However, once a project deviates from standard work flows, custom scripts are needed. Results: We present HTSeq, a Python library to facilitate the rapid development of such scripts. HTSeq offers parsers for many common data formats in HTS projects, as well as classes to represent data such as genomic coordinates, sequences, sequencing reads, alignments, gene model information, variant calls, and provides data structures that allow for querying via genomic coordinates. We also present htseq-count, a tool developed with HTSeq that preprocesses RNA-Seq data for differential expression analysis by counting the overlap of reads with genes. Availability: HTSeq is released as open-source software under the GNU General Public Licence and available from http://www-huber.embl.de/HTSeq or from the Python Package Index, https://pypi.python.org/pypi/HTSeq

scholarly journals Identifying inaccuracies in gene expression estimates from unstranded RNA-seq data

2019 ◽  
Vol 9 (1) ◽  
Author(s):  
Mikhail Pomaznoy ◽  
Ashu Sethi ◽  
Jason Greenbaum ◽  
Bjoern Peters
Keyword(s):  
Gene Expression ◽  
Differential Expression Analysis ◽  
Cell Types ◽  
Library Preparation ◽  
Rna Seq ◽  
Protein Coding ◽  
Protein Coding Genes ◽  
Machine Learning Model ◽  
Specific Manner ◽  
Library Preparation Protocol

Abstract RNA-seq methods are widely utilized for transcriptomic profiling of biological samples. However, there are known caveats of this technology which can skew the gene expression estimates. Specifically, if the library preparation protocol does not retain RNA strand information then some genes can be erroneously quantitated. Although strand-specific protocols have been established, a significant portion of RNA-seq data is generated in non-strand-specific manner. We used a comprehensive stranded RNA-seq dataset of 15 blood cell types to identify genes for which expression would be erroneously estimated if strand information was not available. We found that about 10% of all genes and 2.5% of protein coding genes have a two-fold or higher difference in estimated expression when strand information of the reads was ignored. We used parameters of read alignments of these genes to construct a machine learning model that can identify which genes in an unstranded dataset might have incorrect expression estimates and which ones do not. We also show that differential expression analysis of genes with biased expression estimates in unstranded read data can be recovered by limiting the reads considered to those which span exonic boundaries. The resulting approach is implemented as a package available at https://github.com/mikpom/uslcount.

scholarly journals Comparative evaluation of RNA-Seq library preparation methods for strand-specificity and low input

2019 ◽  
Vol 9 (1) ◽  
Author(s):  
Dimitra Sarantopoulou ◽  
Soon Yew Tang ◽  
Emanuela Ricciotti ◽  
Nicholas F. Lahens ◽  
Damien Lekkas ◽  
...  
Keyword(s):  
Comparative Evaluation ◽  
Differential Expression Analysis ◽  
Functional Results ◽  
Starting Material ◽  
The Other ◽  
Library Preparation ◽  
Rna Seq ◽  
Experimental Conditions ◽  
Preparation Methods ◽  
Low Input

Abstract Library preparation is a key step in sequencing. For RNA sequencing there are advantages to both strand specificity and working with minute starting material, yet until recently there was no kit available enabling both. The Illumina TruSeq stranded mRNA Sample Preparation kit (TruSeq) requires abundant starting material while the Takara Bio SMART-Seq v4 Ultra Low Input RNA kit (V4) sacrifices strand specificity. The SMARTer Stranded Total RNA-Seq Kit v2 - Pico Input Mammalian (Pico) by Takara Bio claims to overcome these limitations. Comparative evaluation of these kits is important for selecting the appropriate protocol. We compared the three kits in a realistic differential expression analysis. We prepared and sequenced samples from two experimental conditions of biological interest with each of the three kits. We report differences between the kits at the level of differential gene expression; for example, the Pico kit results in 55% fewer differentially expressed genes than TruSeq. Nevertheless, the agreement of the observed enriched pathways suggests that comparable functional results can be obtained. In summary we conclude that the Pico kit sufficiently reproduces the results of the other kits at the level of pathway analysis while providing a combination of options that is not available in the other kits.

scholarly journals WIND (Workflow for pIRNAs aNd beyonD): a strategy for in-depth analysis of small RNA-seq data

F1000Research ◽  
2021 ◽  
Vol 10 ◽  
pp. 1
Author(s):  
Konstantinos Geles ◽  
Domenico Palumbo ◽  
Assunta Sellitto ◽  
Giorgio Giurato ◽  
Eleonora Cianflone ◽  
...  
Keyword(s):  
Small Rna ◽  
Differential Expression Analysis ◽  
Small Rna Sequencing ◽  
Rna Seq ◽  
Sequencing Data ◽  
Transcript Quantification ◽  
Annotation Track ◽  
Depth Analysis ◽  
Exploratory Data ◽  
Downstream Analysis

Current bioinformatics workflows for PIWI-interacting RNA (piRNA) analysis focus primarily on germline-derived piRNAs and piRNA-clusters. Frequently, they suffer from outdated piRNA databases, questionable quantification methods, and lack of reproducibility. Often, pipelines specific to miRNA analysis are used for the piRNA research in silico. Furthermore, the absence of a well-established database for piRNA annotation, as for miRNA, leads to uniformity issues between studies and generates confusion for data analysts and biologists. For these reasons, we have developed WIND (Workflow for pIRNAs aNd beyonD), a bioinformatics workflow that addresses the crucial issue of piRNA annotation, thereby allowing a reliable analysis of small RNA sequencing data for the identification of piRNAs and other small non-coding RNAs (sncRNAs) that in the past have been incorrectly classified as piRNAs. WIND allows the creation of a comprehensive annotation track of sncRNAs combining information available in RNAcentral, with piRNA sequences from piRNABank, the first database dedicated to piRNA annotation. WIND was built with Docker containers for reproducibility and integrates widely used bioinformatics tools for sequence alignment and quantification. In addition, it includes Bioconductor packages for exploratory data and differential expression analysis. Moreover, WIND implements a "dual" approach for the evaluation of sncRNAs expression level quantifying the aligned reads to the annotated genome and carrying out an alignment-free transcript quantification using reads mapped to the transcriptome. Therefore, a broader range of piRNAs can be annotated, improving their quantification and easing the subsequent downstream analysis. WIND performance has been tested with several small RNA-seq datasets, demonstrating how our approach can be a useful and comprehensive resource to analyse piRNAs and other classes of sncRNAs.

scholarly journals A relative comparison between Hidden Markov- and Log-Linear-based models for differential expression analysis in a real time course RNA sequencing data

2018 ◽  
Author(s):  
Fatemeh Gholizadeh ◽  
Zahra Salehi ◽  
Ali Mohammad banaei-Moghaddam ◽  
Abbas Rahimi Foroushani ◽  
Kaveh kavousi
Keyword(s):  
Real Time ◽  
Differential Expression ◽  
Rna Sequencing ◽  
Time Course ◽  
Hidden Markov ◽  
Differential Expression Analysis ◽  
Rna Seq ◽  
Sequencing Data ◽  
Normalization Methods ◽  
Log Linear

AbstractWith the advent of the Next Generation Sequencing technologies, RNA-seq has become known as an optimal approach for studying gene expression profiling. Particularly, time course RNA-seq differential expression analysis has been used in many studies to identify candidate genes. However, applying a statistical method to efficiently identify differentially expressed genes (DEGs) in time course studies is challenging due to inherent characteristics of such data including correlation and dependencies over time. Here we aim to relatively compare EBSeq-HMM, a Hidden Markov-based model, with multiDE, a Log-Linear-based model, in a real time course RNA sequencing data. In order to conduct the comparison, common DEGs detected by edgeR, DESeq2 and Voom (referred to as Benchmark DEGs) were utilized as a measure. Each of the two models were compared using different normalization methods. The findings revealed that multiDE identified more Benchmark DEGs and showed a higher agreement with them than EBSeq-HMM. Furthermore, multiDE and EBSeq-HMM displayed their best performance using TMM and Upper-Quartile normalization methods, respectively.

scholarly journals WIND (Workflow for pIRNAs aNd beyonD): a strategy for in-depth analysis of small RNA-seq data

F1000Research ◽  
2021 ◽  
Vol 10 ◽  
pp. 1
Author(s):  
Konstantinos Geles ◽  
Domenico Palumbo ◽  
Assunta Sellitto ◽  
Giorgio Giurato ◽  
Eleonora Cianflone ◽  
...  
Keyword(s):  
Small Rna ◽  
Differential Expression Analysis ◽  
Small Rna Sequencing ◽  
Rna Seq ◽  
Sequencing Data ◽  
Transcript Quantification ◽  
Annotation Track ◽  
Depth Analysis ◽  
Exploratory Data ◽  
Downstream Analysis

Current bioinformatics workflows for PIWI-interacting RNA (piRNA) analysis focus primarily on germline-derived piRNAs and piRNA-clusters. Frequently, they suffer from outdated piRNA databases, questionable quantification methods, and lack of reproducibility. Often, pipelines specific to miRNA analysis are used for the piRNA research in silico. Furthermore, the absence of a well-established database for piRNA annotation, as for miRNA, leads to uniformity issues between studies and generates confusion for data analysts and biologists. For these reasons, we have developed WIND (Workflow for pIRNAs aNd beyonD), a bioinformatics workflow that addresses the crucial issue of piRNA annotation, thereby allowing a reliable analysis of small RNA sequencing data for the identification of piRNAs and other small non-coding RNAs (sncRNAs) that in the past have been incorrectly classified as piRNAs. WIND allows the creation of a comprehensive annotation track of sncRNAs combining information available in RNAcentral, with piRNA sequences from piRNABank, the first database dedicated to piRNA annotation. WIND was built with Docker containers for reproducibility and integrates widely used bioinformatics tools for sequence alignment and quantification. In addition, it includes Bioconductor packages for exploratory data and differential expression analysis. Moreover, WIND implements a "dual" approach for the evaluation of sncRNAs expression level quantifying the aligned reads to the annotated genome and carrying out an alignment-free transcript quantification using reads mapped to the transcriptome. Therefore, a broader range of piRNAs can be annotated, improving their quantification and easing the subsequent downstream analysis. WIND performance has been tested with several small RNA-seq datasets, demonstrating how our approach can be a useful and comprehensive resource to analyse piRNAs and other classes of sncRNAs.

Transcriptome Characterization of Panax Quinquefolius L. seedlings in Response to Arbuscular Mycorrhizal Fungi using RNA-seq

2021 ◽  
Author(s):  
Zhifang Ran ◽  
Xiaotong Yang ◽  
Yongqing Zhang ◽  
Jie Zhou
Keyword(s):  
Arbuscular Mycorrhizal Fungi ◽  
Stress Responses ◽  
Mycorrhizal Fungi ◽  
Differential Expression Analysis ◽  
Enrichment Analysis ◽  
Arbuscular Mycorrhizal ◽  
Transcriptome Profiling ◽  
Underlying Mechanism ◽  
Panax Quinquefolius ◽  
Rna Seq

Abstract Panax quinquefolius L. has been considered as an important traditional Chinese medicine with a history of more than 300 years in China. Ginsenoside is the main bioactive component. Our research group has found that the accumulation of ginsenoside could be affected by arbuscular mycorrhizal fungi (AMF). However the underlying mechanism how AMF affected the biosynthesis of ginsenoside in P. quinquefolius is still unclear. In this study, the RNA-seq analysis was used to evaluate the effects of AMF (Rhizophagus intraradices, R. intraradices) on the expression of ginsenoside synthesis related genes in P. quinquefolius root. The results indicated that a symbiotic relationship between R. intraradices and P. quinquefolius was established. RNA-seq achieved approximately 48.62 G reads of all samples. Assembly of all the reads involved in all samples produced 63420 transcripts and 24137 unigenes. Differential expression analysis was performed between the control and AMF group. A total of 111 differentially expressed genes (DEGs) in response to AMF vs control were identified, 78 and 33 transcripts were upregulated and downregulated, respectively. Based on the functional analysis, Gene ontology (GO) analysis revealed that most DEGs were related to stress responses and cellular metabolic processes. The Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis identified transduction, plant hormone signal transduction and terpenoids and polyketides biosynthesis pathways. Furthermore, the expression of glycolysis-related genes and ginsenoside synthesis related genes was largely induced by AMF. In conclusion, our results comprehensively elucidated the molecular mechanism how AMF affected the biosynthesis of ginsenoside in P.quinquefolius by transcriptome profiling.

scholarly journals Differential Expression Analysis for RNA-Seq: An Overview of Statistical Methods and Computational Software

2015 ◽  
Vol 14s1 ◽  
pp. CIN.S21631 ◽  
Author(s):  
Huei-Chung Huang ◽  
Yi Niu ◽  
Li-Xuan Qin
Keyword(s):  
Differential Expression ◽  
Statistical Methods ◽  
Differential Expression Analysis ◽  
Rna Seq ◽  
Sequencing Data ◽  
Computational Tools ◽  
Estimation Algorithms ◽  
Testing Strategies ◽  
Discrete Nature ◽  
Statistical Background

Deep sequencing has recently emerged as a powerful alternative to microarrays for the high-throughput profiling of gene expression. In order to account for the discrete nature of RNA sequencing data, new statistical methods and computational tools have been developed for the analysis of differential expression to identify genes that are relevant to a disease such as cancer. In this paper, it is thus timely to provide an overview of these analysis methods and tools. For readers with statistical background, we also review the parameter estimation algorithms and hypothesis testing strategies used in these methods.

scholarly journals scBatch: batch-effect correction of RNA-seq data through sample distance matrix adjustment

2020 ◽  
Vol 36 (10) ◽  
pp. 3115-3123 ◽  
Author(s):  
Teng Fei ◽  
Tianwei Yu
Keyword(s):  
Single Cell ◽  
Differential Expression Analysis ◽  
Distance Matrix ◽  
Real Data ◽  
R Package ◽  
Batch Effect ◽  
Supplementary Information ◽  
Rna Seq ◽  
Sequencing Data ◽  
Gene Differential Expression

Abstract Motivation Batch effect is a frequent challenge in deep sequencing data analysis that can lead to misleading conclusions. Existing methods do not correct batch effects satisfactorily, especially with single-cell RNA sequencing (RNA-seq) data. Results We present scBatch, a numerical algorithm for batch-effect correction on bulk and single-cell RNA-seq data with emphasis on improving both clustering and gene differential expression analysis. scBatch is not restricted by assumptions on the mechanism of batch-effect generation. As shown in simulations and real data analyses, scBatch outperforms benchmark batch-effect correction methods. Availability and implementation The R package is available at github.com/tengfei-emory/scBatch. The code to generate results and figures in this article is available at github.com/tengfei-emory/scBatch-paper-scripts. Supplementary information Supplementary data are available at Bioinformatics online.

Sign in / Sign up

Export Citation Format

Share Document