adding go go workshop 3-6 august 2010. 1. goanna results and goanna2ga 2. gene association files 3....

22
Adding GO GO Workshop 3-6 August 2010

Upload: duane-carroll

Post on 13-Jan-2016

213 views

Category:

Documents


0 download

TRANSCRIPT

Page 1: Adding GO GO Workshop 3-6 August 2010. 1. GOanna results and GOanna2ga 2. gene association files 3. getting GO for your dataset 4. adding more GO (introduction)

Adding GO

GO Workshop3-6 August 2010

Page 2: Adding GO GO Workshop 3-6 August 2010. 1. GOanna results and GOanna2ga 2. gene association files 3. getting GO for your dataset 4. adding more GO (introduction)

1. GOanna results and GOanna2ga2. gene association files3. getting GO for your dataset4. adding more GO (introduction)5. requesting GO

Page 3: Adding GO GO Workshop 3-6 August 2010. 1. GOanna results and GOanna2ga 2. gene association files 3. getting GO for your dataset 4. adding more GO (introduction)

only returns existing GO only accepts limited accession types

GOanna does a Blast search against existing GO annotated products.

allows you to quickly transfer GO to gene products where they have similar sequences

accepts fasta files

Page 4: Adding GO GO Workshop 3-6 August 2010. 1. GOanna results and GOanna2ga 2. gene association files 3. getting GO for your dataset 4. adding more GO (introduction)

Incorrect email address – you will not receive your results!

Contact AgBase if you have not received results after 24-48h.

Page 5: Adding GO GO Workshop 3-6 August 2010. 1. GOanna results and GOanna2ga 2. gene association files 3. getting GO for your dataset 4. adding more GO (introduction)
Page 6: Adding GO GO Workshop 3-6 August 2010. 1. GOanna results and GOanna2ga 2. gene association files 3. getting GO for your dataset 4. adding more GO (introduction)
Page 7: Adding GO GO Workshop 3-6 August 2010. 1. GOanna results and GOanna2ga 2. gene association files 3. getting GO for your dataset 4. adding more GO (introduction)
Page 8: Adding GO GO Workshop 3-6 August 2010. 1. GOanna results and GOanna2ga 2. gene association files 3. getting GO for your dataset 4. adding more GO (introduction)

Standard BLAST parameters

Page 9: Adding GO GO Workshop 3-6 August 2010. 1. GOanna results and GOanna2ga 2. gene association files 3. getting GO for your dataset 4. adding more GO (introduction)

default setting

Page 10: Adding GO GO Workshop 3-6 August 2010. 1. GOanna results and GOanna2ga 2. gene association files 3. getting GO for your dataset 4. adding more GO (introduction)
Page 11: Adding GO GO Workshop 3-6 August 2010. 1. GOanna results and GOanna2ga 2. gene association files 3. getting GO for your dataset 4. adding more GO (introduction)

GOanna ResultsIf you enter an incorrect email address – you will not receive your results!

Contact AgBase if you have not received results after 24-48h.

Page 12: Adding GO GO Workshop 3-6 August 2010. 1. GOanna results and GOanna2ga 2. gene association files 3. getting GO for your dataset 4. adding more GO (introduction)

query IDs are hyperlinked to BLAST data(files must be in the same directory)

Page 13: Adding GO GO Workshop 3-6 August 2010. 1. GOanna results and GOanna2ga 2. gene association files 3. getting GO for your dataset 4. adding more GO (introduction)

*WHAT IS A GOOD ALIGNMENT?

1. Manually inspect alignments and delete any lines where there is not a good alignment*.2. Add this additional annotation to the annotations from GORetriever.

Page 14: Adding GO GO Workshop 3-6 August 2010. 1. GOanna results and GOanna2ga 2. gene association files 3. getting GO for your dataset 4. adding more GO (introduction)

GOanna2gaNew to AgBase: an online script to convert your GOanna file to a gene association file format.

• add manually checked GOanna annotations to a GORetriever file

Page 15: Adding GO GO Workshop 3-6 August 2010. 1. GOanna results and GOanna2ga 2. gene association files 3. getting GO for your dataset 4. adding more GO (introduction)

Tutorial 4: adding GO using GOannaFor this example, we will use a list of NCBI Cassava proteins that have no GO annotation. 4.1 Example 1: Cassava

4.1.1 GOanna 4.1.2 GOanna2ga

4.2 Example 2: Working on your own data Add GO to a sample from your own data set. Alternatively, use GOanna to add GO to the file

of accessions from tutorial 2 that had no GO. We will continue to work on this tomorrow.

Page 16: Adding GO GO Workshop 3-6 August 2010. 1. GOanna results and GOanna2ga 2. gene association files 3. getting GO for your dataset 4. adding more GO (introduction)

Requesting GO

Page 17: Adding GO GO Workshop 3-6 August 2010. 1. GOanna results and GOanna2ga 2. gene association files 3. getting GO for your dataset 4. adding more GO (introduction)

Some limitations of GOanna:

• BLAST analysis is slow – results emailed• limit to 1,000 sequences/submission• limit to 3 jobs submitted/user at one time• limit of <3 databases selected for any one job

How do I do to get GO for my 50,000 RNA-Seq dataset?

50 x GOanna submissions + manual interpretation of results – impractical and slow!!

ALTERNATIVELY: Contact AgBase we use internal GO annotation pipelines GO can be kept private and released after

publication.

Page 18: Adding GO GO Workshop 3-6 August 2010. 1. GOanna results and GOanna2ga 2. gene association files 3. getting GO for your dataset 4. adding more GO (introduction)
Page 19: Adding GO GO Workshop 3-6 August 2010. 1. GOanna results and GOanna2ga 2. gene association files 3. getting GO for your dataset 4. adding more GO (introduction)
Page 20: Adding GO GO Workshop 3-6 August 2010. 1. GOanna results and GOanna2ga 2. gene association files 3. getting GO for your dataset 4. adding more GO (introduction)
Page 21: Adding GO GO Workshop 3-6 August 2010. 1. GOanna results and GOanna2ga 2. gene association files 3. getting GO for your dataset 4. adding more GO (introduction)
Page 22: Adding GO GO Workshop 3-6 August 2010. 1. GOanna results and GOanna2ga 2. gene association files 3. getting GO for your dataset 4. adding more GO (introduction)

Tutorial 4: adding GO using GOannaFor this example, we will use a list of NCBI Cassava proteins that have no GO annotation. 4.1 Example 1: Cassava

4.1.1 GOanna 4.1.2 GOanna2ga

4.2 Example 2: Working on your own data Add GO to a sample from your own data set. Alternatively, use GOanna to add GO to the file

of accessions from tutorial 2 that had no GO. We will continue to work on this tomorrow.