Background choice for GO analysis of RNA-Seq
1
0
Entering edit mode
3 months ago
Júlia • 0

Hello everyone! I have a question about which Background to use for the Gene Ontology Analysis of the DEGs obtained by RNA-Seq. According to the instructions on the GO website:

" The reference list should be the list of all the genes from which your smaller analysis list was selected. For example, in a list of differentially expressed genes, the reference list should only contain genes that were detected at all in the experiment, and thus potentially could have been on a list of genes derived from the experiment."

In this context, I'm using the list of genes obtained by the exactTest() function of edgeR, which gives a list based on "pseudo - counts". Since this "pseudo-counts" are for internal use of edgeR pipeline, I don´t know if it is adequate in this case. Any thoughts?

Thank you!

Analysis Background RNA-Seq GO • 350 views
ADD COMMENT
1
Entering edit mode

Did you read the paper mentioned in this post?

ADD REPLY
3
Entering edit mode
3 months ago

Run filterByExpr prior to testing and roll with the remaining genes as your background.

ADD COMMENT
0
Entering edit mode

Hi! We actually use another filtration:

y <- DGEList(counts=count, group=group) 
#Normalization
y <- calcNormFactors(y) 
#Filtering
keep <- rowSums(cpm(y)>1) >= 3 
y <- y[keep, , keep.lib.sizes=FALSE]
# Calculate dispersions 
y <- estimateCommonDisp(y)
y <- estimateTrendedDisp(y)
y <- estimateTagwiseDisp(y)
# Testing for DE genes 
et <- exactTest(y, pair=c("control", "case"))

This approach works as well? We have four samples in case group and six samples in controls, with only one factor. Thanks in advance!

ADD REPLY
1
Entering edit mode

Well, you're free to take whatever approach you like. Just use whatever passes your filtering criteria as the background.

ADD REPLY

Login before adding your answer.

Traffic: 756 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