Question: Snords and MiRs in RNA-Seq data
gravatar for lhawthorn
4.6 years ago by
lhawthorn0 wrote:

I have RNA-Seq data from a large number of samples analyzed using CuffDiff. They are 50bp pe reads and libraries prepared using Truseq (so size selected). The most highly up and down regulated transcripts are snords and miRNAs. They have unusually high fold changes but the q-values are high (~1). But in many cases the q values are high because the mir/snord is absent in one of the samples. Are these artifacts of alignment and would you suggest ignoring them? Or should they be included in my analyses? Thanks

rna-seq alignment • 1.0k views
ADD COMMENTlink modified 4.5 years ago by Biostar ♦♦ 20 • written 4.6 years ago by lhawthorn0

I'd recommend including them. You might be able to use the GTF file from cufflinks (or stringTie or whatever), use that with featureCounts to get counts and then put them in edgeR. The results from that tend to be a little more reliable (you could also try DESeq2, but it seems edgeR performs a bit better when you randomly have a sample with a dropout). Note that it's good to confirm results with an orthogonal technology regardless of what's used.

ADD REPLYlink written 4.5 years ago by Devon Ryan98k
Please log in to add an answer.


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