This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Dataset normalization before gene ontology analysis

I peformed GO analysis from a list of Genes using the Goseq package from bioconductor. After plotting the results, I could see that the bigger the gene list was, most counts it had from each category, there are some way for normalize this by the size of each gene list?

gene ontology normalization rna-seq bioconductor
• 2,572 views • link
updated by svlachavas svlachavas • written by tiago211287 tiago211287

1 answer

Dear Tiago211287,

I believe that you get this result from ploting, because generally in RNA-seq the length of one gene is crusial regarding the levels of its expression (which in turn is associated with power). Thus, one way to possibly adjust for this when performing a GO analysis with RNA-seq data, is to use prior the function ?nullp:

nullp(DEgenes, genome, id, bias.data=NULL,plot.fit=TRUE)

This will produce a set of relative weights which are "somehow proportional" to how "big" are your input genes.

Then, you can feed it directly to goseq()

Hope that helps,
Efstathios

• 0 views • link
updated by Ram Ram • written by svlachavas svlachavas

I did that in goseq, generating a pwf(Probability Weighting Function)

• 0 views • link
written by tiago211287 tiago211287

Well then, excuse me but I misunderstood your question. So, did you meant that you used more than one gene lists ? If so, (without being an expert on RNA-seq analysis) why do you want to normalize for the size of each list?

• 0 views • link
updated by Ram Ram • written by svlachavas svlachavas

Log in to answer this question.

More posts like this