[Bioc-sig-seq] 3'-5' position of Affymetrix probe sets

Christos Hatzis christos.hatzis at nuverabio.com
Wed Mar 2 17:49:06 CET 2011


Hello,

For many genes, as in the example below, there are more than one Affy probe
sets that interrogate the same transcript.  Is there a way to rank those
probe sets relative to their proximity to the 3' end of the gene?  I wasn't
sure whether there are annotation packages that provide the genomic
coordinates of the probes in a probe set.

> testids <- c("206527_at","209459_s_at","209460_at")
> library(biomaRt)
> ensembl <- useMart("ensembl", dataset = "hsapiens_gene_ensembl")
Checking attributes ... ok
Checking filters ... ok
> probAttrs <- c("affy_hg_u133a", "hgnc_symbol", "chromosome_name", 
 	"start_position", "end_position", "strand") 
> testres <- getBM(attributes = probAttrs, filters = "affy_hg_u133a", values
= 
	testids, mart = ensembl) 
> testres
  affy_hg_u133a hgnc_symbol chromosome_name start_position end_position
strand
1   209459_s_at        ABAT              16        8768422      8878432
1
2     209460_at        ABAT              16        8768422      8878432
1
3     206527_at        ABAT              16        8768422      8878432
1

Thank you.   
-Christos

Christos Hatzis, Ph.D.
Nuvera Biosciences, Inc.
400 West Cummings Park, Suite 5350
Woburn, MA 01801



More information about the Bioc-sig-sequencing mailing list