On 3/5/2013 4:45 PM, Naxerova, Kamila wrote:
> Dear all,
> I am analyzing a set of Affymetrix Mouse Gene 2.0 ST arrays. I am
quite familiar with 3'-biased chips, but this is my first time looking
at data from WT arrays. I have a few general questions -- any advice
would be appreciated to speed up my learning process.
> 1) I have already read on this mailing list that the good old affy
package does not work well with WT arrays (can anybody point me to any
literature on why that is?). So I have installed the oligo and xps
packages -- what are the advantages/disadvantages for each? Any
opinions on which one is the right "starter kit"?
The affy package was never intended to work with these arrays - it was
designed specifically for the 3' biased arrays, which had pre-defined
probesets, and which didn't share probes between probesets. In
the makecdfenv package is designed to work with the old style CDF
packages, and Affy has never released a CDF for these new chips that
they are willing to support in any meaningful way.
There were some changes made to the affy package in order to
the fact that probes could be shared between probesets, and it is
possible to use functions in affxparser to re-create conventional CDF
packages using the newer pgf and clf files. So hypothetically you
still use the affy package (and hypothetically you could still use an
Apple IIe for all your computing needs, but that's crazy, so let's
I don't think you will find much difference between oligo and xps,
than the fact that xps requires the additional installation of ROOT.
might play around with both and see which suits you better.
I should throw in my obligatory cautionary statement about summarizing
Gene ST data at the probeset (as compared to the transcript) level. If
you look at the number of probes/probeset, there are a huge number
< 4 probes. So hypothetically you can do this, but I wouldn't.
> 2) I see with some dread that there seems to be no annotation
package for the 2.0 array yet. I have never built my own... any quick
bullet points on how I would go about doing that for a WT array?
No dread should be required. All you need to do is get the
transcript-level annotation file from Affy
and then the AnnotationForge, mouse.db0, and org.Mm.eg.db packages.
outputDir = ".",
manufacturer = "Affymetrix",
chipName = "Human Gene 2.1 ST Array",
manufacturerUrl = "http://www.affymetrix.com
author = "Kamila Naxerova",
maintainer = "Kamila Naxerova <naxerova at="" fas.harvard.edu="">")
should do the trick. You can then install directly from within R by
> 3) It seems that RMA is also used for normalization of WT arrays, so
that part I am comfortable with. But are there any differences in
preprocessing between 3' and WT arrays that I should watch out for?
Not really. I don't use xps, so cannot say for certain how you do
with that package, but with oligo it's a simple
abatch <- read.celfiles(list.celfiles())
eset <- rma(abatch)
To normalize and summarize at the transcript level. Note however that
the annotation for the resulting ExpressionSet will be the
pd.mogene.2.0.st.v1 package, and if you use annotation(eset) in any
further calls to do gene annotation, it won't work out. You need to
annotation(eset) <- "mogene20sttranscriptcluster.db"
One further note: the intronic controls (especially) have an
habit of popping up in lists of differentially expressed genes. This
IMO likely due to mRNA that has not been fully processed to excise the
introns, but regardless, these probesets tend to have no annotation at
all, so are not useful without extra work to figure out what they are
supposed to be measuring. My usual MO is to just summarily excise them
after e.g., the eBayes() step of an analysis using limma. If you are
interested, there is a function in the affycoretools package called
getMainProbes() that will do this for you.
> Thanks so much!
> Bioconductor mailing list
> Bioconductor at r-project.org
> Search the archives:
James W. MacDonald, M.S.
University of Washington
Environmental and Occupational Health Sciences
4225 Roosevelt Way NE, # 100
Seattle WA 98105-6099