Gene Ontology Mapping
Mapping is the process of retrieving GO terms associated with the Hits obtained by the BLAST search. OmicsBox performs four different mappings steps:
BLAST result accessions are used to retrieve gene names or symbols making use of two mapping files provided by the NCBI (gene_info, gene2accession). Identified gene names are then searched in the species-specific entries of the gene-product table of the GO database.
GeneBank identifiers (gi), the primary blast Hit ids, are used to retrieve UniProt IDs making use of a mapping file from PIR (Non-redundant Reference Protein Database) including PSD, UniProt, Swiss-Prot, TrEMBL, RefSeq, GenPept and PDB.
Accessions are searched directly in the dbxref table of the GO database.
BLAST result accessions are searched directly in the gene-product table of the GO database.
Gotz S., Garcia-Gomez JM., Terol J., Williams TD., Nagaraj SH., Nueda MJ., Robles M., Talon M., Dopazo J. and Conesa A. (2008). High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic acids research, 36(10), 3420-35.
Figure 1: Mapping options
Run GO Mapping
The GO Mapping functionality can be found in the Side Panel → GO Mapping after loading a fasta file.
The GO Mapping options are:
GO Mapping. Mapping will start using online mapping database.
Remove Mapping. Delete Mapping results for the selected sequences.
Charts. Generate Mapping statistics charts.
The mapping step needs protein ids to run. Make sure you ran blast against a protein database.
blastx - if one has nucleotide sequences
blastp - if one has protein sequences
Show Mapping Results
For each sequence, it is possible to see the mapping results individually.
Show Mapping Results. A new table will be displayed (see figure 3). The resulting table shows the GO mapping results for a particular sequence. See Table section to manipulate/extract the results from this table.
Show GO Descriptions. GO ID, description, type, and definition are given for all GO terms associated with the selected sequence. The GO ID is linked to the AmiGO browser at the Gene Ontology site while the show option displays the DAG representation of the GO term.
Annotate Sequence. This function allows changing annotation parameters for the selected sequence and re-running automatic annotation.
Change Annotation and Description. This function edits the annotation of the selected and allows typing and deleting of annotation or sequence description. A manual annotation check-box (see figure 5 in Gene Ontology Annotation section) is available for marking sequences with manual annotation. The sequence will get the pink label on the Main Sequence Table.
Make Graph of GO-Mapping-Results with Annotation Score. Displays a DAG with all GO terms related to one sequence. Shows all the GOs from the mapping step as well as final annotations (highlighted). The wizard (figure 4 allows filtering the hits which will be taken into account (see Gene Ontology Graphs section for more details about visualization in OmicsBox)
Hit Filter. Nodes can be filtered out by a number of hits: only nodes with more than a given number of BLAST-Hits will be shown in the graph.
HSP-Hit Coverage CutOff: Includes only those hits which are overage with the HSP for a given percentage.
For Mapping statistic charts see the Charts and Statistics page of this user manual.
Figure 2: Show Mapping Results
Figure 3: Mapping Results for sequence C02006A02
Figure 4: Single Graph Drawing Configuration
Three different charts are available to summarise the mapping step:
GO Mapping Distribution: This shows the distribution of the number of Gene Ontology candidate terms assigned to each sequence during the GO Mapping step.
EC Distribution for Sequences: This chart shows the distribution of GO evidence codes for the functional terms obtained during the mapping step. It gives an idea about how many annotations derive from automatic/computational annotations or manually curated ones.
EC Distribution for Blast Hits: Evidence Codes associated with the obtained GO pool.
Figure 5: Mapping Statistics
Figure6: GO Mapping Distribution
Export Mapping Results
A tab separator text file can be exported with the corresponding mapping results (Side Panel > Export > Export Mapping Results).