Question: Oligo: GeneChip miRNA 4.0 analysis only returns control probes
1
gravatar for Keiron O'Shea
19 months 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 • 1.3k views
ADD COMMENTlink modified 7 weeks ago by anand0 • written 19 months 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 19 months 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 19 months ago by iraia.munoa100
1

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

ADD REPLYlink written 19 months ago by regmkbl66k

Hi Keiron,

Thanks for the detailed response. I am working on Affymetrix miRNA.3.0 and have normalized my data using 'oligo' followed by differential expression analysis using limma package.

I was wondering how can I retrieve miRNA expression p-values after oligo normalization?

Anand

ADD REPLYlink written 7 weeks ago by anand0

Hi Anand. please use limma to perform the differential expression analysis. see the toptable() function.

ADD REPLYlink written 7 weeks ago by regmkbl66k
0
gravatar for regmkbl
19 months ago by
regmkbl66k
regmkbl66k 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 19 months ago • written 19 months ago by regmkbl66k
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: 2026 users visited in the last hour