Abstract
After mapping, RNA-Seq data can be summarized by a sequence of read counts commonly modeled as Poisson variables with constant rates along each transcript, which actually fit data poorly. We suggest using variable rates for different positions, and propose two models to predict these rates based on local sequences. These models explain more than 50% of the variations and can lead to improved estimates of gene and isoform expressions for both Illumina and Applied Biosystems data.
Publication types
-
Research Support, N.I.H., Extramural
-
Research Support, U.S. Gov't, Non-P.H.S.
MeSH terms
-
Animals
-
Apolipoproteins E / genetics
-
Base Sequence
-
Databases, Nucleic Acid*
-
Embryo, Mammalian / metabolism
-
Exons / genetics
-
Gene Expression Profiling
-
Gene Expression Regulation
-
Humans
-
Linear Models
-
Mice
-
Models, Genetic*
-
Poisson Distribution
-
Protein Isoforms / genetics
-
Protein Isoforms / metabolism
-
RNA / genetics*
-
Sequence Analysis, RNA / methods*
-
Statistics, Nonparametric
Substances
-
Apolipoproteins E
-
Protein Isoforms
-
RNA