GO FEAT: a rapid web-based functional annotation tool for genomic and transcriptomic data

<jats:title>Abstract</jats:title><jats:p>Downstream analysis of genomic and transcriptomic sequence data is often executed by functional annotation that can be performed by various bioinformatics tools and biological databases. However, a full fast integrated tool is not available for such analysis. Besides, the current available software is not able to produce analytic lists of annotations and graphs to help users in evaluating the output results. Therefore, we present the Gene Ontology Functional Enrichment Annotation Tool (GO FEAT), a free web platform for functional annotation and enrichment of genomic and transcriptomic data based on sequence homology search. The analysis can be customized and visualized as per users’ needs and specifications. GO FEAT is freely available at <jats:ext-link xmlns:xlink="http://www.w3.org/1999/xlink" ext-link-type="uri" xlink:href="http://computationalbiology.ufpa.br/gofeat/">http://computationalbiology.ufpa.br/gofeat/</jats:ext-link> and its source code is hosted at <jats:ext-link xmlns:xlink="http://www.w3.org/1999/xlink" ext-link-type="uri" xlink:href="https://github.com/fabriciopa/gofeat">https://github.com/fabriciopa/gofeat</jats:ext-link>.</jats:p>

section_content: UniProt public API and SEED database via SEED public API. After the integration, the results are processed and displayed in graphs, charts, and tables to simplify the analysis. 

Since the EBI servers restrict the number of request to 30 at time, a queue control parameter was developed to optimize the server's resources. For projects with 100 or less sequences, resources are allocated dynamically for maximum of 10 users simultaneously (3 requests for each project). If resources are available, the projects can receive more than 3 requests. Projects with more than 100 sequences are put in a queue for local alignment using DIAMOND that process batches of 500 sequences at a time. This allows the server's resources usage to be optimized and more sequences can be processed at the same time. 

To compare the results from GO FEAT with other tools, we performed the functional annotation in six different scenarios: a random sequence with 500 bp from Escherichia coli; the full genome of Escherichia coli K-12 MG1655 (4140 CDS and average size of 321 bp) [RefSeq NC_000913. 3]; the full genome of Drosophila melanogaster BDGP6 (30482 CDS and CDS average size of 668 bp) [Assembly GCA_000001215. 4]; the full genome of Nostoc sp. PCC 7107 (5237 CDS and CDS average size of 330 bp) [RefSeq NC_019676. 1]; the transcriptomic data from E. coli response to five different perturbations (4092 CDS and CDS average size of 326 bp); 19 and the transcriptomic data from M. tuberculosis response to macrophages (4076 CDS and CDS average size of 332 bp) 20. The results of this comparison are shown in the next section.

section_content: Interface. GO FEAT was developed to be executed in any modern internet browser. Also, it has a clean and easy-to-use graphic interface. It's not required any kind of installation of any tool or software and users can execute projects without previous registration.
section_name: Results

section_content: formed by registered users. In the project manager, it's possible to check the project's progress, export data to several formats and share projects to other people users to avoid running the same project multiple times. 

Reports, charts and graphs. GO FEAT allows different ways for result visualization: spreadsheet reports present sequences in tables corresponding to its Blast result, which are integrated to several databases, to perform searches and export results; it's also possible to view the results in graphs and charts, which are divided by molecular function, cellular component and biological process. On each one, it is possible to view all GO terms of each category together with the sequences identification. Finally, the user can view the GO terms with its acyclic graph, downloaded through the Quick GO API. 

Benchmarking. For a 500 bp random sequence chosen from Escherichia coli's genome, GO FEAT takes around 4 minutes for full functional annotation and enrichment while Blast2GO takes around 14 minutes for the same sequence. Direct Blast to NCBI website takes around 2 minutes, however, the mapping between the blast result and other databases are not automatically made. At UniProt, the function annotation and enrichment takes around 2 minutes. Since NCBI's Blast does not perform a full functional annotation and UniProt website has limitations regarding the number of sequences, they will not be included in further analysis. For complete genomes of model organisms, GO FEAT needs around 5 hours for Escherichia coli and 30 hours for Drosophila melanogaster. For transcriptomic data, 5 hours were required to peform the functional annotation described in the Jozefczuk's paper and 5 hours to perform the functional annotation described in Rohde's paper. Blast2GO was unable to perform the full annotation and enrichment of any complete genome or transcriptomic data analysis in less than 10 days. For non-model organism such as Nostoc sp. PCC 7107, around 4 hours is required to finish the processing in GO FEAT. The time varies depending on server loads of the remote APIs. At full load, around 11 hours was necessary to process 10 projects, each one with 1000 different sequences from Drosophila melanogaster and CDS average size of 603 bp. Regarding functionalities, GO FEAT presents useful features in comparison to other functional annotation tools (Table 1 ) and rapidly process the input sequence and generates the results. 

Limitations. GO FEAT was developed to perform functional annotations on previously predicted genes, coding DNA sequences (CDS), open reading frames (ORF) or transcripts. Thus, large sequences, such as full genomes or contigs, are not suitable to be used as inputs in GO FEAT due to size limitation of alignment softwares.
section_name: Project manager. GO FEAT provides a project manager to facilitate the categorization of each analysis per-

section_content: Functional characterization of biological sequences is a required step in the analysis of biological data. GO FEAT is an annotation platform integrated with several databases which can be used for different datasets, such as: coding sequences identified after gene prediction and sequences produced after new sequence assembly of next-generation sequencing data. The user can share results with collaborators through graphic interface and can export the results to many formats. Since the tool uses API to access various databases, the annotations are based on most recent and updated data from those databases. 

We are committed to maintain GO FEAT for at least 2 years and expect to improve its performance as our computational infrastructure grows. For future works, we plan on adding a prediction step before the functional annotation so users can input large sequences, exporting the predicted sequences.
section_name: Conclusions

section_name: Acknowledgements

section_content: This work has been supported by the CNPq ( Conselho Nacional de Pesquisa Científica ) grant # 421528/2016-8, CAPES ( Coordenação de Aperfeiçoamento de Pessoal de Nível Superior) and PROPESP/UFPA (Pró-Reitoria de Pesquisa e Pós-Graduação/Universidade Federal do Pará ).
section_name: Acknowledgements

section_name: Author Contributions

section_content: Rommel Thiago Juca Ramos and Fabricio Almeida Araujo conceived the program's idea and developed it. Debmalya Barh, Artur Silva and Luis Guimarães evaluated the biological informations and defined the databases to be integrated. All authors reviewed the manuscript.
section_name: Author Contributions

section_content: Competing Interests: The authors declare that they have no competing interests. 

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
section_name: Additional Information