Using DESEQ2 with selective gene vs all genes
1
0
Entering edit mode
Shail • 0
@dba6855a
Last seen 6 weeks ago
Hong Kong

Hi,

I am new to Deseq analysis. Recently i conducted metatranscriptomics analysis of bacterial community in two different conditions. My aim to explore the overall significent transcript expression and also some selective transcripts (related to specific function). I used raw featurecounts output to run Deseq and created heatmap on VSD transformed data on to 20 features.

But, due to vast number of transcripts, this does not cover my transcripts of interest. Is it ok to extract raw counts of all transcripts related to a particular function and run Deseq?

Also, the heatmap generated via vsd transformed data shows scale as -1, -0.5, 0.5 and 1.0. Are this suppossed to be the Log2F change?

DESeq2 • 150 views
ADD COMMENT
0
Entering edit mode
@mikelove
Last seen 4 hours ago
United States

This has been asked before on the support site (even within the past month) and we do not recommend running DESeq2 on a subset of transcripts or genes.

For details on the transformations, see the DESeq2 workflow.

https://www.bioconductor.org/packages/devel/workflows/vignettes/rnaseqGene/inst/doc/rnaseqGene.html#the-variance-stabilizing-transformation-and-the-rlog

For genes with high counts, both the VST and the rlog will give similar result to the ordinary log2 transformation of normalized counts

ADD COMMENT
0
Entering edit mode

Thanks for the quck reply.

ADD REPLY

Login before adding your answer.

Traffic: 310 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6