Seqanswers Leaderboard Ad

**fanli** · 10-17-2016, 08:20 AM

I don't see why that'd be a problem...

Quoting from this benchmark paper:

A comprehensive evaluation of ensembl, RefSeq, and UCSC annotations in the context of RNA-seq read mapping and gene quantification - BMC Genomics

http://bmcgenomics.biomedcentral.com/articles/10.1186/s12864-015-1308-8

Background RNA-Seq has become increasingly popular in transcriptome profiling. One aspect of transcriptome research is to quantify the expression levels of genomic elements, such as genes, their transcripts and exons. Acquiring a transcriptome expression profile requires genomic elements to be defined in the context of the genome. Multiple human genome annotation databases exist, including RefGene (RefSeq Gene), Ensembl, and the UCSC annotation database. The impact of the choice of an annotation on estimating gene expression remains insufficiently investigated. Results In this paper, we systematically characterized the impact of genome annotation choice on read mapping and transcriptome quantification by analyzing a RNA-Seq dataset generated by the Human Body Map 2.0 Project. The impact of a gene model on mapping of non-junction reads is different from junction reads. For the RNA-Seq dataset with a read length of 75 bp, on average, 95% of non-junction reads were mapped to exactly the same genomic location regardless of which gene models was used. By contrast, this percentage dropped to 53% for junction reads. In addition, about 30% of junction reads failed to align without the assistance of a gene model, while 10–15% mapped alternatively. There are 21,958 common genes among RefGene, Ensembl, and UCSC annotations. When we compared the gene quantification results in RefGene and Ensembl annotations, 20% of genes are not expressed, and thus have a zero count in both annotations. Surprisingly, identical gene quantification results were obtained for only 16.3% (about one sixth) of genes. Approximately 28.1% of genes’ expression levels differed by 5% or higher, and of those, the relative expression levels for 9.3% of genes (equivalent to 2038) differed by 50% or greater. The case studies revealed that the gene definition differences in gene models frequently result in inconsistency in gene quantification. Conclusions We demonstrated that the choice of a gene model has a dramatic effect on both gene quantification and differential analysis. Our research will help RNA-Seq data analysts to make an informed choice of gene model in practical RNA-Seq data analysis.

...when conducting research that emphasizes reproducible and robust gene expression estimates, a less complex genome annotation, such as RefGene, might be preferred. When conducting more exploratory research, a more complex genome annotation, such as Ensembl, should be chosen.

**graceqy** · 10-19-2016, 04:45 AM

Got it. Thanks!

Topics	Statistics	Last Post
Artificial Intelligence Shows Promise in Automating Functional Genomics Analysis by seqadmin Started by seqadmin, 12-02-2024, 09:29 AM	0 responses 131 views 0 likes	Last Post by seqadmin 12-02-2024, 09:29 AM
Rethinking Epigenetic Clocks for Accurate Aging Measurements by seqadmin Started by seqadmin, 12-02-2024, 09:06 AM	0 responses 48 views 0 likes	Last Post by seqadmin 12-02-2024, 09:06 AM
Researchers Capture Ribosome-mRNA Interactions in Atomic Detail by seqadmin Started by seqadmin, 12-02-2024, 08:03 AM	0 responses 38 views 0 likes	Last Post by seqadmin 12-02-2024, 08:03 AM
Single-Cell Sequencing Links Early Genetic Mutations to Breast Cancer Development by seqadmin Started by seqadmin, 11-22-2024, 07:36 AM	0 responses 68 views 0 likes	Last Post by seqadmin 11-22-2024, 07:36 AM

Seqanswers Leaderboard Ad

Announcement

Kallisto with refseq fasta file

Comment

Comment

Latest Articles

ad_right_rmr

News