Reading Affymetrix data from pd.hugene.2.0.st using oligo for RNA degradation plot
1
0
Entering edit mode
5.3 years ago
Vasu ▴ 540

Hi,

I have some Affymetrix CEL files that use the HuGene 2.0 ST design.

I'm using oligo package to read CEL files. My question is how can I get RNA degradation plot using oligo? 

Thanks

oligo rnadegradation pd.hugene.2.0.st affymetrix • 2.4k views
ADD COMMENT
0
Entering edit mode
5.3 years ago

You can get probe sequences and coordinates (usually) from the Affy website:

You'll want to extract the probe-level summarization and you can then plot intensity as a function of position with your genes of interest.

With respect to the probe-level summarization from the oligo package, the code should look something like this:

library(oligo)

celFiles = list.celfiles("/path/to/folder", full.names=T)
rawData = read.celfiles(celFiles)
probesetData = rma(rawData, target="probeset")
probeset.expression = exprs(probesetData)
probesets = rownames(probesetData)

However, selecting the optimal genes may be tricky. I've tried using the validated genes from this paper, but they didn't seem to show much of a position degradation effect (in the array data that I have checked so far):

ADD COMMENT
0
Entering edit mode

Hi Charles,

Thanks for the reply. I have used "xps" package for this.

And I got the following plot. Help Needed to Interpret RNA degradation plot.

Thank you in Advance

ADD REPLY
0
Entering edit mode

I'm not exactly sure how to interpret the figure - I would have expected more of a 5' degradation bias. At least the intensities are relatively consistent within each gene, so I would say the signal isn't too noisy.

ADD REPLY

Login before adding your answer.

Traffic: 1218 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6