Question: Oligo: GeneChip miRNA 4.0 analysis only returns control probes
1
gravatar for Keiron O'Shea
6 weeks ago by
Aberystwyth University
Keiron O'Shea10 wrote:

This is the first time I've ever posted on biostars, and I'm still not 100% sure of the etiquette - please feel free to berate me if I've forgotten to mention something as I'm here to learn.

I've been trying to understand this problem for the best part of 3 days now, and would like some help.

I've been given a number of CEL files that were produced using a GeneChip™ miRNA 4.0 Array. The samples have come from human tissue. In an attempt to obtain the expression data, I have written the following R script:

library(oligo)
# I've seen some issues with this failing to load on oligo::read.celfiles.
library(pd.mirna.4.0)

celFiles <- oligoClasses::list.celfiles("/home/uid/path/to/data", full.names=TRUE)
rawData <- oligo::read.celfiles(celFiles)

expressionSet <- oligo::rma(rawData)

This returns a "Large ExpressionSet" of 2.5Mb in size (from the original 17.9Mb loaded in as rawData). Sadly, when I want to take a look at the probe identifiers, only control probes seem to be returned (as in everything ending with "_st").

Biobase::featureNames(expressionSet)

I'm not 100% sure where I've gone wrong here as I've analysed quite a lot of miRNA in the past - although this is my first attempt at analysing GeneChip miRNA 4.0 Array data.

Any help would be massively appreciated.

Best wishes,

Keiron.

mirna oligo R • 145 views
ADD COMMENTlink modified 26 days ago by iraia.munoa60 • written 6 weeks ago by Keiron O'Shea10

I think I've confused myself a bit here. I would like some clarification, but would I be right in saying that st is just appended to the end of the probeset identifier to affirm that it's a Affymetrix annotation?

ADD REPLYlink written 6 weeks ago by Keiron O'Shea10

Hi! I am in the same point. Is it posible when doing the rma() for normalization to obtain more columns added in the output file, I mean other information apart from the ID name of the miRNA, such as the specie, transcript ID, etc. Or the only way to know this info is to take the output and do an intersect with the annotation table you provide in the previous link?

Thanks in advance,

Iraia

ADD REPLYlink written 26 days ago by iraia.munoa60
1

These will be contained in the NatAffx annotation file, which you can obtain from Affymetrix / ThermoFisher's website.

ADD REPLYlink written 26 days ago by Kevin Blighe41k
0
gravatar for Kevin Blighe
6 weeks ago by
Kevin Blighe41k
London, England
Kevin Blighe41k wrote:

Hey Keiron,

Yes, you are correct: The probes will have a _st appended to them. This indicates that they are ST array probes. You can see (and download) an annotation table for this array here: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GPL19117

Note that you can access the background corrected, quantile normalised, log (base 2) transformed expression values with exprs(expressionSet).

If you next need to perform differential expression analysis, you can try to use limma.

Note that you should not have to explicitly reference functions with :: - this is only important where namespace issues exist, i.e., where a function has the same name across 2 or more packages.

Kevin

ADD COMMENTlink modified 6 weeks ago • written 6 weeks ago by Kevin Blighe41k
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 2.3.0
Traffic: 823 users visited in the last hour