how can i convert gene symbol to affymetrix IDs ?

Question: how can i convert gene symbol to affymetrix IDs ?

6.0 years ago by

Mo • 920

Mo • 920 wrote:

I have a list of Gene IDs from human, I would like to convert it to Affymetrix IDs,

examples gene IDs

GNAI1
HLA-DQA1
WDR1
LTBP2
MXRA5
MMS19
SEC13
C6
IL16
FAM82A1
GBE1
TUBAL3
CYB5R3
ANXA1
APCS
COL14A1

Do you know how I can do it in R and based on website?

I have tried the following ID convertors online but none of them worked

http://biodbnet.abcc.ncifcrf.gov/db/db2db.php

https://gemdbj.nibio.go.jp/dgdb/ConvertOperation.do

David only convert to three Affymetrix such as 3PRIME_IVT_ID, EXON_gene_ID and SNP_ID

http://david.abcc.ncifcrf.gov/conversion.jsp

I also tried Ensemble Biomart like below but did not convert anything

1- first I web to website

http://www.ensembl.org/biomart/martview/b3cd9cdf1c6bc6862d64b0e60934ef0e

2- select " Ensemble Genes 78" and "homo sapiens genes (GRCh38)"

3- then click on "filter" ------> "gene" -------> and then paste the example list of genes there

4- then "Attribute" ------> "External" ------> and select "Affy HG U133A probeset" Then click on "Results" but nothing showed up

affymetrix microarray R • 9.3k views

ADD COMMENT • link •

modified 6.0 years ago by komal.rathi • 3.7k • written 6.0 years ago by Mo • 920

If nothing showed up then you did something wrong. I used your example gene symbols (not IDs) at Ensembl BioMart and the conversion works fine (click "Results" after loading this URL). Perhaps you forgot to specify "HGNC symbol(s)" as the ID type in Filters?

This can also be done using R/biomaRt. Please search this site for numerous answers to similar questions.

ADD REPLY • link modified 6.0 years ago • written 6.0 years ago by Neilfws ♦ 49k

Which organism and which array? There are annotation packages in Bioconductor for most common ones that can facilitate the conversion.

ADD REPLY • link modified 6.0 years ago • written 6.0 years ago by Devon Ryan ♦ 98k

@Devon Ryan human cell lines like promyelocytic leukemia , breast adenocarcinoma. HG-U133A microarray platform

ADD REPLY • link modified 6.0 years ago • written 6.0 years ago by Mo • 920

6.0 years ago by

komal.rathi • 3.7k

Children's Hospital of Philadelphia, Philadelphia, PA

komal.rathi • 3.7k wrote:

This will give you Affy Probe ID mappings to the Gene Symbols for HG-U133A microarray platform:

source("http://bioconductor.org/biocLite.R")
biocLite("hgu133a.db")
library(hgu133a.db)
library(annotate)
x <- hgu133aSYMBOL
# Get the probe identifiers - gene symbol mappings
mapped_probes <- mappedkeys(x)
# Convert to a dataframe
genesym.probeid <- as.data.frame(x[mapped_probes])
head(genesym.probeid)
   probe_id symbol
1   1053_at   RFC2
2    117_at  HSPA6
3    121_at   PAX8
4 1255_g_at GUCA1A
5   1316_at   THRA
6   1320_at PTPN21

Then you can subset this dataframe using your gene list.

ADD COMMENT • link modified 6.0 years ago • written 6.0 years ago by komal.rathi • 3.7k

@komal.rathi Thanks for your message Lets assume your solution works. look the list is empty

test <- structure(list(V1 = structure(c(8L, 9L, 16L, 11L, 13L, 12L, 14L, 
3L, 10L, 6L, 7L, 15L, 5L, 1L, 2L, 4L), .Label = c("ANXA1", "APCS", 
"C6", "COL14A1", "CYB5R3", "FAM82A1", "GBE1", "GNAI1", "HLA-DQA1", 
"IL16", "LTBP2", "MMS19", "MXRA5", "SEC13", "TUBAL3", "WDR1"), class = "factor")), .Names = "V1", class = "data.frame", row.names = c(NA, 
-16L))

x <- test
# Get the probe identifiers - gene symbol mappings
mapped_probes <- mappedkeys(x)
# Convert to a dataframe
genesym.probeid <- as.data.frame(x[mapped_probes])
head(genesym.probeid)

# data frame with 0 columns and 6 rows

ADD REPLY • link modified 6.0 years ago • written 6.0 years ago by Mo • 920

hgu133aSYMBOL maps the symbols to probe IDs, what I would do is keep everything as @komal.rathi and then do:

genes <- c("ANXA1", "APCS", "C6", "COL14A1", "CYB5R3", "FAM82A1", "GBE1", "GNAI1", "HLA-DQA1", "IL16", "LTBP2", "MMS19", "MXRA5", "SEC13", "TUBAL3", "WDR1")

results <- genesym.probeid[which(genesym.probeid$symbol %in% genes),]

that should make it

ADD REPLY • link modified 6.0 years ago • written 6.0 years ago by TriS • 4.3k

Mo test, as defined by you, is a dataframe. x, as defined by me, is an object of class ProbeAnnDbBimap. Hope you understand why the list looks empty.

TriS That's what the last line in my answer implies.

ADD REPLY • link modified 6.0 years ago • written 6.0 years ago by komal.rathi • 3.7k

@komal.rathi thanks in general you first make the list in genesym.probeid then I should extract my list based on that

@TriS thanks

ADD REPLY • link modified 6.0 years ago • written 6.0 years ago by Mo • 920

6.0 years ago by

TriS • 4.3k

United States, Buffalo

TriS • 4.3k wrote:

Biomart is probably what you look for

-- edit --

I saw you used it from ensembl, I believe this is easier to use.

using the list you provided it worked for me for Affy hgu133 a2

ADD COMMENT • link modified 6.0 years ago • written 6.0 years ago by TriS • 4.3k

@TriS Yes this is easier to do, but many IDs are missing. do you know a way to do it in R ? The order of the list changes and there are many overlapped names as well

ADD REPLY • link modified 6.0 years ago • written 6.0 years ago by Mo • 920

6.0 years ago by

geek_y ♦ 11k

Barcelona

geek_y ♦ 11k wrote:

David conversion tool ? http://david.abcc.ncifcrf.gov/conversion.jsp

ADD COMMENT • link written 6.0 years ago by geek_y ♦ 11k

@Geek_y . there are quite some limitation with david. e.g. limited to 1000 genes, not up dated etc. you can see above that I updated my question

ADD REPLY • link modified 6.0 years ago • written 6.0 years ago by Mo • 920

Please log in to add an answer.

Similar posts • Search »