Relationships in the GO

Alex Warwick Vesztrocy, March 2016

For some analyses, it is possible to only use the is_a definitions given in the Gene Ontology.

However, it is important to remember that this isn't always the case. As such, GOATOOLS includes the option to load the relationship definitions also.

Loading GO graph with the relationship tags

This is possible by using the optional_attrs argument, upon instantiating a GODag. This can be done with either the full or the basic version of the ontology. Here, the full version shall be used.



In [1]:

    
from goatools.obo_parser import GODag
import wget

go_fn = wget.download('http://geneontology.org/ontology/go.obo')
go = GODag(go_fn, optional_attrs=['relationship'])









    



100% [........................................................................] 34671316 / 34671316go.obo: format-version(1.2) data-version(releases/2016-04-27)





    



load obo file go.obo
46518





    









    



 nodes imported

Viewing relationships in the GO graph

So now, when looking at an individual term (which has relationships defined in the GO) these are listed in a nested manner. As an example, look at GO:1901990 which has a single has_part relationship, as well as a regulates one.



In [2]:

    
eg_term = go['GO:1901990']



In [3]:

    
eg_term









    Out[3]:





GOTerm('GO:1901990'):
  name:regulation of mitotic cell cycle phase transition
  relationship: 2 items
    regulates: 1 items
      GO:0044772	level-05	depth-05	mitotic cell cycle phase transition [biological_process] 
    has_part: 1 items
      GO:0051437	level-07	depth-12	positive regulation of ubiquitin-protein ligase activity involved in regulation of mitotic cell cycle transition [biological_process] 
  level:6
  is_obsolete:False
  namespace:biological_process
  id:GO:1901990
  depth:7
  parents: 2 items
    GO:0007346	level-05	depth-05	regulation of mitotic cell cycle [biological_process] 
    GO:1901987	level-06	depth-06	regulation of cell cycle phase transition [biological_process] 
  children: 7 items
    GO:1901991	level-07	depth-08	negative regulation of mitotic cell cycle phase transition [biological_process] 
    GO:1901992	level-07	depth-08	positive regulation of mitotic cell cycle phase transition [biological_process] 
    GO:2000045	level-07	depth-08	regulation of G1/S transition of mitotic cell cycle [biological_process] 
    GO:0010389	level-07	depth-08	regulation of G2/M transition of mitotic cell cycle [biological_process] 
    GO:0007096	level-07	depth-08	regulation of exit from mitosis [biological_process] 
    GO:0030071	level-07	depth-10	regulation of mitotic metaphase/anaphase transition [biological_process] 
    GO:0060564	level-07	depth-13	negative regulation of APC-Cdc20 complex activity [biological_process] 
  _parents: 2 items
    GO:0007346
    GO:1901987
  alt_ids: 0 items

These different relationship types are stored as a dictionary within the relationship attribute on a GO term.



In [4]:

    
print(eg_term.relationship.keys())









    



['regulates', 'has_part']



In [5]:

    
print(eg_term.relationship['regulates'])









    



set([GOTerm('GO:0044772'):
  name:mitotic cell cycle phase transition
  relationship: 1 items
    part_of: 1 items
      GO:0000278	level-04	depth-04	mitotic cell cycle [biological_process] 
  level:5
  is_obsolete:False
  namespace:biological_process
  id:GO:0044772
  depth:5
  parents: 2 items
    GO:0044770	level-04	depth-04	cell cycle phase transition [biological_process] 
    GO:1903047	level-04	depth-04	mitotic cell cycle process [biological_process] 
  children: 4 items
    GO:0000086	level-06	depth-06	G2/M transition of mitotic cell cycle [biological_process] 
    GO:0000082	level-06	depth-06	G1/S transition of mitotic cell cycle [biological_process] 
    GO:0010458	level-06	depth-06	exit from mitosis [biological_process] 
    GO:0007091	level-06	depth-10	metaphase/anaphase transition of mitotic cell cycle [biological_process] 
  _parents: 2 items
    GO:0044770
    GO:1903047
  alt_ids: 0 items])

Example use case

One example use case for the relationship terms, would be to look for all functions which regulate pseudohyphal growth (GO:0007124). That is:

A pattern of cell growth that occurs in conditions of nitrogen limitation and abundant fermentable carbon source. Cells become elongated, switch to a unipolar budding pattern, remain physically attached to each other, and invade the growth substrate.

Source: https://www.ebi.ac.uk/QuickGO/GTerm?id=GO:0007124#term=info&info=1



In [6]:

    
term_of_interest = go['GO:0007124']

First, find the relationship types which contain "regulates":



In [7]:

    
regulates = frozenset([typedef 
                       for typedef in go.typedefs.keys() 
                       if 'regulates' in typedef])
print(regulates)









    



frozenset(['regulates', 'negatively_regulates', 'positively_regulates'])

Now, search through the terms in the tree for those with a relationship in this list and add them to a dictionary dependent on the type of regulation.



In [8]:

    
from collections import defaultdict

regulating_terms = defaultdict(list)

for t in go.values():
    if hasattr(t, 'relationship'):
        for typedef in regulates.intersection(t.relationship.keys()):
            if term_of_interest in t.relationship[typedef]:
                regulating_terms['{:s}d_by'.format(typedef[:-1])].append(t)

Now regulating_terms contains the GO terms which relate to regulating protein localisation to the nucleolus.



In [9]:

    
print('{:s} ({:s}) is:'.format(term_of_interest.name, term_of_interest.id))
for (k, v) in regulating_terms.items():
    print('\n  - {:s}:'.format(k))
    for t in v:
        print('    -- {:s} {:s}'.format(t.id, t.name))









    



pseudohyphal growth (GO:0007124) is:

  - negatively_regulated_by:
    -- GO:2000221 negative regulation of pseudohyphal growth
    -- GO:0100042 negative regulation of pseudohyphal growth by transcription from RNA polymerase II promoter
    -- GO:1900462 negative regulation of pseudohyphal growth by negative regulation of transcription from RNA polymerase II promoter

  - regulated_by:
    -- GO:2000220 regulation of pseudohyphal growth

  - positively_regulated_by:
    -- GO:2000222 positive regulation of pseudohyphal growth
    -- GO:0100041 positive regulation of pseudohyphal growth by transcription from RNA polymerase II promoter
    -- GO:1900461 positive regulation of pseudohyphal growth by positive regulation of transcription from RNA polymerase II promoter



In [ ]: