Question: What people do with the probes that doesn't have any GENE SYMBOL in annotation databases?
gravatar for arronar
2.5 years ago by
arronar230 wrote:


I'm doing some microarray analysis, and I figured out that some probe ids have not any GENE SYMBOL. Am I suppose to delete those probes or let them exist ?

The way I apply annotation info for an Affymetrix HGU133plus2 array after RMA application, is the following:

Symbols = unlist(mget(probes, hgu133plus2SYMBOL, ifnotfound=NA))
Entrez_IDs = unlist(mget(probes, hgu133plus2ENTREZID, ifnotfound=NA))

Is something wrong with my code or such behavior is expected ? What do you guys do with these NA genes ?

annotation microoarrays • 1.3k views
ADD COMMENTlink modified 2.5 years ago by Maxime Lamontagne2.2k • written 2.5 years ago by arronar230
gravatar for Maxime Lamontagne
2.5 years ago by
Maxime Lamontagne2.2k wrote:

It’s not unusual in older array. Due to changes in gene annotation and reference sequence, some probe set may be outdated (probes not mapping to any gene in the new build). To find a gene for your probe set, you need to blat all the probes of the probe set and look where they are mapping. ​

ADD COMMENTlink written 2.5 years ago by Maxime Lamontagne2.2k

It would be helpful enough if you could explain the procedure with more details because I'm not so experienced with such tasks. Thank you very much.

ADD REPLYlink written 2.5 years ago by arronar230

You could get the probe sequences from the available Affymetrix annotation files, and then paste the sequences into a BLAT tool like the one at UCSC here to find where the probes map to the latest genome build.

ADD REPLYlink written 2.5 years ago by Ahill1.8k

I think this paper is doing exactly what you want.

Framework for reanalysis of publicly available Affymetrix® GeneChip® data sets based on functional regions of interest

ADD REPLYlink written 2.5 years ago by Maxime Lamontagne2.2k

Those probe sets that didn't map in gene symbols, were about 12000 . So I needed a more automatic way to retrieve a gene symbol for them and not one by one . So I found this site in which you configure the input as Affyid and the output as gene symbol. Then, I wrote these lines in R to retrieve the results (for anyone else that will need such a thing):

for ( probe in missed_probes ){
  q <- paste(sep="","",probe,"&outputs=genesymbol&taxonId=9606&format=row")
  results <- rbind(results,fromJSON(txt=q))

The thing now is that for many probe sets, there are more than one gene symbols. :s And maybe I have to decide randomly which one to keep :-p

EDIT: Also I found these sites that do this work.

ADD REPLYlink modified 2.5 years ago • written 2.5 years ago by arronar230
gravatar for Hussain Ather
2.5 years ago by
Hussain Ather940
National Institutes of Health, Bethesda, MD
Hussain Ather940 wrote:

You probably want to get the current annotated reference (maybe look here on Bioconductor) and map the sequence of each probe to it. It's possible that the gene list you were using is out of date.

ADD COMMENTlink written 2.5 years ago by Hussain Ather940

I'm already using the lateset version of that package.

Package hgu133plus2.db version 3.2.3

ADD REPLYlink written 2.5 years ago by arronar230
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: 2115 users visited in the last hour