SPCC548.03c from "GO annotation dataset" (plus question, how to query the "with" filed?)

[ValWood]

There are 2 genes represented in pombemine for
SPCC548.03c.01
SPCC548.03c.02

Neither @kimrutherford or I can figure out where these originate (they are isoform IDs but not separate genes).
We can't see where we export these. Are they coming from another source?

thanks
v

[ValWood]

I should mention that KIms digging says this comes from:
"GO Annotation data set"

but this is not a gene name in our annotation, it's an isoform identifier. ...and I can't see that we have used it in the PomBase GO annotations.

[ValWood]

@kimrutherford

[kimrutherford]

There are 2 genes represented in pombemine for
SPCC548.03c.01
SPCC548.03c.02

That should be:
SPCC548.03c.1
SPCC548.03c.2

One case where we use these IDs in a place that is mostly gene IDs is the "with" column of the GAF file. For example:

PomBase SPCC1906.03     wtf19           GO:0005737      PMID:32032353   ISS     PomBase:SPCC548.03c.1   C       wtf meiotic drive antidote Wtf19                protein taxon:4896      20200914      PomBase part_of(CL:0000607)
PomBase SPCC1906.03     wtf19           GO:0005737      PMID:32032353   ISS     PomBase:SPCC548.03c.2   C       wtf meiotic drive antidote Wtf19                protein taxon:4896      20200914      PomBase part_of(CL:0000607)

Maybe they are being misunderstood as gene identifiers in that context?

[ValWood]

Right that makes sense.Hmm this is a real edge case. We can infer the location of the different specific versions of this protein (poison and antidote) , and in this case we have specified the isoform(alternative transcript) ID in the with column.

I checked the docs http://geneontology.org/docs/go-annotation-file-gaf-format-2.1/#with-or-from-column-8 to see if this field is restricted to "gene" and it isn't but isoform is not documented:

I suspect if we discussed this the format would be the same as an allele, so it would be
DB:gene_symbol[isoform_symbol]

I will check this with GO

[ValWood]

GO ticket
geneontology/helpdesk#394

@danielabutano I have taken this ticket and I 'll report back.
It might be possible to improve how InterMIne handles this field if IDs can be typed.
Note that "protein complex" identifiers can also be used in this field (I am not sure how?)

[ValWood]

@danielabutano one thing I did wonder was about the value of adding the genes from the "with" field.
The genes of interest will be loaded already from other routes.

[ValWood]

OK I have a response from GO.
geneontology/helpdesk#394
basically it isn't safe to assume that the IDs in the "with" field refer to genes.

But I think that is OK, we don't need to use these "with field" entries in any queries. They are really arbitrary sources of support for an annotation, but they aren't useful for querying, and therefore probably shouldn't be loaded as independent genes (as long as the string is visible (Prefix plus ID) people can look up the sources if they want to validate a specific annotation).

I wanted to see what the with field output looks like but I can't get a query to output this column. Where are the instructions for this?

[danielabutano]

Hi @ValWood, this query shows the genes created from the with column

[danielabutano]

this query is more precise

below the XML if you want to import the it:
<query model="genomic" view="GOEvidence.withText GOEvidence.with.name GOEvidence.with.primaryIdentifier GOEvidence.with.dataSets.name" constraintLogic="(A)" sortOrder=""> <constraint path="GOEvidence.with.dataSets.name" op="NONE OF" code="A"><value>Pombe-gene data set</value><value>cerevisiae-orthologs data set</value><value>human-orthologs data set</value><value>BioGRID interaction data set</value><value>PomBase disease data set</value><value>PomBase phenotypes data set</value></constraint> </query>

[ValWood]

Got it. I forgot I need to switch to "GO evidence"