[BioC] Definition of probes in lumiHumanAll annotation database

Mark Dunning mark.dunning at gmail.com
Mon Apr 2 12:11:31 CEST 2012


Hi Javier,

I'm afraid I'm not clear about what you mean by probe definitions.
Rather than being Refseq-centric, the illuminaHumanv1.db,
illuminaHumanv2.db, illuminaHumanv3.db, illuminaHumanv4.db packages
provide an interface to individual probe sequences and quality
information, if this is helpful?

probes <- unlist(mget("A1CF", revmap(illuminaHumanv4SYMBOL)))

> probes
         A1CF1          A1CF2          A1CF3
"ILMN_1779670" "ILMN_1806310" "ILMN_2383229"

>  mget(probes, illuminaHumanv4PROBESEQUENCE)
$ILMN_1779670
[1] "GGCACATGCCCAGAGCCAGAAGCGAGCATGAGCACAGCAATTCCTGGCCT"

$ILMN_1806310
[1] "GAGGTCTACCCAACTTTTGCAGTGACTGCCCGAGGGGATGGATATGGCAC"

$ILMN_2383229
[1] "TGCTGTCCCTAATGCAACTGCACCCGTGTCTGCAGCCCAGCTCAAGCAAG"

>  mget(probes, illuminaHumanv4GENOMICLOCATION)
$ILMN_1779670
[1] "chr10:52610480:52610529:-"

$ILMN_1806310
[1] "chr10:52566496:52566545:-"

$ILMN_2383229
[1] "chr10:52566587:52566636:-"

Regards,

Mark

2012/4/2 Javier Pérez Florido <jpflorido at gmail.com>:
> Dear list,
> I've checked all the annotation elements given by lumiHumanAll.db and, as
> far as I know, none of them provides the definition of probes. If I'm not
> wrong, this definition can distinguishes isoforms and other factors.
> How can I get such information from the lumiHumanAll annotation file? For
> example, gene A1CF has three different probes, corresponding each one to a
> different transcript variant, but I cannot obtain such information from the
> lumiHumanAll database.
>
> Thanks,
> Javier
>
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at r-project.org
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives:
> http://news.gmane.org/gmane.science.biology.informatics.conductor



More information about the Bioconductor mailing list