Adding information to the DESeq2 output file (gene name and expression levels for each sample(

Question: Adding information to the DESeq2 output file (gene name and expression levels for each sample(

19 months ago by

amyfm • 0

Ireland

amyfm • 0 wrote:

Hi,

After doing DESeq2 with RNA seq data, the output I obtain is an excel file with ensembl transcript ID, baseMean, log2FoldChange, lfcSE, stat, pvalue, and padj.

I would like to know 2 things:

-If you know if there is any code that I can use in order to obtain the gene name or description in the output file instead of the ensembl transcript ID

-If there is any code to obtain more columns in the output excel with the expression levels of each of my samples

Thank you

ADD COMMENT • link •

modified 19 months ago by Michael Love ♦ 12k • written 19 months ago by amyfm • 0

19 months ago by

Michael Love ♦ 12k

United States

Michael Love ♦ 12k wrote:

See here:

http://www.bioconductor.org/help/workflows/rnaseqGene/#annotate

ADD COMMENT • link written 19 months ago by Michael Love ♦ 12k

The normalized counts for each sample can be obtained with

counts(dds, normalized=TRUE)

and you can add to the results table this with cbind ("column bind").

ADD REPLY • link written 19 months ago by Michael Love ♦ 12k

thank you very much for your answer

ADD REPLY • link modified 18 months ago • written 19 months ago by amyfm • 0

Hi,

I can't get the information I want (gene names on my excel output) as I get some error messages: Error in .testForValidKeys(x, keys, keytype, fks) :
None of the keys entered are valid keys for 'ENSEMBLTRANS'. Please use the keys method to see a listing of valid arguments.

I am using cows as a model, and my excel output is a columns with the ensembl transcripts IDs (appearing like this for example: transcript:ENSBTAT00000019061

I would appreciate if someone give me the exact code for my case, as I am not very familiar with R and differential expression analysis but I need to have this information.

Thank you

ADD REPLY • link written 18 months ago by amyfm • 0

Annotation tasks like this typically take some amount of exploration/trial and error. If I came upon this error myself, I would compare the keys I was trying to use to query the database with the keys that are in the database, for example:

k <- keys(org.Xx.eg.db, keytype="ENSEMBLTRANS")
head(k)

Where Xx would be your organism. You can take a look at the AnnotationDbi vignette for more guidance.

ADD REPLY • link written 18 months ago by Michael Love ♦ 12k

Similar posts • Search »