Question: Use of slice to estimate rpkm/fpkm
gravatar for a.morris
10 months ago by
a.morris40 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 10 months ago by Jennifer Hillman Jackson25k • written 10 months ago by a.morris40
gravatar for Jennifer Hillman Jackson
10 months ago by
United States
Jennifer Hillman Jackson25k 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 10 months ago by Jennifer Hillman Jackson25k
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: 172 users visited in the last hour