IMG: the Integrated Microbial Genomes database and comparative analysis system

Microbial life in a fjord: metagenomic analysis of a microbial mat in Chilean patagonia

The current study describes the taxonomic and functional composition of metagenomic sequences obtained from a filamentous microbial mat isolated from the Comau fjord, located in the northernmost part of the Chilean Patagonia. The taxonomic composition of the microbial community showed a high proportion of members of the Gammaproteobacteria, including a high number of sequences that were recruited to the genomes of Moritella marina MP-1 and Colwelliapsycherythraea 34H, suggesting the presence of populations related to these two psychrophilic bacterial species. Functional analysis of the community indicated a high proportion of genes coding for the transport and metabolism of amino acids, as well as in energy production. Among the energy production functions, we found protein-coding genes for sulfate and nitrate reduction, both processes associated with Gammaproteobacteria-related sequences. This report provides the first examination of the taxonomic composition and genetic diversity associated with these conspicuous microbial mat communities and provides a framework for future microbial studies in the Comau fjord.

Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences

Profiling phylogenetic marker genes, such as the 16S rRNA gene, is a key tool for studies of microbial communities but does not provide direct evidence of a community’s functional capabilities. Here we describe PICRUSt (Phylogenetic Investigation of Communities by Reconstruction of Unobserved States), a computational approach to predict the functional composition of a metagenome using marker gene data and a database of reference genomes. PICRUSt uses an extended ancestral-state reconstruction algorithm to predict which gene families are present and then combines gene families to estimate the composite metagenome. Using 16S information, PICRUSt recaptures key findings from the Human Microbiome Project and accurately predicts the abundance of gene families in host-associated and environmental communities, with quantifiable uncertainty. Our results demonstrate that phylogeny and function are sufficiently linked that this ‘predictive metagenomic’ approach should provide useful insights into the thousands of uncultivated microbial communities for which only marker gene surveys are currently available.

Integrating metagenomic and amplicon databases to resolve the phylogenetic and ecological diversity of the Chlamydiae

In the era of metagenomics and amplicon sequencing, comprehensive analyses of available sequence data remain a challenge. Here we describe an approach exploiting metagenomic and amplicon data sets from public databases to elucidate phylogenetic diversity of defined microbial taxa. We investigated the phylum Chlamydiae whose known members are obligate intracellular bacteria that represent important pathogens of humans and animals, as well as symbionts of protists. Despite their medical relevance, our knowledge about chlamydial diversity is still scarce. Most of the nine known families are represented by only a few isolates, while previous clone library-based surveys suggested the existence of yet uncharacterized members of this phylum. Here we identified more than 22,000 high quality, non-redundant chlamydial 16S rRNA gene sequences in diverse databases, as well as 1900 putative chlamydial protein-encoding genes. Even when applying the most conservative approach, clustering of chlamydial 16S rRNA gene sequences into operational taxonomic units revealed an unexpectedly high species, genus and family-level diversity within the Chlamydiae, including 181 putative families. These in silico findings were verified experimentally in one Antarctic sample, which contained a high diversity of novel Chlamydiae. In our analysis, the Rhabdochlamydiaceae, whose known members infect arthropods, represents the most diverse and species-rich chlamydial family, followed by the protist-associated Parachlamydiaceae, and a putative new family (PCF8) with unknown host specificity. Available information on the origin of metagenomic samples indicated that marine environments contain the majority of the newly discovered chlamydial lineages, highlighting this environment as an important chlamydial reservoir.

Structure and evolution of chlorate reduction composite transposons

The genes for chlorate reduction in six bacterial strains were analyzed in order to gain insight into the metabolism. A newly isolated chlorate-reducing bacterium (Shewanella algae ACDC) and three previously isolated strains (Ideonella dechloratans, Pseudomonas sp. strain PK, and Dechloromarinus chlorophilus NSS) were genome sequenced and compared to published sequences (Alicycliphilus denitrificans BC plasmid pALIDE01 and Pseudomonas chloritidismutans AW-1). De novo assembly of genomes failed to join regions adjacent to genes involved in chlorate reduction, suggesting the presence of repeat regions. Using a bioinformatics approach and finishing PCRs to connect fragmented contigs, we discovered that chlorate reduction genes are flanked by insertion sequences, forming composite transposons in all four newly sequenced strains. These insertion sequences delineate regions with the potential to move horizontally and define a set of genes that may be important for chlorate reduction. In addition to core metabolic components, we have highlighted several such genes through comparative analysis and visualization. Phylogenetic analysis places chlorate reductase within a functionally diverse clade of type II dimethyl sulfoxide (DMSO) reductases, part of a larger family of enzymes with reactivity toward chlorate. Nucleotide-level forensics of regions surrounding chlorite dismutase (cld), as well as its phylogenetic clustering in a betaproteobacterial Cld clade, indicate that cld has been mobilized at least once from a perchlorate reducer to build chlorate respiration.

IMPORTANCE

Genome sequencing has identified, for the first time, chlorate reduction composite transposons. These transposons are constructed with flanking insertion sequences that differ in type and orientation between organisms, indicating that this mobile element has formed multiple times and is important for dissemination. Apart from core metabolic enzymes, very little is known about the genetic factors involved in chlorate reduction. Comparative analysis has identified several genes that may also be important, but the relative absence of accessory genes suggests that this mobile metabolism relies on host systems for electron transport, regulation, and cofactor synthesis. Phylogenetic analysis of Cld and ClrA provides support for the hypothesis that chlorate reduction was built multiple times from type II dimethyl sulfoxide (DMSO) reductases and cld. In at least one case, cld has been coopted from a perchlorate reduction island for this purpose. This work is a significant step toward understanding the genetics and evolution of chlorate reduction.

Towards a predictive systems-level model of the human microbiome: progress, challenges, and opportunities

The human microbiome represents a vastly complex ecosystem that is tightly linked to our development, physiology, and health. Our increased capacity to generate multiple channels of omic data from this system, brought about by recent advances in high throughput molecular technologies, calls for the development of systems-level methods and models that take into account not only the composition of genes and species in a microbiome but also the interactions between these components. Such models should aim to study the microbiome as a community of species whose metabolisms are tightly intertwined with each other and with that of the host, and should be developed with a view towards an integrated, comprehensive, and predictive modeling framework. Here, we review recent work specifically in metabolic modeling of the human microbiome, highlighting both novel methodologies and pressing challenges. We discuss various modeling approaches that lay the foundation for a full-scale predictive model, focusing on models of interactions between microbial species, metagenome-scale models of community-level metabolism, and models of the interaction between the microbiome and the host. Continued development of such models and of their integration into a multi-scale model of the microbiome will lead to a deeper mechanistic understanding of how variation in the microbiome impacts the host, and will promote the discovery of clinically relevant and ecologically relevant insights from the rich trove of data now available.

Proteogenomic analysis of a thermophilic bacterial consortium adapted to deconstruct switchgrass

Thermophilic bacteria are a potential source of enzymes for the deconstruction of lignocellulosic biomass. However, the complement of proteins used to deconstruct biomass and the specific roles of different microbial groups in thermophilic biomass deconstruction are not well-explored. Here we report on the metagenomic and proteogenomic analyses of a compost-derived bacterial consortium adapted to switchgrass at elevated temperature with high levels of glycoside hydrolase activities. Near-complete genomes were reconstructed for the most abundant populations, which included composite genomes for populations closely related to sequenced strains of Thermus thermophilus and Rhodothermus marinus, and for novel populations that are related to thermophilic Paenibacilli and an uncultivated subdivision of the little-studied Gemmatimonadetes phylum. Partial genomes were also reconstructed for a number of lower abundance thermophilic Chloroflexi populations. Identification of genes for lignocellulose processing and metabolic reconstructions suggested Rhodothermus, Paenibacillus and Gemmatimonadetes as key groups for deconstructing biomass, and Thermus as a group that may primarily metabolize low molecular weight compounds. Mass spectrometry-based proteomic analysis of the consortium was used to identify >3000 proteins in fractionated samples from the cultures, and confirmed the importance of Paenibacillus and Gemmatimonadetes to biomass deconstruction. These studies also indicate that there are unexplored proteins with important roles in bacterial lignocellulose deconstruction.

Metabolic modeling of species interaction in the human microbiome elucidates community-level assembly rules

The human microbiome plays a key role in human health and is associated with numerous diseases. Metagenomic-based studies are now generating valuable information about the composition of the microbiome in health and in disease, demonstrating nonneutral assembly processes and complex co-occurrence patterns. However, the underlying ecological forces that structure the microbiome are still unclear. Specifically, compositional studies alone with no information about mechanisms of interaction, potential competition, or syntrophy, cannot clearly distinguish habitat-filtering and species assortment assembly processes. To address this challenge, we introduce a computational framework, integrating metagenomic-based compositional data with genome-scale metabolic modeling of species interaction. We use in silico metabolic network models to predict levels of competition and complementarity among 154 microbiome species and compare predicted interaction measures to species co-occurrence. Applying this approach to two large-scale datasets describing the composition of the gut microbiome, we find that species tend to co-occur across individuals more frequently with species with which they strongly compete, suggesting that microbiome assembly is dominated by habitat filtering. Moreover, species' partners and excluders exhibit distinct metabolic interaction levels. Importantly, we show that these trends cannot be explained by phylogeny alone and hold across multiple taxonomic levels. Interestingly, controlling for host health does not change the observed patterns, indicating that the axes along which species are filtered are not fully defined by macroecological host states. The approach presented here lays the foundation for a reverse-ecology framework for addressing key questions concerning the assembly of host-associated communities and for informing clinical efforts to manipulate the microbiome.

Single-cell genomics reveals complex carbohydrate degradation patterns in poribacterial symbionts of marine sponges

Many marine sponges are hosts to dense and phylogenetically diverse microbial communities that are located in the extracellular matrix of the animal. The candidate phylum Poribacteria is a predominant member of the sponge microbiome and its representatives are nearly exclusively found in sponges. Here we used single-cell genomics to obtain comprehensive insights into the metabolic potential of individual poribacterial cells representing three distinct phylogenetic groups within Poribacteria. Genome sizes were up to 5.4 Mbp and genome coverage was as high as 98.5%. Common features of the poribacterial genomes indicated that heterotrophy is likely to be of importance for this bacterial candidate phylum. Carbohydrate-active enzyme database screening and further detailed analysis of carbohydrate metabolism suggested the ability to degrade diverse carbohydrate sources likely originating from seawater and from the host itself. The presence of uronic acid degradation pathways as well as several specific sulfatases provides strong support that Poribacteria degrade glycosaminoglycan chains of proteoglycans, which are important components of the sponge host matrix. Dominant glycoside hydrolase families further suggest degradation of other glycoproteins in the host matrix. We therefore propose that Poribacteria are well adapted to an existence in the sponge extracellular matrix. Poribacteria may be viewed as efficient scavengers and recyclers of a particular suite of carbon compounds that are unique to sponges as microbial ecosystems.

A cell wall recycling shortcut that bypasses peptidoglycan de novo biosynthesis

We report a salvage pathway in Gram-negative bacteria that bypasses de novo biosynthesis of UDP N-acetylmuramic acid (UDP-MurNAc), the first committed peptidoglycan precursor, and thus provides a rationale for intrinsic fosfomycin resistance. The anomeric sugar kinase AmgK and the MurNAc α-1-phosphate uridylyl transferase MurU, defining this new cell wall sugar-recycling route in Pseudomonas putida, were characterized and engineered into Escherichia coli, channeling external MurNAc directly to peptidoglycan biosynthesis.

Quantification of endospore-forming firmicutes by quantitative PCR with the functional gene spo0A

Bacterial endospores are highly specialized cellular forms that allow endospore-forming Firmicutes (EFF) to tolerate harsh environmental conditions. EFF are considered ubiquitous in natural environments, in particular, those subjected to stress conditions. In addition to natural habitats, EFF are often the cause of contamination problems in anthropogenic environments, such as industrial production plants or hospitals. It is therefore desirable to assess their prevalence in environmental and industrial fields. To this end, a high-sensitivity detection method is still needed. The aim of this study was to develop and evaluate an approach based on quantitative PCR (qPCR). For this, the suitability of functional genes specific for and common to all EFF were evaluated. Seven genes were considered, but only spo0A was retained to identify conserved regions for qPCR primer design. An approach based on multivariate analysis was developed for primer design. Two primer sets were obtained and evaluated with 16 pure cultures, including representatives of the genera Bacillus, Paenibacillus, Brevibacillus, Geobacillus, Alicyclobacillus, Sulfobacillus, Clostridium, and Desulfotomaculum, as well as with environmental samples. The primer sets developed gave a reliable quantification when tested on laboratory strains, with the exception of Sulfobacillus and Desulfotomaculum. A test using sediment samples with a diverse EFF community also gave a reliable quantification compared to 16S rRNA gene pyrosequencing. A detection limit of about 10(4) cells (or spores) per gram of initial material was calculated, indicating this method has a promising potential for the detection of EFF over a wide range of applications.

An ICEBs1-like element may be associated with the extreme radiation and desiccation resistance of Bacillus pumilus SAFR-032 spores

Comparisons of the genomes of Bacillus pumilus SAFR-032 and the closely related type strain, B. pumilus ATCC7061(T), exposed an extended region of non-homologous genes. A detailed examination of this region revealed the presence of an ICEBs1-like integrative conjugative element in SAFR-032. A similar element was subsequently located elsewhere in the ATCC7061(T) genome. A detailed comparison of these elements and the ICEBs1 of B. subtilis revealed extremely rapid flux in gene content, genome organization and sequence similarity. It is not clear if the B. pumilus elements as they are currently structured are functional. However, it is clear that the past involvement of these elements has brought multiple genes of unknown function to the SAFR-032 genome and these genes may be responsible for the rapid evolution that led to the extreme radiation and desiccation resistance of this organism's spores.

Draft Genome Sequence for Ralstonia sp. Strain OR214, a Bacterium with Potential for Bioremediation

Ralstonia sp. strain OR214 belongs to the class Betaproteobacteria and was isolated from subsurface sediments in Oak Ridge, TN. A member of this genus has been described as a potential bioremediation agent. Strain OR214 is tolerant to various heavy metals, such as uranium, nickel, cobalt, and cadmium. We present its draft genome sequence here.

Draft Genome Sequence for Caulobacter sp. Strain OR37, a Bacterium Tolerant to Heavy Metals

Caulobacter sp. strain OR37 belongs to the class Alphaproteobacteria and was isolated from subsurface sediments in Oak Ridge, TN. Strain OR37 is noteworthy due to its tolerance to high concentrations of heavy metals, such as uranium, nickel, cobalt, and cadmium, and we present its draft genome sequence here.

MS/MS networking guided analysis of molecule and gene cluster families

The ability to correlate the production of specialized metabolites to the genetic capacity of the organism that produces such molecules has become an invaluable tool in aiding the discovery of biotechnologically applicable molecules. Here, we accomplish this task by matching molecular families with gene cluster families, making these correlations to 60 microbes at one time instead of connecting one molecule to one organism at a time, such as how it is traditionally done. We can correlate these families through the use of nanospray desorption electrospray ionization MS/MS, an ambient pressure MS technique, in conjunction with MS/MS networking and peptidogenomics. We matched the molecular families of peptide natural products produced by 42 bacilli and 18 pseudomonads through the generation of amino acid sequence tags from MS/MS data of specific clusters found in the MS/MS network. These sequence tags were then linked to biosynthetic gene clusters in publicly accessible genomes, providing us with the ability to link particular molecules with the genes that produced them. As an example of its use, this approach was applied to two unsequenced Pseudoalteromonas species, leading to the discovery of the gene cluster for a molecular family, the bromoalterochromides, in the previously sequenced strain P. piscicida JCM 20779(T). The approach itself is not limited to 60 related strains, because spectral networking can be readily adopted to look at molecular family-gene cluster families of hundreds or more diverse organisms in one single MS/MS network.

The abundance and variety of carbohydrate-active enzymes in the human gut microbiota

Descriptions of the microbial communities that live on and in the human body have progressed at a spectacular rate over the past 5 years, fuelled primarily by highly parallel DNA-sequencing technologies and associated advances in bioinformatics, and by the expectation that understanding how to manipulate the structure and functions of our microbiota will allow us to affect health and prevent or treat diseases. Among the myriad of genes that have been identified in the human gut microbiome, those that encode carbohydrate-active enzymes (CAZymes) are of particular interest, as these enzymes are required to digest most of our complex repertoire of dietary polysaccharides. In this Analysis article, we examine the carbohydrate-digestive capacity of a simplified but representative mini-microbiome in order to highlight the abundance and variety of bacterial CAZymes that are represented in the human gut microbiota.

Permanent draft genome sequence of Comamonas testosteroni KF-1

Comamonas testosteroni KF-1 is a model organism for the elucidation of the novel biochemical degradation pathways for xenobiotic 4-sulfophenylcarboxylates (SPC) formed during biodegradation of synthetic 4-sulfophenylalkane surfactants (linear alkylbenzenesulfonates, LAS) by bacterial communities. Here we describe the features of this organism, together with the complete genome sequence and annotation. The 6,026,527 bp long chromosome (one sequencing gap) exhibits an average G+C content of 61.79% and is predicted to encode 5,492 protein-coding genes and 114 RNA genes.

Contrasting genomic properties of free-living and particle-attached microbial assemblages within a coastal ecosystem

The Columbia River (CR) is a powerful economic and environmental driver in the US Pacific Northwest. Microbial communities in the water column were analyzed from four diverse habitats: (1) an estuarine turbidity maximum (ETM), (2) a chlorophyll maximum of the river plume, (3) an upwelling-associated hypoxic zone, and (4) the deep ocean bottom. Three size fractions, 0.1-0.8, 0.8-3, and 3-200 μm were collected for each habitat in August 2007, and used for DNA isolation and 454 sequencing, resulting in 12 metagenomes of >5 million reads (>1.6 Gbp). To characterize the dominant microorganisms and metabolisms contributing to coastal biogeochemistry, we used predicted peptide and rRNA data. The 3- and 0.8-μm metagenomes, representing particulate fractions, were taxonomically diverse across habitats. The 3-μm size fractions contained a high abundance of eukaryota with diatoms dominating the hypoxic water and plume, while cryptophytes were more abundant in the ETM. The 0.1-μm metagenomes represented mainly free-living bacteria and archaea. The most abundant archaeal hits were observed in the deep ocean and hypoxic water (19% of prokaryotic peptides in the 0.1-μm metagenomes), and were homologous to Nitrosopumilus maritimus (ammonia-oxidizing Thaumarchaeota). Bacteria dominated metagenomes of all samples. In the euphotic zone (estuary, plume and hypoxic ocean), the most abundant bacterial taxa (≥40% of prokaryotic peptides) represented aerobic photoheterotrophs. In contrast, the low-oxygen, deep water metagenome was enriched with sequences for strict and facultative anaerobes. Interestingly, many of the same anaerobic bacterial families were enriched in the 3-μm size fraction of the ETM (2-10X more abundant relative to the 0.1-μm metagenome), indicating possible formation of anoxic microniches within particles. Results from this study provide a metagenome perspective on ecosystem-scale metabolism in an upwelling-influenced river-dominated coastal margin.

Draft Genome Sequence of the Cellulolytic, Mesophilic, Anaerobic Bacterium Clostridium termitidis Strain CT1112 (DSM 5398)

Here, we report the draft genome sequence of Clostridium termitidis strain CT1112 (DSM 5398), a mesophilic, cellulolytic bacterium that can utilize a variety of sugars, as well as pure cellulose, as a sole carbon source; it also synthesizes fermentation end products with potential industrial applications.

Genome Sequence of Pseudomonas aeruginosa PA45, a Highly Virulent Strain Isolated from a Patient with Bloodstream Infection

Pseudomonas aeruginosa is a ubiquitous opportunistic pathogen causing a broad range of infections in humans. We provide the draft genome sequence of the recently identified and highly virulent P. aeruginosa PA45 strain. Its 6.6-Mb genome contains 6,822 genes, including an unparalleled number of virulence genes, which might explain its aggressive phenotype.

Computational meta'omics for microbial community studies

Complex microbial communities are an integral part of the Earth's ecosystem and of our bodies in health and disease. In the last two decades, culture-independent approaches have provided new insights into their structure and function, with the exponentially decreasing cost of high-throughput sequencing resulting in broadly available tools for microbial surveys. However, the field remains far from reaching a technological plateau, as both computational techniques and nucleotide sequencing platforms for microbial genomic and transcriptional content continue to improve. Current microbiome analyses are thus starting to adopt multiple and complementary meta'omic approaches, leading to unprecedented opportunities to comprehensively and accurately characterize microbial communities and their interactions with their environments and hosts. This diversity of available assays, analysis methods, and public data is in turn beginning to enable microbiome-based predictive and modeling tools. We thus review here the technological and computational meta'omics approaches that are already available, those that are under active development, their success in biological discovery, and several outstanding challenges.

Identification of new protein coding sequences and signal peptidase cleavage sites of Helicobacter pylori strain 26695 by proteogenomics

Correct annotation of protein coding genes is the basis of conventional data analysis in proteomic studies. Nevertheless, most protein sequence databases almost exclusively rely on gene finding software and inevitably also miss protein annotations or possess errors. Proteogenomics tries to overcome these issues by matching MS data directly against a genome sequence database. Here we report an in-depth proteogenomics study of Helicobacter pylori strain 26695. MS data was searched against a combined database of the NCBI annotations and a six-frame translation of the genome. Database searches with Mascot and X! Tandem revealed 1115 proteins identified by at least two peptides with a peptide false discovery rate below 1%. This represents 71% of the predicted proteome. So far this is the most extensive proteome study of Helicobacter pylori. Our proteogenomic approach unambiguously identified four previously missed annotations and furthermore allowed us to correct sequences of six annotated proteins. Since secreted proteins are often involved in pathogenic processes we further investigated signal peptidase cleavage sites. By applying a database search that accommodates the identification of semi-specific cleaved peptides, 63 previously unknown signal peptides were detected. The motif LXA showed to be the predominant recognition sequence for signal peptidases.

BIOLOGICAL SIGNIFICANCE

The results of MS-based proteomic studies highly rely on correct annotation of protein coding genes which is the basis of conventional data analysis. However, the annotation of protein coding sequences in genomic data is usually based on gene finding software. These tools are limited in their prediction accuracy such as the problematic determination of exact gene boundaries. Thus, protein databases own partly erroneous or incomplete sequences. Additionally, some protein sequences might also be missing in the databases. Proteogenomics, a combination of proteomic and genomic data analyses, is well suited to detect previously not annotated proteins and to correct erroneous sequences. For this purpose, the existing database of the investigated species is typically supplemented with a six-frame translation of the genome. Here, we studied the proteome of the major human pathogen Helicobacter pylori that is responsible for many gastric diseases such as duodenal ulcers and gastric cancer. Our in-depth proteomic study highly reliably identified 1115 proteins (FDR<0.01%) by at least two peptides (FDR<1%) which represent 71% of the predicted proteome deposited at NCBI. The proteogenomic data analysis of our data set resulted in the unambiguous identification of four previously missed annotations, the correction of six annotated proteins as well as the detection of 63 previously unknown signal peptides. We have annotated proteins of particular biological interest like the ferrous iron transport protein A, the coiled-coil-rich protein HP0058 and the lipopolysaccharide biosynthesis protein HP0619. For instance, the protein HP0619 could be a drug target for the inhibition of the LPS synthesis pathway. Furthermore it has been proven that the motif "LXA" is the predominant recognition sequence for the signal peptidase I of H. pylori. Signal peptidases are essential enzymes for the viability of bacterial cells and are involved in pathogenesis. Therefore signal peptidases could be novel targets for antibiotics. The inclusion of the corrected and new annotated proteins as well as the information of signal peptide cleavage sites will help in the study of biological pathways involved in pathogenesis or drug response of H. pylori.

The YNP Metagenome Project: Environmental Parameters Responsible for Microbial Distribution in the Yellowstone Geothermal Ecosystem

The Yellowstone geothermal complex contains over 10,000 diverse geothermal features that host numerous phylogenetically deeply rooted and poorly understood archaea, bacteria, and viruses. Microbial communities in high-temperature environments are generally less diverse than soil, marine, sediment, or lake habitats and therefore offer a tremendous opportunity for studying the structure and function of different model microbial communities using environmental metagenomics. One of the broader goals of this study was to establish linkages among microbial distribution, metabolic potential, and environmental variables. Twenty geochemically distinct geothermal ecosystems representing a broad spectrum of Yellowstone hot-spring environments were used for metagenomic and geochemical analysis and included approximately equal numbers of: (1) phototrophic mats, (2) "filamentous streamer" communities, and (3) archaeal-dominated sediments. The metagenomes were analyzed using a suite of complementary and integrative bioinformatic tools, including phylogenetic and functional analysis of both individual sequence reads and assemblies of predominant phylotypes. This volume identifies major environmental determinants of a large number of thermophilic microbial lineages, many of which have not been fully described in the literature nor previously cultivated to enable functional and genomic analyses. Moreover, protein family abundance comparisons and in-depth analyses of specific genes and metabolic pathways relevant to these hot-spring environments reveal hallmark signatures of metabolic capabilities that parallel the distribution of phylotypes across specific types of geochemical environments.

Poles apart: Arctic and Antarctic Octadecabacter strains share high genome plasticity and a new type of xanthorhodopsin

The genus Octadecabacter is a member of the ubiquitous marine Roseobacter clade. The two described species of this genus, Octadecabacter arcticus and Octadecabacter antarcticus, are psychrophilic and display a bipolar distribution. Here we provide the manually annotated and finished genome sequences of the type strains O. arcticus 238 and O. antarcticus 307, isolated from sea ice of the Arctic and Antarctic, respectively. Both genomes exhibit a high genome plasticity caused by an unusually high density and diversity of transposable elements. This could explain the discrepancy between the low genome synteny and high 16S rRNA gene sequence similarity between both strains. Numerous characteristic features were identified in the Octadecabacter genomes, which show indications of horizontal gene transfer and may represent specific adaptations to the habitats of the strains. These include a gene cluster encoding the synthesis and degradation of cyanophycin in O. arcticus 238, which is absent in O. antarcticus 307 and unique among the Roseobacter clade. Furthermore, genes representing a new subgroup of xanthorhodopsins as an adaptation to icy environments are present in both Octadecabacter strains. This new xanthorhodopsin subgroup differs from the previously characterized xanthorhodopsins of Salinibacter ruber and Gloeobacter violaceus in phylogeny, biogeography and the potential to bind 4-keto-carotenoids. Biochemical characterization of the Octadecabacter xanthorhodopsins revealed that they function as light-driven proton pumps.

CRISPR-Cas systems target a diverse collection of invasive mobile genetic elements in human microbiomes

Background

Bacteria and archaea develop immunity against invading genomes by incorporating pieces of the invaders' sequences, called spacers, into a clustered regularly interspaced short palindromic repeats (CRISPR) locus between repeats, forming arrays of repeat-spacer units. When spacers are expressed, they direct CRISPR-associated (Cas) proteins to silence complementary invading DNA. In order to characterize the invaders of human microbiomes, we use spacers from CRISPR arrays that we had previously assembled from shotgun metagenomic datasets, and identify contigs that contain these spacers' targets.

Results

We discover 95,000 contigs that are putative invasive mobile genetic elements, some targeted by hundreds of CRISPR spacers. We find that oral sites in healthy human populations have a much greater variety of mobile genetic elements than stool samples. Mobile genetic elements carry genes encoding diverse functions: only 7% of the mobile genetic elements are similar to known phages or plasmids, although a much greater proportion contain phage- or plasmid-related genes. A small number of contigs share similarity with known integrative and conjugative elements, providing the first examples of CRISPR defenses against this class of element. We provide detailed analyses of a few large mobile genetic elements of various types, and a relative abundance analysis of mobile genetic elements and putative hosts, exploring the dynamic activities of mobile genetic elements in human microbiomes. A joint analysis of mobile genetic elements and CRISPRs shows that protospacer-adjacent motifs drive their interaction network; however, some CRISPR-Cas systems target mobile genetic elements lacking motifs.

Conclusions

We identify a large collection of invasive mobile genetic elements in human microbiomes, an important resource for further study of the interaction between the CRISPR-Cas immune system and invaders.

Phylogeny of bacterial and archaeal genomes using conserved genes: supertrees and supermatrices

Over 3000 microbial (bacterial and archaeal) genomes have been made publically available to date, providing an unprecedented opportunity to examine evolutionary genomic trends and offering valuable reference data for a variety of other studies such as metagenomics. The utility of these genome sequences is greatly enhanced when we have an understanding of how they are phylogenetically related to each other. Therefore, we here describe our efforts to reconstruct the phylogeny of all available bacterial and archaeal genomes. We identified 24, single-copy, ubiquitous genes suitable for this phylogenetic analysis. We used two approaches to combine the data for the 24 genes. First, we concatenated alignments of all genes into a single alignment from which a Maximum Likelihood (ML) tree was inferred using RAxML. Second, we used a relatively new approach to combining gene data, Bayesian Concordance Analysis (BCA), as implemented in the BUCKy software, in which the results of 24 single-gene phylogenetic analyses are used to generate a "primary concordance" tree. A comparison of the concatenated ML tree and the primary concordance (BUCKy) tree reveals that the two approaches give similar results, relative to a phylogenetic tree inferred from the 16S rRNA gene. After comparing the results and the methods used, we conclude that the current best approach for generating a single phylogenetic tree, suitable for use as a reference phylogeny for comparative analyses, is to perform a maximum likelihood analysis of a concatenated alignment of conserved, single-copy genes.

Complete genome sequence of Halorhodospira halophila SL1

Halorhodospira halophila is among the most halophilic organisms known. It is an obligately photosynthetic and anaerobic purple sulfur bacterium that exhibits autotrophic growth up to saturated NaCl concentrations. The type strain H. halophila SL1 was isolated from a hypersaline lake in Oregon. Here we report the determination of its entire genome in a single contig. This is the first genome of a phototrophic extreme halophile. The genome consists of 2,678,452 bp, encoding 2,493 predicted genes as determined by automated genome annotation. Of the 2,407 predicted proteins, 1,905 were assigned to a putative function. Future detailed analysis of this genome promises to yield insights into the halophilic adaptations of this organism, its ability for photoautotrophic growth under extreme conditions, and its characteristic sulfur metabolism.

Draft Genome Sequence for Desulfovibrio africanus Strain PCS

Desulfovibrio africanus strain PCS is an anaerobic sulfate-reducing bacterium (SRB) isolated from sediment from Paleta Creek, San Diego, CA. Strain PCS is capable of reducing metals such as Fe(III) and Cr(VI), has a cell cycle, and is predicted to produce methylmercury. We present the D. africanus PCS genome sequence.

Shifts in microbial community composition and function in the acidification of a lead/zinc mine tailings

In an attempt to link the microbial community composition and function in mine tailings to the generation of acid mine drainage, we simultaneously explored the geochemistry and microbiology of six tailings collected from a lead/zinc mine, i.e. primary tailings (T1), slightly acidic tailings (T2), extremely acidic tailings (T3, T4 and T5) and orange-coloured oxidized tailings (T6). Geochemical results showed that the six tailings (from T1 to T6) likely represented sequential stages of the acidification process of the mine tailings. 16S rRNA pyrosequencing revealed a contrasting microbial composition between the six tailings: Proteobacteria-related sequences dominated T1-T3 with relative abundance ranging from 56 to 93%, whereas Ferroplasma-related sequences dominated T4-T6 with relative abundance ranging from 28 to 58%. Furthermore, metagenomic analysis of the microbial communities of T2 and T6 indicated that the genes encoding key enzymes for microbial carbon fixation, nitrogen fixation and sulfur oxidation in T2 were largely from Thiobacillus and Acidithiobacillus, Methylococcus capsulatus, and Thiobacillus denitrificans respectively; while those in T6 were mostly identified in Acidithiobacillus and Leptospirillum, Acidithiobacillus and Leptospirillum, and Acidithiobacillus respectively. The microbial communities in T2 and T6 harboured more genes suggesting diverse metabolic capacities for sulfur oxidation/heavy metal detoxification and tolerating low pH respectively.

MGcV: the microbial genomic context viewer for comparative genome analysis

Background

Conserved gene context is used in many types of comparative genome analyses. It is used to provide leads on gene function, to guide the discovery of regulatory sequences, but also to aid in the reconstruction of metabolic networks. We present the Microbial Genomic context Viewer (MGcV), an interactive, web-based application tailored to strengthen the practice of manual comparative genome context analysis for bacteria.

Results

MGcV is a versatile, easy-to-use tool that renders a visualization of the genomic context of any set of selected genes, genes within a phylogenetic tree, genomic segments, or regulatory elements. It is tailored to facilitate laborious tasks such as the interactive annotation of gene function, the discovery of regulatory elements, or the sequence-based reconstruction of gene regulatory networks. We illustrate that MGcV can be used in gene function annotation by visually integrating information on prokaryotic genes, like their annotation as available from NCBI with other annotation data such as Pfam domains, sub-cellular location predictions and gene-sequence characteristics such as GC content. We also illustrate the usefulness of the interactive features that allow the graphical selection of genes to facilitate data gathering (e.g. upstream regions, ID’s or annotation), in the analysis and reconstruction of transcription regulation. Moreover, putative regulatory elements and their corresponding scores or data from RNA-seq and microarray experiments can be uploaded, visualized and interpreted in (ranked-) comparative context maps. The ranked maps allow the interpretation of predicted regulatory elements and experimental data in light of each other.

Conclusion

MGcV advances the manual comparative analysis of genes and regulatory elements by providing fast and flexible integration of gene related data combined with straightforward data retrieval. MGcV is available at http://mgcv.cmbi.ru.nl.

Aspects of nitrogen-fixing Actinobacteria, in particular free-living and symbiotic Frankia

Studies of nitrogen-fixing properties among the Gram-positive Actinobacteria revealed that some species of Arthrobacter, Agromyces, Corynebacterium, Mycobacterium, Micromonospora, Propionibacteria and Streptomyces have nitrogen-fixing capacity. This is also valid for Frankia that fix nitrogen both in free-living and in symbiotic conditions. Frankia symbiosis results from interaction between the Frankia bacteria and dicotyledonous plants, that is, actinorhiza. These plants, which are important in forestry and agroforestry, form, together with the legumes (Fabales), a single nitrogen-fixing clade. It has been shown that a receptor-like kinase gene, SymRK, is necessary for nodulation in actinorhizal plants as well as in legumes and arbuscular mycorrhizal fungi. Recently, the involvement of isoflavonoids as signal molecules during nodulation of an actinorhizal plant was shown. The genome sizes of three Frankia species, Frankia EANpec, ACN14a and CcI3, are different, revealing a relationship between genome size and geographical distribution. Recent genomic sequencing data of Frankia represent genomes from cluster I to IV, indicating that the genome of DgI is one of the smallest genomes in Frankia. In addition, nonsymbiotic Frankiales such as Acidothermus cellulolyticus, Blastococcus saxoobsidens, Geodermatophilus obscurus and Modestobacter marinus have a variety of genome sizes ranging from 2.4 to 5.57 Mb.

Multiple single-cell genomes provide insight into functions of uncultured Deltaproteobacteria in the human oral cavity

Despite a long history of investigation, many bacteria associated with the human oral cavity have yet to be cultured. Studies that correlate the presence or abundance of uncultured species with oral health or disease highlight the importance of these community members. Thus, we sequenced several single-cell genomic amplicons from Desulfobulbus and Desulfovibrio (class Deltaproteobacteria) to better understand their function within the human oral community and their association with periodontitis, as well as other systemic diseases. Genomic data from oral Desulfobulbus and Desulfovibrio species were compared to other available deltaproteobacterial genomes, including from a subset of host-associated species. While both groups share a large number of genes with other environmental Deltaproteobacteria genomes, they encode a wide array of unique genes that appear to function in survival in a host environment. Many of these genes are similar to virulence and host adaptation factors of known human pathogens, suggesting that the oral Deltaproteobacteria have the potential to play a role in the etiology of periodontal disease.

Metabolic model for the filamentous 'Candidatus Microthrix parvicella' based on genomic and metagenomic analyses

'Candidatus Microthrix parvicella' is a lipid-accumulating, filamentous bacterium so far found only in activated sludge wastewater treatment plants, where it is a common causative agent of sludge separation problems. Despite attracting considerable interest, its detailed physiology is still unclear. In this study, the genome of the RN1 strain was sequenced and annotated, which facilitated the construction of a theoretical metabolic model based on available in situ and axenic experimental data. This model proposes that under anaerobic conditions, this organism accumulates preferentially long-chain fatty acids as triacylglycerols. Utilisation of trehalose and/or polyphosphate stores or partial oxidation of long-chain fatty acids may supply the energy required for anaerobic lipid uptake and storage. Comparing the genome sequence of this isolate with metagenomes from two full-scale wastewater treatment plants with enhanced biological phosphorus removal reveals high similarity, with few metabolic differences between the axenic and the dominant community 'Ca. M. parvicella' strains. Hence, the metabolic model presented in this paper could be considered generally applicable to strains in full-scale treatment systems. The genomic information obtained here will provide the basis for future research into in situ gene expression and regulation. Such information will give substantial insight into the ecophysiology of this unusual and biotechnologically important filamentous bacterium.

Genome mining for methanobactins

Background

Methanobactins (Mbns) are a family of copper-binding natural products involved in copper uptake by methanotrophic bacteria. The few Mbns that have been structurally characterized feature copper coordination by two nitrogen-containing heterocycles next to thioamide groups embedded in a peptidic backbone of varying composition. Mbns are proposed to derive from post-translational modification of ribosomally synthesized peptides, but only a few genes encoding potential precursor peptides have been identified. Moreover, the relevance of neighboring genes in these genomes has been unclear.

Results

The potential for Mbn production in a wider range of bacterial species was assessed by mining microbial genomes. Operons encoding Mbn-like precursor peptides, MbnAs, were identified in 16 new species, including both methanotrophs and, surprisingly, non-methanotrophs. Along with MbnA, the core of the operon is formed by two putative biosynthetic genes denoted MbnB and MbnC. The species can be divided into five groups on the basis of their MbnA and MbnB sequences and their operon compositions. Additional biosynthetic proteins, including aminotransferases, sulfotransferases and flavin adenine dinucleotide (FAD)-dependent oxidoreductases were also identified in some families. Beyond biosynthetic machinery, a conserved set of transporters was identified, including MATE multidrug exporters and TonB-dependent transporters. Additional proteins of interest include a di-heme cytochrome c peroxidase and a partner protein, the roles of which remain a mystery.

Conclusions

This study indicates that Mbn-like compounds may be more widespread than previously thought, but are not present in all methanotrophs. This distribution of species suggests a broader role in metal homeostasis. These data provide a link between precursor peptide sequence and Mbn structure, facilitating predictions of new Mbn structures and supporting a post-translational modification biosynthetic pathway. In addition, testable models for Mbn transport and for methanotrophic copper regulation have emerged. Given the unusual modifications observed in Mbns characterized thus far, understanding the roles of the putative biosynthetic proteins is likely to reveal novel pathways and chemistry.

A metagenomic insight into freshwater methane-utilizing communities and evidence for cooperation between the Methylococcaceae and the Methylophilaceae

We investigated microbial communities active in methane oxidation in lake sediment at different oxygen tensions and their response to the addition of nitrate, via stable isotope probing combined with deep metagenomic sequencing. Communities from a total of four manipulated microcosms were analyzed, supplied with ¹³C-methane in, respectively, ambient air, ambient air with the addition of nitrate, nitrogen atmosphere and nitrogen atmosphere with the addition of nitrate, and these were compared to the community from an unamended sediment sample. We found that the major group involved in methane oxidation in both aerobic and microaerobic conditions were members of the family Methylococcaceae, dominated by species of the genus Methylobacter, and these were stimulated by nitrate in aerobic but not microaerobic conditions. In aerobic conditions, we also noted a pronounced response to both methane and nitrate by members of the family Methylophilaceae that are non-methane-oxidizing methylotrophs, and predominantly by the members of the genus Methylotenera. The relevant abundances of the Methylococcaceae and the Methylophilaceae and their coordinated response to methane and nitrate suggest that these species may be engaged in cooperative behavior, the nature of which remains unknown.

De novo prediction of the genomic components and capabilities for microbial plant biomass degradation from (meta-)genomes

BACKGROUND

Understanding the biological mechanisms used by microorganisms for plant biomass degradation is of considerable biotechnological interest. Despite of the growing number of sequenced (meta)genomes of plant biomass-degrading microbes, there is currently no technique for the systematic determination of the genomic components of this process from these data.

RESULTS

We describe a computational method for the discovery of the protein domains and CAZy families involved in microbial plant biomass degradation. Our method furthermore accurately predicts the capability to degrade plant biomass for microbial species from their genome sequences. Application to a large, manually curated data set of microbial degraders and non-degraders identified gene families of enzymes known by physiological and biochemical tests to be implicated in cellulose degradation, such as GH5 and GH6. Additionally, genes of enzymes that degrade other plant polysaccharides, such as hemicellulose, pectins and oligosaccharides, were found, as well as gene families which have not previously been related to the process. For draft genomes reconstructed from a cow rumen metagenome our method predicted Bacteroidetes-affiliated species and a relative to a known plant biomass degrader to be plant biomass degraders. This was supported by the presence of genes encoding enzymatically active glycoside hydrolases in these genomes.

CONCLUSIONS

Our results show the potential of the method for generating novel insights into microbial plant biomass degradation from (meta-)genome data, where there is an increasing production of genome assemblages for uncultured microbes.

Improving microbial genome annotations in an integrated database context

Effective comparative analysis of microbial genomes requires a consistent and complete view of biological data. Consistency regards the biological coherence of annotations, while completeness regards the extent and coverage of functional characterization for genomes. We have developed tools that allow scientists to assess and improve the consistency and completeness of microbial genome annotations in the context of the Integrated Microbial Genomes (IMG) family of systems. All publicly available microbial genomes are characterized in IMG using different functional annotation and pathway resources, thus providing a comprehensive framework for identifying and resolving annotation discrepancies. A rule based system for predicting phenotypes in IMG provides a powerful mechanism for validating functional annotations, whereby the phenotypic traits of an organism are inferred based on the presence of certain metabolic reactions and pathways and compared to experimentally observed phenotypes. The IMG family of systems are available at http://img.jgi.doe.gov/.

Simultaneous quantification of multiple bacteria by the BactoChip microarray designed to target species-specific marker genes

Bacteria are ubiquitous throughout the environment, the most abundant inhabitants of the healthy human microbiome, and causal pathogens in a variety of diseases. Their identification in disease is often an essential step in rapid diagnosis and targeted intervention, particularly in clinical settings. At present, clinical bacterial detection and discrimination is primarily culture-based, requiring both time and microbiological expertise, especially for bacteria that are not easily cultivated. Higher-throughput molecular methods based on PCR amplification or, recently, microarrays are reaching the clinic as well. However, these methods are currently restricted to a small set of microbes or based on conserved phylogenetic markers such as the 16S rRNA gene, which are difficult to resolve at the species or strain levels. Here, we designed and experimentally validated the BactoChip, an oligonucleotide microarray for bacterial detection and quantification. The chip allows the culture-independent identification of bacterial species, also determining their relative abundances in complex communities as occur in the commensal microbiota or in clinical settings. The microarray successfully distinguished among bacterial species from 21 different genera using 60-mer probes targeting a novel set of in silico identified high-resolution marker genes. The BactoChip additionally proved accurate in determining species-level relative abundances over a 100-fold dynamic range in complex bacterial communities and with a low limit of detection (0.1%). In combination with the continually increasing number of sequenced bacterial genomes, future iterations of the technology could enable to highly accurate clinically-oriented tools for rapid assessment of bacterial community composition and relative abundances.

Use of degradation tags to control protein levels in the Cyanobacterium Synechocystis sp. Strain PCC 6803

We generated a collection of ssrA-based C-terminal protein degradation tags with different degradation strengths. The steady-state fluorescence levels of different enhanced yellow fluorescent protein (eYFP) tag variants in a Synechocystis sp. indicated a tunable range from 1% to 50% of untagged eYFP.

Comparative metagenomic and rRNA microbial diversity characterization using archaeal and bacterial synthetic communities

Next-generation sequencing has dramatically changed the landscape of microbial ecology, large-scale and in-depth diversity studies being now widely accessible. However, determining the accuracy of taxonomic and quantitative inferences and comparing results obtained with different approaches are complicated by incongruence of experimental and computational data types and also by lack of knowledge of the true ecological diversity. Here we used highly diverse bacterial and archaeal synthetic communities assembled from pure genomic DNAs to compare inferences from metagenomic and SSU rRNA amplicon sequencing. Both Illumina and 454 metagenomic data outperformed amplicon sequencing in quantifying the community composition, but the outcome was dependent on analysis parameters and platform. New approaches in processing and classifying amplicons can reconstruct the taxonomic composition of the community with high reproducibility within primer sets, but all tested primers sets lead to significant taxon-specific biases. Controlled synthetic communities assembled to broadly mimic the phylogenetic richness in target environments can provide important validation for fine-tuning experimental and computational parameters used to characterize natural communities.

Complete Genome Sequence of Bacillus thuringiensis Strain 407 Cry-

Bacillus thuringiensis is an insect pathogen that has been used widely as a biopesticide. Here, we report the genome sequence of strain 407 Cry-, which is used to study the genetic determinants of pathogenicity. The genome consists of a 5.5-Mb chromosome and nine plasmids, including a novel 502-kb megaplasmid.

Draft Genome Sequence of Clostridium sp. Maddingley, Isolated from Coal-Seam Gas Formation Water

Clostridium sp. Maddingley was isolated as an axenic culture from a brown coal-seam formation water sample collected from Victoria, Australia. It lacks the solventogenesis genes found in closely related clostridial strains. Metabolic reconstructions suggest that volatile fatty acids are the main fermentation end products.

Draft Genome Sequence of Methanobacterium sp. Maddingley, Reconstructed from Metagenomic Sequencing of a Methanogenic Microbial Consortium Enriched from Coal-Seam Gas Formation Water

The draft genome of Methanobacterium sp. Maddingley was reconstructed from metagenomic sequencing of a methanogenic microbial consortium enriched from coal-seam gas formation water. It is a hydrogenotrophic methanogen predicted to grow using hydrogen and carbon dioxide.

Twenty-one genome sequences from Pseudomonas species and 19 genome sequences from diverse bacteria isolated from the rhizosphere and endosphere of Populus deltoides

To aid in the investigation of the Populus deltoides microbiome, we generated draft genome sequences for 21 Pseudomonas strains and 19 other diverse bacteria isolated from Populus deltoides roots. Genome sequences for isolates similar to Acidovorax, Bradyrhizobium, Brevibacillus, Caulobacter, Chryseobacterium, Flavobacterium, Herbaspirillum, Novosphingobium, Pantoea, Phyllobacterium, Polaromonas, Rhizobium, Sphingobium, and Variovorax were generated.