deFUME: Dynamic exploration of functional metagenomic sequencing data
© van der Helm et al. 2015
Received: 19 May 2015
Accepted: 15 July 2015
Published: 31 July 2015
Functional metagenomic selections represent a powerful technique that is widely applied for identification of novel genes from complex metagenomic sources. However, whereas hundreds to thousands of clones can be easily generated and sequenced over a few days of experiments, analyzing the data is time consuming and constitutes a major bottleneck for experimental researchers in the field.
Here we present the deFUME web server, an easy-to-use web-based interface for processing, annotation and visualization of functional metagenomics sequencing data, tailored to meet the requirements of non-bioinformaticians. The web-server integrates multiple analysis steps into one single workflow: read assembly, open reading frame prediction, and annotation with BLAST, InterPro and GO classifiers. Analysis results are visualized in an online dynamic web-interface.
The deFUME webserver provides a fast track from raw sequence to a comprehensive visual data overview that facilitates effortless inspection of gene function, clustering and distribution. The webserver is available at cbs.dtu.dk/services/deFUME/and the source code is distributed at github.com/EvdH0/deFUME.
KeywordsFunctional metagenomics Web services Visualization Sequence analysis Gene ontology
Functional selection represents a powerful technique for discovery and functionally validated annotation of genes [1, 2]. The technique relies on the expression of randomly cloned genomic or metagenomic DNA, typically as short (1–3 kb) fragments of DNA into an expression vector. The expression library is subsequently transformed into a suitable host where the desired functionality can be selected. DNA inserts from clones exhibiting the desired phenotype can be sequenced, enabling functional isolation of novel genes. The approach has been applied for identification of genes from complex metagenomic sources and examples include DNA polymerases , antibiotics resistance genes [1, 2, 4, 5], xenobiotic degradative enzymes , and more .
Most commonly, inserts of 1–3 kb are sequenced either by next generation sequencing [4, 8] or by conventional bi-directional Sanger sequencing, followed by base-calling, quality trimming of reads and assembly of reads into contigs. Finally data analysis is performed including BLAST searches and other functional annotations. Whereas hundreds to thousands of clones can be easily generated over a few days of experiments, analyzing the data is time consuming and constitutes a major bottleneck for experimental researchers in the field. To address this challenge we developed deFUME; an easy-to-use web server that automatically processes and annotates large amounts of sequencing data obtained through functional selections.
As input, deFUME accepts nucleotide sequences generated from next generation sequencing technologies as well as raw reads from Sanger sequencing. Sequences can be uploaded via the submission page either as preassembled projects, when using next generation sequencing data in. FASTA format, or as raw Sanger sequencing chromatograms (.ab1 format). For Sanger sequencing chromatograms, base-calling and assembly of the resulting reads is carried out with Phred  (default parameters) and Phrap  (using the non-default parameters -minscore 25, -trim_score 20, and -min_match 20 to reflect medium stringency cutoffs in order to ensure good quality assemblies of Sanger reads) respectively. Prior to assembly the vector sequence is optionally masked in the reads using Cross_match (using -minmatch 12 -minscore 20)  (Additional file 1: Figure S1).
Open reading frames (ORFs) are predicted using MetaGeneMark  (default parameters) from the resulting assembly (generated by Phrap or user input). The translated ORFs are aligned to the nr protein database using BLASTp (reporting only the 25 most significant hits with a minimal E value of 0.001)  and submitted to InterPro  (default parameters). The InterPro database contains signatures of known proteins families that can be queried to functionally characterize new sequences. To ensure using the most recent databases, InterPro is accessed using the simple object access protocol (SOAP) via InterProScan 5 .
Results and discussion
To demonstrate and test the deFUME web server, we analyzed Sanger sequencing data derived from a functional metagenomics selection for genes conferring tolerance of E. coli to high levels of lysine. Briefly, metagenomic DNA derived from cow fecal matter was mechanically sheared and subsequently cloned and expressed in E. coli. The resulting cell library was subjected to high lysine levels on Luria–Bertani agar plates (Additional file 1). The inserts of 80 individual colonies tolerant to high lysine levels were sequenced. The resulting 160 raw Sanger sequencing chromatograms were submitted as.ab1 files to the deFUME web server. In less than 2 h, all reads were trimmed and assembled, resulting in 69 unique contigs, 117 predicted ORFs, 134 GO annotations and 622 InterPro functional predictions. A screenshot of the output page is shown in Fig. 1, displaying one of the 69 inserts. In this particular insert, deFUME predicted three ORFs with a sequence coverage (calculated by dividing the protein sequence length of the predicted ORF by the length of the best BLASTp hit) of 12, 100 and 90%. Together with an E value of 4−52 this indicates that ORF 2 likely encodes the mechanism for lysine tolerance. The translated ORF is annotated by BLASTp as an “acetyltransferase GNAT family” protein catalyzing the acetylation of the nitrogen group of lysine, thus providing tolerance to high levels of lysine. The example demonstrates the ability of deFUME to accelerate the overall process of going from experimental raw data to functional annotation.
deFUME is the first web server to integrate all steps from sequencing read assembly to a comprehensive visual output for functional annotation of metagenomic insert libraries. Additionally, deFUME reduces the hands-on time required for analysis compared to packages like CLC Main Workbench (CLC Bio, Aarhus, Denmark) and Mobyle , where the user has to transfer the intermediate data from one tool to the other. Furthermore, the data integration and visualization of deFUME is substantially more advanced compared to the current state-of-art by providing interactive exploration of heterogeneous data.
Availability and requirements
Project name: deFUME.
Project home page: http://www.cbs.dtu.dk/services/deFUME.
Operating system(s): Platform independent.
Other requirements: Browser supporting HTML5.
License: Creative Commons BY 2.0
Any restrictions to use by non-academics: no.
EvdH, HJG, HMGH and MOAS conceived the study. EvdH and HMGH wrote the software. SM performed the experimental work. All authors read and approved the final manuscript.
The authors thank Dionísio S. Paiva for providing the cow fecal library and Hans H. Stærfeldt for advice on the web server implementation.
Compliance with ethical guidelines
Competing interests The authors declare that they have no competing interests.
Funding This work was supported by the European Union Seventh Framework Programme—ITN [FP7/2012/317058 to EvdH]; Novozymes A/S [HJG, HMGH]; the European Union Seventh Framework Programme (FP7-KBBE-2013-7-single-stage) under grant agreement no. 613745, Promys [MOAS, HJG, EvdH] and the Novo Nordisk Foundation.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Riesenfeld CS, Goodman RM, Handelsman J (2004) Uncultured soil bacteria are a reservoir of new antibiotic resistance genes. Environ Microbiol 6:981–989PubMedView ArticleGoogle Scholar
- Sommer MOA, Dantas G, Church GM (2009) Functional characterization of the antibiotic resistance reservoir in the human microflora. Science 325:1128–1131PubMedView ArticleGoogle Scholar
- Simon C, Herath J, Rockstroh S, Daniel R (2009) Rapid identification of genes encoding DNA polymerases by function-based screening of metagenomic libraries derived from glacial ice. Appl Environ Microbiol 75:2964–2968PubMed CentralPubMedView ArticleGoogle Scholar
- Forsberg KJ, Reyes A, Wang B, Selleck EM, Sommer MOA, Dantas G (2012) The shared antibiotic resistome of soil bacteria and human pathogens. Science 337(6098):1107–1111PubMed CentralPubMedView ArticleGoogle Scholar
- Forsberg KJ, Patel S, Gibson MK, Lauber CL, Knight R, Fierer N et al (2014) Bacterial phylogeny structures soil resistomes across habitats. Nature 509:612–616PubMed CentralPubMedView ArticleGoogle Scholar
- Ono A, Miyazaki R, Sota M, Ohtsubo Y, Nagata Y, Tsuda M (2007) Isolation and characterization of naphthalene-catabolic genes and plasmids from oil-contaminated soil by using two cultivation-independent approaches. Appl Microbiol Biotechnol 74:501–510PubMedView ArticleGoogle Scholar
- Uchiyama T, Miyazaki K (2009) Functional metagenomics for enzyme discovery: challenges to efficient screening. Curr Opin Biotechnol 20:616–622PubMedView ArticleGoogle Scholar
- Udikovic-Kolic N, Wichmann F, Broderick NA, Handelsman J (2014) Bloom of resident antibiotic-resistant bacteria in soil following manure fertilization. Proc Natl Acad Sci 111:15202–15207PubMed CentralPubMedView ArticleGoogle Scholar
- Ewing B, Green P (1998) Base-calling of automated sequencer traces using Phred. II. Error probabilities. Genome Res 8:175–185PubMedView ArticleGoogle Scholar
- Gordon D, Abajian C, Green P (1998) Consed: a graphical tool for sequence finishing. Genome Res 8:195–202PubMedView ArticleGoogle Scholar
- Zhu W, Lomsadze A, Borodovsky M (2010) Ab initio gene identification in metagenomic sequences. Nucleic Acids Res 38:e132PubMed CentralPubMedView ArticleGoogle Scholar
- Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K et al (2009) BLAST+ : architecture and applications. BMC Bioinform 10:421View ArticleGoogle Scholar
- Hunter S, Jones P, Mitchell A, Apweiler R, Attwood TK, Bateman A et al (2012) InterPro in 2011: New developments in the family and domain prediction database. Nucleic Acids Res 40(November 2011):306–312View ArticleGoogle Scholar
- Jones P, Binns D, Chang HY, Fraser M, Li W, McAnulla C et al (2014) InterProScan 5: Genome-scale protein function classification. Bioinformatics 30:1236–1240PubMed CentralPubMedView ArticleGoogle Scholar
- Néron B, Ménager H, Maufrais C, Joly N, Maupetit J, Letort S et al (2009) Mobyle: a new full web bioinformatics framework. Bioinformatics 25:3005–3011PubMed CentralPubMedView ArticleGoogle Scholar