Correction To: P_RNA_scaffolder: a Fast and Accurate Genome Scaffolder Using Paired-End RNAsequencing Reads

Bai-Han Zhu,Jun Xiao,Wei Xue,Gui-Cai Xu,Ming-Yuan Sun,Jiong-Tang Li

DOI: https://doi.org/10.1186/s12864-019-5855-2

IF: 4.547

2019-01-01

BMC Genomics

Abstract:Following the publication of this article

What problem does this paper attempt to address?

P_RNA_scaffolder: a Fast and Accurate Genome Scaffolder Using Paired-End RNA-sequencing Reads

Bai-Han Zhu,Jun Xiao,Wei Xue,Gui-Cai Xu,Ming-Yuan Sun,Jiong-Tang Li

DOI: https://doi.org/10.1186/s12864-018-4567-3

IF: 4.547

2018-01-01

BMC Genomics

Abstract:Background Obtaining complete gene structures is one major goal of genome assembly. Some gene regions are fragmented in low quality and high-quality assemblies. Therefore, new approaches are needed to recover gene regions. Genomes are widely transcribed, generating messenger and non-coding RNAs. These widespread transcripts can be used to scaffold genomes and complete transcribed regions. Results We present P_RNA_scaffolder, a fast and accurate tool using paired-end RNA-sequencing reads to scaffold genomes. This tool aims to improve the completeness of both protein-coding and non-coding genes. After this tool was applied to scaffolding human contigs, the structures of both protein-coding genes and circular RNAs were almost completely recovered and equivalent to those in a complete genome, especially for long proteins and long circular RNAs. Tested in various species, P_RNA_scaffolder exhibited higher speed and efficiency than the existing state-of-the-art scaffolders. This tool also improved the contiguity of genome assemblies generated by current mate-pair scaffolding and third-generation single-molecule sequencing assembly. Conclusions The P_RNA_scaffolder can improve the contiguity of genome assembly and benefit gene prediction. This tool is available at http://www.fishbrowser.org/software/P_RNA_scaffolder .
Correction To: L_RNA_scaffolder: Scaffolding Genomes with Transcripts

Wei Xue,Jiong-Tang Li,Yaping Zhu,Guangyuan Hou,Xue Kong,Youyi Kuang,Xiaohui Sun

DOI: https://doi.org/10.1186/s12864-019-5856-1

IF: 4.547

2019-01-01

BMC Genomics

Abstract:Following the publication of this article [1], the authors reported that the link to the software described in the article is no longer valid.
Pseudo-Sanger Sequencing: Massively Parallel Production of Long and Near Error-Free Reads Using NGS Technology

Jue Ruan,Lan Jiang,Zechen Chong,Qiang Gong,Heng Li,Chunyan Li,Yong Tao,Caihong Zheng,Weiwei Zhai,David Turissini,Charles H Cannon,Xuemei Lu,Chung-I Wu

DOI: https://doi.org/10.1186/1471-2164-14-711

IF: 4.547

2013-01-01

BMC Genomics

Abstract:Background Usually, next generation sequencing (NGS) technology has the property of ultra-high throughput but the read length is remarkably short compared to conventional Sanger sequencing. Paired-end NGS could computationally extend the read length but with a lot of practical inconvenience because of the inherent gaps. Now that Illumina paired-end sequencing has the ability of read both ends from 600 bp or even 800 bp DNA fragments, how to fill in the gaps between paired ends to produce accurate long reads is intriguing but challenging. Results We have developed a new technology, referred to as pseudo-Sanger (PS) sequencing. It tries to fill in the gaps between paired ends and could generate near error-free sequences equivalent to the conventional Sanger reads in length but with the high throughput of the Next Generation Sequencing. The major novelty of PS method lies on that the gap filling is based on local assembly of paired-end reads which have overlaps with at either end. Thus, we are able to fill in the gaps in repetitive genomic region correctly. The PS sequencing starts with short reads from NGS platforms, using a series of paired-end libraries of stepwise decreasing insert sizes. A computational method is introduced to transform these special paired-end reads into long and near error-free PS sequences, which correspond in length to those with the largest insert sizes. The PS construction has 3 advantages over untransformed reads: gap filling, error correction and heterozygote tolerance. Among the many applications of the PS construction is de novo genome assembly, which we tested in this study. Assembly of PS reads from a non-isogenic strain of Drosophila melanogaster yields an N50 contig of 190 kb, a 5 fold improvement over the existing de novo assembly methods and a 3 fold advantage over the assembly of long reads from 454 sequencing. Conclusions Our method generated near error-free long reads from NGS paired-end sequencing. We demonstrated that de novo assembly could benefit a lot from these Sanger-like reads. Besides, the characteristic of the long reads could be applied to such applications as structural variations detection and metagenomics.
L_RNA_scaffolder: Scaffolding Genomes with Transcripts (vol 14, 604, 2013)

Wei Xue,Jiong-Tang Li,Ya-Ping Zhu,Guang-Yuan Hou,Xiang-Fei Kong,You-Yi Kuang,Xiao-Wen Sun

DOI: https://doi.org/10.1186/s12864-019-5856-1

IF: 4.547

2019-01-01

BMC Genomics

Abstract:Following the publication of this article [1], the authors reported that the link to the software described in the article is no longer valid.
L_RNA_scaffolder: scaffolding genomes with transcripts

Wei Xue,Jiong-Tang Li,Ya-Ping Zhu,Guang-Yuan Hou,Xiang-Fei Kong,You-Yi Kuang,Xiao-Wen Sun

DOI: https://doi.org/10.1186/1471-2164-14-604

IF: 4.547

2013-01-01

BMC Genomics

Abstract:Background Generation of large mate-pair libraries is necessary for de novo genome assembly but the procedure is complex and time-consuming. Furthermore, in some complex genomes, it is hard to increase the N50 length even with large mate-pair libraries, which leads to low transcript coverage. Thus, it is necessary to develop other simple scaffolding approaches, to at least solve the elongation of transcribed fragments. Results We describe L_RNA_scaffolder, a novel genome scaffolding method that uses long transcriptome reads to order, orient and combine genomic fragments into larger sequences. To demonstrate the accuracy of the method, the zebrafish genome was scaffolded. With expanded human transcriptome data, the N50 of human genome was doubled and L_RNA_scaffolder out-performed most scaffolding results by existing scaffolders which employ mate-pair libraries. In these two examples, the transcript coverage was almost complete, especially for long transcripts. We applied L_RNA_scaffolder to the highly polymorphic pearl oyster draft genome and the gene model length significantly increased. Conclusions The simplicity and high-throughput of RNA-seq data makes this approach suitable for genome scaffolding. L_RNA_scaffolder is available at http://www.fishbrowser.org/software/L_RNA_scaffolder .
Corrigendum: Long-range RNA Pairings Contribute to Mutually Exclusive Splicing

Yue Yuan,Yun Yang,Lanzhi Dai,Guang Cao,Ran Chen,Weiling Hong,Baoping Liu,Yang Shi,Yijun Meng,Feng Shi,Mu Xiao,Yongfeng Jin

DOI: https://doi.org/10.1261/rna.058008.116

2016-01-01

RNA

Abstract:In the Discussion section (first paragraph, third sentence) on page 105 of the above-noted article, the sentence “The RNA architectures identified in our study were not predicted by in silico approaches (Pervouchine 2014), nor were they detected using recent high-throughput RNA structure sequencing techniques (Daines et al. 2011; Li et al. 2012)” was not clear and should instead read as follows: “The RNA bidirectional structural architectures identified in our study were not detected using recent high-throughput RNA structure sequencing techniques (Daines et al. 2011; Li et al. 2012), nor were they discovered by in silico approaches, although stem II of the Drosophila srp pre-mRNA was predicted (Raker et al. 2009; Pervouchine 2014).” Accordingly, the following reference should be included in the References section: Raker VA, Mironov AA, Gelfand MS, Pervouchine DD. 2009. Modulation of alternative splicing by long-range RNA structures in Drosophila. Nucleic Acids Res 37: 4533–4544. The authors apologize for any confusion this may have caused. This article has been corrected in both the PDF and full-text HTML files online.
GeneScissors: a comprehensive approach to detecting and correcting spurious transcriptome inference owing to RNA-seq reads misalignment.

Zhaojun Zhang,Shunping Huang,Jack Wang,Xiang Zhang,Fernando Pardo Manuel de Villena,Leonard McMillan,Wei Wang

DOI: https://doi.org/10.1093/bioinformatics/btt216

IF: 5.8

2013-01-01

Bioinformatics

Abstract:Motivation: RNA-seq techniques provide an unparalleled means for exploring a transcriptome with deep coverage and base pair level resolution. Various analysis tools have been developed to align and assemble RNA-seq data, such as the widely used TopHat/Cufflinks pipeline. A common observation is that a sizable fraction of the fragments/reads align to multiple locations of the genome. These multiple alignments pose substantial challenges to existing RNA-seq analysis tools. Inappropriate treatment may result in reporting spurious expressed genes (false positives) and missing the real expressed genes (false negatives). Such errors impact the subsequent analysis, such as differential expression analysis. In our study, we observe that similar to 3.5% of transcripts reported by TopHat/Cufflinks pipeline correspond to annotated nonfunctional pseudogenes. Moreover, similar to 10.0% of reported transcripts are not annotated in the Ensembl database. These genes could be either novel expressed genes or false discoveries. Results: We examine the underlying genomic features that lead to multiple alignments and investigate how they generate systematic errors in RNA-seq analysis. We develop a general tool, GeneScissors, which exploits machine learning techniques guided by biological knowledge to detect and correct spurious transcriptome inference by existing RNA-seq analysis methods. In our simulated study, GeneScissors can predict spurious transcriptome calls owing to misalignment with an accuracy close to 90%. It provides substantial improvement over the widely used TopHat/Cufflinks or MapSplice/Cufflinks pipelines in both precision and F-measurement. On real data, GeneScissors reports 53.6% less pseudogenes and 0.97% more expressed and annotated transcripts, when compared with the TopHat/Cufflinks pipeline. In addition, among the 10.0% unannotated transcripts reported by TopHat/Cufflinks, GeneScissors finds that > 16.3% of them are false positives.
PEP_scaffolder: Using (homologous) Proteins to Scaffold Genomes

Bai-Han Zhu,Ying-Nan Song,Wei Xue,Gui-Cai Xu,Jun Xiao,Ming-Yuan Sun,Xiao-Wen Sun,Jiong-Tang Li

DOI: https://doi.org/10.1093/bioinformatics/btw378

IF: 5.8

2016-01-01

Bioinformatics

Abstract:Abstract Motivation: Recovering the gene structures is one of the important goals of genome assembly. In low-quality assemblies, and even some high-quality assemblies, certain gene regions are still incomplete; thus, novel scaffolding approaches are required to complete gene regions. Results: We developed an efficient and fast genome scaffolding method called PEP_scaffolder, using proteins to scaffold genomes. The pipeline aims to recover protein-coding gene structures. We tested the method on human contigs; using human UniProt proteins as guides, the improvement on N50 size was 17% increase with an accuracy of ∼97%. PEP_scaffolder improved the proportion of fully covered proteins among all proteins, which was close to the proportion in the finished genome. The method provided a high accuracy of 91% using orthologs of distant species. Tested on simulated fly contigs, PEP_scaffolder outperformed other scaffolders, with the shortest running time and the highest accuracy. Availability and Implementation: The software is freely available at http://www.fishbrowser.org/software/PEP_scaffolder/ Contact: lijt@cafs.ac.cn Supplementary information: Supplementary data are available at Bioinformatics online.
PEP_scaffolder: using (homologous) proteins to scaffold genomes.

Bai-Han Zhu,Ying-Nan Song,Wei Xue,Gui-Cai Xu,Jun Xiao,Ming-Yuan Sun,Xiaowen Sun,Jiongtang Li

DOI: https://doi.org/10.1093/bioinformatics/btw378

IF: 5.8

2016-01-01

Bioinformatics

Abstract:Motivation: Recovering the gene structures is one of the important goals of genome assembly. In low-quality assemblies, and even some high-quality assemblies, certain gene regions are still incomplete; thus, novel scaffolding approaches are required to complete gene regions. Results: We developed an efficient and fast genome scaffolding method called PEP_scaffolder, using proteins to scaffold genomes. The pipeline aims to recover protein-coding gene structures. We tested the method on human contigs; using human UniProt proteins as guides, the improvement on N50 size was 17% increase with an accuracy of similar to 97%. PEP_scaffolder improved the proportion of fully covered proteins among all proteins, which was close to the proportion in the finished genome. The method provided a high accuracy of 91% using orthologs of distant species. Tested on simulated fly contigs, PEP_scaffolder outperformed other scaffolders, with the shortest running time and the highest accuracy.
ResSeq: Enhancing Short-Read Sequencing Alignment by Rescuing Error-Containing Reads

Weixing Feng,Peichao Sang,Deyuan Lian,Yansheng Dong,Fengfei Song,Meng Li,Bo He,Fenglin Cao,Yunlong Liu

DOI: https://doi.org/10.1109/tcbb.2014.2366103

2015-01-01

IEEE/ACM Transactions on Computational Biology and Bioinformatics

Abstract:Next-generation short-read sequencing is widely utilized in genomic studies. Biological applications require an alignment step to map sequencing reads to the reference genome, before acquiring expected genomic information. This requirement makes alignment accuracy a key factor for effective biological interpretation. Normally, when accounting for measurement errors and single nucleotide polymorphisms, short read mappings with a few mismatches are generally considered acceptable. However, to further improve the efficiency of short-read sequencing alignment, we propose a method to retrieve additional reliably aligned reads (reads with more than a pre-defined number of mismatches), using a Bayesian-based approach. In this method, we first retrieve the sequence context around the mismatched nucleotides within the already aligned reads; these loci contain the genomic features where sequencing errors occur. Then, using the derived pattern, we evaluate the remaining (typically discarded) reads with more than the allowed number of mismatches, and calculate a score that represents the probability that a specific alignment is correct. This strategy allows the extraction of more reliably aligned reads, therefore improving alignment sensitivity.Implementation: The source code of our tool, ResSeq, can be downloaded from: https://github.com/hrbeubiocenter/Resseq.
Efficient assembly of nanopore reads via highly accurate and intact error correction

Ying Chen,Fan Nie,Shang-Qian Xie,Ying-Feng Zheng,Qi Dai,Thomas Bray,Yao-Xin Wang,Jian-Feng Xing,Zhi-Jian Huang,De-Peng Wang,Li-Juan He,Feng Luo,Jian-Xin Wang,Yi-Zhi Liu,Chuan-Le Xiao

DOI: https://doi.org/10.1038/s41467-020-20236-7

IF: 16.6

2021-01-04

Nature Communications

Abstract:Abstract Long nanopore reads are advantageous in de novo genome assembly. However, nanopore reads usually have broad error distribution and high-error-rate subsequences. Existing error correction tools cannot correct nanopore reads efficiently and effectively. Most methods trim high-error-rate subsequences during error correction, which reduces both the length of the reads and contiguity of the final assembly. Here, we develop an error correction, and de novo assembly tool designed to overcome complex errors in nanopore reads. We propose an adaptive read selection and two-step progressive method to quickly correct nanopore reads to high accuracy. We introduce a two-stage assembler to utilize the full length of nanopore reads. Our tool achieves superior performance in both error correction and de novo assembling nanopore reads. It requires only 8122 hours to assemble a 35X coverage human genome and achieves a 2.47-fold improvement in NG50. Furthermore, our assembly of the human WERI cell line shows an NG50 of 22 Mbp. The high-quality assembly of nanopore reads can significantly reduce false positives in structure variation detection.

multidisciplinary sciences
SLR-superscaffolder: a de novo scaffolding tool for synthetic long reads using a top-to-bottom scheme

Lidong Guo,Mengyang Xu,Wenchao Wang,Shengqiang Gu,Xia Zhao,Fang Chen,Ou Wang,Xun Xu,Inge Seim,Guangyi Fan,Li Deng,Xin Liu

DOI: https://doi.org/10.1186/s12859-021-04081-z

IF: 3.307

2021-03-25

BMC Bioinformatics

Abstract:Abstract Background Synthetic long reads (SLR) with long-range co-barcoding information are now widely applied in genomics research. Although several tools have been developed for each specific SLR technique, a robust standalone scaffolder with high efficiency is warranted for hybrid genome assembly. Results In this work, we developed a standalone scaffolding tool, SLR-superscaffolder, to link together contigs in draft assemblies using co-barcoding and paired-end read information. Our top-to-bottom scheme first builds a global scaffold graph based on Jaccard Similarity to determine the order and orientation of contigs, and then locally improves the scaffolds with the aid of paired-end information. We also exploited a screening algorithm to reduce the negative effect of misassembled contigs in the input assembly. We applied SLR-superscaffolder to a human single tube long fragment read sequencing dataset and increased the scaffold NG50 of its corresponding draft assembly 1349 fold. Moreover, benchmarking on different input contigs showed that this approach overall outperformed existing SLR scaffolders, providing longer contiguity and fewer misassemblies, especially for short contigs assembled by next-generation sequencing data. The open-source code of SLR-superscaffolder is available at https://github.com/BGI-Qingdao/SLR-superscaffolder . Conclusions SLR-superscaffolder can dramatically improve the contiguity of a draft assembly by integrating a hybrid assembly strategy.

biochemical research methods,biotechnology & applied microbiology,mathematical & computational biology
Improving transcriptome assembly through error correction of high-throughput sequence reads

Matthew D. MacManes,Michael B. Eisen

DOI: https://doi.org/10.7717/peerj.113

IF: 3.061

2013-07-23

PeerJ

Abstract:The study of functional genomics, particularly in non-model organisms, has been dramatically improved over the last few years by the use of transcriptomes and RNAseq. While these studies are potentially extremely powerful, a computationally intensive procedure, the de novo construction of a reference transcriptome must be completed as a prerequisite to further analyses. The accurate reference is critically important as all downstream steps, including estimating transcript abundance are critically dependent on the construction of an accurate reference. Though a substantial amount of research has been done on assembly, only recently have the pre-assembly procedures been studied in detail. Specifically, several stand-alone error correction modules have been reported on and, while they have shown to be effective in reducing errors at the level of sequencing reads, how error correction impacts assembly accuracy is largely unknown. Here, we show via use of a simulated and empiric dataset, that applying error correction to sequencing reads has significant positive effects on assembly accuracy, and should be applied to all datasets. A complete collection of commands which will allow for the production of Reptile corrected reads is available at https://github.com/macmanes/error_correction/tree/master/scripts and as File S1.

multidisciplinary sciences
Publisher Correction: Integrated Analysis of a Compendium of RNA-Seq Datasets for Splicing Factors

Peng Yu,Jin Li,Su-Ping Deng,Feiran Zhang,Petar N. Grozdanov,Eunice W. M. Chin,Sheree D. Martin,Laurent Vergnes,M. Saharul Islam,Deqiang Sun,Janine M. LaSalle,Sean L. McGee,Eyleen Goh,Clinton C. MacDonald,Peng Jin

DOI: https://doi.org/10.1038/s41597-020-00607-x

2020-01-01

Scientific Data

Abstract:An amendment to this paper has been published and can be accessed via a link at the top of the paper.
Correction: De Novo Sequencing and Assembly Analysis of the Pseudostellaria Heterophylla Transcriptome

Jun Li,Wei Zheng,Long Dengkai,Ling Ding,Anhui Gong,Chao Xiao,Weike Jiang,Xiaoqing Liu,Tao Zhou,Luqi Huang

DOI: https://doi.org/10.1371/journal.pone.0170134

IF: 3.7

2017-01-01

PLoS ONE

Abstract:[This corrects the article DOI: 10.1371/journal.pone.0164235.].
Efficient Approach to Correct Read Alignment for Pseudogene Abundance Estimates.

Chelsea J. -T. Ju,Zhuangtian Zhao,Wei Wang

DOI: https://doi.org/10.1109/tcbb.2016.2591533

2017-01-01

IEEE/ACM Transactions on Computational Biology and Bioinformatics

Abstract:RNA-Sequencing has been the leading technology to quantify expression of thousands of genes simultaneously. The data analysis of an RNA-Seq experiment starts from aligning short reads to the reference genome/transcriptome or reconstructed transcriptome. However, current aligners lack the sensitivity to distinguish reads that come from homologous regions of an genome. One group of these homologies is the paralog pseudogenes. Pseudogenes arise from duplication of a set of protein coding genes, and have been considered as degraded paralogs in the genome due to their lost of functionality. Recent studies have provided evidence to support their novel regulatory roles in biological processes. With the growing interests in quantifying the expression level of pseudogenes at different tissues or cell lines, it is critical to have a sensitive method that can correctly align ambiguous reads and accurately estimate the expression level among homologous genes. Previously in PseudoLasso, we proposed a linear regression approach to learn read alignment behaviors, and to leverage this knowledge for abundance estimation and alignment correction. In this paper, we extend the work of PseudoLasso by grouping the homologous genomic regions into different communities using a community detection algorithm, followed by building a linear regression model separately for each community. The results show that this approach is able to retain the same accuracy as PseudoLasso. By breaking the genome into smaller homologous communities, the running time is improved from quadratic growth to linear with respect to the number of genes.
PseudoLasso: leveraging read alignment in homologous regions to correct pseudogene expression estimates via RNASeq.

Chelsea J.-T. Ju,Zhuangtian Zhao,Wei Wang

DOI: https://doi.org/10.1145/2649387.2649447

2014-01-01

Abstract:ABSTRACTPseudogenes have long been considered to be nonfunctional segments in the genome, but recent studies have provided evidence to support their novel regulatory roles in biological processes. With the growing interests in pseudogene research, scientists rely on RNA sequencing technology to estimate expression level of pseudogenes at different tissues or cell lines. The major challenge of RNASeq on pseudogene quantification falls in the high sequence similarity between pseudogenes and their homologous parents. Reads can be ambiguously aligned to multiple homologous regions. In this article, we present PseudoLasso, a genome-wide approach to accurately estimate the abundance of pseudogenes and their parents, and correctly align reads to their origins. Our approach focuses on learning read alignment behaviors, and leveraging this knowledge for abundance estimation and alignment correction. Compared to the read count estimates reported by TopHat2, PseudoLasso is able to provide estimates with a reduced error rate of 10-fold.
Author Correction: Determination of isoform-specific RNA structure with nanopore long reads

Jong Ghut Ashley Aw,Shaun W. Lim,Jia Xu Wang,Finnlay R. P. Lambert,Wen Ting Tan,Yang Shen,Yu Zhang,Pornchai Kaewsapsak,Chenhao Li,Sarah B. Ng,Leah A. Vardy,Meng How Tan,Niranjan Nagarajan,Yue Wan

DOI: https://doi.org/10.1038/s41587-020-00755-w

IF: 46.9

2020-11-12

Nature Biotechnology

Abstract:An amendment to this paper has been published and can be accessed via a link at the top of the paper.

biotechnology & applied microbiology
STAR: ultrafast universal RNA-seq aligner

Alexander Dobin,Carrie A. Davis,Felix Schlesinger,Jorg Drenkow,Chris Zaleski,Sonali Jha,Philippe Batut,Mark Chaisson,Thomas R. Gingeras

DOI: https://doi.org/10.1093/bioinformatics/bts635

IF: 5.8

2012-10-25

Bioinformatics

Abstract:MOTIVATION: Accurate alignment of high-throughput RNA-seq data is a challenging and yet unsolved problem because of the non-contiguous transcript structure, relatively short read lengths and constantly increasing throughput of the sequencing technologies. Currently available RNA-seq aligners suffer from high mapping error rates, low mapping speed, read length limitation and mapping biases.RESULTS: To align our large (>80 billon reads) ENCODE Transcriptome RNA-seq dataset, we developed the Spliced Transcripts Alignment to a Reference (STAR) software based on a previously undescribed RNA-seq alignment algorithm that uses sequential maximum mappable seed search in uncompressed suffix arrays followed by seed clustering and stitching procedure. STAR outperforms other aligners by a factor of >50 in mapping speed, aligning to the human genome 550 million 2 × 76 bp paired-end reads per hour on a modest 12-core server, while at the same time improving alignment sensitivity and precision. In addition to unbiased de novo detection of canonical junctions, STAR can discover non-canonical splices and chimeric (fusion) transcripts, and is also capable of mapping full-length RNA sequences. Using Roche 454 sequencing of reverse transcription polymerase chain reaction amplicons, we experimentally validated 1960 novel intergenic splice junctions with an 80-90% success rate, corroborating the high precision of the STAR mapping strategy.AVAILABILITY AND IMPLEMENTATION: STAR is implemented as a standalone C++ code. STAR is free open source software distributed under GPLv3 license and can be downloaded from http://code.google.com/p/rna-star/.

biochemical research methods,biotechnology & applied microbiology,mathematical & computational biology
Comparative analysis of RNA-Seq alignment algorithms and the RNA-Seq unified mapper (RUM)

Gregory R. Grant,Michael H. Farkas,Angel D. Pizarro,Nicholas F. Lahens,Jonathan Schug,Brian P. Brunk,Christian J. Stoeckert,John B. Hogenesch,Eric A. Pierce

DOI: https://doi.org/10.1093/bioinformatics/btr427

IF: 5.8

2011-07-19

Bioinformatics

Abstract:MOTIVATION: A critical task in high-throughput sequencing is aligning millions of short reads to a reference genome. Alignment is especially complicated for RNA sequencing (RNA-Seq) because of RNA splicing. A number of RNA-Seq algorithms are available, and claim to align reads with high accuracy and efficiency while detecting splice junctions. RNA-Seq data are discrete in nature; therefore, with reasonable gene models and comparative metrics RNA-Seq data can be simulated to sufficient accuracy to enable meaningful benchmarking of alignment algorithms. The exercise to rigorously compare all viable published RNA-Seq algorithms has not been performed previously.RESULTS: We developed an RNA-Seq simulator that models the main impediments to RNA alignment, including alternative splicing, insertions, deletions, substitutions, sequencing errors and intron signal. We used this simulator to measure the accuracy and robustness of available algorithms at the base and junction levels. Additionally, we used reverse transcription-polymerase chain reaction (RT-PCR) and Sanger sequencing to validate the ability of the algorithms to detect novel transcript features such as novel exons and alternative splicing in RNA-Seq data from mouse retina. A pipeline based on BLAT was developed to explore the performance of established tools for this problem, and to compare it to the recently developed methods. This pipeline, the RNA-Seq Unified Mapper (RUM), performs comparably to the best current aligners and provides an advantageous combination of accuracy, speed and usability.AVAILABILITY: The RUM pipeline is distributed via the Amazon Cloud and for computing clusters using the Sun Grid Engine (http://cbil.upenn.edu/RUM).CONTACT: ggrant@pcbi.upenn.edu; epierce@mail.med.upenn.eduSUPPLEMENTARY INFORMATION: The RNA-Seq sequence reads described in the article are deposited at GEO, accession GSE26248.

biochemical research methods,biotechnology & applied microbiology,mathematical & computational biology

Correction To: P_RNA_scaffolder: a Fast and Accurate Genome Scaffolder Using Paired-End RNAsequencing Reads

P_RNA_scaffolder: a Fast and Accurate Genome Scaffolder Using Paired-End RNA-sequencing Reads

Correction To: L_RNA_scaffolder: Scaffolding Genomes with Transcripts

Pseudo-Sanger Sequencing: Massively Parallel Production of Long and Near Error-Free Reads Using NGS Technology

L_RNA_scaffolder: Scaffolding Genomes with Transcripts (vol 14, 604, 2013)

L_RNA_scaffolder: scaffolding genomes with transcripts

Corrigendum: Long-range RNA Pairings Contribute to Mutually Exclusive Splicing

GeneScissors: a comprehensive approach to detecting and correcting spurious transcriptome inference owing to RNA-seq reads misalignment.

PEP_scaffolder: Using (homologous) Proteins to Scaffold Genomes

PEP_scaffolder: using (homologous) proteins to scaffold genomes.

ResSeq: Enhancing Short-Read Sequencing Alignment by Rescuing Error-Containing Reads

Efficient assembly of nanopore reads via highly accurate and intact error correction

SLR-superscaffolder: a de novo scaffolding tool for synthetic long reads using a top-to-bottom scheme

Improving transcriptome assembly through error correction of high-throughput sequence reads

Publisher Correction: Integrated Analysis of a Compendium of RNA-Seq Datasets for Splicing Factors

Correction: De Novo Sequencing and Assembly Analysis of the Pseudostellaria Heterophylla Transcriptome

Efficient Approach to Correct Read Alignment for Pseudogene Abundance Estimates.

PseudoLasso: leveraging read alignment in homologous regions to correct pseudogene expression estimates via RNASeq.

Author Correction: Determination of isoform-specific RNA structure with nanopore long reads

STAR: ultrafast universal RNA-seq aligner

Comparative analysis of RNA-Seq alignment algorithms and the RNA-Seq unified mapper (RUM)