Get the located metabolic subpathways

Share:

Description

Locate metabolic subpathways mediated by miRNAs.

Usage

1
2
getLocSubGraph(moleculeList, graphList, type="gene_miRNA", 
   n=1, s=10, method = "shortestPaths")

Arguments

moleculeList

A character vector. Such as differentially expressed miRNAs and/or genes under disease phenotypes.

graphList

A graph list. There nodes must be represented by genes.

type

A character string. Should be one of "gene", "miRNA" or "gene_miRNA".

n

An integer. The maximum acceptable quantities of non-signature node at the shortest path between each two differential molecules.

s

An integer. The minimum acceptable quantities of nodes in located subpathways.

method

A character string. In which the shorest path algorithms will be used. See the function get.shortest.paths.

Details

We apply lenient distance similarity method to locate metabolic subpathways mediated by miRNAs. We first map user interested differentially expressed miRNAs and/or genes to pathways as signatures. For a given pathway, we compute the shortest path between any two signatures. In shortest path, if the number of non-signature nodes between two signatures is no more than n, then these two signature nodes and other nodes at the shortest path are added into the same node set. We extract the corresponding subgraph in the pathway graph according to each node set. We finally define these subgraphs with node number >= s as the subpathway regions of the pathway. The argument n is maximum number of permitted non-signature nodes at the shortest path between signature nodes. The default parameter n=1. The argument s is used to filter subpathways in which the number of nodes are less than the parameter s. The default parameter s=10. The argument method determines which shortest path algorithms will be used. We set the default value as "get.shortest.paths".

Value

A list of graphs.

Author(s)

Li Feng, Chunquan Li and Xia Li

See Also

identifyGraph, get.shortest.paths

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
## Not run: 

### Integrate miRNAs to KEGG pathway graphs ###

## get hsa-specificd miRNA-target interactions ##
 expMir2Tar <- GetK2riData(K2riData="expMir2Tar")
 row1 <- which(expMir2Tar[["LowTHExps"]]=="YES")
 row2 <- which(expMir2Tar[["Species"]]=="hsa")
 relations <- unique(expMir2Tar[intersect(row1,row2),c(2:3)])

# get user-interested miRNAs and genes sets.
 moleculeList <- c(getBackground(type="gene")[1:1000],
                 getBackground(type="miRNA")[1:2000])

## get direct KEGG metabolic pathway graphs ##
 graphList <- GetK2riData(K2riData="MetabolicGEGEEMGraph")
# get reconstructed pathway graph list.
 InteGraphList <- getInteGraphList(graphList, relations) 
# get locate subpathways.
 subGraphList <- getLocSubGraph(moleculeList,InteGraphList,
                 type="gene_miRNA",n=1,s=10)
# visualize the located subpathways.
 plotGraph(subGraphList[[1]],layout=layout.random)


## get undirect KEGG metabolic pathway graphs ##
 graphList <- GetK2riData(K2riData="MetabolicGEGEUEMGraph")
# get reconstructed pathway graph list.
 InteGraphList <- getInteGraphList(graphList, relations) 
# get locate subpathways.
 subGraphList <- getLocSubGraph(moleculeList,InteGraphList,
                 type="gene_miRNA",n=1,s=10)
# visualize the located subpathways.
 plotGraph(subGraphList[[1]],layout=layout.random)

## End(Not run)