[BioC] unmapped keys in hugene10stprobeset.db
Paul Shannon
pshannon at systemsbiology.org
Mon Aug 16 23:17:54 CEST 2010
Here's an annotation question someone might be able to help me out with. I'll be grateful.
Affymetrix describes their 'GeneChip Human Gene 1.0 ST Array':
Each of the 28,869 genes is represented on the array by approximately 26 probes spread across the full length of the gene, providing a more complete and more accurate picture of gene expression than 3’ based expression array designs. ... The Gene 1.0 ST Array uses a subset of probes from the Human Exon 1.0 ST Array and covers only well-annotated content.
This sounds to me like affy started with sequence from exons of ~29k genes and created probes.
But when I look at the bioc annotation for this chip (hugene10stprobeset.db, hugene10sttranscriptcluster.db), I find that about 7% of the probes are NOT annotated to geneIDs. The sibling array, hugene10sttranscriptclusterENTREZID, though smaller, has a higher proportion of unmapped probes.
library (hugene10stprobeset.db)
library (hugene10sttranscriptcluster.db)
bm = hugene10stprobesetENTREZID
length (keys (bm)) # 257022
count.mappedkeys (bm) # 238141
# unmapped: 18881
cm = hugene10sttranscriptclusterENTREZID
length (keys (cm)); # 33257
count.mappedkeys (cm) # 21787
The same proportions (unmapped/mapped seems to hold true for hugene10stprobesetENSEMBL as well.
Can anyone suggest where I can get entrez geneID annotations for these unmapped probes? Or otherwise clear up my confusion?
Thanks!
- Paul
More information about the Bioconductor
mailing list