Functional enrichment and visualization of metatranscriptomics data

Question

[shail.nair05]

Hi, I am new to clusterProfiler. Is it possible to use clusterProfiler for metatranscriptomics DEG data (deseq2)? As i read from http://yulab-smu.top/clusterProfiler-book/chapter3.html#input-data for unsupported organisms, one need to create TERM2GENE and TERM2NAME file, which according to the guide requires a organism specfic gmt file.

My sample is of mixed bacterial transcript expression (metatranscriptomics). I ran deseq2 and wanted to do functional gene enrichment analysis.

[andres.firrincieli]

TERM2NAME is optional. TERM2GENE is a data.frame with first column of term ID (wpid)and second column of corresponding mapped gene:

  wpid  gene
1 WP23  4690
2 WP23  5781
3 WP23 11184
4 WP23  6195
5 WP23  4609
6 WP23   998

In the example above the term ID correspond to the wikipathway IDs but can be any annotation term (e.g. Gene Ontology ID). For metagenomic/metatranscriptomic data you can use eggNOG to annotate your transcripts with GO

[shail.nair05]

Thanks. is it possible to use kegg for metagenomic/metatranscriptomic?

[andres.firrincieli]

Yes, you can use KEGG pathways id instead of GO. eggNOG will also annotate your transcripts with the KEGG annotation terms.

[shail.nair05]

will try with kegg. Thanks :)