Question: Use of slice to estimate rpkm/fpkm
gravatar for a.morris
5 days ago by
a.morris10 wrote:

When estimating transcript levels from RNA Seq data that have been normalized using rpkm or fpkm, would the slice SAMtool (which manually sets a list of coordinates and extracts the data say for the few genes selected out of the whole transcriptome) differ from the results obtained from the (unsliced) whole transcriptome based RNA-Seq rpkm/fpkm?

ADD COMMENTlink modified 4 days ago by Jennifer Hillman Jackson23k • written 5 days ago by a.morris10
gravatar for Jennifer Hillman Jackson
4 days ago by
United States
Jennifer Hillman Jackson23k wrote:


Yes. As far as I know, it would be better to use the complete hit distribution for the input BAM/SAM datasets when performing the differential expression testing. After that is done, filter results for just those target genes of interest. Using subsets of data is fine sometimes, but those are special use-cases.

To see if this is true for your data, you could run both data though and then compare to find out exactly how different these results are. The DE calculations for R/FPKM are all relative to the specific tool run/inputs and so cannot be directly compared to different runs/different inputs, but you could compare overall counts for significant over/under expression, the genes and other features associated with those, and possibly a few target gene/transcript regions in a visualization that contains both results (in a genome browser, etc).

If I have misunderstood your question, please explain what you are doing step-by-step to clarify.

Galaxy tutorials:

Thanks! Jen, Galaxy team

ADD COMMENTlink written 4 days ago by Jennifer Hillman Jackson23k
Please log in to add an answer.


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