Question: StringTie reports 0 coverage for all reference transcripts
0
gravatar for PoGibas
21 months ago by
PoGibas4.7k
Vilnius
PoGibas4.7k wrote:

I have used StringTie on mice RNA-seq data. StringTie was run with parameters given in Nature paper. Reference gtf annotation was downloaded from the Gencode Genes.
Most of the transcripts in the reference annotation have very low expression values (> 95% transcripts have 0 fpkm; for one sample all reference transcripts have 0 fpkm). However, half of newly merged transcripts have non-zero expression values. My questions are:

  1. Is it normal for all the reference transcripts to have expression values equal to zero?
  2. After running StringTie with default parameters I got many (maybe too many) newly merged transcripts (140K vs 25K in the reference annotation). This looks kinda strange and I would prefer to use transcripts from the reference, but as I mentioned above - they have 0 expression values.
rna-seq stringtie • 939 views
ADD COMMENTlink modified 20 months ago by k.mooo10 • written 21 months ago by PoGibas4.7k
1
gravatar for k.mooo
20 months ago by
k.mooo10
k.mooo10 wrote:

0 expression values could indicate an incongruous reference genome between alignment and annotation. The same version, as well as the same source (for example, Ensembl or UCSC), must be used when aligning reads and annotating reads. This is mentioned at the bottom of this blog with respect to reference genome builds. https://github.com/griffithlab/rnaseq_tutorial/wiki/Annotation

ADD COMMENTlink written 20 months ago by k.mooo10
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 644 users visited in the last hour