Question: The current annotation for Affymetrix HTA 2.0 array is problematic
1
gravatar for tangboyun
4.3 years ago by
tangboyun30
China
tangboyun30 wrote:

It seems that Affy is silently changing their procedure for the probeset annotation of HTA 2.0.

cut -d"," -f 1,10,13 HTA-2_0.na33.hg19.probeset.csv

"probeset_id","psr_id","probeset_type"
"2824546_st","2824546_st","normgene->exon"
"2824549_st","2824549_st","normgene->exon"
"2824551_st","2824551_st","normgene->exon"
"2824554_st","2824554_st","normgene->exon"
"2827992_st","2827992_st","normgene->exon"
cut -d"," -f 1,10,13 HTA-2_0.na34.hg19.probeset.csv

"probeset_id","psr_id","probeset_type"
"18670005","2824546_st","normgene->exon"
"18670007","2824549_st","normgene->exon"
"18670009","2824551_st","normgene->exon"
"18670011","2824554_st","normgene->exon"
"18670020","2827992_st","normgene->exon"
"18670022","2827995_st","normgene->exon"

cut -d"," -f 1,10,13 HTA-2_0.na35.hg19.probeset.csv

"probeset_id","psr_id","probeset_type"
"18670000","---","Exonic normalization control (Positive Control)"
"18670001","---","Exonic normalization control (Positive Control)"
"18670002","---","Exonic normalization control (Positive Control)"
"18670003","---","Exonic normalization control (Positive Control)"
"18670004","---","Exonic normalization control (Positive Control)"

The probesets in the "main" category are not affected by this changing, as their probeset_name are the same as probeset_id. But all the ctrl probes lost their probeset_type annotations.

> all.probes <- getProbeInfo(dat,field=c('type'), target='probeset')
> unique(all.probes$type)
[1] "main"                       NA 

In na35 probeset annotations, the probeset_type is also changed:

na35, previous
"Positive Control", "pos_control"
"ERCC (External RNA Controls Consortium) step control","control->affx->ercc->step"
"Antigenomic background control", "control->bgp->antigenomic"
"Intronic normalization control (Negative Control)", "normgene->intron"
 "Exonic normalization control (Positive Control)","normgene->exon"

Uers could build na35 annotation db by manually doing the above conversion on "HTA-2_0.na35.hg19.probeset.csv". Currently, Affy hasn't changed the probeset annotation for MTA 1.0 and RTA 1.0 yet ;). 

 

oligo pdinfobuilder pd.hta.2.0 • 1.4k views
ADD COMMENTlink modified 4.3 years ago • written 4.3 years ago by tangboyun30
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 225 users visited in the last hour