I am having trouble grabbing specific pathway info using databases in R. I have RNAseq results and I want to remove immune related genes from the current list that I have.

With a vector of gene names/symbols from RNAseq , I queried org.hs.eg.db to grab entrezIDs keys etc.

mapIds(org.Hs.eg.db,keys=mygenes,keytype ="ENTREZID",column="PATHID")

I thought about trying to use pathids from this database which I think are KEGG pathway ids which I can subsequently filter out the appropriate ones but my list goes from 19,000 to 4,000 before I drop immune related genes due to a lot of NAs. I tried reactomedb, but the pathnames don't really let me filter by immune related or not.

Any advice is appreciated.

mapIds Bioconductor https://bioconductor.org/packages/devel/bioc/vignettes/AnnotationDbi/inst/doc/IntroToAnnotationPackages.pdf

KEGG https://www.genome.jp/kegg/ "KEGG is a database resource for understanding high-level functions and utilities of the biological system, such as the cell, the organism and the ecosystem, from molecular-level information, especially large-scale molecular datasets generated by genome sequencing and other high-throughput experimental technologies. See Release notes (April 1, 2022) for new and updated features."

  • $\begingroup$ Hi James, those pathID come from the KEGG database and aren't updated since 2016 due to a license change. There isn't an easy solution to this problem (what is the threshold for a gene to be considered immune related?) and most people either use REACTOME, wikipathways or MSigDB (from the broad institute). $\endgroup$
    – llrs
    Jun 15, 2020 at 9:53
  • $\begingroup$ I was a bit confused on how to query REACTOME, as it was hard to visually see which terms were immune related when I looked at the pathname keys. Would a threshold use some sort of ranked pvalue (in which case probably something low)? Thanks for the suggestions on the others! $\endgroup$
    – James
    Jun 15, 2020 at 12:11
  • $\begingroup$ Well the problem is that the "olfactory pathway" on your epithelial cells might indicate an immune response and there isn't a test to prove that. People use some kind of network/system-biology approach or the databases available. Afaik there isn't much progress on this field. I'm attempting to tackle this problem by creating new methods about pathways but my thesis is in other topics so I can't invest much time on this... $\endgroup$
    – llrs
    Jun 15, 2020 at 13:17


Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Browse other questions tagged or ask your own question.