PARRoT- a homology-based strategy to quantify and compare RNA-sequencing from non-model organisms

Abstract

BackgroundNext-generation sequencing promises the de novo genomic and transcriptomic analysis of samples of interests. However, there are only a few organisms having reference genomic sequences and even fewer having well-defined or curated annotations. For transcriptome studies focusing on organisms lacking proper reference genomes, the common strategy is de novo assembly followed by functional annotation. However, things become even more complicated when multiple transcriptomes are compared.ResultsHere, we propose a new analysis strategy and quantification methods for quantifying expression level which not only generate a virtual reference from sequencing data, but also provide comparisons between transcriptomes. First, all reads from the transcriptome datasets are pooled together for de novo assembly. The assembled contigs are searched against NCBI NR databases to find potential homolog sequences. Based on the searched result, a set of virtual transcripts are generated and served as a reference transcriptome. By using the same reference, normalized quantification values including RC (read counts), eRPKM (estimated RPKM) and eTPM (estimated TPM) can be obtained that are comparable across transcriptome datasets. In order to demonstrate the feasibility of our strategy, we implement it in the web service PARRoT. PARRoT stands for Pipeline for Analyzing RNA Reads of Transcriptomes. It analyzes gene expression profiles for two transcriptome sequencing datasets. For better understanding of the biological meaning from the comparison among transcriptomes, PARRoT further provides linkage between these virtual transcripts and their potential function through showing best hits in SwissProt, NR database, assigning GO terms. Our demo datasets showed that PARRoT can analyze two paired-end transcriptomic datasets of approximately 100 million reads within just three hours.ConclusionsIn this study, we proposed and implemented a strategy to analyze transcriptomes from non-reference organisms which offers the opportunity to quantify and compare transcriptome profiles through a homolog based virtual transcriptome reference. By using the homolog based reference, our strategy effectively avoids the problems that may cause from inconsistencies among transcriptomes. This strategy will shed lights on the field of comparative genomics for non-model organism. We have implemented PARRoT as a web service which is freely available at http://parrot.cgu.edu.tw.

Long non-coding RNA ZFAS1 sponges miR-486 to promote osteosarcoma cells progression and metastasis in vitro and vivo

“…In the past decades, accompany with the fast innovation and development of sequencing technique, many important progression have been achieved [ 12 , 13 ]. High-throughput sequencing and bioinformatics analysis powerfully assist the non-coding RNA research on the epigenetic modification, discovering huge amounts of new identified functional lncRNAs in tumorigenesis [ 14 ].…”

Comparative transcriptome analysis of shoot and root tissue of Bacopa monnieri identifies potential genes related to triterpenoid saponin biosynthesis

“…Efforts were also made to identify potential transcription factors, CYP450s and UGTs that play a key role in regulation and diversification of secondary metabolites. Paired-end sequencing strategy extends the mapped fragment length to 200–500 bp so it is expected to be useful in obtaining longer reads, allows for detecting alternate splice junctions, deletions, insertions, and is useful for de novo transcriptome assembly [ 28 ]. Transcripts generated after assembly were functionally elucidated in different gene ontology, some of the secondary metabolic pathway related transcripts were selected for further validation.…”

Bioinformatics and systems biology research update from the 15th International Conference on Bioinformatics (InCoB2016)

“…The latter requires less memory than BBb as the partial digestion problem is solved in two stages in breadth-first mode. Gan et al [5] have developed a comparative transcriptomic analysis webserver, PARRoT, for non-reference organisms, using a homologue-based virtual transcriptome reference, which provided homologues from sequence databases as well as gene ontology (GO) terms. For entire microbial communities, Xie et al [6] have developed RiboTagger, for fast and accurate recovery of small subunit ribosomal RNA sequences, from all three domains of life.…”