Processes KEGG Pathways data extracted from a KGML file. Joins the entries and relations into a single data frame and translates the Gene Symbols to UniProt IDs.
Arguments
- entries
A data frames with entries extracted from a KGML file by
kegg_pathway_download.- relations
A data frames with relations extracted from a KGML file by
kegg_pathway_download.- max_expansion
Numeric: the maximum number of relations derived from a single relation record. As one entry might represent more than one molecular entities, one relation might yield a large number of relations in the processing. This happens in a combinatorial way, e.g. if the two entries represent 3 and 4 entities, that results 12 relations. If
NULL, all relations will be expanded.- simplify
Logical: remove KEGG's internal identifiers and the pathway annotations, keep only unique interactions with direction and effect sign.
Value
A data frame (tibble) of interactions. In rare cases when a
pathway doesn't contain any relation, returns NULL.
Examples
hsa04350 <- kegg_pathway_download('hsa04350', process = FALSE)
tgf_pathway <- kegg_process(hsa04350$entries, hsa04350$relations)
tgf_pathway
#> # A tibble: 55 × 12
#> source target type effect arrow relation_id kegg_id_source genesymbol_source
#> <chr> <chr> <chr> <chr> <chr> <chr> <chr> <chr>
#> 1 51 49 PPrel activ… --> hsa04350:1 hsa:7040 hsa:… TGFB1
#> 2 57 55 PPrel activ… --> hsa04350:2 hsa:151449 hs… GDF7
#> 3 34 32 PPrel activ… --> hsa04350:3 hsa:3624 hsa:… INHBA
#> 4 20 17 PPrel activ… --> hsa04350:4 hsa:4838 NODAL
#> 5 60 46 PPrel activ… --> hsa04350:5 hsa:4086 hsa:… SMAD1
#> 6 43 41 PPrel activ… --> hsa04350:6 hsa:4087 hsa:… SMAD2
#> 7 22 16 PPrel activ… --> hsa04350:7 hsa:4087 hsa:… SMAD2
#> 8 19 15 PPrel activ… --> hsa04350:8 hsa:4087 hsa:… SMAD2
#> 9 27 26 PPrel inhib… --| hsa04350:11 hsa:4609 MYC
#> 10 47 43 PPrel inhib… --| hsa04350:13 hsa:4091 hsa:… SMAD6
#> # ℹ 45 more rows
#> # ℹ 4 more variables: uniprot_source <chr>, kegg_id_target <chr>,
#> # genesymbol_target <chr>, uniprot_target <chr>
# # A tibble: 50 x 12
# source target type effect arrow relation_id kegg_id_source
# <chr> <chr> <chr> <chr> <chr> <chr> <chr>
# 1 51 49 PPrel activ. --> hsa04350:1 hsa:7040 hsa:.
# 2 57 55 PPrel activ. --> hsa04350:2 hsa:151449 hs.
# 3 34 32 PPrel activ. --> hsa04350:3 hsa:3624 hsa:.
# 4 20 17 PPrel activ. --> hsa04350:4 hsa:4838
# 5 60 46 PPrel activ. --> hsa04350:5 hsa:4086 hsa:.
# # . with 45 more rows, and 5 more variables: genesymbol_source <chr>,
# # uniprot_source <chr>, kegg_id_target <chr>,
# # genesymbol_target <chr>, uniprot_target <chr>