how i can get the name of genes in pathway?
1
0
Entering edit mode
chavaleab • 0
@330334d9
Last seen 20 months ago
Israel

Enter the body of text here

Code should be placed in three backticks as shown below



peakonly_DMSO2 <- readPeakFile("path_to_peak_bed_DMSO")
geneonly_DMSO <- seq2gene(peakonly_DMSO2, tssRegion = c(-1000, 1000), flankDistance = 3000, TxDb=txdb)
pathway2only_DMSO <- enrichPathway(geneonly_DMSO)
dotplot(pathway2only_DMSO)

i get dotplot graph and in the first annotation i get that there is 130 gene that involve, and i want to know what is this gene, how can i extract it? if i write:

pathway2only_DMSO[1]$geneID
gene_list <- strsplit(pathway2only_DMSO[1]$geneID, "/")[[1]]

i get: this list: gene_list [1] "2475" "842" "3065" "5525" "80351" "1937" "572" "84447" "4041" "9846" "8322"
[12] "867" "7450" "4040" "3845" "8496" "2065" "9612" "2100" "207" "83464" "5604"
[23] "8312" "7249" "930" "5432" "9611" "6774" "5296" "7311" "5714" "6868" "10018" [34] "182" "6238" "7529" "22916" "2261" "55109" "9542" "5159" "57670" "11160" "11064" [45] "79109" "5293" "1185" "3932" "5906" "89780" "818" "10718" "8036" "6934" "5682"
[56] "408" "2120" "51411" "5908" "4254" "283455" "1024" "8660" "5529" "5663" "9321"
[67] "4140" "4088" "8826" "23164" "399687" "4763" "5717" "6777" "10197" "10014" "3690"
[78] "5573" "2885" "57674" "4089" "5605" "208" "238" "6654" "27436" "27248" "7039"
[89] "6772" "2066" "598" "5594" "613" "9978" "10533" "1499" "22906" "79718" "5290"
[100] "152831" "5156" "27123" "2247" "55534" "6500" "8841" "2246" "815" "10817" "5528"
[111] "2309" "2534" "2099" "9734" "1956" "1523" "27044" "8805" "80853" "673" "8454"
[122] "8658" "3084" "23189" "6907" "22866" "55869" "5519" "7314" "471" so how can i know the full name gene from this R names?

thank you very nutch for any answer!!!!

chipseq ChIPseeker • 1.4k views
ADD COMMENT
1
Entering edit mode
@james-w-macdonald-5106
Last seen 14 hours ago
United States

Those aren't 'R names'. Those are NCBI Gene IDs, which you can map using the correct orgDb package. I assume you are using human data.

> library(org.Hs.eg.db)
> mapIds(org.Hs.eg.db, c("2475", "842", "3065", "5525", "80351", "1937", "572", "84447", "4041", "9846", "8322"), "SYMBOL","ENTREZID")
'select()' returned 1:1 mapping between keys and columns
     2475       842      3065      5525     80351      1937       572     84447 
   "MTOR"   "CASP9"   "HDAC1" "PPP2R5A"   "TNKS2"   "EEF1G"     "BAD"   "SYVN1" 
     4041      9846      8322 
   "LRP5"    "GAB2"    "FZD4" 

## OR
> select(org.Hs.eg.db, c("2475", "842", "3065", "5525", "80351", "1937", "572", "84447", "4041", "9846", "8322"), "SYMBOL","ENTREZID")
'select()' returned 1:1 mapping between keys and columns
   ENTREZID  SYMBOL
1      2475    MTOR
2       842   CASP9
3      3065   HDAC1
4      5525 PPP2R5A
5     80351   TNKS2
6      1937   EEF1G
7       572     BAD
8     84447   SYVN1
9      4041    LRP5
10     9846    GAB2
11     8322    FZD4
ADD COMMENT
0
Entering edit mode

thanks!!!!!!!!!!!!!!!

ADD REPLY

Login before adding your answer.

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