IMG: the Integrated Microbial Genomes database and comparative analysis system

Coupling FACS and genomic methods for the characterization of uncultivated symbionts

Symbioses between microbes are likely widespread and functionally relevant in diverse biological systems; however, they are difficult to discover. Most microbes remain uncultivated, symbioses can be relatively rare or dynamic, and intercellular connections can be delicate. Thus, traditional methods such as microscopy are inadequate for efficient discovery and precise characterization of novel interactions, their metabolic basis, and the species involved. High-throughput metagenomic sequencing of entire microbial communities has revolutionized the field of microbial ecology; however, genomic signals from symbionts can get buried in sequences from abundant organisms and evidence for direct links between microbial species cannot be gained from bulk samples. Thus, a specialized approach to the characterization of symbioses between naturally occurring microbes is required. This chapter presents methods for combining fluorescence-activated cell sorting to isolate and separate uncultivated symbionts with molecular biology techniques for DNA amplification in order to characterize uncultivated symbionts through genomic and metagenomic techniques.

Integrating grant-funded research into the undergraduate biology curriculum using IMG-ACT

It has become clear in current scientific pedagogy that the emersion of students in the scientific process in terms of designing, implementing, and analyzing experiments is imperative for their education; as such, it has been our goal to model this active learning process in the classroom and laboratory in the context of a genuine scientific question. Toward this objective, the National Science Foundation funded a collaborative research grant between a primarily undergraduate institution and a research-intensive institution to study the chemotactic responses of the bacterium Pseudomonas putida F1. As part of the project, a new Bioinformatics course was developed in which undergraduates annotate relevant regions of the P. putida F1 genome using Integrated Microbial Genomes Annotation Collaboration Toolkit, a bioinformatics interface specifically developed for undergraduate programs by the Department of Energy Joint Genome Institute. Based on annotations of putative chemotaxis genes in P. putida F1 and comparative genomics studies, undergraduate students from both institutions developed functional genomics research projects that evolved from the annotations. The purpose of this study is to describe the nature of the NSF grant, the development of the Bioinformatics lecture and wet laboratory course, and how undergraduate student involvement in the project that was initiated in the classroom has served as a springboard for independent undergraduate research projects.

A cyclic form of N6-threonylcarbamoyladenosine as a widely distributed tRNA hypermodification

N(6)-threonylcarbamoyladenosine (t(6)A) is a universally conserved, essential modified nucleoside found in transfer RNAs (tRNAs) responsible for ANN codons in all three domains of life. t(6)A has a crucial role in maintaining decoding accuracy during protein synthesis. The presence of t(6)A in cellular tRNAs has been well documented for more than four decades. However, under conditions optimized for nucleoside preparation, we detected little t(6)A in tRNAs from Escherichia coli. Instead, we identified a new modified base named 'cyclic t(6)A' (ct(6)A), which is a cyclized active ester with an oxazolone ring. An E1-like enzyme, CsdL (renamed as TcdA), which catalyzes ATP-dependent dehydration of t(6)A to form ct(6)A, was also identified. Two yeast homologs of tcdA, YHR003C (TCD1) and YKL027W (TCD2), were required for ct(6)A formation and respiratory cell growth. ct(6)A was involved in promoting decoding efficiency. Structural modeling suggests that ct(6)A recognizes the first adenine base of ANN codon at the ribosomal A site.

Evaluation method for the potential functionome harbored in the genome and metagenome

Background

One of the main goals of genomic analysis is to elucidate the comprehensive functions (functionome) in individual organisms or a whole community in various environments. However, a standard evaluation method for discerning the functional potentials harbored within the genome or metagenome has not yet been established. We have developed a new evaluation method for the potential functionome, based on the completion ratio of Kyoto Encyclopedia of Genes and Genomes (KEGG) functional modules.

Results

Distribution of the completion ratio of the KEGG functional modules in 768 prokaryotic species varied greatly with the kind of module, and all modules primarily fell into 4 patterns (universal, restricted, diversified and non-prokaryotic modules), indicating the universal and unique nature of each module, and also the versatility of the KEGG Orthology (KO) identifiers mapped to each one. The module completion ratio in 8 phenotypically different bacilli revealed that some modules were shared only in phenotypically similar species. Metagenomes of human gut microbiomes from 13 healthy individuals previously determined by the Sanger method were analyzed based on the module completion ratio. Results led to new discoveries in the nutritional preferences of gut microbes, believed to be one of the mutualistic representations of gut microbiomes to avoid nutritional competition with the host.

Conclusions

The method developed in this study could characterize the functionome harbored in genomes and metagenomes. As this method also provided taxonomical information from KEGG modules as well as the gene hosts constructing the modules, interpretation of completion profiles was simplified and we could identify the complementarity between biochemical functions in human hosts and the nutritional preferences in human gut microbiomes. Thus, our method has the potential to be a powerful tool for comparative functional analysis in genomics and metagenomics, able to target unknown environments containing various uncultivable microbes within unidentified phyla.

Using Bioinformatics to Develop and Test Hypotheses: E. coli-Specific Virulence Determinants

Bioinformatics, the use of computer resources to understand biological information, is an important tool in research, and can be easily integrated into the curriculum of undergraduate courses. Such an example is provided in this series of four activities that introduces students to the field of bioinformatics as they design PCR based tests for pathogenic E. coli strains. A variety of computer tools are used including BLAST searches at NCBI, bacterial genome searches at the Integrated Microbial Genomes (IMG) database, protein analysis at Pfam and literature research at PubMed. In the process, students also learn about virulence factors, enzyme function and horizontal gene transfer. Some or all of the four activities can be incorporated into microbiology or general biology courses taken by students at a variety of levels, ranging from high school through college. The activities build on one another as they teach and reinforce knowledge and skills, promote critical thinking, and provide for student collaboration and presentation. The computer-based activities can be done either in class or outside of class, thus are appropriate for inclusion in online or blended learning formats. Assessment data showed that students learned general microbiology concepts related to pathogenesis and enzyme function, gained skills in using tools of bioinformatics and molecular biology, and successfully developed and tested a scientific hypothesis.

By their genes ye shall know them: genomic signatures of predatory bacteria

Predatory bacteria are taxonomically disparate, exhibit diverse predatory strategies and are widely distributed in varied environments. To date, their predatory phenotypes cannot be discerned in genome sequence data thereby limiting our understanding of bacterial predation, and of its impact in nature. Here, we define the 'predatome,' that is, sets of protein families that reflect the phenotypes of predatory bacteria. The proteomes of all sequenced 11 predatory bacteria, including two de novo sequenced genomes, and 19 non-predatory bacteria from across the phylogenetic and ecological landscapes were compared. Protein families discriminating between the two groups were identified and quantified, demonstrating that differences in the proteomes of predatory and non-predatory bacteria are large and significant. This analysis allows predictions to be made, as we show by confirming from genome data an over-looked bacterial predator. The predatome exhibits deficiencies in riboflavin and amino acids biosynthesis, suggesting that predators obtain them from their prey. In contrast, these genomes are highly enriched in adhesins, proteases and particular metabolic proteins, used for binding to, processing and consuming prey, respectively. Strikingly, predators and non-predators differ in isoprenoid biosynthesis: predators use the mevalonate pathway, whereas non-predators, like almost all bacteria, use the DOXP pathway. By defining predatory signatures in bacterial genomes, the predatory potential they encode can be uncovered, filling an essential gap for measuring bacterial predation in nature. Moreover, we suggest that full-genome proteomic comparisons are applicable to other ecological interactions between microbes, and provide a convenient and rational tool for the functional classification of bacteria.

OrtholugeDB: a bacterial and archaeal orthology resource for improved comparative genomic analysis

Prediction of orthologs (homologous genes that diverged because of speciation) is an integral component of many comparative genomics methods. Although orthologs are more likely to have similar function versus paralogs (genes that diverged because of duplication), recent studies have shown that their degree of functional conservation is variable. Also, there are inherent problems with several large-scale ortholog prediction approaches. To address these issues, we previously developed Ortholuge, which uses phylogenetic distance ratios to provide more precise ortholog assessments for a set of predicted orthologs. However, the original version of Ortholuge required manual intervention and was not easily accessible; therefore, we now report the development of OrtholugeDB, available online at http://www.pathogenomics.sfu.ca/ortholugedb. OrtholugeDB provides ortholog predictions for completely sequenced bacterial and archaeal genomes from NCBI based on reciprocal best Basic Local Alignment Search Tool hits, supplemented with further evaluation by the more precise Ortholuge method. The OrtholugeDB web interface facilitates user-friendly and flexible ortholog analysis, from single genes to genomes, plus flexible data download options. We compare Ortholuge with similar methods, showing how it may more consistently identify orthologs with conserved features across a wide range of taxonomic distances. OrtholugeDB facilitates rapid, and more accurate, bacterial and archaeal comparative genomic analysis and large-scale ortholog predictions.

Comparative genomic and phylogenetic approaches to characterize the role of genetic recombination in mycobacterial evolution

The genus Mycobacterium encompasses over one hundred named species of environmental and pathogenic organisms, including the causative agents of devastating human diseases such as tuberculosis and leprosy. The success of these human pathogens is due in part to their ability to rapidly adapt to their changing environment and host. Recombination is the fastest way for bacterial genomes to acquire genetic material, but conflicting results about the extent of recombination in the genus Mycobacterium have been reported. We examined a data set comprising 18 distinct strains from 13 named species for evidence of recombination. Genomic regions common to all strains (accounting for 10% to 22% of the full genomes of all examined species) were aligned and concatenated in the chromosomal order of one mycobacterial reference species. The concatenated sequence was screened for evidence of recombination using a variety of statistical methods, with each proposed event evaluated by comparing maximum-likelihood phylogenies of the recombinant section with the non-recombinant portion of the dataset. Incongruent phylogenies were identified by comparing the site-wise log-likelihoods of each tree using multiple tests. We also used a phylogenomic approach to identify genes that may have been acquired through horizontal transfer from non-mycobacterial sources. The most frequent associated lineages (and potential gene transfer partners) in the Mycobacterium lineage-restricted gene trees are other members of suborder Corynebacterinae, but more-distant partners were identified as well. In two examined cases of potentially frequent and habitat-directed transfer (M. abscessus to Segniliparus and M. smegmatis to Streptomyces), observed sequence distances were small and consistent with a hypothesis of transfer, while in a third case (M. vanbaalenii to Streptomyces) distances were larger. The analyses described here indicate that whereas evidence of recombination in core regions within the genus is relatively sparse, the acquisition of genes from non-mycobacterial lineages is a significant feature of mycobacterial evolution.

Functional Annotation Analytics of Bacillus Genomes Reveals Stress Responsive Acetate Utilization and Sulfate Uptake in the Biotechnologically Relevant Bacillus megaterium

Bacillus species form an heterogeneous group of Gram-positive bacteria that include members that are disease-causing, biotechnologically-relevant, and can serve as biological research tools. A common feature of Bacillus species is their ability to survive in harsh environmental conditions by formation of resistant endospores. Genes encoding the universal stress protein (USP) domain confer cellular and organismal survival during unfavorable conditions such as nutrient depletion. As of February 2012, the genome sequences and a variety of functional annotations for at least 123 Bacillus isolates including 45 Bacillus cereus isolates were available in public domain bioinformatics resources. Additionally, the genome sequencing status of 10 of the B. cereus isolates were annotated as finished with each genome encoded 3 USP genes. The conservation of gene neighborhood of the 140 aa universal stress protein in the B. cereus genomes led to the identification of a predicted plasmid-encoded transcriptional unit that includes a USP gene and a sulfate uptake gene in the soil-inhabiting Bacillus megaterium. Gene neighborhood analysis combined with visual analytics of chemical ligand binding sites data provided knowledge-building biological insights on possible cellular functions of B. megaterium universal stress proteins. These functions include sulfate and potassium uptake, acid extrusion, cellular energy-level sensing, survival in high oxygen conditions and acetate utilization. Of particular interest was a two-gene transcriptional unit that consisted of genes for a universal stress protein and a sirtuin Sir2 (deacetylase enzyme for NAD+-dependent acetate utilization). The predicted transcriptional units for stress responsive inorganic sulfate uptake and acetate utilization could explain biological mechanisms for survival of soil-inhabiting Bacillus species in sulfate and acetate limiting conditions. Considering the key role of sirtuins in mammalian physiology additional research on the USP-Sir2 transcriptional unit of B. megaterium could help explain mammalian acetate metabolism in glucose-limiting conditions such as caloric restriction. Finally, the deep-rooted position of B. megaterium in the phylogeny of Bacillus species makes the investigation of the functional coupling acetate utilization and stress response compelling.

Lateral gene transfer of an ABC transporter complex between major constituents of the human gut microbiome

Background

Several links have been established between the human gut microbiome and conditions such as obesity and inflammatory bowel syndrome. This highlights the importance of understanding what properties of the gut microbiome can affect the health of the human host. Studies have been undertaken to determine the species composition of this microbiome and infer functional profiles associated with such host properties. However, lateral gene transfer (LGT) between community members may result in misleading taxonomic attributions for the recipient organisms, thus making species-function links difficult to establish.

Results

We identified a peptides/nickel transport complex whose components differed in abundance based upon levels of host obesity, and assigned the encoded proteins to members of the microbial community. Each protein was assigned to several distinct taxonomic groups, with moderate levels of agreement observed among different proteins in the complex. Phylogenetic trees of these proteins produced clusters that differed greatly from taxonomic attributions and indicated that habitat-directed LGT of this complex is likely to have occurred, though not always between the same partners.

Conclusions

These findings demonstrate that certain membrane transport systems may be an important factor within an obese-associated gut microbiome and that such complexes may be acquired several times by different strains of the same species. Additionally, an example of individual proteins from different organisms being transferred into one operon was observed, potentially demonstrating a functional complex despite the donors of the subunits being taxonomically disparate. Our results also highlight the potential impact of habitat-directed LGT on the resident microbiota.

MBGD update 2013: the microbial genome database for exploring the diversity of microbial world

The microbial genome database for comparative analysis (MBGD, available at http://mbgd.genome.ad.jp/) is a platform for microbial genome comparison based on orthology analysis. As its unique feature, MBGD allows users to conduct orthology analysis among any specified set of organisms; this flexibility allows MBGD to adapt to a variety of microbial genomic study. Reflecting the huge diversity of microbial world, the number of microbial genome projects now becomes several thousands. To efficiently explore the diversity of the entire microbial genomic data, MBGD now provides summary pages for pre-calculated ortholog tables among various taxonomic groups. For some closely related taxa, MBGD also provides the conserved synteny information (core genome alignment) pre-calculated using the CoreAligner program. In addition, efficient incremental updating procedure can create extended ortholog table by adding additional genomes to the default ortholog table generated from the representative set of genomes. Combining with the functionalities of the dynamic orthology calculation of any specified set of organisms, MBGD is an efficient and flexible tool for exploring the microbial genome diversity.

Resuscitation promoting factor (Rpf) from Tomitella biformata AHU 1821(T) promotes growth and resuscitates non-dividing cells

Functional variation of Rpf, a growth factor found exclusively in Actinobacteria, is differentiated by its source and amino acid sequences. Only purified Rpf proteins from three species have been studied so far. To seek new Rpfs for use in future studies to understand their role in Actinobacteria, the objective of this study was to identify rpf gene homologs in Tomitella biformata AHU 1821^T, a novel Actinobacteria isolated from permafrost ice wedge. Amplification using degenerate primers targeting the essential Rpf domain led to the discovery of a new rpf gene in T. biformata. Gene structure and the deduced Rpf domain amino acid sequence indicated that this rpf gene was not identical to previously studied Rpf. Phylogenetic analysis placed T. biformata Rpf in a monophyletic branch in the RpfB subfamily. The deduced amino acid sequence was 44.9% identical to RpfB in Mycobacterium tuberculosis, the closest functionally tested Rpf. The gene was cloned and expressed in Escherichia coli; the recombinant Rpf protein (rRpf) promoted the growth of dividing cells and resuscitated non-dividing cells of T. biformata. Compared to other studies, this Rpf was required at higher concentrations to promote its growth and to resuscitate itself from a non-dividing state. The resuscitation function was likely due to the highly conserved Rpf domain. This study provides evidence that a genetically unique but functional Rpf can be found in novel members of Actinobacteria and can lead to a better understanding of bacterial cytokines in this phylum.

Marine proteomics: a critical assessment of an emerging technology

The application of proteomics to marine sciences has increased in recent years because the proteome represents the interface between genotypic and phenotypic variability and, thus, corresponds to the broadest possible biomarker for eco-physiological responses and adaptations. Likewise, proteomics can provide important functional information regarding biosynthetic pathways, as well as insights into mechanism of action, of novel marine natural products. The goal of this review is to (1) explore the application of proteomics methodologies to marine systems, (2) assess the technical approaches that have been used, and (3) evaluate the pros and cons of this proteomic research, with the intent of providing a critical analysis of its future roles in marine sciences. To date, proteomics techniques have been utilized to investigate marine microbe, plant, invertebrate, and vertebrate physiology, developmental biology, seafood safety, susceptibility to disease, and responses to environmental change. However, marine proteomics studies often suffer from poor experimental design, sample processing/optimization difficulties, and data analysis/interpretation issues. Moreover, a major limitation is the lack of available annotated genomes and proteomes for most marine organisms, including several "model species". Even with these challenges in mind, there is no doubt that marine proteomics is a rapidly expanding and powerful integrative molecular research tool from which our knowledge of the marine environment, and the natural products from this resource, will be significantly expanded.

Sifting through genomes with iterative-sequence clustering produces a large, phylogenetically diverse protein-family resource

BACKGROUND

New computational resources are needed to manage the increasing volume of biological data from genome sequencing projects. One fundamental challenge is the ability to maintain a complete and current catalog of protein diversity. We developed a new approach for the identification of protein families that focuses on the rapid discovery of homologous protein sequences.

RESULTS

We implemented fully automated and high-throughput procedures to de novo cluster proteins into families based upon global alignment similarity. Our approach employs an iterative clustering strategy in which homologs of known families are sifted out of the search for new families. The resulting reduction in computational complexity enables us to rapidly identify novel protein families found in new genomes and to perform efficient, automated updates that keep pace with genome sequencing. We refer to protein families identified through this approach as "Sifting Families," or SFams. Our analysis of ~10.5 million protein sequences from 2,928 genomes identified 436,360 SFams, many of which are not represented in other protein family databases. We validated the quality of SFam clustering through statistical as well as network topology-based analyses.

CONCLUSIONS

We describe the rapid identification of SFams and demonstrate how they can be used to annotate genomes and metagenomes. The SFam database catalogs protein-family quality metrics, multiple sequence alignments, hidden Markov models, and phylogenetic trees. Our source code and database are publicly available and will be subject to frequent updates (http://edhar.genomecenter.ucdavis.edu/sifting_families/).

Genomic and physiological characterization of the chromate-reducing, aquifer-derived Firmicute Pelosinus sp. strain HCF1

Pelosinus spp. are fermentative firmicutes that were recently reported to be prominent members of microbial communities at contaminated subsurface sites in multiple locations. Here we report metabolic characteristics and their putative genetic basis in Pelosinus sp. strain HCF1, an isolate that predominated anaerobic, Cr(VI)-reducing columns constructed with aquifer sediment. Strain HCF1 ferments lactate to propionate and acetate (the methylmalonyl-coenzyme A [CoA] pathway was identified in the genome), and its genome encodes two [NiFe]- and four [FeFe]-hydrogenases for H(2) cycling. The reduction of Cr(VI) and Fe(III) may be catalyzed by a flavoprotein with 42 to 51% sequence identity to both ChrR and FerB. This bacterium has unexpected capabilities and gene content associated with reduction of nitrogen oxides, including dissimilatory reduction of nitrate to ammonium (two copies of NrfH and NrfA were identified along with NarGHI) and a nitric oxide reductase (NorCB). In this strain, either H(2) or lactate can act as a sole electron donor for nitrate, Cr(VI), and Fe(III) reduction. Transcriptional studies demonstrated differential expression of hydrogenases and nitrate and nitrite reductases. Overall, the unexpected metabolic capabilities and gene content reported here broaden our perspective on what biogeochemical and ecological roles this species might play as a prominent member of microbial communities in subsurface environments.

Sequence, structure, and evolution of cellulases in glycoside hydrolase family 48

Currently, the cost of cellulase enzymes remains a key economic impediment to commercialization of biofuels. Enzymes from glycoside hydrolase family 48 (GH48) are a critical component of numerous natural lignocellulose-degrading systems. Although computational mining of large genomic data sets is a promising new approach for identifying novel cellulolytic activities, current computational methods are unable to distinguish between cellulases and enzymes with different substrate specificities that belong to the same protein family. We show that by using a robust computational approach supported by experimental studies, cellulases and non-cellulases can be effectively identified within a given protein family. Phylogenetic analysis of GH48 showed non-monophyletic distribution, an indication of horizontal gene transfer. Enzymatic function of GH48 proteins coded by horizontally transferred genes was verified experimentally, which confirmed that these proteins are cellulases. Computational and structural studies of GH48 enzymes identified structural elements that define cellulases and can be used to computationally distinguish them from non-cellulases. We propose that the structural element that can be used for in silico discrimination between cellulases and non-cellulases belonging to GH48 is an ω-loop located on the surface of the molecule and characterized by highly conserved rare amino acids. These markers were used to screen metagenomics data for "true" cellulases.

Engineering ecosystems and synthetic ecologies

Microbial ecosystems play an important role in nature. Engineering these systems for industrial, medical, or biotechnological purposes are important pursuits for synthetic biologists and biological engineers moving forward. Here we provide a review of recent progress in engineering natural and synthetic microbial ecosystems. We highlight important forward engineering design principles, theoretical and quantitative models, new experimental and manipulation tools, and possible applications of microbial ecosystem engineering. We argue that simply engineering individual microbes will lead to fragile homogenous populations that are difficult to sustain, especially in highly heterogeneous and unpredictable environments. Instead, engineered microbial ecosystems are likely to be more robust and able to achieve complex tasks at the spatial and temporal resolution needed for truly programmable biology.

Uncovering the co-evolutionary network among prokaryotic genes

MOTIVATION

Correlated events of gains and losses enable inference of co-evolution relations. The reconstruction of the co-evolutionary interactions network in prokaryotic species may elucidate functional associations among genes.

RESULTS

We developed a novel probabilistic methodology for the detection of co-evolutionary interactions between pairs of genes. Using this method we inferred the co-evolutionary network among 4593 Clusters of Orthologous Genes (COGs). The number of co-evolutionary interactions substantially differed among COGs. Over 40% were found to co-evolve with at least one partner. We partitioned the network of co-evolutionary relations into clusters and uncovered multiple modular assemblies of genes with clearly defined functions. Finally, we measured the extent to which co-evolutionary relations coincide with other cellular relations such as genomic proximity, gene fusion propensity, co-expression, protein-protein interactions and metabolic connections. Our results show that co-evolutionary relations only partially overlap with these other types of networks. Our results suggest that the inferred co-evolutionary network in prokaryotes is highly informative towards revealing functional relations among genes, often showing signals that cannot be extracted from other network types.

AVAILABILITY AND IMPLEMENTATION

Available under GPL license as open source.

CONTACT

talp@post.tau.ac.il.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

The Human Microbiome Project: a community resource for the healthy human microbiome

This manuscript describes the NIH Human Microbiome Project, including a brief review of human microbiome research, a history of the project, and a comprehensive overview of the consortium's recent collection of publications analyzing the human microbiome.

Metagenomics and its connection to microbial community organization

Microbes dominate most global biogeochemical cycles, and microbial metagenomics (studying the collective microbial genomes) provides invaluable new insights into microbial systems, independent of cultivation. Metagenomic approaches targeting specific genes, e.g. small subunit (ssu) ribosomal RNA (rRNA), can be used to investigate microbial community organization by efficiently showing which taxa of organisms are present, while shotgun approaches show all genes and can indicate what functions the organisms are capable of. But collecting and organizing comprehensive shotgun data is extremely challenging and costly, and, in theory, predicting functionalities from microbial identities alone would save immense effort. However, we don’t yet know to what extent such predictions are applicable.

Sequencing of seven haloarchaeal genomes reveals patterns of genomic flux

We report the sequencing of seven genomes from two haloarchaeal genera, Haloferax and Haloarcula. Ease of cultivation and the existence of well-developed genetic and biochemical tools for several diverse haloarchaeal species make haloarchaea a model group for the study of archaeal biology. The unique physiological properties of these organisms also make them good candidates for novel enzyme discovery for biotechnological applications. Seven genomes were sequenced to ∼20×coverage and assembled to an average of 50 contigs (range 5 scaffolds-168 contigs). Comparisons of protein-coding gene compliments revealed large-scale differences in COG functional group enrichment between these genera. Analysis of genes encoding machinery for DNA metabolism reveals genera-specific expansions of the general transcription factor TATA binding protein as well as a history of extensive duplication and horizontal transfer of the proliferating cell nuclear antigen. Insights gained from this study emphasize the importance of haloarchaea for investigation of archaeal biology.

Inability to catabolize galactose leads to increased ability to compete for nodule occupancy in Sinorhizobium meliloti

A mutant unable to utilize galactose was isolated in Sinorhizobium meliloti strain Rm1021. The mutation was found to be in a gene annotated dgoK1, a putative 2-keto-3-deoxygalactonokinase. The genetic region was isolated on a complementing cosmid and subsequently characterized. Based on genetic and bioinformatic evidence, the locus encodes all five enzymes (galD, dgoK, dgoA, SMc00883, and ilvD1) involved in the De Ley-Doudoroff pathway for galactose catabolism. Although all five genes are present, genetic analysis suggests that the galactonase (SMc00883) and the dehydratase (ilvD1) are dispensable with respect to the ability to catabolize galactose. In addition, we show that the transport of galactose is partially facilitated by the arabinose transporter (AraABC) and that both glucose and galactose compete with arabinose for transport. Quantitative reverse transcription-PCR (qRT-PCR) data show that in a dgoK background, the galactose locus is constitutively expressed, and the induction of the ara locus seems to be enhanced. Assays of competition for nodule occupancy show that the inability to catabolize galactose is correlated with an increased ability to compete for nodule occupancy.

The Biological Observation Matrix (BIOM) format or: how I learned to stop worrying and love the ome-ome

BACKGROUND

We present the Biological Observation Matrix (BIOM, pronounced "biome") format: a JSON-based file format for representing arbitrary observation by sample contingency tables with associated sample and observation metadata. As the number of categories of comparative omics data types (collectively, the "ome-ome") grows rapidly, a general format to represent and archive this data will facilitate the interoperability of existing bioinformatics tools and future meta-analyses.

FINDINGS

The BIOM file format is supported by an independent open-source software project (the biom-format project), which initially contains Python objects that support the use and manipulation of BIOM data in Python programs, and is intended to be an open development effort where developers can submit implementations of these objects in other programming languages.

CONCLUSIONS

The BIOM file format and the biom-format project are steps toward reducing the "bioinformatics bottleneck" that is currently being experienced in diverse areas of biological sciences, and will help us move toward the next phase of comparative omics where basic science is translated into clinical and environmental applications. The BIOM file format is currently recognized as an Earth Microbiome Project Standard, and as a Candidate Standard by the Genomic Standards Consortium.

Genome sequence of the mesophilic Thermotogales bacterium Mesotoga prima MesG1.Ag.4.2 reveals the largest Thermotogales genome to date

Here we describe the genome of Mesotoga prima MesG1.Ag4.2, the first genome of a mesophilic Thermotogales bacterium. Mesotoga prima was isolated from a polychlorinated biphenyl (PCB)-dechlorinating enrichment culture from Baltimore Harbor sediments. Its 2.97 Mb genome is considerably larger than any previously sequenced Thermotogales genomes, which range between 1.86 and 2.30 Mb. This larger size is due to both higher numbers of protein-coding genes and larger intergenic regions. In particular, the M. prima genome contains more genes for proteins involved in regulatory functions, for instance those involved in regulation of transcription. Together with its closest relative, Kosmotoga olearia, it also encodes different types of proteins involved in environmental and cell-cell interactions as compared with other Thermotogales bacteria. Amino acid composition analysis of M. prima proteins implies that this lineage has inhabited low-temperature environments for a long time. A large fraction of the M. prima genome has been acquired by lateral gene transfer (LGT): a DarkHorse analysis suggests that 766 (32%) of predicted protein-coding genes have been involved in LGT after Mesotoga diverged from the other Thermotogales lineages. A notable example of a lineage-specific LGT event is a reductive dehalogenase gene-a key enzyme in dehalorespiration, indicating M. prima may have a more active role in PCB dechlorination than was previously assumed.

IMG/M-HMP: a metagenome comparative analysis system for the Human Microbiome Project

The Integrated Microbial Genomes and Metagenomes (IMG/M) resource is a data management system that supports the analysis of sequence data from microbial communities in the integrated context of all publicly available draft and complete genomes from the three domains of life as well as a large number of plasmids and viruses. IMG/M currently contains thousands of genomes and metagenome samples with billions of genes. IMG/M-HMP is an IMG/M data mart serving the US National Institutes of Health (NIH) Human Microbiome Project (HMP), focussed on HMP generated metagenome datasets, and is one of the central resources provided from the HMP Data Analysis and Coordination Center (DACC). IMG/M-HMP is available at http://www.hmpdacc-resources.org/imgm_hmp/.

Next-generation and whole-genome sequencing in the diagnostic clinical microbiology laboratory

The identification and/or prediction of the antimicrobial resistance of microorganisms in clinical materials solely by molecular means in the diagnostic microbiology laboratory is not novel. However, the ability to sequence multitudes of bacterial genomes and deliver and interpret the resultant sequence information in near "real-time" is the basis of next-generation sequencing (NGS) technologies. There have been numerous applications and successes of NGS applications in the clinical and public health domain. However, none have, as yet, delivered perhaps the most sought after application, i.e., the generation of microbial sequence data for "real-time" patient management. In this review, we discuss the use of NGS and whole-genome sequencing (WGS) of microorganisms as a logical next step for the routine diagnosis of infection and the prediction of antimicrobial susceptibility in the clinical microbiology laboratory.

Unlocking the potential of metagenomics through replicated experimental design

Metagenomics holds enormous promise for discovering novel enzymes and organisms that are biomarkers or drivers of processes relevant to disease, industry and the environment. In the past two years, we have seen a paradigm shift in metagenomics to the application of cross-sectional and longitudinal studies enabled by advances in DNA sequencing and high-performance computing. These technologies now make it possible to broadly assess microbial diversity and function, allowing systematic investigation of the largely unexplored frontier of microbial life. To achieve this aim, the global scientific community must collaborate and agree upon common objectives and data standards to enable comparative research across the Earth's microbiome. Improvements in comparability of data will facilitate the study of biotechnologically relevant processes, such as bioprospecting for new glycoside hydrolases or identifying novel energy sources.

Phenotypic and transcriptomic response of auxotrophic Mycobacterium avium subsp. paratuberculosis leuD mutant under environmental stress

Mycobacterium avium subsp. paratuberculosis (MAP) is the causative agent of severe gastroenteritis in cattle. To gain a better understanding of MAP virulence, we investigated the role of leuD gene in MAP metabolism and stress response. For this, we have constructed an auxotrophic strain of MAP by deleting the leuD gene using allelic exchange. The wildtype and mutant strains were then compared for metabolic phenotypic changes using Biolog phenotype microarrays. The responses of both strains to physiologically relevant stress conditions were assessed using DNA microarrays. Transcriptomic data was then analyzed in the context of cellular metabolic pathways and gene networks. Our results showed that deletion of leuD gene has a global effect on both MAP phenotypic and transcriptome response. At the metabolic level, the mutant strain lost the ability to utilize most of the carbon, nitrogen, sulphur, phosphorus and nutrient supplements as energy source. At the transcriptome level, more than 100 genes were differentially expressed in each of the stress condition tested. Systems level network analysis revealed that the differentially expressed genes were distributed throughout the gene network, thus explaining the global impact of leuD deletion in metabolic phenotype. Further, we find that leuD deletion impacted metabolic pathways associated with fatty acids. We verified this by experimentally estimating the total fatty acid content of both mutant and wildtype. The mutant strain had 30% less fatty acid content when compared to wildtype, thus supporting the results from transcriptional and computational analyses. Our results therefore reveal the intricate connection between the metabolism and virulence in MAP.

Caldicellulosiruptor core and pangenomes reveal determinants for noncellulosomal thermophilic deconstruction of plant biomass

Extremely thermophilic bacteria of the genus Caldicellulosiruptor utilize carbohydrate components of plant cell walls, including cellulose and hemicellulose, facilitated by a diverse set of glycoside hydrolases (GHs). From a biofuel perspective, this capability is crucial for deconstruction of plant biomass into fermentable sugars. While all species from the genus grow on xylan and acid-pretreated switchgrass, growth on crystalline cellulose is variable. The basis for this variability was examined using microbiological, genomic, and proteomic analyses of eight globally diverse Caldicellulosiruptor species. The open Caldicellulosiruptor pangenome (4,009 open reading frames [ORFs]) encodes 106 GHs, representing 43 GH families, but only 26 GHs from 17 families are included in the core (noncellulosic) genome (1,543 ORFs). Differentiating the strongly cellulolytic Caldicellulosiruptor species from the others is a specific genomic locus that encodes multidomain cellulases from GH families 9 and 48, which are associated with cellulose-binding modules. This locus also encodes a novel adhesin associated with type IV pili, which was identified in the exoproteome bound to crystalline cellulose. Taking into account the core genomes, pangenomes, and individual genomes, the ancestral Caldicellulosiruptor was likely cellulolytic and evolved, in some cases, into species that lost the ability to degrade crystalline cellulose while maintaining the capacity to hydrolyze amorphous cellulose and hemicellulose.

MicrobeDB: a locally maintainable database of microbial genomic sequences

Summary: Analysis of microbial genomes often requires the general organization and comparison of tens to thousands of genomes both from public repositories and unpublished sources. MicrobeDB provides a foundation for such projects by the automation of downloading published, completed bacterial and archaeal genomes from key sources, parsing annotations of all genomes (both public and private) into a local database, and allowing interaction with the database through an easy to use programming interface. MicrobeDB creates a simple to use, easy to maintain, centralized local resource for various large-scale comparative genomic analyses and a back-end for future microbial application design.

Availability: MicrobeDB is freely available under the GNU-GPL at: http://github.com/mlangill/microbedb/

Contact:morgan.g.i.langille@gmail.com

Informatic strategies for the discovery of polyketides and nonribosomal peptides

A modern challenge and opportunity exists for in the ability to link genomic and metabolomic data, using novel informatic methods to find new bioactive natural products.