The metagenomics RAST server – a public resource for the automatic phylogenetic and functional analysis of metagenomes
© Meyer et al; licensee BioMed Central Ltd. 2008
Received: 14 February 2008
Accepted: 19 September 2008
Published: 19 September 2008
Random community genomes (metagenomes) are now commonly used to study microbes in different environments. Over the past few years, the major challenge associated with metagenomics shifted from generating to analyzing sequences. High-throughput, low-cost next-generation sequencing has provided access to metagenomics to a wide range of researchers.
A high-throughput pipeline has been constructed to provide high-performance computing to all researchers interested in using metagenomics. The pipeline produces automated functional assignments of sequences in the metagenome by comparing both protein and nucleotide databases. Phylogenetic and functional summaries of the metagenomes are generated, and tools for comparative metagenomics are incorporated into the standard views. User access is controlled to ensure data privacy, but the collaborative environment underpinning the service provides a framework for sharing datasets between multiple users. In the metagenomics RAST, all users retain full control of their data, and everything is available for download in a variety of formats.
The open-source metagenomics RAST service provides a new paradigm for the annotation and analysis of metagenomes. With built-in support for multiple data sources and a back end that houses abstract data types, the metagenomics RAST is stable, extensible, and freely available to all researchers. This service has removed one of the primary bottlenecks in metagenome sequence analysis – the availability of high-performance computing for annotating the data.
The genomic revolution of the 1990s has yielded almost a thousand sequenced microbial genomes. More recently, the explosion of random community genomics, or metagenomics, where DNA is sequenced directly from environmental samples has provided insights into microbial communities. Currently, two approaches to sequencing metagenome samples are commonly used. In the traditional approach, DNA is cloned into BACs, or small plasmids, and dideoxy chain termination sequencing ("Sanger sequencing") is used to determine the sequences [1, 2]. In the alternative approach, DNA is sequenced without cloning, using one of the so-called next-generation sequencing techniques, usually pyrosequencing. Both approaches have advantages and disadvantages. For example, Sanger sequencing generates longer sequence reads but has inherent biases due to the cloning. Pyrosequencing has much higher throughput and a lower error rate per base sequenced compared to Sanger sequencing, but those errors are biased toward certain mistakes .
Regardless of the sequencing approach used to generate the data, the first steps in analysis of any metagenome involve comparing those sequences to known sequence databases. This computationally intensive task provides the basic data types for many subsequent analyses, including phylogenetic comparisons, functional annotations, binning of sequences, phylogenomic profiling, and metabolic reconstructions.
Here we describe the development of a freely available, fully automated open source system for processing metagenome sequence data to generate these basic elements. A public implementation of this system has been provided for all researchers to analyze their metagenomes. Our service, the metagenomics RAST server (mg-RAST for short), is available over the web to all researchers, and access is not limited to specific groups or data types. Almost 500 metagenomes have been processed through the beta version of the pipeline so far.
Implementation and results
User Registration and Management
The user registration serves two functions: to limit access to each data set to the user and their colleagues and to secure a valid email address in case correspondence is required, for example if a data-processing problem occurs. Once logged in, users can view their own metagenomes, those to which the owner has granted them rights, and the default set of publicly available metagenomes. The system supports delegation of authorization so that users can allow others to access one or more of their metagenomes. In addition, data owners can release their metagenomes to the public at any point, allowing all users of the system to view their data.
The pipeline accepts data in a number of formats: 454 reads may be uploaded directly in the format delivered by 454 , and fasta files typical of Sanger-sequences and used by other platforms may also be uploaded. The pipeline will also accept assembled sequences in fasta format. Sequence data may be compressed by one of several common computer programs to speed upload.
Users may choose to upload raw unassembled reads or assembled contigs. As discussed below, each approach has advantages and disadvantages. Users with a limited number of larger contigs, where the average contig length exceeds 40 kb, should consider using the RAST server for the analysis of complete Bacterial and Archaeal genomes .
The Genomics Standards Consortium has proposed a minimal set of data, called the Minimum Information about a Genome Sequence (MIGS) , that should be collected with every metagenome sequence. Although this is an evolving standard, the metagenomics-RAST server is MIGS-compliant. Metadata, accessory data about the metagenome (e.g., date and location where the sample was collected), is requested from the user at the time of sequence submission. This data is stored with the user's data and can be provided to the GSC genome catalogue, and other archives, when the sequence data is ready for public release.
Implementation and Core Analyses
The pipeline is implemented in Perl by using a number of open source components, including the SEED framework , NCBI BLAST , SQLite, and Sun Grid Engine  as components. The system also uses the publicly available SEED subsystems, SEED nr, and FIGfam protein families (see http://www.theSEED.org).
The distinct steps are implemented to provide a flexible, extensible processing pipeline. The steps incrementally add data to a self-contained "job directory" that contains all job-relevant data in flat file and SQLite  format. Relational database technology is used to efficiently provide a mapping of sequences in a metagenome to both organisms and metabolic functions and at the same time allow the user to change the parameters for the underlying sequence matches. The user interface enables the download of the user's job directories, and a future version of the software will allow uploading of user-created directories into the server.
In the second step, the sequences are screened for potential protein encoding genes (PEGs) via a BLASTX  search against the SEED comprehensive nonredundant database sourced from the INSDC databases, sequencing centers, and other sources . An expect value (E) cut-off of 0.01 is used to pick up potentially coding elements. (This was chosen empirically to increase the number of potentially coding elements while not being overwhelming for data analysis.) In parallel with the BLASTX searches, the sequence data is compared to all accessory databases by using the appropriate algorithms and significance selection criteria. These databases include several rDNA databases, including GREENGENES , RDP-II , and the European 16S RNA database , and boutique databases such as the chloroplast database, mitochondrial database, and ACLAME database of mobile elements . The search criteria are specific for each database. For example, screens for ribosomal RNA genes are performed by using BLASTN against the rDNA databases, but much more stringent selection criteria are used to identify candidate RNA genes than for identifying protein-encoding genes (by default, the similarity must exceed 50 bp in length and have an expect value less than 1 × 10-5).
In the third step, these matches to external databases are used to compute the derived data. First, a phylogenomic reconstruction of the sample is computed by using both the phylogenetic information contained in the SEED nr database and the similarities to the ribosomal RNA database. Functional classifications of the PEGs are computed by projecting against SEED FIGfams  and subsystems based on these similarity searches . These functional assignments become the raw input to an automatically generated initial metabolic reconstruction of the sample, providing suggestions for metabolic fluxes and flows, reactions, and enzymes.
The abundance of comparative metagenomics tools is central to the utility of the mg-RAST platform. Various tools have been built into the framework, allowing users to compare their data against other metagenomes or complete genomes taken from the SEED  environment. The subsystems heat map and the taxonomic heat map provide comparative metagenomics summaries that encapsulate the differences between samples.
The subsystem comparison tools identify the number of pegs in each metagenome that are connected to a subsystem via protein level similarity. Based on these connections, each subsystem present in a sample is scored by counting the number of sequences that are similar to a protein in each subsystem. This score is divided by the total number of sequences from the sample that are similar to any protein in a subsystem, to give a fraction of sequences in subsystems that are in a given subsystem. This approach allows comparisons between samples that have different numbers of sequences. Since the fractions tend to be small (a few sequences hit each subsystem, but there are now over 600 subsystems in the SEED), the scores can be factored for display purposes. Furthermore, a nonquantitative approach is provided to group the subsystem scores, emphasizing those subsystems that are most different between the samples. Moreover, the display can be limited to specific areas of metabolism, or other subsystem groups, as desired by the user.
The taxonomic heat map works in an analogous fashion but highlights the different taxonomic profiles in each sample, as determined by the phylogenetic or phylogenomic approaches selected by the end user (e.g., 16S comparisons, phylogenomics from BLAST results). Again, samples may be grouped in a nonquantitative fashion to rapidly highlight particular phylogenetic groups that predominate in different samples.
A completely new public metagenome annotation system has been developed and released. The process is the result of several years of planning and engineering. Designed to leverage the SEED microbial genome annotation platform, the mg-RAST platform provides seamless integration of metagenome data, microbial genomics, and manually curated annotations. Each metagenome project has its own requirements for stringency, datasets to be analyzed, and output format for results. The metagenomics SEED pipeline was designed to allow alterations to the parameters for the sequence matches underlying both the phylogenetic and metabolic reconstructions to restrict matches. It has been built by using an extensible format allowing the integration of new datasets and algorithms without a need for recomputation of existing results.
The mg-RAST service handles both assembled and unassembled data. Each approach has advantages that should be considered when comparing metagenomes. For example, if one is carrying out comparative metagenomics or if statistics are being used to compare samples [18, 19], the sequences cannot be assembled, since the assembly process loses the frequency information critical for determining differences between samples. In contrast, assembled sequences tend to be longer and therefore more likely to accurately identify gene function or phylogenetic source from binning .
The analytical methods integrated into the pipeline provide core annotation and analysis tools to compare and contrast a diverse set of metagenomes [21–24]. The approach underlying the subsystems-based functional analysis of metagenomes has been validated with 90 different samples from nine major biomes. The analysis demonstrated that the biomes could clearly be separated by their functional composition . All of the metagenomes present in that study are included in the publicly available datasets visible in the mg-RAST server.
Although the service contains core functionality for the annotation and analysis of metagenomes, many of the techniques traditionally used for genome analysis (e.g., approaches for the prediction of coding sequences) either do not work with metagenomes or show a significant performance degradation . Many of the differences between complete genome annotation and metagenome annotation are reminiscent of those encountered previously with the analysis of expressed sequenced tags . Therefore, new analytical methods are needed to fully understand metagenomics data. The most obvious problem is with the large number of unknown sequences in any sample. Depending on the specific sample processed, as few as 10% of the sequences or as many as 98% of the sequences may have no known similarity to anything in the database . We and others are developing new binning, clustering, and coding region prediction tools to handle these unknown sequences, and effective tools will be incorporated into the pipeline when available. Another problem is that the rapid pace with which sequence data is being generated outpaces increases in computational speed, and therefore improvements in common search algorithms are required to ensure that sequence space can be accurately and efficiently searched. A third problem, common to all annotation platforms, is that metabolic reconstructions and analyses are dependent on the underlying quality of the data. The SEED has the most consistent and accurate microbial genome annotations of any publicly available source because of the subsystems approach to annotation. However, the SEED subsystems are necessarily focused on core metabolism and pathogenesis of a select few organisms. Comprehensive subsystem coverage of secondary metabolism, and especially of metabolism specific to diverse environments, is required to truly comprehend those data sets.
We have provided a free, public resource for the analysis of metagenome sequence data. Our service does not require a specific type of sequence data and has no requirement for release or control of the data. All sequence data remains protected by a password mechanism and is visible only to permitted users. This metagenomics annotation pipeline was specifically developed to handle pyrosequencing data and accommodate some of the nuances associated with that data. However, the tools and approaches we have developed are applicable and available for any metagenome project, regardless of sequence type. This service has removed one of the primary bottlenecks in metagenome sequence analysis – the availability of high-performance computing for annotating the data.
The service is available to all users after a simple registration process. In addition to being available through the integrated SEED-Viewer  interface, all results are available for download in variety of formats, including GFF3, GenBank, and flat text formats (e.g., tab-separated text for use in spreadsheets). The server is made available on a best-effort basis, and all underlying data and software are open source (please see ftp://ftp.theseed.org).
Plans include the development of novel tools to allow systematic data mining in the samples and improved support for in-depth analysis of 16S-based metagenome data sets.
Rapid Annotation using Subsystems Technology
RAST for metagenomes.
Part of this project has been funded with federal funds from the National Institute of Allergy and Infectious Diseases, National Institutes of Health, Department of Health and Human Services, under Contract No. HHSN266200400042C. Argonne National Laboratory's work was supported under U.S. Department of Energy contract DE-AC02-06CH11357. We thank Forest Rohwer and all members of his group for helpful advice while developing the pipeline. We thank Veronika Vonstein and Ross Overbeek and the members of their group for advice and comments in using subsystem technology.
The submitted manuscript has been created by UChicago Argonne, LLC, operator of Argonne National Laboratory ("Argonne"). Argonne, a U.S. Department of Energy Office of Science laboratory, is operated under Contract No. DE-AC02-06CH11357. The U.S. Government retains for itself, and others acting on its behalf, a paid-up, nonexclusive, irrevocable worldwide license in said article to reproduce, prepare derivative works, distribute copies to the public, and perform publicly and display publicly, by or on behalf of the Government.
- Venter JC, Remington K, Heidelberg JF, Halpern AL, Rusch D, Eisen JA, Wu D, Paulsen I, Nelson KE, Nelson W, et al.: Environmental genome shotgun sequencing of the Sargasso Sea. Science 2004, 304(5667):66–74. 10.1126/science.1093857View ArticlePubMedGoogle Scholar
- Tyson GW, Chapman J, Hugenholtz P, Allen EE, Ram RJ, Richardson PM, Solovyev VV, Rubin EM, Rokhsar DS, Banfield JF: Community structure and metabolism through reconstruction of microbial genomes from the environment. Nature 2004, 428(6978):37–43. 10.1038/nature02340View ArticlePubMedGoogle Scholar
- Huse SM, Huber JA, Morrison HG, Sogin ML, Welch DM: Accuracy and quality of massively parallel DNA pyrosequencing. Genome biology 2007, 8(7):R143. 10.1186/gb-2007-8-7-r143PubMed CentralView ArticlePubMedGoogle Scholar
- Overbeek R, Begley T, Butler RM, Choudhuri JV, Diaz N, Chuang H-Y, Cohoon M, de Crécy-Lagard V, Disz T, Edwards R, et al.: The Subsystems Approach to Genome Annotation and its Use in the Project to Annotate 1000 Genomes. Nucleic Acids Res 2005., 33(17):Google Scholar
- McNeil LK, Reich C, Aziz RK, Bartels D, Cohoon M, Disz T, Edwards RA, Gerdes S, Hwang K, Kubal M, et al.: The National Microbial Pathogen Database Resource (NMPDR): a genomics platform based on subsystem annotation. Nucleic Acids Res 2007, (35 Database):D347–353. 10.1093/nar/gkl947Google Scholar
- Margulies M, Egholm M, Altman WE, Attiya S, Bader JS, Bemben LA, Berka J, Braverman MS, Chen YJ, Chen Z, et al.: Genome sequencing in microfabricated high-density picolitre reactors. Nature 2005, 437(7057):376–380.PubMed CentralPubMedGoogle Scholar
- Aziz RK, Bartels D, Best AA, DeJongh M, Disz T, Edwards RA, Formsma K, Gerdes S, Glass EM, Kubal M, et al.: The RAST Server: Rapid Annotations using Subsystems Technology. BMC Genomics 2008, in press.Google Scholar
- Field D, Morrison N, Selengut J, Sterk P: Meeting report: eGenomics: Cataloguing our Complete Genome Collection II. Omics 2006, 10(2):100–104. 10.1089/omi.2006.10.100View ArticlePubMedGoogle Scholar
- Altschul SFl, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 1997, 25: 3389–3402. 10.1093/nar/25.17.3389PubMed CentralView ArticlePubMedGoogle Scholar
- gridengine – Project home[http://gridengine.sunsource.net/]
- Jarvie T: personal communication. 2006.Google Scholar
- DeSantis TZ, Hugenholtz P, Larsen N, Rojas M, Brodie EL, Keller K, Huber T, Dalevil D, Hu P, Andersen GL: Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl Environ Microbiol 2006, 72(7):5069–5072. 10.1128/AEM.03006-05PubMed CentralView ArticlePubMedGoogle Scholar
- Cole JR, Chai B, Farris RJ, Wang Q, Kulam-Syed-Mohideen AS, McGarrell DM, Bandela AM, Cardenas E, Garrity GM, Tiedje JM: The ribosomal database project (RDP-II): introducing myRDP space and quality controlled public data. Nucleic Acids Res 2007, (35 Database):D169–172. 10.1093/nar/gkl889Google Scholar
- Wuyts J, Peer Y, Winkelmans T, De Wachter R: The European database on small subunit ribosomal RNA. Nucleic Acids Res 2002, 30(1):183–185. 10.1093/nar/30.1.183PubMed CentralView ArticlePubMedGoogle Scholar
- Leplae R, Hebrant A, Wodak SJ, Toussaint A: ACLAME: a CLAssification of Mobile genetic Elements. Nucleic Acids Res 2004, (32 Database):D45–49. 10.1093/nar/gkh084Google Scholar
- Meyer F, Overbeek R, Rodriquez A: FIGfams – Yet another protein family collection. 2008, in press.Google Scholar
- Tringe SG, von Mering C, Kobayashi A, Salamov AA, Chen K, Chang HW, Podar M, Short JM, Mathur EJ, Detter JC, et al.: Comparative metagenomics of microbial communities. Science 2005, 308(5721):554–557. 10.1126/science.1107851View ArticlePubMedGoogle Scholar
- Rodriguez-Brito B, Rohwer F, Edwards RA: An application of statistics to comparative metagenomics. BMC bioinformatics [electronic resource] 2006, 7: 162. 10.1186/1471-2105-7-162View ArticleGoogle Scholar
- McHardy AC, Martin HG, Tsirigos A, Hugenholtz P, Rigoutsos I: Accurate phylogenetic classification of variable-length DNA fragments. Nature methods 2007, 4(1):63–72. 10.1038/nmeth976View ArticlePubMedGoogle Scholar
- Edwards RA, Rodriguez-Brito B, Wegley L, Haynes M, Breitbart M, Peterson DM, Saar MO, Alexander S, Alexander EC Jr, Rohwer F: Using pyrosequencing to shed light on deep mine microbial ecology. BMC Genomics 2006, 7: 57. 10.1186/1471-2164-7-57PubMed CentralView ArticlePubMedGoogle Scholar
- Fierer N, Breitbart M, Nulton J, Salamon P, Lozupone C, Jones R, Robeson M, Edwards RA, Felts B, Rayhawk S, et al.: Metagenomic and small-subunit rRNA analyses reveal the genetic diversity of bacteria archaea, fungi, and viruses in soil. Appl Environ Microbiol 2007, 73(21):7059–7066. 10.1128/AEM.00358-07PubMed CentralView ArticlePubMedGoogle Scholar
- Wegley L, Edwards R, Rodriguez-Brito B, Liu H, Rohwer F: Metagenomic analysis of the microbial community associated with the coral Porites astreoides. Environmental microbiology 2007, 9(11):2707–2719. 10.1111/j.1462-2920.2007.01383.xView ArticlePubMedGoogle Scholar
- Mou XSS, Edwards RA, Hodson RE, Moran MA: Generalist Species Dominate Bacterial Carbon Processing in the Coastal Ocean. Nature 2008.Google Scholar
- Dinsdale EA, Edwards RA, Hall D, Angly F, Breitbart M, Brulc JM, Furlan M, Desnues C, Haynes M, Li LM, et al.: Functional Metagenomic Profiling of Nine Biomes. Nature 2008.Google Scholar
- Krause L, Diaz NN, Bartels D, Edwards RA, Puhler A, Rohwer F, Meyer F, Stoye J: Finding novel genes in bacterial communities isolated from the environment. Bioinformatics 2006, 22(14):e281–289. 10.1093/bioinformatics/btl247View ArticlePubMedGoogle Scholar
- Liang F, Holt I, Pertea G, Karamycheva S, Salzberg SL, Quackenbush J: An optimized protocol for analysis of EST sequences. Nucleic Acids Res 2000, 28(18):3657–3665. 10.1093/nar/28.18.3657PubMed CentralView ArticlePubMedGoogle Scholar
- Rohwer F: personal communication. 2007.Google Scholar
- Paarmann D, Paczian T, Meyer F: SEED-Viewer – a web user interface for studying genomics datasets. unpublished
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.