Question: Pathway Graphs for gene set from RNA-seq
gravatar for a.james
4.6 years ago by
a.james180 wrote:

Hello All,

I would like to display the data from results in a graphical way rfom the functional pathway analysis.

I have tried several R packages like KEGGgraph , PAthview, But they have there own input data format.

My data looks like this, where first col is genes and second column is the pathway enriched .


ALDH7A1; BPGM; hsa00010: Glycolysis / Gluconeogenesis
ALDH7A1; hsa00053: Ascorbate and aldarate metabolism
ACSL6; ALDH7A1; hsa00071: Fatty acid metabolism
SOAT2; hsa00100: Steroid biosynthesis
PDE7B; GMPR; PDE4B; hsa00230: Purine metabolism
UPP1; hsa00240: Pyrimidine metabolism
ALAS2; ALDH7A1; BPGM; hsa00260: Glycine
TST; hsa00270: Cysteine and methionine metabolism
ALDH7A1; BCAT1; hsa00280: Valine
BCAT1; hsa00290: Valine
ALDH7A1; hsa00300: Lysine biosynthesis
ALDH7A1; hsa00310: Lysine degradation

So in there are many pathways were more than one gene is involved .. For this kind of data is there any available method for data display I want it in pathway graphs not as bar plots.

Thank you so much for help


rna-seq • 2.9k views
ADD COMMENTlink modified 4.6 years ago • written 4.6 years ago by a.james180

>> But they have their own input data format.

A big part of the job is restructuring or re-organizing data. That's unrelated to graph displays.

ADD REPLYlink written 4.6 years ago by karl.stamm3.4k

Thank you ppl. I will look into all the solutions

ADD REPLYlink written 4.6 years ago by a.james180

By the way, what tools did you use to do pathway enrichment analysis? Have you done GO enrichment?


ADD REPLYlink written 4.6 years ago by liupfskygre190

You may try GAGE for both pathway and GO analysis if you know R. It supports ~3000 KEGG species, and 19 GO species.:

It works with well with Pathview.

ADD REPLYlink written 4.6 years ago by bigmawen310

OK, thanks!

ADD REPLYlink written 4.5 years ago by liupfskygre190

If you want to find GO or Pathway enrichment for your set of genes you can try this. 


Gene Set Clustering based on Functional annotation (GeneSCF)

ADD REPLYlink modified 2.2 years ago • written 4.5 years ago by EagleEye6.2k
gravatar for bigmawen
4.6 years ago by
United States
bigmawen310 wrote:

Pathview is actually quite flexible. It accepts different data formats (matrices, vectors), different data attributes (discrete, continuous, numeric, characters). It works with tens of different gene/compound IDs and ~3000 species. Please check the package tutorial and documentation:


If you supply your data as it is, propably no program can figure out what that is. Most programs want you do basic data preparation and provide some information about the data.

ADD COMMENTlink written 4.6 years ago by bigmawen310
gravatar for karl.stamm
4.6 years ago by
United States
karl.stamm3.4k wrote:

If you want to see nodes and edges on a graphical display, you need a way to define those connections.

Cytoscape is a good tool for viewing and manipulating networks. A plugin for Cytoscape called GenMapp-CS can automatically take networks from and draw your gene sets on them. Its nice if you have a gene list with differential expression, because cytoscape can color each gene correspondingly.

I guess you want to find HSA00010 as a diagram, then highlight ALDH7A1 and BPGM on it. Iteration would make a diagram for every line of the file.

If you don't want to review thousands of pictures; what were you trying to do with this dataset?

ADD COMMENTlink written 4.6 years ago by karl.stamm3.4k
gravatar for EagleEye
4.6 years ago by
EagleEye6.2k wrote:
Hi, try Cytoscape... It has various plugins to play with your data.
ADD COMMENTlink written 4.6 years ago by EagleEye6.2k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1075 users visited in the last hour