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

BMC Bioinformatics

Software

Ting-Wen Chen¹  Po-Jung Huang¹  Chi-Ching Lee¹  Yuan-Ming Yeh¹  Petrus Tang²  Ruei-Chi Gan³  Hsien-Da Huang⁴  Timothy H. Wu⁵  Cheng-Hsun Chiu⁶ 

[1] Bioinformatics Center, Molecular Medicine Research Center, Chang Gung University, Taoyuan, Taiwan;Bioinformatics Center, Molecular Medicine Research Center, Chang Gung University, Taoyuan, Taiwan;Molecular Infectious Diseases Research Center, Chang Gung Memorial Hospital, Taoyuan, Taiwan;Molecular Regulation & Bioinformatics Laboratory, Chang Gung University, Taoyuan, Taiwan;Department of Biological Science and Technology, National Chiao Tung University, 300, Hsin-Chu, Taiwan;Bioinformatics Center, Molecular Medicine Research Center, Chang Gung University, Taoyuan, Taiwan;Department of Biological Science and Technology, National Chiao Tung University, 300, Hsin-Chu, Taiwan;Institute of Bioinformatics and Systems Biology, National Chiao Tung University, 300, Hsin-Chu, Taiwan;Institute of Biomedical Informatics, National Yang-Ming University, Taipei City, Taiwan;Molecular Infectious Diseases Research Center, Chang Gung Memorial Hospital, Taoyuan, Taiwan;

关键词: Comparative transcriptome;    Transcriptome quantification;    De novo transcriptome assembly;    Non-model transcriptome;    Web service;   

DOI  :  10.1186/s12859-016-1366-1

来源: Springer

【 摘 要 】

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.

【 授权许可】

CC BY   
© The Author(s). 2016