This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Pipeline for analyzing Microarray & RNA-seq GSE files from NCBI GEO

0 down vote favorite I have very limited experience in R but would like to know if anyone can share their pipeline for analyzing GSE files from GEO, both for microarray and/or RNA-seq. The eventual goal would be to look at the differentially expressed genes.

For example, I would like to analyze GSE113590 which is a RNA-seq data and GSE47045 which is a microarray data.

The general consensus seems to be that you download the data using this:

source("http://bioconductor.org/biocLite.R")
biocLite("GEOquery")
library("GEOquery")
gset <- getGEO("GSE113590", GSEMatrix =TRUE)

But I'm not sure how to move forward from here, and there seems to be a different pipeline depending on whether it is a microarray / RNA-seq.

Thanks for your help.

rna-seq microarray r pipline
• 5,870 views • link
updated by ewre ewre • written by jsl • jsl

Understood, my apologies. I have taken down the post on stackoverflow.

• 0 views • link
written by jsl • jsl

2 answers

The distributions of RNAseq counts and array intensities is very different, hence the need for different packages. limma is the go-to package for microarray analysis, for RNAseq counts, the main options are edgeR, DESeq2 and limma, after using the voom transformtion on the counts.

• 0 views • link
written by h.mon h.mon

I see. But I cant seem to extract the counts on the gset. I know that for microarray, one could use 

gset <- getGEO("GSE47045", GSEMatrix =TRUE)
if (length(gset) > 1) idx <- grep("GPL6246", attr(gset, "names")) else idx <- 1
gset <- gset[[idx]]

str(exprs(gset))

I get 

 num [1:34760, 1:24] 12.85 11.2 7.58 13.42 6.72 ...
 - attr(*, "dimnames")=List of 2
  ..$ : chr [1:34760] "10338001" "10338003" "10338004" "10338017" ...
  ..$ : chr [1:24] "GSM1143711" "GSM1143712" "GSM1143713" "GSM1143714" ..

But when it comes to RNA-seq Illumina data, the "if" command line generated NULL counts..

• 0 views • link
written by jsl • jsl

Can you provide an example of a RNAseq accession which causes trouble?

• 0 views • link
written by h.mon h.mon

this is the one im trying to analyze GSE113590

• 0 views • link
written by jsl • jsl

Not sure if this one is helpful for you cause it deals with ArrayExpress which is EBI instead of NCBI. As I know, most of the datasets in GEO can be found in ArrayExpress.

• 1 views • link
written by ewre ewre

Log in to answer this question.

More posts like this