Up and down pathways are the same?
1
1
Entering edit mode
6.4 years ago
Tania ▴ 180

Hello everyone

I am using gage to get the pathways of some affymatrix sets.

The set of up and down pathways are the same but with different order in one one pathway only when I set same.dr = true.

When I set same.dr = false I got the same set of pathways. This does not change whether or not I set paired or unpaired. Does this mean the pathways are both up and down? And what does this really mean? Or something is wrong?

Here is a snapshot from my code:

res$symbol <- mapIds(org.Hs.eg.db,
                     keys=row.names(res),
                     column="SYMBOL",
                     keytype="SYMBOL",
                     multiVals="first")
res$entrez <- mapIds(org.Hs.eg.db,
                     keys=row.names(res),
                     column="ENTREZID",
                     keytype="SYMBOL",
                     multiVals="first")
foldchanges = res$logfc
names(foldchanges) = res$entrez
#---------------Use Kegg and gage to get upregulated and downregulated pathways
datakegg.gs)
keggres = gage(foldchanges, gsets =kegg.gs, same.dir =,True,  compare="unpaired" )

lapply(keggres, head)
keggrespathwaysup = data.frame(id=rownames(keggres$greater), keggres$greater) %>%
  tbl_df() %>%
  filter(row_number()<=5) %>%
  .$id %>%
  as.character()
keggrespathwaysdn = data.frame(id=rownames(keggres$less), keggres$less) %>%
  tbl_df() %>%
  filter(row_number()<=5) %>%
  .$id %>%
  as.character()
write.csv(keggrespathwaysup, "keggspathsup.csv")
write.csv(keggrespathwaysdn, "keggspathsdn.csv")
#-------------------------------------------------------
keggresidsup = substr(keggrespathwaysup, start=1, stop=8)
keggresidsup
keggresidsdn = substr(keggrespathwaysdn, start=1, stop=8)
gobpres = gage(foldchanges, gsets=kegg.gs, same.dir =TRUE, compare ="unpaired")
lapply(gobpres, head)
plot_pathway = function(pid) pathview(gene.data=foldchanges,gene.idtype="ENTREZID", pathway.id=pid, 
species="hsa", new.signature=FALSE)

#------------------------------plot multiple pathways ups and downs
 tmpup = sapply(keggresidsup, function(pid) pathview(gene.data=foldchanges,gene.idtype="ENTREZID", pathway.id=pid, species="hsa"))
tmpdn = sapply(keggresidsdn, function(pid) pathview(gene.data=foldchanges,gene.idtype="ENTREZID", pathway.id=pid, species="hsa"))

Thanks

pathways gage • 2.3k views
ADD COMMENT
3
Entering edit mode
6.4 years ago

Why should it be surprising ? If a pathway is strongly affected one could expect that positive regulators will go the opposite way of negative regulators. Check the pathway genes in both situations. If they are the same then you may have a problem, otherwise think whether this makes sense in the biological context of your experiment.

ADD COMMENT
0
Entering edit mode

Thank you.

"If a pathway is strongly affected one could expect that positive regulators will go the opposite way of negative regulators"

Does this means for two sets of genes, one set up-regulate a pathway and the other set down-regulate the same pathway? I feel there is something I am still missing

ADD REPLY
3
Entering edit mode

I agree entirely with Jean-Karim. This is actually a common occurrence in in silico analyses. It does not necessarily mean that the entire pathway is being both up- and down-regulated. It could indicate that certain parts of the pathway are adversely and differently affected. Biology is pretty complex and, right now, even with all information that we have produced, cellular pathways are still very much unknown entities to us. Many different pathways are undoubtedly also inter-connected.

ADD REPLY
0
Entering edit mode

Thanks Jean-Karim and Kevin. I think it is much clear now !

ADD REPLY
0
Entering edit mode

Did the data validate? Were those gene entities non-overlapping (to a large extent)?

ADD REPLY
0
Entering edit mode

I am trying to print the genes but can't

ADD REPLY
2
Entering edit mode

A pathway often contains activators and inhibitors. Imagine that a cell treatment leads to activation of a signaling pathway. To activate the pathway, the cell may need to express some members of this pathway (e.g. positive regulators) while at the same time removing inhibitors (e.g. negative regulators part of some feedback loop) of this pathway. For the purpose of enrichment analysis, pathways are simply lists of genes participating in the same biological process. The way biological processes and pathways are defined is arbitrary and different resources/databases will have different definitions (for example, compare the DNA replication pathway in Panther and in KEGG).

ADD REPLY
0
Entering edit mode

Thanks Jean-Karim, very much appreciated !

ADD REPLY

Login before adding your answer.

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