How to annotate edgeR results?
1
0
Entering edit mode
4.3 years ago

I have got log2fold change by using edgeR.

geneID        logFC         logCPM       PValue
AG1IA_00007   0.457824277   6.711014893  8.44E-53
AG1IA_00008   1.514709777   4.376560603  3.28E-111
AG1IA_00009   -0.409884055  5.357685435  1.39E-19
AG1IA_00010   0.36806713    4.554970223  5.85E-10
AG1IA_00011   -0.338575977  1.144437033  0.095242116

I needed to functional annotate this table. This belongs to fungi (Rhizoctonia solani), I also have gtf and gff files.. can anyone help me by giving commands to annotate this.

R rna-seq next-gen annotation • 1.4k views
ADD COMMENT
1
Entering edit mode

What do you mean by "functionally annotate"? What information are you trying to add? Can you give an example of your desired output?

ADD REPLY
0
Entering edit mode

geneID logFC logCPM PValue product

AG1IA_00007 0.457824277 6.711014893 8.44E-53 vesicular-fusion protein SEC18

AG1IA_00008 1.514709777 4.376560603 3.28E-111 hypothetical protein AG1IA_00008

AG1IA_00009 -0.409884055 5.357685435 1.39E-19 ran binding protein 11

AG1IA_00010 0.36806713 4.554970223 5.85E-10 peroxisomal biogenesis factor 2

AG1IA_00011 -0.338575977 1.144437033 0.095242116 hypothetical protein AG1IA_00011

I needed to functional annotate this last row(product).

I have got these names by individually copy-pasting the gene names in NCBI. There are 10538 genes, I can't do this manually ... so I need some help Jared.andrews.

ADD REPLY
1
Entering edit mode

There is a genome sequence and gene annotation for that fungus available here. Maybe redo your analysis to use those data. Alternatively, you can do ORF prediction/translation on your transcripts and annotate using HMMer, blastp, or interproscan.

ADD REPLY
1
Entering edit mode
4.3 years ago
h.mon 35k

The annotation is available at the Rice Sheath Blight Pathogen Genomics Database download page (here). You then have to format the gene set of interest into some format the functional enrichment software you want to use - for example, here is a list of formats Broad GSEA accepts: Data formats.

ADD COMMENT

Login before adding your answer.

Traffic: 2822 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