MG-RAST, a Metagenomics Service for Analysis of Microbial Community Structure and Function.

Approaches in molecular biology, particularly those that deal with high-throughput sequencing of entire microbial communities (the field of metagenomics), are rapidly advancing our understanding of the composition and functional content of microbial communities involved in climate change, environmental pollution, human health, biotechnology, etc. Metagenomics provides researchers with the most complete picture of the taxonomic (i.e., what organisms are there) and functional (i.e., what are those organisms doing) composition of natively sampled microbial communities, making it possible to perform investigations that include organisms that were previously intractable to laboratory-controlled culturing; currently, these constitute the vast majority of all microbes on the planet. All organisms contained in environmental samples are sequenced in a culture-independent manner, most often with 16S ribosomal amplicon methods to investigate the taxonomic or whole-genome shotgun-based methods to investigate the functional content of sampled communities. Metagenomics allows researchers to characterize the community composition and functional content of microbial communities, but it cannot show which functional processes are active; however, near parallel developments in transcriptomics promise a dramatic increase in our knowledge in this area as well. Since 2008, MG-RAST (Meyer et al., BMC Bioinformatics 9:386, 2008) has served as a public resource for annotation and analysis of metagenomic sequence data, providing a repository that currently houses more than 150,000 data sets (containing 60+ tera-base-pairs) with more than 23,000 publically available. MG-RAST, or the metagenomics RAST (rapid annotation using subsystems technology) server makes it possible for users to upload raw metagenomic sequence data in (preferably) fastq or fasta format. Assessments of sequence quality, annotation with respect to multiple reference databases, are performed automatically with minimal input from the user (see Subheading 4 at the end of this chapter for more details). Post-annotation analysis and visualization are also possible, directly through the web interface, or with tools like matR (metagenomic analysis tools for R, covered later in this chapter) that utilize the MG-RAST API ( http://api.metagenomics.anl.gov/api.html ) to easily download data from any stage in the MG-RAST processing pipeline. Over the years, MG-RAST has undergone substantial revisions to keep pace with the dramatic growth in the number, size, and types of sequence data that accompany constantly evolving developments in metagenomics and related -omic sciences (e.g., metatranscriptomics).

A RESTful API for accessing microbial community data for MG-RAST.

Metagenomic sequencing has produced significant amounts of data in recent years. For example, as of summer 2013, MG-RAST has been used to annotate over 110,000 data sets totaling over 43 Terabases. With metagenomic sequencing finding even wider adoption in the scientific community, the existing web-based analysis tools and infrastructure in MG-RAST provide limited capability for data retrieval and analysis, such as comparative analysis between multiple data sets. Moreover, although the system provides many analysis tools, it is not comprehensive. By opening MG-RAST up via a web services API (application programmers interface) we have greatly expanded access to MG-RAST data, as well as provided a mechanism for the use of third-party analysis tools with MG-RAST data. This RESTful API makes all data and data objects created by the MG-RAST pipeline accessible as JSON objects. As part of the DOE Systems Biology Knowledgebase project (KBase, http://kbase.us) we have implemented a web services API for MG-RAST. This API complements the existing MG-RAST web interface and constitutes the basis of KBase's microbial community capabilities. In addition, the API exposes a comprehensive collection of data to programmers. This API, which uses a RESTful (Representational State Transfer) implementation, is compatible with most programming environments and should be easy to use for end users and third parties. It provides comprehensive access to sequence data, quality control results, annotations, and many other data types. Where feasible, we have used standards to expose data and metadata. Code examples are provided in a number of languages both to show the versatility of the API and to provide a starting point for users. We present an API that exposes the data in MG-RAST for consumption by our users, greatly enhancing the utility of the MG-RAST service.

Target selection for structural genomics of infectious diseases.

This chapter describes the protocols used to identify, filter, and annotate potential protein targets from an organism associated with infectious diseases. Protocols often combine computational approaches for mining information in public databases or for checking whether the protein has already been targeted for structure determination, with manual strategies that examine the literature for information on the biological role of the protein or the experimental strategies that explore the effects of knocking out the protein. Publicly available computational tools have been cited as much as possible. Where these do not exist, the concepts underlying in-house tools developed for the Center for Structural Genomics of Infectious Diseases have been described.

The complete genome sequence for putative H2- and S-oxidizer Candidatus Sulfuricurvum sp., assembled de novo from an aquifer-derived metagenome.

We reconstructed the complete 2.4 Mb-long genome of a previously uncultivated epsilonproteobacterium, Candidatus Sulfuricurvum sp. RIFRC-1, via assembly of short-read shotgun metagenomic data using a complexity reduction approach. Genome-based comparisons indicate the bacterium is a novel species within the Sulfuricurvum genus, which contains one cultivated representative, S. kujiense. Divergence between the species appears due in part to extensive genomic rearrangements, gene loss and chromosomal versus plasmid encoding of certain (respiratory) genes by RIFRC-1. Deoxyribonucleic acid for the genome was obtained from terrestrial aquifer sediment, in which RIFRC-1 comprised ∼ 47% of the bacterial community. Genomic evidence suggests RIFRC-1 is a chemolithoautotrophic diazotroph capable of deriving energy for growth by microaerobic or nitrate-/nitric oxide-dependent oxidation of S°, sulfide or sulfite or H2oxidation. Carbon may be fixed via the reductive tricarboxylic acid cycle. Consistent with these physiological attributes, the local aquifer was microoxic with small concentrations of available nitrate, small but elevated concentrations of reduced sulfur and NH(4)(+) /NH3-limited. Additionally, various mechanisms for heavy metal and metalloid tolerance and virulence point to a lifestyle well-adapted for metal(loid)-rich environments and a shared evolutionary past with pathogenic Epsilonproteobacteria. Results expand upon recent findings highlighting the potential importance of sulfur and hydrogen metabolism in the terrestrial subsurface.

A metagenomics portal for a democratized sequencing world.

The democratized world of sequencing is leading to numerous data analysis challenges; MG-RAST addresses many of these challenges for diverse datasets, including amplicon datasets, shotgun metagenomes, and metatranscriptomes. The changes from version 2 to version 3 include the addition of a dedicated gene calling stage using FragGenescan, clustering of predicted proteins at 90% identity, and the use of BLAT for the computation of similarities. Together with changes in the underlying software infrastructure, this has enabled the dramatic scaling up of pipeline throughput while remaining on a limited hardware budget. The Web-based service allows upload, fully automated analysis, and visualization of results. As a result of the plummeting cost of sequencing and the readily available analytical power of MG-RAST, over 78,000 metagenomic datasets have been analyzed, with over 12,000 of them publicly available in MG-RAST.

The M5nr: a novel non-redundant database containing protein sequences and annotations from multiple sources and associated tools.

BACKGROUND: Computing of sequence similarity results is becoming a limiting factor in metagenome analysis. Sequence similarity search results encoded in an open, exchangeable format have the potential to limit the needs for computational reanalysis of these data sets. A prerequisite for sharing of similarity results is a common reference. DESCRIPTION: We introduce a mechanism for automatically maintaining a comprehensive, non-redundant protein database and for creating a quarterly release of this resource. In addition, we present tools for translating similarity searches into many annotation namespaces, e.g. KEGG or NCBI's GenBank. CONCLUSIONS: The data and tools we present allow the creation of multiple result sets using a single computation, permitting computational results to be shared between groups for large sequence data sets.

Identification of chromosomal genes in Yersinia pestis that influence type III secretion and delivery of Yops into target cells.

Pathogenic Yersinia species possess a type III secretion system, which is required for the delivery of effector Yop proteins into target cells during infection. Genes encoding the type III secretion machinery, its substrates, and several regulatory proteins all reside on a 70-Kb virulence plasmid. Genes encoded in the chromosome of yersiniae are thought to play important roles in bacterial perception of host environments and in the coordinated activation of the type III secretion pathway. Here, we investigate the contribution of chromosomal genes to the complex regulatory process controlling type III secretion in Yersinia pestis. Using transposon mutagenesis, we identified five chromosomal genes required for expression or secretion of Yops in laboratory media. Four out of the five chromosomal mutants were defective to various extents at injecting Yops into tissue culture cells. Interestingly, we found one mutant that was not able to secrete in vitro but was fully competent for injecting Yops into host cells, suggesting independent mechanisms for activation of the secretion apparatus. When tested in a mouse model of plague disease, three mutants were avirulent, whereas two strains were severely attenuated. Together these results demonstrate the importance of Y. pestis chromosomal genes in the proper function of type III secretion and in the pathogenesis of plague.

From genomics to metagenomics.

Next-generation sequencing has changed metagenomics. However, sequencing DNA is no longer the bottleneck, rather, the bottleneck is computational analysis and also interpretation. Computational cost is the obvious issue, as is tool limitations, considering most of the tools we routinely use have been built for clonal genomics or are being adapted to microbial communities. The current trend in metagenomics analysis is toward reducing computational costs through improved algorithms and through analysis strategies. Data sharing and interoperability between tools are critical, since computation for metagenomic datasets is very high.

Meeting report: the 2 annual argonne soils workshop, argonne national laboratory, chicago illinois, USA, october 6-8, 2010.

This report summarizes the proceedings of the 2(nd) Annual Argonne Soils Workshop held at Argonne National Laboratory October 6-8, 2010. The workshop assembled a diverse group of soil ecologists, microbiologists, molecular biologists, and computational scientists to discuss the challenges and opportunities related to implementation of metagenomics approaches in soil microbial ecology. The overarching theme of the workshop was "designing ecologically meaningful soil metagenomics research", which encouraged presentations on both ecological and computational topics. The workshop fostered valuable cross-discipline communication and delivered the message that soil metagenomics research must be based on an iterative process between biological inquiry and bioinformatics tools.

Generating a collection of insertion mutations in the Staphylococcus aureus genome using bursa aurealis.

Staphylococcus aureus is the leading cause of wound and hospital-acquired infections. The emergence of strains with resistance to all antibiotics has created a serious public health problem. Transposon-based mutagenesis can be used to generate libraries of mutants and to query genomes for factors involved in nonessential pathways, such as virulence and antibiotic resistance. Ideally, such studies should employ defined and complete sets of isogenic mutants and should be conducted so as to permit acquisition and comparison of the complete data sets. Such systematic knowledge can reveal entire pathways and can be exploited for the rational design of therapies. The mariner-based transposon, bursa aurealis, can be used to generate random libraries of mutants in laboratory strains and clinical isolates of S. aureus. This chapter describes a procedure for isolating mutants and mapping the insertion sites on the chromosome.

The RAST Server: rapid annotations using subsystems technology.

BACKGROUND: The number of prokaryotic genome sequences becoming available is growing steadily and is growing faster than our ability to accurately annotate them. DESCRIPTION: We describe a fully automated service for annotating bacterial and archaeal genomes. The service identifies protein-encoding, rRNA and tRNA genes, assigns functions to the genes, predicts which subsystems are represented in the genome, uses this information to reconstruct the metabolic network and makes the output easily downloadable for the user. In addition, the annotated genome can be browsed in an environment that supports comparative analysis with the annotated genomes maintained in the SEED environment. The service normally makes the annotated genome available within 12-24 hours of submission, but ultimately the quality of such a service will be judged in terms of accuracy, consistency, and completeness of the produced annotations. We summarize our attempts to address these issues and discuss plans for incrementally enhancing the service. CONCLUSION: By providing accurate, rapid annotation freely to the community we have created an important community resource. The service has now been utilized by over 120 external users annotating over 350 distinct genomes.

Identification of Francisella tularensis Himar1-based transposon mutants defective for replication in macrophages.

Francisella tularensis, the etiologic agent of tularemia in humans, is a potential biological threat due to its low infectious dose and multiple routes of entry. F. tularensis replicates within several cell types, eventually causing cell death by inducing apoptosis. In this study, a modified Himar1 transposon (HimarFT) was used to mutagenize F. tularensis LVS. Approximately 7,000 Km(r) clones were screened using J774A.1 macrophages for reduction in cytopathogenicity based on retention of the cell monolayer. A total of 441 candidates with significant host cell retention compared to the parent were identified following screening in a high-throughput format. Retesting at a defined multiplicity of infection followed by in vitro growth analyses resulted in identification of approximately 70 candidates representing 26 unique loci involved in macrophage replication and/or cytotoxicity. Mutants carrying insertions in seven hypothetical genes were screened in a mouse model of infection, and all strains tested appeared to be attenuated, which validated the initial in vitro results obtained with cultured macrophages. Complementation and reverse transcription-PCR experiments suggested that the expression of genes adjacent to the HimarFT insertion may be affected depending on the orientation of the constitutive groEL promoter region used to ensure transcription of the selective marker in the transposon. A hypothetical gene, FTL_0706, postulated to be important for lipopolysaccharide biosynthesis, was confirmed to be a gene involved in O-antigen expression in F. tularensis LVS and Schu S4. These and other studies demonstrate that therapeutic targets, vaccine candidates, or virulence-related genes may be discovered utilizing classical genetic approaches in Francisella.

Sentra: a database of signal transduction proteins for comparative genome analysis.

Sentra (http://compbio.mcs.anl.gov/sentra), a database of signal transduction proteins encoded in completely sequenced prokaryotic genomes, has been updated to reflect recent advances in understanding signal transduction events on a whole-genome scale. Sentra consists of two principal components, a manually curated list of signal transduction proteins in 202 completely sequenced prokaryotic genomes and an automatically generated listing of predicted signaling proteins in 235 sequenced genomes that are awaiting manual curation. In addition to two-component histidine kinases and response regulators, the database now lists manually curated Ser/Thr/Tyr protein kinases and protein phosphatases, as well as adenylate and diguanylate cyclases and c-di-GMP phosphodiesterases, as defined in several recent reviews. All entries in Sentra are extensively annotated with relevant information from public databases (e.g. UniProt, KEGG, PDB and NCBI). Sentra's infrastructure was redesigned to support interactive cross-genome comparisons of signal transduction capabilities of prokaryotic organisms from a taxonomic and phenotypic perspective and in the framework of signal transduction pathways from KEGG. Sentra leverages the PUMA2 system to support interactive analysis and annotation of signal transduction proteins by the users.

Transposon mutagenesis of Bacillus anthracis strain Sterne using Bursa aurealis.

Bacillus anthracis, a spore forming Gram-positive microbe, is the causative agent of anthrax. Although plasmid encoded factors such as lethal toxin (LeTx), edema toxin (EdTx), and gamma-poly-d-glutamic acid (PGA) capsule are known to be required for disease pathogenesis, B. anthracis genes that enable spore invasion, phagosomal escape and macrophage replication are still unknown. To establish transposon mutagenesis as a tool for the characterization of anthrax genes, we employed the mariner-based mini-transposon Bursa aurealis in B. anthracis strain Sterne 7702. B. aurealis carrying an erythromycin resistance cassette and its cognate transposase were delivered by transformation of two plasmids. B. aurealis transposition can be selected for by temperature shift to prevent plasmid replication and by screening colonies for erythromycin resistance. Using inverse polymerase chain reaction, DNA fragments of 129 random erythromycin-resistant transposon mutants were amplified and submitted to DNA sequence analysis. These studies demonstrate that B. aurealis inserts randomly into the genome of B. anthracis and can therefore be employed for finding genes involved in virulence.

The subsystems approach to genome annotation and its use in the project to annotate 1000 genomes.

The release of the 1000th complete microbial genome will occur in the next two to three years. In anticipation of this milestone, the Fellowship for Interpretation of Genomes (FIG) launched the Project to Annotate 1000 Genomes. The project is built around the principle that the key to improved accuracy in high-throughput annotation technology is to have experts annotate single subsystems over the complete collection of genomes, rather than having an annotation expert attempt to annotate all of the genes in a single genome. Using the subsystems approach, all of the genes implementing the subsystem are analyzed by an expert in that subsystem. An annotation environment was created where populated subsystems are curated and projected to new genomes. A portable notion of a populated subsystem was defined, and tools developed for exchanging and curating these objects. Tools were also developed to resolve conflicts between populated subsystems. The SEED is the first annotation environment that supports this model of annotation. Here, we describe the subsystem approach, and offer the first release of our growing library of populated subsystems. The initial release of data includes 180 177 distinct proteins with 2133 distinct functional roles. This data comes from 173 subsystems and 383 different organisms.

Bacillus anthracis sortase A (SrtA) anchors LPXTG motif-containing surface proteins to the cell wall envelope.

Cell wall-anchored surface proteins of gram-positive pathogens play important roles during the establishment of many infectious diseases, but the contributions of surface proteins to the pathogenesis of anthrax have not yet been revealed. Cell wall anchoring in Staphylococcus aureus occurs by a transpeptidation mechanism requiring surface proteins with C-terminal sorting signals as well as sortase enzymes. The genome sequence of Bacillus anthracis encodes three sortase genes and eleven surface proteins with different types of cell wall sorting signals. Purified B. anthracis sortase A cleaved peptides encompassing LPXTG motif-type sorting signals between the threonine (T) and the glycine (G) residues in vitro. Sortase A activity could be inhibited by thiol-reactive reagents, similar to staphylococcal sortases. B. anthracis parent strain Sterne 34F(2), but not variants lacking the srtA gene, anchored the collagen-binding MSCRAMM (microbial surface components recognizing adhesive matrix molecules) BasC (BA5258/BAS4884) to the bacterial cell wall. These results suggest that B. anthracis SrtA anchors surface proteins bearing LPXTG motif sorting signals to the cell wall envelope of vegetative bacilli.

Molecular mechanisms involved in robustness of yeast central metabolism against null mutations.

Adaptive strategies employed by the yeast Saccharomyces cerevisiae provide robustness and adaptability of its central metabolism. Since central metabolism in yeast has been well studied at the enzymatic and genetic levels, it represents an excellent system for evaluating the relative roles of duplicate genes and alternative metabolic pathways as possible mechanisms for the stability of central metabolism against null mutations. Yeast appears to employ a variety of mechanisms to ensure functional robustness of its central metabolism. Uninterrupted flow of energy and precursor metabolites through the pathways of central metabolism via glycolysis (EMP), pentose phosphate shunt (PPS), and the tricarboxylic acid (TCA) cycle are ensured by a variety of adaptive mechanisms. One of the most significant mechanisms appears to be gene duplication events that have produced a number of isozymes functioning under variable environmental and physiological conditions. Alternative pathways represent another important mechanism for increasing the robustness of the system. The robustness of the pathways of central metabolism is apparently higher than that of the other parts of metabolism, because of its exceptional importance to the organism's vitality. The proportion of duplicated viable genes also is substantially larger in central metabolism than that in a pool of other metabolic genes.