error in topGO?
1
0
Entering edit mode
7 months ago

Hi,

I am encountering a strange occurrence with topGO and I don't understand some of the output I am getting.

Using the following relatively straight forward code

myGOdata <- new("topGOdata", description="My project", ontology="BP", allGenes=geneList,  annot = annFUN.gene2GO, gene2GO = geneID2GO)

resultFisherweight <- runTest(myGOdata, algorithm="weight01", statistic="fisher")
resultFisherclassic <- runTest(myGOdata, algorithm="classic", statistic="fisher")
resultFisherelim <- runTest(myGOdata, algorithm="elim", statistic="fisher")

sig.tab <- GenTable(myGOdata, pFisweight = resultFisherweight,pFisClassic=resultFisherclassic, pelimFisher=resultFisherelim, topNodes = length(go), numChar=1000)##for all of the GO terms
goid<- "GO:0050890"
termStat(myGOdata, goid)

I generated output that looks like this. Note the 2 significant genes

                Annotated Significant Expected
GO:0050890       164           2     1.34

From some work out of R* and topGO, i find that the 2 significant genes (out of 88 possible ones) are: Ccalc.v3.00391 Ccalc.v3.14166

Now, look at the associated GO terms for each of these genes (below). The problem is that GO:0050890 is missing from Ccalc.v3.00391, but is present in Ccalc.v3.14166. But, topGO seems to detect GO:0050890 in Ccalc.v3.00391 even though it is not there! Why? Is this an error on topGO?

There are similar GO terms within Ccalc.v3.00391, like GO:0050896, but topGO shouldn't be detecting this, right?

Ccalc.v3.00391 GO:0000139, GO:0000166, GO:0001664, GO:0002027, GO:0002031, GO:0003674, GO:0003676, GO:0003677, GO:0003697, GO:0003700, GO:0003723, GO:0003724, GO:0003735, GO:0003887, GO:0003924, GO:0004190, GO:0004386, GO:0004519, GO:0004527, GO:0004674, GO:0005085, GO:0005102, GO:0005198, GO:0005509, GO:0005515, GO:0005524, GO:0005576, GO:0005634, GO:0005667, GO:0005730, GO:0005737, GO:0005743, GO:0005770, GO:0005773, GO:0005774, GO:0005783, GO:0005789, GO:0005791, GO:0005794, GO:0005829, GO:0005834, GO:0005840, GO:0005886, GO:0005887, GO:0005975, GO:0006260, GO:0006261, GO:0006269, GO:0006270, GO:0006310, GO:0006325, GO:0006342, GO:0006355, GO:0006364, GO:0006412, GO:0006417, GO:0006468, GO:0006508, GO:0006865, GO:0006914, GO:0007007, GO:0007154, GO:0007165, GO:0007186, GO:0007188, GO:0007219, GO:0007220, GO:0007411, GO:0007412, GO:0007417, GO:0007422, GO:0007480, GO:0007614, GO:0008053, GO:0008083, GO:0008152, GO:0008168, GO:0008233, GO:0008306, GO:0008408, GO:0008593, GO:0008833, GO:0009295, GO:0015074, GO:0015677, GO:0016020, GO:0016021, GO:0016032, GO:0016055, GO:0016209, GO:0016301, GO:0016310, GO:0016319, GO:0016324, GO:0016485, GO:0016491, GO:0016539, GO:0016740, GO:0016757, GO:0016779, GO:0016787, GO:0016788, GO:0016798, GO:0016813, GO:0016874, GO:0018024, GO:0018142, GO:0019001, GO:0019028, GO:0019033, GO:0019050, GO:0019062, GO:0019068, GO:0019076, GO:0019835, GO:0020002, GO:0030036, GO:0030424, GO:0030426, GO:0030430, GO:0030683, GO:0031151, GO:0031234, GO:0031532, GO:0031648, GO:0031683, GO:0032259, GO:0032974, GO:0033644, GO:0034729, GO:0035025, GO:0035333, GO:0035556, GO:0035641, GO:0039573, GO:0039664, GO:0039678, GO:0039693, GO:0042025, GO:0042407, GO:0042803, GO:0042981, GO:0043409, GO:0044659, GO:0045746, GO:0045751, GO:0045879, GO:0045887, GO:0046718, GO:0046729, GO:0046797, GO:0046798, GO:0046872, GO:0046983, GO:0048190, GO:0048260, GO:0050728, GO:0050770, GO:0050790, GO:0050896, GO:0051726, GO:0051920, GO:0055036, GO:0055037, GO:0055085, GO:0060158, GO:0060189, GO:0061617, GO:0070765, GO:0071627, GO:0071772, GO:0071897, GO:0075521, GO:0075606, GO:0090305, GO:0098003, GO:0098015, GO:0098024, GO:0098671, GO:0098869, GO:0098994, GO:0099001, GO:0099002, GO:0120177

Ccalc.v3.14166 GO:0000166, GO:0003677, GO:0003697, GO:0003723, GO:0003824, GO:0004045, GO:0004176, GO:0004252, GO:0004519, GO:0004617, GO:0005198, GO:0005388, GO:0005524, GO:0005576, GO:0005634, GO:0005737, GO:0005768, GO:0005769, GO:0005798, GO:0005886, GO:0005905, GO:0006210, GO:0006212, GO:0006259, GO:0006260, GO:0006281, GO:0006310, GO:0006508, GO:0006564, GO:0006816, GO:0006974, GO:0007399, GO:0008094, GO:0008233, GO:0008652, GO:0009267, GO:0009432, GO:0009986, GO:0016020, GO:0016021, GO:0016491, GO:0016616, GO:0016627, GO:0016740, GO:0016787, GO:0016788, GO:0017113, GO:0030426, GO:0031410, GO:0042995, GO:0043204, GO:0046872, GO:0050485, GO:0050890, GO:0051287, GO:0051536, GO:0055037, GO:0070403, GO:0070588, GO:0090305, GO:0140603

gene set topGO significant • 425 views
ADD COMMENT
0
Entering edit mode
7 months ago

https://amigo.geneontology.org/amigo/term/GO:0050890

GO:0050890 is a really high-level GO-term, 'cognition'. My hunch is that Ccalc.v3.00391 has one or more of the child-terms of GO:0050890 and topGO then chose to 'walk up the graph' of the GO-term children.

For starters, Ccalc.v3.00391 has the GO-term GO:0007614 short-term memory which is one of the child-terms of GO:0050890, see https://amigo.geneontology.org/amigo/term/GO:0007614#display-lineage-tab The other gene does not have this term, so TopGO used the highest-connected GO-term instead.

ADD COMMENT

Login before adding your answer.

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