Question: nearby genes expression
0
gravatar for ta_awwad
2.5 years ago by
ta_awwad280
Frankfurt am Main
ta_awwad280 wrote:

Hello Guy, I have deseq2 table with differentially expressed genes. I would like to extract the log2FC as well as raw counts for 5 genes That are flanking my gene of interest to use them for heatmap .. can anyone help?

rna-seq deseq2 • 643 views
ADD COMMENTlink modified 2.5 years ago • written 2.5 years ago by ta_awwad280
2

Sure.

  1. Get the coordinates of your gene of interest from GTF.
  2. use bedtools to find what are the nearest genes using closest feature
  3. Take the gene IDs from step.2 and grep for log2FC from DESeq2 results

or

  1. Go to UCSC genome browser and enter your gene of interest
  2. Zoom out to find what are the flanking genes and note down
  3. Take the gene IDs from step.2 and grep for log2FC from DESeq2 results
ADD REPLYlink written 2.5 years ago by venu6.6k

thanks much venu ... but there is no shortcut to do this task?

ADD REPLYlink written 2.5 years ago by ta_awwad280
1

Both of them takes <5 min to complete the task. I can't think of any other shortcut.

ADD REPLYlink written 2.5 years ago by venu6.6k

If you know the names of those genes then fish them out of the results object (which can be written to a file).

ADD REPLYlink written 2.5 years ago by genomax85k

Thanks Genomax .. I was just thinking of direct way to fetch the names of those genes then extract the expression values .. I will do it manually ... thanks once more

ADD REPLYlink written 2.5 years ago by ta_awwad280
0
gravatar for arta
2.5 years ago by
arta580
Sweden
arta580 wrote:
dds <- DESeq(dds)
res <- results(dds)

res is the data frame where you can find log2FC in second column res[,2] and rownames are genes.

dds <- DESeqDataSetFromMatrix(countData = cts,
                          colData = coldata,
                          design = ~ condition)

The counts are stored in countData.

ADD COMMENTlink modified 2.5 years ago • written 2.5 years ago by arta580

thanks arta for the contribution ... but this not what I want actually

ADD REPLYlink written 2.5 years ago by ta_awwad280
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: 1969 users visited in the last hour