Community-wide analysis of microbial genome sequence signatures

Single-Cell-Genomics-Facilitated Read Binning of Candidate Phylum EM19 Genomes from Geothermal Spring Metagenomes

The vast majority of microbial life remains uncatalogued due to the inability to cultivate these organisms in the laboratory. This "microbial dark matter" represents a substantial portion of the tree of life and of the populations that contribute to chemical cycling in many ecosystems. In this work, we leveraged an existing single-cell genomic data set representing the candidate bacterial phylum "Calescamantes" (EM19) to calibrate machine learning algorithms and define metagenomic bins directly from pyrosequencing reads derived from Great Boiling Spring in the U.S. Great Basin. Compared to other assembly-based methods, taxonomic binning with a read-based machine learning approach yielded final assemblies with the highest predicted genome completeness of any method tested. Read-first binning subsequently was used to extract Calescamantes bins from all metagenomes with abundant Calescamantes populations, including metagenomes from Octopus Spring and Bison Pool in Yellowstone National Park and Gongxiaoshe Spring in Yunnan Province, China. Metabolic reconstruction suggests that Calescamantes are heterotrophic, facultative anaerobes, which can utilize oxidized nitrogen sources as terminal electron acceptors for respiration in the absence of oxygen and use proteins as their primary carbon source. Despite their phylogenetic divergence, the geographically separate Calescamantes populations were highly similar in their predicted metabolic capabilities and core gene content, respiring O2, or oxidized nitrogen species for energy conservation in distant but chemically similar hot springs.

Genomic and transcriptomic evidence for scavenging of diverse organic compounds by widespread deep-sea archaea

Microbial activity is one of the most important processes to mediate the flux of organic carbon from the ocean surface to the seafloor. However, little is known about the microorganisms that underpin this key step of the global carbon cycle in the deep oceans. Here we present genomic and transcriptomic evidence that five ubiquitous archaeal groups actively use proteins, carbohydrates, fatty acids and lipids as sources of carbon and energy at depths ranging from 800 to 4,950 m in hydrothermal vent plumes and pelagic background seawater across three different ocean basins. Genome-enabled metabolic reconstructions and gene expression patterns show that these marine archaea are motile heterotrophs with extensive mechanisms for scavenging organic matter. Our results shed light on the ecological and physiological properties of ubiquitous marine archaea and highlight their versatile metabolic strategies in deep oceans that might play a critical role in global carbon cycling.

Metagenome and Metatranscriptome Revealed a Highly Active and Intensive Sulfur Cycle in an Oil-Immersed Hydrothermal Chimney in Guaymas Basin

The hydrothermal vent system is a typical chemosynthetic ecosystem in which microorganisms play essential roles in the geobiochemical cycling. Although it has been well-recognized that the inorganic sulfur compounds are abundant and actively converted through chemosynthetic pathways, the sulfur budget in a hydrothermal vent is poorly characterized due to the complexity of microbial sulfur cycling resulting from the numerous parties involved in the processes. In this study, we performed an integrated metagenomic and metatranscriptomic analysis on a chimney sample from Guaymas Basin to achieve a comprehensive study of each sulfur metabolic pathway and its hosting microorganisms and constructed the microbial sulfur cycle that occurs in the site. Our results clearly illustrated the stratified sulfur oxidation and sulfate reduction at the chimney wall. Besides, sulfur metabolizing is closely interacting with carbon cycles, especially the hydrocarbon degradation process in Guaymas Basin. This work supports that the internal sulfur cycling is intensive and the net sulfur budget is low in the hydrothermal ecosystem.

Experimental Design and Bioinformatics Analysis for the Application of Metagenomics in Environmental Sciences and Biotechnology

Recent advances in DNA sequencing technologies have prompted the widespread application of metagenomics for the investigation of novel bioresources (e.g., industrial enzymes and bioactive molecules) and unknown biohazards (e.g., pathogens and antibiotic resistance genes) in natural and engineered microbial systems across multiple disciplines. This review discusses the rigorous experimental design and sample preparation in the context of applying metagenomics in environmental sciences and biotechnology. Moreover, this review summarizes the principles, methodologies, and state-of-the-art bioinformatics procedures, tools and database resources for metagenomics applications and discusses two popular strategies (analysis of unassembled reads versus assembled contigs/draft genomes) for quantitative or qualitative insights of microbial community structure and functions. Overall, this review aims to facilitate more extensive application of metagenomics in the investigation of uncultured microorganisms, novel enzymes, microbe-environment interactions, and biohazards in biotechnological applications where microbial communities are engineered for bioenergy production, wastewater treatment, and bioremediation.

MaxBin 2.0: an automated binning algorithm to recover genomes from multiple metagenomic datasets

The recovery of genomes from metagenomic datasets is a critical step to defining the functional roles of the underlying uncultivated populations. We previously developed MaxBin, an automated binning approach for high-throughput recovery of microbial genomes from metagenomes. Here we present an expanded binning algorithm, MaxBin 2.0, which recovers genomes from co-assembly of a collection of metagenomic datasets. Tests on simulated datasets revealed that MaxBin 2.0 is highly accurate in recovering individual genomes, and the application of MaxBin 2.0 to several metagenomes from environmental samples demonstrated that it could achieve two complementary goals: recovering more bacterial genomes compared to binning a single sample as well as comparing the microbial community composition between different sampling environments.

AVAILABILITY AND IMPLEMENTATION

MaxBin 2.0 is freely available at http://sourceforge.net/projects/maxbin/ under BSD license.

Identification and Resolution of Microdiversity through Metagenomic Sequencing of Parallel Consortia

To gain a predictive understanding of the interspecies interactions within microbial communities that govern community function, the genomic complement of every member population must be determined. Although metagenomic sequencing has enabled the de novo reconstruction of some microbial genomes from environmental communities, microdiversity confounds current genome reconstruction techniques. To overcome this issue, we performed short-read metagenomic sequencing on parallel consortia, defined as consortia cultivated under the same conditions from the same natural community with overlapping species composition. The differences in species abundance between the two consortia allowed reconstruction of near-complete (at an estimated >85% of gene complement) genome sequences for 17 of the 20 detected member species. Two Halomonas spp. indistinguishable by amplicon analysis were found to be present within the community. In addition, comparison of metagenomic reads against the consensus scaffolds revealed within-species variation for one of the Halomonas populations, one of the Rhodobacteraceae populations, and the Rhizobiales population. Genomic comparison of these representative instances of inter- and intraspecies microdiversity suggests differences in functional potential that may result in the expression of distinct roles in the community. In addition, isolation and complete genome sequence determination of six member species allowed an investigation into the sensitivity and specificity of genome reconstruction processes, demonstrating robustness across a wide range of sequence coverage (9× to 2,700×) within the metagenomic data set.

Development of self-compressing BLSOM for comprehensive analysis of big sequence data

With the remarkable increase in genomic sequence data from various organisms, novel tools are needed for comprehensive analyses of available big sequence data. We previously developed a Batch-Learning Self-Organizing Map (BLSOM), which can cluster genomic fragment sequences according to phylotype solely dependent on oligonucleotide composition and applied to genome and metagenomic studies. BLSOM is suitable for high-performance parallel-computing and can analyze big data simultaneously, but a large-scale BLSOM needs a large computational resource. We have developed Self-Compressing BLSOM (SC-BLSOM) for reduction of computation time, which allows us to carry out comprehensive analysis of big sequence data without the use of high-performance supercomputers. The strategy of SC-BLSOM is to hierarchically construct BLSOMs according to data class, such as phylotype. The first-layer BLSOM was constructed with each of the divided input data pieces that represents the data subclass, such as phylotype division, resulting in compression of the number of data pieces. The second BLSOM was constructed with a total of weight vectors obtained in the first-layer BLSOMs. We compared SC-BLSOM with the conventional BLSOM by analyzing bacterial genome sequences. SC-BLSOM could be constructed faster than BLSOM and cluster the sequences according to phylotype with high accuracy, showing the method's suitability for efficient knowledge discovery from big sequence data.

Identification and analysis of the bacterial endosymbiont specialized for production of the chemotherapeutic natural product ET-743

Ecteinascidin 743 (ET-743, Yondelis) is a clinically approved chemotherapeutic natural product isolated from the Caribbean mangrove tunicate Ecteinascidia turbinata. Researchers have long suspected that a microorganism may be the true producer of the anticancer drug, but its genome has remained elusive due to our inability to culture the bacterium in the laboratory using standard techniques. Here, we sequenced and assembled the complete genome of the ET-743 producer, Candidatus Endoecteinascidia frumentensis, directly from metagenomic DNA isolated from the tunicate. Analysis of the ∼ 631 kb microbial genome revealed strong evidence of an endosymbiotic lifestyle and extreme genome reduction. Phylogenetic analysis suggested that the producer of the anti-cancer drug is taxonomically distinct from other sequenced microorganisms and could represent a new family of Gammaproteobacteria. The complete genome has also greatly expanded our understanding of ET-743 production and revealed new biosynthetic genes dispersed across more than 173 kb of the small genome. The gene cluster's architecture and its preservation demonstrate that the drug is likely essential to the interactions of the microorganism with its mangrove tunicate host. Taken together, these studies elucidate the lifestyle of a unique, and pharmaceutically important microorganism and highlight the wide diversity of bacteria capable of making potent natural products.

Next-generation sequencing (NGS) for assessment of microbial water quality: current progress, challenges, and future opportunities

Water quality is an emergent property of a complex system comprised of interacting microbial populations and introduced microbial and chemical contaminants. Studies leveraging next-generation sequencing (NGS) technologies are providing new insights into the ecology of microbially mediated processes that influence fresh water quality such as algal blooms, contaminant biodegradation, and pathogen dissemination. In addition, sequencing methods targeting small subunit (SSU) rRNA hypervariable regions have allowed identification of signature microbial species that serve as bioindicators for sewage contamination in these environments. Beyond amplicon sequencing, metagenomic and metatranscriptomic analyses of microbial communities in fresh water environments reveal the genetic capabilities and interplay of waterborne microorganisms, shedding light on the mechanisms for production and biodegradation of toxins and other contaminants. This review discusses the challenges and benefits of applying NGS-based methods to water quality research and assessment. We will consider the suitability and biases inherent in the application of NGS as a screening tool for assessment of biological risks and discuss the potential and limitations for direct quantitative interpretation of NGS data. Secondly, we will examine case studies from recent literature where NGS based methods have been applied to topics in water quality assessment, including development of bioindicators for sewage pollution and microbial source tracking, characterizing the distribution of toxin and antibiotic resistance genes in water samples, and investigating mechanisms of biodegradation of harmful pollutants that threaten water quality. Finally, we provide a short review of emerging NGS platforms and their potential applications to the next generation of water quality assessment tools.

Refining the phylum Chlorobi by resolving the phylogeny and metabolic potential of the representative of a deeply branching, uncultivated lineage

Recent studies have expanded the phylum Chlorobi, demonstrating that the green sulfur bacteria (GSB), the original cultured representatives of the phylum, are a part of a broader lineage whose members have more diverse metabolic capabilities that overlap with members of the phylum Bacteroidetes. The 16S rRNA gene of an uncultivated clone, OPB56, distantly related to the phyla Chlorobi and Bacteroidetes, was recovered from Obsidian Pool in Yellowstone National Park; however, the detailed phylogeny and function of OPB56 and related clones have remained unknown. Culturing of thermophilic bacterial consortia from compost by adaptation to grow on ionic-liquid pretreated switchgrass provided a consortium in which one of the most abundant members, NICIL-2, clustered with OPB56-related clones. Phylogenetic analysis using the full-length 16S rRNA gene from NICIL-2 demonstrated that it was part of a monophyletic clade, referred to as OPB56, distinct from the Bacteroidetes and Chlorobi. A near complete draft genome (>95% complete) was recovered from metagenomic data from the culture adapted to grow on ionic-liquid pretreated switchgrass using an automated binning algorithm, and this genome was used for marker gene-based phylogenetic analysis and metabolic reconstruction. Six additional genomes related to NICIL-2 were reconstructed from metagenomic data sets obtained from thermal springs at Yellowstone National Park and Nevada Great Boiling Spring. In contrast to the 16S rRNA gene phylogenetic analysis, protein phylogenetic analysis was most consistent with the clustering of the Chlorobea, Ignavibacteria and OPB56 into a single phylum level clade. Metabolic reconstruction of NICIL-2 demonstrated a close linkage with the class Ignavibacteria and the family Rhodothermaceae, a deeply branching Bacteroidetes lineage. The combined phylogenetic and functional analysis of the NICIL-2 genome has refined the membership in the phylum Chlorobi and emphasized the close evolutionary and metabolic relationship between the phyla Chlorobi and the Bacteroidetes.

From cultured to uncultured genome sequences: metagenomics and modeling microbial ecosystems

Microorganisms and the viruses that infect them are the most numerous biological entities on Earth and enclose its greatest biodiversity and genetic reservoir. With strength in their numbers, these microscopic organisms are major players in the cycles of energy and matter that sustain all life. Scientists have only scratched the surface of this vast microbial world through culture-dependent methods. Recent developments in generating metagenomes, large random samples of nucleic acid sequences isolated directly from the environment, are providing comprehensive portraits of the composition, structure, and functioning of microbial communities. Moreover, advances in metagenomic analysis have created the possibility of obtaining complete or nearly complete genome sequences from uncultured microorganisms, providing important means to study their biology, ecology, and evolution. Here we review some of the recent developments in the field of metagenomics, focusing on the discovery of genetic novelty and on methods for obtaining uncultured genome sequences, including through the recycling of previously published datasets. Moreover we discuss how metagenomics has become a core scientific tool to characterize eco-evolutionary patterns of microbial ecosystems, thus allowing us to simultaneously discover new microbes and study their natural communities. We conclude by discussing general guidelines and challenges for modeling the interactions between uncultured microorganisms and viruses based on the information contained in their genome sequences. These models will significantly advance our understanding of the functioning of microbial ecosystems and the roles of microbes in the environment.

The reduced genomes of Parcubacteria (OD1) contain signatures of a symbiotic lifestyle

Candidate phylum OD1 bacteria (also referred to as Parcubacteria) have been identified in a broad range of anoxic environments through community survey analysis. Although none of these species have been isolated in the laboratory, several genome sequences have been reconstructed from metagenomic sequence data and single-cell sequencing. The organisms have small (generally <1 Mb) genomes with severely reduced metabolic capabilities. We have reconstructed 8 partial to near-complete OD1 genomes from oxic groundwater samples, and compared them against existing genomic data. The conserved core gene set comprises 202 genes, or ~28% of the genomic complement. “Housekeeping” genes and genes for biosynthesis of peptidoglycan and Type IV pilus production are conserved. Gene sets for biosynthesis of cofactors, amino acids, nucleotides, and fatty acids are absent entirely or greatly reduced. The only aspects of energy metabolism conserved are the non-oxidative branch of the pentose-phosphate shunt and central glycolysis. These organisms also lack some activities conserved in almost all other known bacterial genomes, including signal recognition particle, pseudouridine synthase A, and FAD synthase. Pan-genome analysis indicates a broad genotypic diversity and perhaps a highly fluid gene complement, indicating historical adaptation to a wide range of growth environments and a high degree of specialization. The genomes were examined for signatures suggesting either a free-living, streamlined lifestyle, or a symbiotic lifestyle. The lack of biosynthetic capabilities and DNA repair, along with the presence of potential attachment and adhesion proteins suggest that the Parcubacteria are ectosymbionts or parasites of other organisms. The wide diversity of genes that potentially mediate cell-cell contact suggests a broad range of partner/prey organisms across the phylum.

CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes

Large-scale recovery of genomes from isolates, single cells, and metagenomic data has been made possible by advances in computational methods and substantial reductions in sequencing costs. Although this increasing breadth of draft genomes is providing key information regarding the evolutionary and functional diversity of microbial life, it has become impractical to finish all available reference genomes. Making robust biological inferences from draft genomes requires accurate estimates of their completeness and contamination. Current methods for assessing genome quality are ad hoc and generally make use of a limited number of "marker" genes conserved across all bacterial or archaeal genomes. Here we introduce CheckM, an automated method for assessing the quality of a genome using a broader set of marker genes specific to the position of a genome within a reference genome tree and information about the collocation of these genes. We demonstrate the effectiveness of CheckM using synthetic data and a wide range of isolate-, single-cell-, and metagenome-derived genomes. CheckM is shown to provide accurate estimates of genome completeness and contamination and to outperform existing approaches. Using CheckM, we identify a diverse range of errors currently impacting publicly available isolate genomes and demonstrate that genomes obtained from single cells and metagenomic data vary substantially in quality. In order to facilitate the use of draft genomes, we propose an objective measure of genome quality that can be used to select genomes suitable for specific gene- and genome-centric analyses of microbial communities.

Phylogeny and physiology of candidate phylum 'Atribacteria' (OP9/JS1) inferred from cultivation-independent genomics

The 'Atribacteria' is a candidate phylum in the Bacteria recently proposed to include members of the OP9 and JS1 lineages. OP9 and JS1 are globally distributed, and in some cases abundant, in anaerobic marine sediments, geothermal environments, anaerobic digesters and reactors and petroleum reservoirs. However, the monophyly of OP9 and JS1 has been questioned and their physiology and ecology remain largely enigmatic due to a lack of cultivated representatives. Here cultivation-independent genomic approaches were used to provide a first comprehensive view of the phylogeny, conserved genomic features and metabolic potential of members of this ubiquitous candidate phylum. Previously available and heretofore unpublished OP9 and JS1 single-cell genomic data sets were used as recruitment platforms for the reconstruction of atribacterial metagenome bins from a terephthalate-degrading reactor biofilm and from the monimolimnion of meromictic Sakinaw Lake. The single-cell genomes and metagenome bins together comprise six species- to genus-level groups that represent most major lineages within OP9 and JS1. Phylogenomic analyses of these combined data sets confirmed the monophyly of the 'Atribacteria' inclusive of OP9 and JS1. Additional conserved features within the 'Atribacteria' were identified, including a gene cluster encoding putative bacterial microcompartments that may be involved in aldehyde and sugar metabolism, energy conservation and carbon storage. Comparative analysis of the metabolic potential inferred from these data sets revealed that members of the 'Atribacteria' are likely to be heterotrophic anaerobes that lack respiratory capacity, with some lineages predicted to specialize in either primary fermentation of carbohydrates or secondary fermentation of organic acids, such as propionate.

Unusual biology across a group comprising more than 15% of domain Bacteria

A prominent feature of the bacterial domain is a radiation of major lineages that are defined as candidate phyla because they lack isolated representatives. Bacteria from these phyla occur in diverse environments and are thought to mediate carbon and hydrogen cycles. Genomic analyses of a few representatives suggested that metabolic limitations have prevented their cultivation. Here we reconstructed 8 complete and 789 draft genomes from bacteria representing >35 phyla and documented features that consistently distinguish these organisms from other bacteria. We infer that this group, which may comprise >15% of the bacterial domain, has shared evolutionary history, and describe it as the candidate phyla radiation (CPR). All CPR genomes are small and most lack numerous biosynthetic pathways. Owing to divergent 16S ribosomal RNA (rRNA) gene sequences, 50-100% of organisms sampled from specific phyla would evade detection in typical cultivation-independent surveys. CPR organisms often have self-splicing introns and proteins encoded within their rRNA genes, a feature rarely reported in bacteria. Furthermore, they have unusual ribosome compositions. All are missing a ribosomal protein often absent in symbionts, and specific lineages are missing ribosomal proteins and biogenesis factors considered universal in bacteria. This implies different ribosome structures and biogenesis mechanisms, and underlines unusual biology across a large part of the bacterial domain.

Metagenomic resolution of microbial functions in deep-sea hydrothermal plumes across the Eastern Lau Spreading Center

Microbial processes within deep-sea hydrothermal plumes affect ocean biogeochemistry on global scales. In rising hydrothermal plumes, a combination of microbial metabolism and particle formation processes initiate the transformation of reduced chemicals like hydrogen sulfide, hydrogen, methane, iron, manganese and ammonia that are abundant in hydrothermal vent fluids. Despite the biogeochemical importance of this rising portion of plumes, it is understudied in comparison to neutrally buoyant plumes. Here we use metagenomics and bioenergetic modeling to describe the abundance and genetic potential of microorganisms in relation to available electron donors in five different hydrothermal plumes and three associated background deep-sea waters from the Eastern Lau Spreading Center located in the Western Pacific Ocean. Three hundred and thirty one distinct genomic 'bins' were identified, comprising an estimated 951 genomes of archaea, bacteria, eukarya and viruses. A significant proportion of these genomes is from novel microorganisms and thus reveals insights into the energy metabolism of heretofore unknown microbial groups. Community-wide analyses of genes encoding enzymes that oxidize inorganic energy sources showed that sulfur oxidation was the most abundant and diverse chemolithotrophic microbial metabolism in the community. Genes for sulfur oxidation were commonly present in genomic bins that also contained genes for oxidation of hydrogen and methane, suggesting metabolic versatility in these microbial groups. The relative diversity and abundance of genes encoding hydrogen oxidation was moderate, whereas that of genes for methane and ammonia oxidation was low in comparison to sulfur oxidation. Bioenergetic-thermodynamic modeling supports the metagenomic analyses, showing that oxidation of elemental sulfur with oxygen is the most dominant catabolic reaction in the hydrothermal plumes. We conclude that the energy metabolism of microbial communities inhabiting rising hydrothermal plumes is dictated by the underlying plume chemistry, with a dominant role for sulfur-based chemolithoautotrophy.

Lifestyle evolution in cyanobacterial symbionts of sponges

The "Candidatus Synechococcus spongiarum" group includes different clades of cyanobacteria with high 16S rRNA sequence identity (~99%) and is the most abundant and widespread cyanobacterial symbiont of marine sponges. The first draft genome of a "Ca. Synechococcus spongiarum" group member was recently published, providing evidence of genome reduction by loss of genes involved in several nonessential functions. However, "Ca. Synechococcus spongiarum" includes a variety of clades that may differ widely in genomic repertoire and consequently in physiology and symbiotic function. Here, we present three additional draft genomes of "Ca. Synechococcus spongiarum," each from a different clade. By comparing all four symbiont genomes to those of free-living cyanobacteria, we revealed general adaptations to life inside sponges and specific adaptations of each phylotype. Symbiont genomes shared about half of their total number of coding genes. Common traits of "Ca. Synechococcus spongiarum" members were a high abundance of DNA modification and recombination genes and a reduction in genes involved in inorganic ion transport and metabolism, cell wall biogenesis, and signal transduction mechanisms. Moreover, these symbionts were characterized by a reduced number of antioxidant enzymes and low-weight peptides of photosystem II compared to their free-living relatives. Variability within the "Ca. Synechococcus spongiarum" group was mostly related to immune system features, potential for siderophore-mediated iron transport, and dependency on methionine from external sources. The common absence of genes involved in synthesis of residues, typical of the O antigen of free-living Synechococcus species, suggests a novel mechanism utilized by these symbionts to avoid sponge predation and phage attack.

IMPORTANCE

While the Synechococcus/Prochlorococcus-type cyanobacteria are widely distributed in the world's oceans, a subgroup has established its niche within marine sponge tissues. Recently, the first genome of sponge-associated cyanobacteria, "Candidatus Synechococcus spongiarum," was described. The sequencing of three representatives of different clades within this cyanobacterial group has enabled us to investigate intraspecies diversity, as well as to give a more comprehensive understanding of the common symbiotic features that adapt "Ca. Synechococcus spongiarum" to its life within the sponge host.

Ecological roles of dominant and rare prokaryotes in acid mine drainage revealed by metagenomics and metatranscriptomics

High-throughput sequencing is expanding our knowledge of microbial diversity in the environment. Still, understanding the metabolic potentials and ecological roles of rare and uncultured microbes in natural communities remains a major challenge. To this end, we applied a 'divide and conquer' strategy that partitioned a massive metagenomic data set (>100 Gbp) into subsets based on K-mer frequency in sequence assembly to a low-diversity acid mine drainage (AMD) microbial community and, by integrating with an additional metatranscriptomic assembly, successfully obtained 11 draft genomes most of which represent yet uncultured and/or rare taxa (relative abundance <1%). We report the first genome of a naturally occurring Ferrovum population (relative abundance >90%) and its metabolic potentials and gene expression profile, providing initial molecular insights into the ecological role of these lesser known, but potentially important, microorganisms in the AMD environment. Gene transcriptional analysis of the active taxa revealed major metabolic capabilities executed in situ, including carbon- and nitrogen-related metabolisms associated with syntrophic interactions, iron and sulfur oxidation, which are key in energy conservation and AMD generation, and the mechanisms of adaptation and response to the environmental stresses (heavy metals, low pH and oxidative stress). Remarkably, nitrogen fixation and sulfur oxidation were performed by the rare taxa, indicating their critical roles in the overall functioning and assembly of the AMD community. Our study demonstrates the potential of the 'divide and conquer' strategy in high-throughput sequencing data assembly for genome reconstruction and functional partitioning analysis of both dominant and rare species in natural microbial assemblages.

CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes

Large-scale recovery of genomes from isolates, single cells, and metagenomic data has been made possible by advances in computational methods and substantial reductions in sequencing costs. Although this increasing breadth of draft genomes is providing key information regarding the evolutionary and functional diversity of microbial life, it has become impractical to finish all available reference genomes. Making robust biological inferences from draft genomes requires accurate estimates of their completeness and contamination. Current methods for assessing genome quality are ad hoc and generally make use of a limited number of “marker” genes conserved across all bacterial or archaeal genomes. Here we introduce CheckM, an automated method for assessing the quality of a genome using a broader set of marker genes specific to the position of a genome within a reference genome tree and information about the collocation of these genes. We demonstrate the effectiveness of CheckM using synthetic data and a wide range of isolate-, single-cell-, and metagenome-derived genomes. CheckM is shown to provide accurate estimates of genome completeness and contamination and to outperform existing approaches. Using CheckM, we identify a diverse range of errors currently impacting publicly available isolate genomes and demonstrate that genomes obtained from single cells and metagenomic data vary substantially in quality. In order to facilitate the use of draft genomes, we propose an objective measure of genome quality that can be used to select genomes suitable for specific gene- and genome-centric analyses of microbial communities.

Metagenomics: tools and insights for analyzing next-generation sequencing data derived from biodiversity studies

Advances in next-generation sequencing (NGS) have allowed significant breakthroughs in microbial ecology studies. This has led to the rapid expansion of research in the field and the establishment of "metagenomics", often defined as the analysis of DNA from microbial communities in environmental samples without prior need for culturing. Many metagenomics statistical/computational tools and databases have been developed in order to allow the exploitation of the huge influx of data. In this review article, we provide an overview of the sequencing technologies and how they are uniquely suited to various types of metagenomic studies. We focus on the currently available bioinformatics techniques, tools, and methodologies for performing each individual step of a typical metagenomic dataset analysis. We also provide future trends in the field with respect to tools and technologies currently under development. Moreover, we discuss data management, distribution, and integration tools that are capable of performing comparative metagenomic analyses of multiple datasets using well-established databases, as well as commonly used annotation standards.

Draft Genome Sequence of Ruminoclostridium sp. Ne3, Clostridia from an Enrichment Culture Obtained from Australian Subterranean Termite, Nasutitermes exitiosus

The draft genome sequence of Ruminoclostridium sp. Ne3 was reconstructed from the metagenome of a hydrogenogenic microbial consortium growing on xylan. The organism is likely the primary hemicellulose degrader within the consortium.

Draft Genome Sequence of Clostridium beijerinckii Ne1, Clostridia from an Enrichment Culture Obtained from Australian Subterranean Termite, Nasutitermes exitiosus

The draft genome of Clostridium beijerinckii strain Ne1 was reconstructed from the metagenomic sequence of a mixed-microbial consortium that produced commercially significant quantities of hydrogen from xylan as a sole feedstock. The organism possesses relatively limited hemicellulolytic capacity and likely requires the action of other organisms to completely degrade xylan.

Draft Genome Sequence of Clostridium sp. Ne2, Clostridia from an Enrichment Culture Obtained from Australian Subterranean Termite, Nasutitermes exitiosus

The draft genome sequence of Clostridium sp. Ne2 was reconstructed from a metagenome of a hydrogenogenic microbial consortium. The organism is most closely related to Clostridium magnum and is a strict anaerobe that is predicted to ferment a range of simple sugars.

The changing landscape of microbial biodiversity exploration and its implications for systematics

A vast diversity of Bacteria and Archaea exists in nature that has evaded axenic culture. Advancements in single-cell genomics, metagenomics, and molecular microbial ecology approaches provide ever-improving insight into the biology of this so-called "microbial dark matter"; however, due to the International Code of Nomenclature of Prokaryotes, yet-uncultivated microorganisms are not accommodated in formal taxonomy regardless of the quantity or quality of data. Meanwhile, efforts to calibrate the existing taxonomy with phylogenetic anchors and genomic data are increasingly robust. The current climate provides an exciting opportunity to leverage rapidly expanding single-cell genomics and metagenomics datasets to improve the taxonomy of Bacteria and Archaea. However, this opportunity must be weighted carefully in light of the strengths and limitations of these approaches. We propose to expand the definition of the Candidatus taxonomy to include taxa, from the phylum level to the species level, that are described genomically, particularly when genomic work is coupled with advanced molecular ecology approaches to probe metabolic functions in situ. This system would preserve the rigor and value of traditional microbial systematics while enabling growth of a provisional taxonomic structure to facilitate communication about "dark" lineages on the tree of life.

Accurate, multi-kb reads resolve complex populations and detect rare microorganisms

Accurate evaluation of microbial communities is essential for understanding global biogeochemical processes and can guide bioremediation and medical treatments. Metagenomics is most commonly used to analyze microbial diversity and metabolic potential, but assemblies of the short reads generated by current sequencing platforms may fail to recover heterogeneous strain populations and rare organisms. Here we used short (150-bp) and long (multi-kb) synthetic reads to evaluate strain heterogeneity and study microorganisms at low abundance in complex microbial communities from terrestrial sediments. The long-read data revealed multiple (probably dozens of) closely related species and strains from previously undescribed Deltaproteobacteria and Aminicenantes (candidate phylum OP8). Notably, these are the most abundant organisms in the communities, yet short-read assemblies achieved only partial genome coverage, mostly in the form of short scaffolds (N50 = ∼ 2200 bp). Genome architecture and metabolic potential for these lineages were reconstructed using a new synteny-based method. Analysis of long-read data also revealed thousands of species whose abundances were <0.1% in all samples. Most of the organisms in this "long tail" of rare organisms belong to phyla that are also represented by abundant organisms. Genes encoding glycosyl hydrolases are significantly more abundant than expected in rare genomes, suggesting that rare species may augment the capability for carbon turnover and confer resilience to changing environmental conditions. Overall, the study showed that a diversity of closely related strains and rare organisms account for a major portion of the communities. These are probably common features of many microbial communities and can be effectively studied using a combination of long and short reads.

Targeted diversity generation by intraterrestrial archaea and archaeal viruses

In the evolutionary arms race between microbes, their parasites, and their neighbours, the capacity for rapid protein diversification is a potent weapon. Diversity-generating retroelements (DGRs) use mutagenic reverse transcription and retrohoming to generate myriad variants of a target gene. Originally discovered in pathogens, these retroelements have been identified in bacteria and their viruses, but never in archaea. Here we report the discovery of intact DGRs in two distinct intraterrestrial archaeal systems: a novel virus that appears to infect archaea in the marine subsurface, and, separately, two uncultivated nanoarchaea from the terrestrial subsurface. The viral DGR system targets putative tail fibre ligand-binding domains, potentially generating >10¹⁸ protein variants. The two single-cell nanoarchaeal genomes each possess ≥4 distinct DGRs. Against an expected background of low genome-wide mutation rates, these results demonstrate a previously unsuspected potential for rapid, targeted sequence diversification in intraterrestrial archaea and their viruses.

Diversity-generating retroelements (DGRs) are genetic elements that introduce sequence variation within target genes in bacteria and their viruses. Here, Paul et al. report the discovery of DGRs in an archaeal virus and in two archaea from marine and terrestrial subsurface environments, respectively.

Counts and sequences, observations that continue to change our understanding of viruses in nature

The discovery of abundant viruses in the oceans and on land has ushered in a quarter century of groundbreaking advancements in our understanding of viruses within ecosystems. Two types of observations from environmental samples--direct counts of viral particles and viral metagenomic sequences--have been critical to these discoveries. Accurate direct counts have established ecosystem-scale trends in the impacts of viral infection on microbial host populations and have shown that viral communities within aquatic and soil environments respond to both short term and seasonal environmental change. Direct counts have been critical for estimating viral production rate, a measurement essential to quantifying the implications of viral infection for the biogeochemical cycling of nutrients within ecosystems. While direct counts have defined the magnitude of viral processes; shotgun sequences of environmental viral DNA--virome sequences--have enabled researchers to estimate the diversity and composition of natural viral communities. Virome-enabled studies have found the virioplankton to contain thousands of viral genotypes in communities where the most dominant viral population accounts for a small fraction of total abundance followed by a long tail of diverse populations. Detailed examination of long virome sequences has led to new understanding of genotype-to-phenotype connections within marine viruses and revealed that viruses carry metabolic genes that are important to maintaining cellular energy during viral replication. Increased access to long virome sequences will undoubtedly reveal more genetic secrets of viruses and enable us to build a genomics rulebook for predicting key biological and ecological features of unknown viruses.

Gut bacteria are rarely shared by co-hospitalized premature infants, regardless of necrotizing enterocolitis development

Premature infants are highly vulnerable to aberrant gastrointestinal tract colonization, a process that may lead to diseases like necrotizing enterocolitis. Thus, spread of potential pathogens among hospitalized infants is of great concern. Here, we reconstructed hundreds of high-quality genomes of microorganisms that colonized co-hospitalized premature infants, assessed their metabolic potential, and tracked them over time to evaluate bacterial strain dispersal among infants. We compared microbial communities in infants who did and did not develop necrotizing enterocolitis. Surprisingly, while potentially pathogenic bacteria of the same species colonized many infants, our genome-resolved analysis revealed that strains colonizing each baby were typically distinct. In particular, no strain was common to all infants who developed necrotizing enterocolitis. The paucity of shared gut colonizers suggests the existence of significant barriers to the spread of bacteria among infants. Importantly, we demonstrate that strain-resolved comprehensive community analysis can be accomplished on potentially medically relevant time scales.

DOI:http://dx.doi.org/10.7554/eLife.05477.001

The spread of potentially harmful bacteria is a major source of disease in patients staying in hospitals. Premature babies—born before 37 weeks of pregnancy—can be particularly vulnerable to these infections because their organs may not yet be fully developed. Also, young babies do not have the fully established populations of beneficial microbes that help to protect us from dangerous bacteria.

Necrotizing enterocolitis—a life-threatening disease that can cause portions of the bowel to die—is mostly seen in extremely premature babies. Although it is not known what causes this serious condition, research has suggested that a contagious microbe may be responsible.

The development of methods that can sequence DNA from whole communities of microbes, known as metagenomics, allows researchers to identify the presence of individual strains of bacteria within these communities. This makes it possible to compare and contrast the strains of bacteria present in both diseased and healthy individuals, to help identify the bacteria responsible for a disease.

Here, Raveh-Sadka et al. used a metagenomics approach to study the communities of microbes present in premature babies in a hospital unit during an outbreak of necrotizing enterocolitis. The study found that very few bacterial strains were present in more than one baby, suggesting that bacterial strains are not readily transferred between the babies while they are in the hospital. Furthermore, Raveh-Sadka et al. reveal that no single bacterial strain was shared among all the babies who developed necrotizing enterocolitis.

These findings indicate that necrotizing enterocolitis is not caused by a single strain of bacterium. Instead, if bacteria do contribute to the disease, it maybe that it is caused by a variety of potentially harmful bacteria colonizing the gut at the cost of beneficial bacteria. In future, better understanding of the barriers that limit the transfer of bacteria between premature babies could help inform efforts to reduce the spread of infections between patients in hospitals.

DOI:http://dx.doi.org/10.7554/eLife.05477.002

Metatranscriptomic analysis of diminutive Thiomargarita-like bacteria ("Candidatus Thiopilula" spp.) from abyssal cold seeps of the Barbados Accretionary Prism

Large sulfur-oxidizing bacteria in the family Beggiatoaceae are important players in the global sulfur cycle. This group contains members of the well-known genera Beggiatoa, Thioploca, and Thiomargarita but also recently identified and relatively unknown candidate taxa, including "Candidatus Thiopilula" spp. and "Ca. Thiophysa" spp. We discovered a population of "Ca. Thiopilula" spp. colonizing cold seeps near Barbados at a ∼4.7-km water depth. The Barbados population consists of spherical cells that are morphologically similar to Thiomargarita spp., with elemental sulfur inclusions and a central vacuole, but have much smaller cell diameters (5 to 40 μm). Metatranscriptomic analysis revealed that when exposed to anoxic sulfidic conditions, Barbados "Ca. Thiopilula" organisms expressed genes for the oxidation of elemental sulfur and the reduction of nitrogenous compounds, consistent with their vacuolated morphology and intracellular sulfur storage capability. Metatranscriptomic analysis further revealed that anaerobic methane-oxidizing and sulfate-reducing organisms were active in the sediment, which likely provided reduced sulfur substrates for "Ca. Thiopilula" and other sulfur-oxidizing microorganisms in the community. The novel observations of "Ca. Thiopilula" and associated organisms reported here expand our knowledge of the globally distributed and ecologically successful Beggiatoaceae group and thus offer insight into the composition and ecology of deep cold seep microbial communities.

Genomic expansion of domain archaea highlights roles for organisms from new phyla in anaerobic carbon cycling

BACKGROUND

Archaea represent a significant fraction of Earth's biodiversity, yet they remain much less well understood than Bacteria. Gene surveys, a few metagenomic studies, and some single-cell sequencing projects have revealed numerous little-studied archaeal phyla. Certain lineages appear to branch deeply and may be part of a major phylum radiation. The structure of this radiation and the physiology of the organisms remain almost unknown.

RESULTS

We used genome-resolved metagenomic analyses to investigate the diversity, genomes sizes, metabolic capacities, and potential roles of Archaea in terrestrial subsurface biogeochemical cycles. We sequenced DNA from complex sediment and planktonic consortia from an aquifer adjacent to the Colorado River (USA) and reconstructed the first complete genomes for Archaea using cultivation-independent methods. To provide taxonomic context, we analyzed an additional 151 newly sampled archaeal sequences. We resolved two new phyla within a major, apparently deep-branching group of phyla (a superphylum). The organisms have small genomes, and metabolic predictions indicate that their primary contributions to Earth's biogeochemical cycles involve carbon and hydrogen metabolism, probably associated with symbiotic and/or fermentation-based lifestyles.

CONCLUSIONS

The results dramatically expand genomic sampling of the domain Archaea and clarify taxonomic designations within a major superphylum. This study, in combination with recently published work on bacterial phyla lacking cultivated representatives, reveals a fascinating phenomenon of major radiations of organisms with small genomes, novel proteome composition, and strong interdependence in both domains.

The effects of variable sample biomass on comparative metagenomics

Longitudinal studies that integrate samples with variable biomass are essential to understand microbial community dynamics across space or time. Shotgun metagenomics is widely used to investigate these communities at the functional level, but little is known about the effects of combining low and high biomass samples on downstream analysis. We investigated the interacting effects of DNA input and library amplification by polymerase chain reaction on comparative metagenomic analysis using dilutions of a single complex template from an Arabidopsis thaliana-associated microbial community. We modified the Illumina Nextera kit to generate high-quality large-insert (680 bp) paired-end libraries using a range of 50 pg to 50 ng of input DNA. Using assembly-based metagenomic analysis, we demonstrate that DNA input level has a significant impact on community structure due to overrepresentation of low-GC genomic regions following library amplification. In our system, these differences were largely superseded by variations between biological replicates, but our results advocate verifying the influence of library amplification on a case-by-case basis. Overall, this study provides recommendations for quality filtering and de-replication prior to analysis, as well as a practical framework to address the issue of low biomass or biomass heterogeneity in longitudinal metagenomic surveys.

Metagenomic analysis reveals adaptations to a cold-adapted lifestyle in a low-temperature acid mine drainage stream

An acid mine drainage (pH 2.5-2.7) stream biofilm situated 250 m below ground in the low-temperature (6-10°C) Kristineberg mine, northern Sweden, contained a microbial community equipped for growth at low temperature and acidic pH. Metagenomic sequencing of the biofilm and planktonic fractions identified the most abundant microorganism to be similar to the psychrotolerant acidophile, Acidithiobacillus ferrivorans. In addition, metagenome contigs were most similar to other Acidithiobacillus species, an Acidobacteria-like species, and a Gallionellaceae-like species. Analyses of the metagenomes indicated functional characteristics previously characterized as related to growth at low temperature including cold-shock proteins, several pathways for the production of compatible solutes and an anti-freeze protein. In addition, genes were predicted to encode functions related to pH homeostasis and metal resistance related to growth in the acidic metal-containing mine water. Metagenome analyses identified microorganisms capable of nitrogen fixation and exhibiting a primarily autotrophic lifestyle driven by the oxidation of the ferrous iron and inorganic sulfur compounds contained in the sulfidic mine waters. The study identified a low diversity of abundant microorganisms adapted to a low-temperature acidic environment as well as identifying some of the strategies the microorganisms employ to grow in this extreme environment.

Metagenomic insights into S(0) precipitation in a terrestrial subsurface lithoautotrophic ecosystem

The Frasassi and Acquasanta Terme cave systems in Italy host isolated lithoautotrophic ecosystems characterized by sulfur-oxidizing biofilms with up to 50% S(0) by mass. The net contributions of microbial taxa in the biofilms to production and consumption of S(0) are poorly understood and have implications for understanding the formation of geological sulfur deposits as well as the ecological niches of sulfur-oxidizing autotrophs. Filamentous Epsilonproteobacteria are among the principal biofilm architects in Frasassi and Acquasanta Terme streams, colonizing high-sulfide, low-oxygen niches relative to other major biofilm-forming populations. Metagenomic sequencing of eight biofilm samples indicated the presence of diverse and abundant Epsilonproteobacteria. Populations of Sulfurovum-like organisms were the most abundant Epsilonproteobacteria regardless of differences in biofilm morphology, temperature, or water chemistry. After assembling and binning the metagenomic data, we retrieved four nearly-complete genomes of Sulfurovum-like organisms as well as a Sulfuricurvum spp. Analyses of the binned and assembled metagenomic data indicate that the Epsilonproteobacteria are autotrophic and therefore provide organic carbon to the isolated subsurface ecosystem. Multiple homologs of sulfide-quinone oxidoreductase (Sqr), together with incomplete or absent Sox pathways, suggest that cave Sulfurovum-like Epsilonproteobacteria oxidize sulfide incompletely to S(0) using either O₂ or nitrate as a terminal electron acceptor, consistent with previous evidence that they are most successful in niches with high dissolved sulfide to oxygen ratios. In contrast, we recovered homologs of the complete complement of Sox proteins affiliated Gammaproteobacteria and with less abundant Sulfuricurvum spp. and Arcobacter spp., suggesting that these populations are capable of the complete oxidation of sulfide to sulfate. These and other genomic data presented here offer new clues into the physiology and genetic potential of the largely uncultivated and ecologically successful cave Sulfurovum-like populations, and suggest that they play an integral role in subsurface S(0) formation.

Genomic resolution of linkages in carbon, nitrogen, and sulfur cycling among widespread estuary sediment bacteria

BACKGROUND

Estuaries are among the most productive habitats on the planet. Bacteria in estuary sediments control the turnover of organic carbon and the cycling of nitrogen and sulfur. These communities are complex and primarily made up of uncultured lineages, thus little is known about how ecological and metabolic processes are partitioned in sediments.

RESULTS

De novo assembly and binning resulted in the reconstruction of 82 bacterial genomes from different redox regimes of estuary sediments. These genomes belong to 23 bacterial groups, including uncultured candidate phyla (for example, KSB1, TA06, and KD3-62) and three newly described phyla (White Oak River (WOR)-1, WOR-2, and WOR-3). The uncultured phyla are generally most abundant in the sulfate-methane transition (SMTZ) and methane-rich zones, and genomic data predict that they mediate essential biogeochemical processes of the estuarine environment, including organic carbon degradation and fermentation. Among the most abundant organisms in the sulfate-rich layer are novel Gammaproteobacteria that have genes for the oxidation of sulfur and the reduction of nitrate and nitrite. Interestingly, the terminal steps of denitrification (NO3 to N2O and then N2O to N2) are present in distinct bacterial populations.

CONCLUSIONS

This dataset extends our knowledge of the metabolic potential of several uncultured phyla. Within the sediments, there is redundancy in the genomic potential in different lineages, often distinct phyla, for essential biogeochemical processes. We were able to chart the flow of carbon and nutrients through the multiple geochemical layers of bacterial processing and reveal potential ecological interactions within the communities.

VizBin - an application for reference-independent visualization and human-augmented binning of metagenomic data

BACKGROUND

Metagenomics is limited in its ability to link distinct microbial populations to genetic potential due to a current lack of representative isolate genome sequences. Reference-independent approaches, which exploit for example inherent genomic signatures for the clustering of metagenomic fragments (binning), offer the prospect to resolve and reconstruct population-level genomic complements without the need for prior knowledge.

RESULTS

We present VizBin, a Java™-based application which offers efficient and intuitive reference-independent visualization of metagenomic datasets from single samples for subsequent human-in-the-loop inspection and binning. The method is based on nonlinear dimension reduction of genomic signatures and exploits the superior pattern recognition capabilities of the human eye-brain system for cluster identification and delineation. We demonstrate the general applicability of VizBin for the analysis of metagenomic sequence data by presenting results from two cellulolytic microbial communities and one human-borne microbial consortium. The superior performance of our application compared to other analogous metagenomic visualization and binning methods is also presented.

CONCLUSIONS

VizBin can be applied de novo for the visualization and subsequent binning of metagenomic datasets from single samples, and it can be used for the post hoc inspection and refinement of automatically generated bins. Due to its computational efficiency, it can be run on common desktop machines and enables the analysis of complex metagenomic datasets in a matter of minutes. The software implementation is available at https://claczny.github.io/VizBin under the BSD License (four-clause) and runs under Microsoft Windows™, Apple Mac OS X™ (10.7 to 10.10), and Linux.

Metagenomic evidence for metabolism of trace atmospheric gases by high-elevation desert Actinobacteria

Previous surveys of very dry Atacama Desert mineral soils have consistently revealed sparse communities of non-photosynthetic microbes. The functional nature of these microorganisms remains debatable given the harshness of the environment and low levels of biomass and diversity. The aim of this study was to gain an understanding of the phylogenetic community structure and metabolic potential of a low-diversity mineral soil metagenome that was collected from a high-elevation Atacama Desert volcano debris field. We pooled DNA extractions from over 15 g of volcanic material, and using whole genome shotgun sequencing, observed only 75-78 total 16S rRNA gene OTUs3%. The phylogenetic structure of this community is significantly under dispersed, with actinobacterial lineages making up 97.9-98.6% of the 16S rRNA genes, suggesting a high degree of environmental selection. Due to this low diversity and uneven community composition, we assembled and analyzed the metabolic pathways of the most abundant genome, a Pseudonocardia sp. (56-72% of total 16S genes). Our assembly and binning efforts yielded almost 4.9 Mb of Pseudonocardia sp. contigs, which accounts for an estimated 99.3% of its non-repetitive genomic content. This genome contains a limited array of carbohydrate catabolic pathways, but encodes for CO2 fixation via the Calvin cycle. The genome also encodes complete pathways for the catabolism of various trace gases (H2, CO and several organic C1 compounds) and the assimilation of ammonia and nitrate. We compared genomic content among related Pseudonocardia spp. and estimated rates of non-synonymous and synonymous nucleic acid substitutions between protein coding homologs. Collectively, these comparative analyses suggest that the community structure and various functional genes have undergone strong selection in the nutrient poor desert mineral soils and high-elevation atmospheric conditions.

Comparison of environmental and isolate Sulfobacillus genomes reveals diverse carbon, sulfur, nitrogen, and hydrogen metabolisms

BACKGROUND

Bacteria of the genus Sulfobacillus are found worldwide as members of microbial communities that accelerate sulfide mineral dissolution in acid mine drainage environments (AMD), acid-rock drainage environments (ARD), as well as in industrial bioleaching operations. Despite their frequent identification in these environments, their role in biogeochemical cycling is poorly understood.

RESULTS

Here we report draft genomes of five species of the Sulfobacillus genus (AMDSBA1-5) reconstructed by cultivation-independent sequencing of biofilms sampled from the Richmond Mine (Iron Mountain, CA). Three of these species (AMDSBA2, AMDSBA3, and AMDSBA4) have no cultured representatives while AMDSBA1 is a strain of S. benefaciens, and AMDSBA5 a strain of S. thermosulfidooxidans. We analyzed the diversity of energy conservation and central carbon metabolisms for these genomes and previously published Sulfobacillus genomes. Pathways of sulfur oxidation vary considerably across the genus, including the number and type of subunits of putative heterodisulfide reductase complexes likely involved in sulfur oxidation. The number and type of nickel-iron hydrogenase proteins varied across the genus, as does the presence of different central carbon pathways. Only the AMDSBA3 genome encodes a dissimilatory nitrate reducatase and only the AMDSBA5 and S. thermosulfidooxidans genomes encode assimilatory nitrate reductases. Within the genus, AMDSBA4 is unusual in that its electron transport chain includes a cytochrome bc type complex, a unique cytochrome c oxidase, and two distinct succinate dehydrogenase complexes.

CONCLUSIONS

Overall, the results significantly expand our understanding of carbon, sulfur, nitrogen, and hydrogen metabolism within the Sulfobacillus genus.

Systems-based approaches to unravel multi-species microbial community functioning

Some of the most transformative discoveries promising to enable the resolution of this century's grand societal challenges will most likely arise from environmental science and particularly environmental microbiology and biotechnology. Understanding how microbes interact in situ, and how microbial communities respond to environmental changes remains an enormous challenge for science. Systems biology offers a powerful experimental strategy to tackle the exciting task of deciphering microbial interactions. In this framework, entire microbial communities are considered as metaorganisms and each level of biological information (DNA, RNA, proteins and metabolites) is investigated along with in situ environmental characteristics. In this way, systems biology can help unravel the interactions between the different parts of an ecosystem ultimately responsible for its emergent properties. Indeed each level of biological information provides a different level of characterisation of the microbial communities. Metagenomics, metatranscriptomics, metaproteomics, metabolomics and SIP-omics can be employed to investigate collectively microbial community structure, potential, function, activity and interactions. Omics approaches are enabled by high-throughput 21st century technologies and this review will discuss how their implementation has revolutionised our understanding of microbial communities.

Evaluation of metatranscriptomic protocols and application to the study of freshwater microbial communities

Metatranscriptomics of environmental samples enables the identification of community activities without a priori knowledge of taxonomic or functional composition. However, several technical challenges associated with the RNA preparation protocols can affect the relative representation of transcripts and data interpretation. Here, seven replicate metatranscriptomes from planktonic freshwater samples (Lake Lanier, USA) were sequenced to evaluate technical and biological reproducibility of different RNA extraction protocols. Organic versus bead-beating extraction showed significant enrichment for low versus high G + C% mRNA populations respectively. The sequencing data were best modelled by a negative binomial distribution to account for the large technical and biological variation observed. Despite the variation, the transcriptional activities of populations that persisted in year-round metagenomes from the same site consistently showed distinct expression patterns, reflecting different ecologic strategies and allowing us to test prevailing models on the contribution of both rare biosphere and abundant members to community activity. For instance, abundant members of the Verrucomicrobia phylum systematically showed low transcriptional activity compared with other abundant taxa. Our results provide a practical guide to the analysis of metatranscriptomes and advance understanding of the activity and ecology of abundant and rare members of temperate freshwater microbial communities.

Microscale sulfur cycling in the phototrophic pink berry consortia of the Sippewissett Salt Marsh

Microbial metabolism is the engine that drives global biogeochemical cycles, yet many key transformations are carried out by microbial consortia over short spatiotemporal scales that elude detection by traditional analytical approaches. We investigate syntrophic sulfur cycling in the 'pink berry' consortia of the Sippewissett Salt Marsh through an integrative study at the microbial scale. The pink berries are macroscopic, photosynthetic microbial aggregates composed primarily of two closely associated species: sulfide-oxidizing purple sulfur bacteria (PB-PSB1) and sulfate-reducing bacteria (PB-SRB1). Using metagenomic sequencing and (34) S-enriched sulfate stable isotope probing coupled with nanoSIMS, we demonstrate interspecies transfer of reduced sulfur metabolites from PB-SRB1 to PB-PSB1. The pink berries catalyse net sulfide oxidation and maintain internal sulfide concentrations of 0-500 μm. Sulfide within the berries, captured on silver wires and analysed using secondary ion mass spectrometer, increased in abundance towards the berry interior, while δ(34) S-sulfide decreased from 6‰ to -31‰ from the exterior to interior of the berry. These values correspond to sulfate-sulfide isotopic fractionations (15-53‰) consistent with either sulfate reduction or a mixture of reductive and oxidative metabolisms. Together this combined metagenomic and high-resolution isotopic analysis demonstrates active sulfur cycling at the microscale within well-structured macroscopic consortia consisting of sulfide-oxidizing anoxygenic phototrophs and sulfate-reducing bacteria.

Deciphering diversity and ecological function from marine metagenomes

Metagenomic sequencing now represents a common, powerful approach for investigating diversity and functional relationships in marine ecosystems. High-throughput datasets generated from random fragments of environmental DNA can provide a less biased view of organismal abundance (versus PCR-based amplicon sequencing) and enable novel exploration of microbial genomes by recovering genome assemblies from uncultured species, identifying ecological functions, and reconstructing metabolic pathways. This review highlights the current state of knowledge in marine metagenomics, focusing on biological insights gained from recent environmental studies and detailing commonly employed methods for data collection and analysis.

Binning metagenomic contigs by coverage and composition

Shotgun sequencing enables the reconstruction of genomes from complex microbial communities, but because assembly does not reconstruct entire genomes, it is necessary to bin genome fragments. Here we present CONCOCT, a new algorithm that combines sequence composition and coverage across multiple samples, to automatically cluster contigs into genomes. We demonstrate high recall and precision on artificial as well as real human gut metagenome data sets.

Coupled reductive and oxidative sulfur cycling in the phototrophic plate of a meromictic lake

Mahoney Lake represents an extreme meromictic model system and is a valuable site for examining the organisms and processes that sustain photic zone euxinia (PZE). A single population of purple sulfur bacteria (PSB) living in a dense phototrophic plate in the chemocline is responsible for most of the primary production in Mahoney Lake. Here, we present metagenomic data from this phototrophic plate--including the genome of the major PSB, as obtained from both a highly enriched culture and from the metagenomic data--as well as evidence for multiple other taxa that contribute to the oxidative sulfur cycle and to sulfate reduction. The planktonic PSB is a member of the Chromatiaceae, here renamed Thiohalocapsa sp. strain ML1. It produces the carotenoid okenone, yet its closest relatives are benthic PSB isolates, a finding that may complicate the use of okenone (okenane) as a biomarker for ancient PZE. Favorable thermodynamics for non-phototrophic sulfide oxidation and sulfate reduction reactions also occur in the plate, and a suite of organisms capable of oxidizing and reducing sulfur is apparent in the metagenome. Fluctuating supplies of both reduced carbon and reduced sulfur to the chemocline may partly account for the diversity of both autotrophic and heterotrophic species. Collectively, the data demonstrate the physiological potential for maintaining complex sulfur and carbon cycles in an anoxic water column, driven by the input of exogenous organic matter. This is consistent with suggestions that high levels of oxygenic primary production maintain episodes of PZE in Earth's history and that such communities should support a diversity of sulfur cycle reactions.

Expanding our view of genomic diversity in Candidatus Accumulibacter clades

Enhanced biological phosphorus removal (EBPR) is an important industrial wastewater treatment process mediated by polyphosphate-accumulating organisms (PAOs). Members of the genus Candidatus Accumulibacter are one of the most extensively studied PAO as they are commonly enriched in lab-scale EBPR reactors. Members of different Accumulibacter clades are often enriched through changes in reactor process conditions; however, the two currently sequenced Accumulibacter genomes show extensive metabolic similarity. Here, we expand our understanding of Accumulibacter genomic diversity through recovery of eight population genomes using deep metagenomics, including seven from phylogenetic clades with no previously sequenced representative. Comparative genomic analysis revealed a core of shared genes involved primarily in carbon and phosphorus metabolism; however, each Accumulibacter genome also encoded a substantial number of unique genes (> 700 genes). A major difference between the Accumulibacter clades was the type of nitrate reductase encoded and the capacity to perform subsequent steps in denitrification. The Accumulibacter clade IIF genomes also contained acetaldehyde dehydrogenase that may allow ethanol to be used as carbon source. These differences in metabolism between Accumulibacter genomes provide a molecular basis for niche differentiation observed in lab-scale reactors and may offer new opportunities for process optimization.

Impact of single-cell genomics and metagenomics on the emerging view of extremophile "microbial dark matter"

Despite >130 years of microbial cultivation studies, many microorganisms remain resistant to traditional cultivation approaches, including numerous candidate phyla of bacteria and archaea. Unraveling the mysteries of these candidate phyla is a grand challenge in microbiology and is especially important in habitats where they are abundant, including some extreme environments and low-energy ecosystems. Over the past decade, parallel advances in DNA amplification, DNA sequencing and computing have enabled rapid progress on this problem, particularly through metagenomics and single-cell genomics. Although each approach suffers limitations, metagenomics and single-cell genomics are particularly powerful when combined synergistically. Studies focused on extreme environments have revealed the first substantial genomic information for several candidate phyla, encompassing putative acidophiles (Parvarchaeota), halophiles (Nanohaloarchaeota), thermophiles (Acetothermia, Aigarchaeota, Atribacteria, Calescamantes, Korarchaeota, and Fervidibacteria), and piezophiles (Gracilibacteria). These data have enabled insights into the biology of these organisms, including catabolic and anabolic potential, molecular adaptations to life in extreme environments, unique genomic features such as stop codon reassignments, and predictions about cell ultrastructure. In addition, the rapid expansion of genomic coverage enabled by these studies continues to yield insights into the early diversification of microbial lineages and the relationships within and between the phyla of Bacteria and Archaea. In the next 5 years, the genomic foliage within the tree of life will continue to grow and the study of yet-uncultivated candidate phyla will firmly transition into the post-genomic era.