Genomic insights to SAR86, an abundant and uncultivated marine bacterial lineage

Global metagenomic survey reveals a new bacterial candidate phylum in geothermal springs

Analysis of the increasing wealth of metagenomic data collected from diverse environments can lead to the discovery of novel branches on the tree of life. Here we analyse 5.2 Tb of metagenomic data collected globally to discover a novel bacterial phylum (‘Candidatus Kryptonia') found exclusively in high-temperature pH-neutral geothermal springs. This lineage had remained hidden as a taxonomic ‘blind spot' because of mismatches in the primers commonly used for ribosomal gene surveys. Genome reconstruction from metagenomic data combined with single-cell genomics results in several high-quality genomes representing four genera from the new phylum. Metabolic reconstruction indicates a heterotrophic lifestyle with conspicuous nutritional deficiencies, suggesting the need for metabolic complementarity with other microbes. Co-occurrence patterns identifies a number of putative partners, including an uncultured Armatimonadetes lineage. The discovery of Kryptonia within previously studied geothermal springs underscores the importance of globally sampled metagenomic data in detection of microbial novelty, and highlights the extraordinary diversity of microbial life still awaiting discovery.

The analysis of existing metagenomic data can lead to discovery of new microorganisms. Here, Eloe-Fadrosh et al. perform a large-scale analysis of global metagenomic data, followed by genome reconstruction and single-cell genomics, to describe a new bacterial phylum that inhabits geothermal springs.

Functional Profiling of Unfamiliar Microbial Communities Using a Validated De Novo Assembly Metatranscriptome Pipeline

BACKGROUND

Metatranscriptomic landscapes can provide insights in functional relationships within natural microbial communities. Analysis of complex metatranscriptome datasets of these communities poses a considerable bioinformatic challenge since they are non-restricted with a varying number of participating strains and species. For RNA-Seq data a standard approach is to align the generated reads to a set of closely related reference genomes. This only works well for microbial communities for which a near complete catalogue of reference genomes is available at a small evolutionary distance. In this study, we focus on the design of a validated de novo metatranscriptome assembly pipeline for single-end Illumina RNA-Seq data to obtain functional and taxonomic profiles of murine microbial communities.

RESULTS

The here developed de novo assembly metatranscriptome pipeline combined rRNA removal, IDBA-UD assembler, functional annotation and taxonomic classification. Different assemblers were tested and validated using RNA-Seq data from an in silico generated mock community and in vivo RNA-Seq data from a restricted microbial community taken from a mouse model colonized with Altered Schaedler Flora (ASF). Precision and recall of resulting gene expression, functional and taxonomic profiles were compared to those obtained with a standard alignment method. The validated pipeline was subsequently used to generate expression profiles from non-restricted cecal communities of four C57BL/6J mice fed on a high-fat high-protein diet spiked with an RNA-Seq data set from a well-characterized human sample. The spike in control was used to estimate precision and recall at assembly, functional and taxonomic level of non-restricted communities.

CONCLUSIONS

A generic de novo assembly pipeline for metatranscriptome data analysis was designed for microbial ecosystems, which can be applied for microbial metatranscriptome analysis in any chosen niche.

Expanding the World of Marine Bacterial and Archaeal Clades

Determining which microbial taxa are out there, where they live, and what they are doing is a driving approach in marine microbial ecology. The importance of these questions is underlined by concerted, large-scale, and global ocean sampling initiatives, for example the International Census of Marine Microbes, Ocean Sampling Day, or Tara Oceans. Given decades of effort, we know that the large majority of marine Bacteria and Archaea belong to about a dozen phyla. In addition to the classically culturable Bacteria and Archaea, at least 50 “clades,” at different taxonomic depths, exist. These account for the majority of marine microbial diversity, but there is still an underexplored and less abundant portion remaining. We refer to these hitherto unrecognized clades as unknown, as their boundaries, names, and classifications are not available. In this work, we were able to characterize up to 92 of these unknown clades found within the bacterial and archaeal phylogenetic diversity currently reported for marine water column environments. We mined the SILVA 16S rRNA gene datasets for sequences originating from the marine water column. Instead of the usual subjective taxa delineation and nomenclature methods, we applied the candidate taxonomic unit (CTU) circumscription system, along with a standardized nomenclature to the sequences in newly constructed phylogenetic trees. With this new phylogenetic and taxonomic framework, we performed an analysis of ICoMM rRNA gene amplicon datasets to gain insights into the global distribution of the new marine clades, their ecology, biogeography, and interaction with oceanographic variables. Most of the new clades we identified were interspersed by known taxa with cultivated members, whose genome sequences are available. This result encouraged us to perform metabolic predictions for the novel marine clades using the PICRUSt approach. Our work also provides an update on the taxonomy of several phyla and widely known marine clades as our CTU approach breaks down these randomly lumped clades into smaller objectively calculated subgroups. Finally, all taxa were classified and named following standards compatible with the Bacteriological Code rules, enhancing their digitization, and comparability with future microbial ecological and taxonomy studies.

A Phylogenomic Analysis of the Bacterial Phylum Fibrobacteres

The Fibrobacteres has been recognized as a bacterial phylum for over a decade, but little is known about the group beyond its environmental distribution, and characterization of its sole cultured representative genus, Fibrobacter, after which the phylum was named. Based on these incomplete data, it is thought that cellulose hydrolysis, anaerobic metabolism, and lack of motility are unifying features of the phylum. There are also contradicting views as to whether an uncultured sister lineage, candidate phylum TG3, should be included in the Fibrobacteres. Recently, chitin-degrading cultured representatives of TG3 were isolated from a hypersaline soda lake, and the genome of one species, Chitinivibrio alkaliphilus, sequenced and described in detail. Here, we performed a comparative analysis of Fibrobacter succinogenes, C. alkaliphilus and eight near or substantially complete Fibrobacteres/TG3 genomes of environmental populations recovered from termite gut, anaerobic digester, and sheep rumen metagenomes. We propose that TG3 should be amalgamated with the Fibrobacteres phylum based on robust monophyly of the two lineages and shared character traits. Polymer hydrolysis, using a distinctive set of glycoside hydrolases and binding domains, appears to be a prominent feature of members of the Fibrobacteres. Not all members of this phylum are strictly anaerobic as some termite gut Fibrobacteres have respiratory chains adapted to the microaerophilic conditions found in this habitat. Contrary to expectations, flagella-based motility is predicted to be an ancestral and common trait in this phylum and has only recently been lost in F. succinogenes and its relatives based on phylogenetic distribution of flagellar genes. Our findings extend current understanding of the Fibrobacteres and provide an improved basis for further investigation of this phylum.

Metagenomic Assembly Reveals Hosts of Antibiotic Resistance Genes and the Shared Resistome in Pig, Chicken, and Human Feces

The risk associated with antibiotic resistance disseminating from animal and human feces is an urgent public issue. In the present study, we sought to establish a pipeline for annotating antibiotic resistance genes (ARGs) based on metagenomic assembly to investigate ARGs and their co-occurrence with associated genetic elements. Genetic elements found on the assembled genomic fragments include mobile genetic elements (MGEs) and metal resistance genes (MRGs). We then explored the hosts of these resistance genes and the shared resistome of pig, chicken and human fecal samples. High levels of tetracycline, multidrug, erythromycin, and aminoglycoside resistance genes were discovered in these fecal samples. In particular, significantly high level of ARGs (7762 ×/Gb) was detected in adult chicken feces, indicating higher ARG contamination level than other fecal samples. Many ARGs arrangements (e.g., macA-macB and tetA-tetR) were discovered shared by chicken, pig and human feces. In addition, MGEs such as the aadA5-dfrA17-carrying class 1 integron were identified on an assembled scaffold of chicken feces, and are carried by human pathogens. Differential coverage binning analysis revealed significant ARG enrichment in adult chicken feces. A draft genome, annotated as multidrug resistant Escherichia coli, was retrieved from chicken feces metagenomes and was determined to carry diverse ARGs (multidrug, acriflavine, and macrolide). The present study demonstrates the determination of ARG hosts and the shared resistome from metagenomic data sets and successfully establishes the relationship between ARGs, hosts, and environments. This ARG annotation pipeline based on metagenomic assembly will help to bridge the knowledge gaps regarding ARG-associated genes and ARG hosts with metagenomic data sets. Moreover, this pipeline will facilitate the evaluation of environmental risks in the genetic context of ARGs.

Genome Reconstruction from Metagenomic Data Sets Reveals Novel Microbes in the Brackish Waters of the Caspian Sea

We present here the findings from a study of the microbiome of the southern basin of the Caspian Sea, the largest water body on Earth disconnected from any ocean and a brackish inland sea. By high-throughput metagenomics, we were able to reconstruct the genomes of representative microbes. The gross community structure (at the phylum level) was different from the structure of typical marine and freshwater communities in temperate open oceans, with the Caspian Sea having freshwater-like amounts of Actinobacteria and Alphaproteobacteria, while Gammaproteobacteria and Betaproteobacteria were present at intermediate levels. We assembled the genomes of several groups and provide detailed descriptions of partial genomes from Actinobacteria, Thaumarchaea, and Alphaproteobacteria. Most belonged to hitherto unknown groups, although they were related to either marine or freshwater groups. The phylogenetic placement of the Caspian genomes indicates that the organisms have multiple and separate phylogenetic origins and that they are related to organisms with both freshwater and marine lineages. Comparative recruitment from global aquatic metagenomes indicated that most Caspian microbes are endemic. However, some Caspian genomes were recruited significantly from either marine water (a member of the Alphaproteobacteria) or freshwater (a member of the Actinobacteria). Reciprocally, some genomes of other origins, such as the marine thaumarchaeon " Candidatus Nitrosopelagicus" or the actinobacterium "Candidatus Actinomarina," were recruited from the Caspian Sea, indicating some degree of overlap with the microbiota of other water bodies. Some of these microbes seem to have a remarkably widespread geographic and environmental distribution.

Trace Metal Acquisition by Marine Heterotrophic Bacterioplankton with Contrasting Trophic Strategies

Heterotrophic bacteria in the SAR11 and Roseobacter lineages shape the marine carbon, nitrogen, phosphorous, and sulfur cycles, yet they do so having adopted divergent ecological strategies. Currently, it is unknown whether these globally significant groups partition into specific niches with respect to micronutrients (e.g., trace metals) and how that may affect marine trace metal cycling. Here, we used comparative genomics to identify diverse iron, cobalt, nickel, copper, and zinc uptake capabilities in SAR11 and Roseobacter genomes and uncover surprising unevenness within and between lineages. The strongest predictors for the extent of the metal uptake gene content are the total number of transporters per genome, genome size, total metal transporters, and GC content, but numerous exceptions exist in both groups. Taken together, our results suggest that SAR11 have strongly minimized their trace metal uptake versatility, with high-affinity zinc uptake being a unique exception. The larger Roseobacter genomes have greater trace metal uptake versatility on average, but they also appear to have greater plasticity, resulting in phylogenetically similar genomes having largely different capabilities. Ultimately, phylogeny is predictive of the diversity and extent of 20 to 33% of all metal uptake systems, suggesting that specialization in metal utilization mostly occurred independently from overall lineage diversification in both SAR11 and Roseobacter. We interpret these results as reflecting relatively recent trace metal niche partitioning in both lineages, suggesting that concentrations and chemical forms of metals in the marine environment are important factors shaping the gene content of marine heterotrophic Alphaproteobacteria of the SAR11 and Roseobacter lineages.

Metagenomic analysis and functional characterization of the biogas microbiome using high throughput shotgun sequencing and a novel binning strategy

BACKGROUND

Biogas production is an economically attractive technology that has gained momentum worldwide over the past years. Biogas is produced by a biologically mediated process, widely known as "anaerobic digestion." This process is performed by a specialized and complex microbial community, in which different members have distinct roles in the establishment of a collective organization. Deciphering the complex microbial community engaged in this process is interesting both for unraveling the network of bacterial interactions and for applicability potential to the derived knowledge.

RESULTS

In this study, we dissect the bioma involved in anaerobic digestion by means of high throughput Illumina sequencing (~51 gigabases of sequence data), disclosing nearly one million genes and extracting 106 microbial genomes by a novel strategy combining two binning processes. Microbial phylogeny and putative taxonomy performed using >400 proteins revealed that the biogas community is a trove of new species. A new approach based on functional properties as per network representation was developed to assign roles to the microbial species. The organization of the anaerobic digestion microbiome is resembled by a funnel concept, in which the microbial consortium presents a progressive functional specialization while reaching the final step of the process (i.e., methanogenesis). Key microbial genomes encoding enzymes involved in specific metabolic pathways, such as carbohydrates utilization, fatty acids degradation, amino acids fermentation, and syntrophic acetate oxidation, were identified. Additionally, the analysis identified a new uncultured archaeon that was putatively related to Methanomassiliicoccales but surprisingly having a methylotrophic methanogenic pathway.

CONCLUSION

This study is a pioneer research on the phylogenetic and functional characterization of the microbial community populating biogas reactors. By applying for the first time high-throughput sequencing and a novel binning strategy, the identified genes were anchored to single genomes providing a clear understanding of their metabolic pathways and highlighting their involvement in anaerobic digestion. The overall research established a reference catalog of biogas microbial genomes that will greatly simplify future genomic studies.

Cellular adhesiveness and cellulolytic capacity in Anaerolineae revealed by omics-based genome interpretation

BACKGROUND

The Anaerolineae lineage of Chloroflexi had been identified as one of the core microbial populations in anaerobic digesters; however, the ecological role of the Anaerolineae remains uncertain due to the scarcity of isolates and annotated genome sequences. Our previous metatranscriptional analysis revealed this prevalent population that showed minimum involvement in the main pathways of cellulose hydrolysis and subsequent methanogenesis in the thermophilic cellulose fermentative consortium (TCF).

RESULTS

In further pursuit, five high-quality curated draft genomes (>98 % completeness) of this population, including two affiliated with the inaccessible lineage of SBR1031, were retrieved by sequence-based multi-dimensional coverage binning. Comparative genomic analyses revealed versatile genetic capabilities for carbohydrate-based fermentative lifestyle including key genes catalyzing cellulose hydrolysis in Anaerolinea phylotypes. However, the low transcriptional activities of carbohydrate-active genes (CAGs) excluded cellulolytic capability as the selective advantage for their prevalence in the community. Instead, a substantially active type VI pili (Tfp) assembly was observed. Expression of the tight adherence protein on the Tfp indicated its function for cellular attachment which was further testified to be more likely related to cell aggregation other than cellulose surface adhesion. Meanwhile, this Tfp structure was found not contributing to syntrophic methanogenesis. Members of the SBR1031 encoded key genes for acetogenic dehydrogenation that may allow ethanol to be used as a carbon source.

CONCLUSION

The common prevalence of Anaerolineae in anaerobic digesters should be originated from advantageous cellular adhesiveness enabled by Tfp assembly other than its potential as cellulose degrader or anaerobic syntrophs.

Complete nitrification by Nitrospira bacteria

Nitrification, the oxidation of ammonia via nitrite to nitrate, has always been considered to be a two-step process catalysed by chemolithoautotrophic microorganisms oxidizing either ammonia or nitrite. No known nitrifier carries out both steps, although complete nitrification should be energetically advantageous. This functional separation has puzzled microbiologists for a century. Here we report on the discovery and cultivation of a completely nitrifying bacterium from the genus Nitrospira, a globally distributed group of nitrite oxidizers. The genome of this chemolithoautotrophic organism encodes the pathways both for ammonia and nitrite oxidation, which are concomitantly activated during growth by ammonia oxidation to nitrate. Genes affiliated with the phylogenetically distinct ammonia monooxygenase and hydroxylamine dehydrogenase genes of Nitrospira are present in many environments and were retrieved on Nitrospira-contigs in new metagenomes from engineered systems. These findings fundamentally change our picture of nitrification and point to completely nitrifying Nitrospira as key components of nitrogen-cycling microbial communities.

Microbial Surface Colonization and Biofilm Development in Marine Environments

Biotic and abiotic surfaces in marine waters are rapidly colonized by microorganisms. Surface colonization and subsequent biofilm formation and development provide numerous advantages to these organisms and support critical ecological and biogeochemical functions in the changing marine environment. Microbial surface association also contributes to deleterious effects such as biofouling, biocorrosion, and the persistence and transmission of harmful or pathogenic microorganisms and their genetic determinants. The processes and mechanisms of colonization as well as key players among the surface-associated microbiota have been studied for several decades. Accumulating evidence indicates that specific cell-surface, cell-cell, and interpopulation interactions shape the composition, structure, spatiotemporal dynamics, and functions of surface-associated microbial communities. Several key microbial processes and mechanisms, including (i) surface, population, and community sensing and signaling, (ii) intraspecies and interspecies communication and interaction, and (iii) the regulatory balance between cooperation and competition, have been identified as critical for the microbial surface association lifestyle. In this review, recent progress in the study of marine microbial surface colonization and biofilm development is synthesized and discussed. Major gaps in our knowledge remain. We pose questions for targeted investigation of surface-specific community-level microbial features, answers to which would advance our understanding of surface-associated microbial community ecology and the biogeochemical functions of these communities at levels from molecular mechanistic details through systems biological integration.

Metagenome-assembled genomes uncover a global brackish microbiome

BACKGROUND

Microbes are main drivers of biogeochemical cycles in oceans and lakes. Although the genome is a foundation for understanding the metabolism, ecology and evolution of an organism, few bacterioplankton genomes have been sequenced, partly due to difficulties in cultivating them.

RESULTS

We use automatic binning to reconstruct a large number of bacterioplankton genomes from a metagenomic time-series from the Baltic Sea, one of world's largest brackish water bodies. These genomes represent novel species within typical freshwater and marine clades, including clades not previously sequenced. The genomes' seasonal dynamics follow phylogenetic patterns, but with fine-grained lineage-specific variations, reflected in gene-content. Signs of streamlining are evident in most genomes, and estimated genome sizes correlate with abundance variation across filter size fractions. Comparing the genomes with globally distributed metagenomes reveals significant fragment recruitment at high sequence identity from brackish waters in North America, but little from lakes or oceans. This suggests the existence of a global brackish metacommunity whose populations diverged from freshwater and marine relatives over 100,000 years ago, long before the Baltic Sea was formed (8000 years ago). This markedly contrasts to most Baltic Sea multicellular organisms, which are locally adapted populations of freshwater or marine counterparts.

CONCLUSIONS

We describe the gene content, temporal dynamics and biogeography of a large set of new bacterioplankton genomes assembled from metagenomes. We propose that brackish environments exert such strong selection that lineages adapted to them flourish globally with limited influence from surrounding aquatic communities.

Arsenic rich Himalayan hot spring metagenomics reveal genetically novel predator-prey genotypes

Bdellovibrio bacteriovorus are small Deltaproteobacteria that invade, kill and assimilate their prey. Metagenomic assembly analysis of the microbial mats of an arsenic rich, hot spring was performed to describe the genotypes of the predator Bdellovibrio and the ecogenetically adapted taxa Enterobacter. The microbial mats were enriched with Bdellovibrio (1.3%) and several Gram-negative bacteria including Bordetella (16%), Enterobacter (6.8%), Burkholderia (4.8%), Acinetobacter (2.3%) and Yersinia (1%). A high-quality (47 contigs, 25X coverage; 3.5 Mbp) draft genome of Bdellovibrio (strain ArHS; Arsenic Hot Spring) was reassembled, which lacked the marker gene Bd0108 associated with the usual method of prey interaction and invasion for this genus, while maintaining genes coding for the hydrolytic enzymes necessary for prey assimilation. By filtering microbial mat samples (< 0.45 μm) to enrich for small predatory cell sizes, we observed Bdellovibrio-like cells attached side-on to E. coli through electron microscopy. Furthermore, a draft pan-genome of the dominant potential host taxon, Enterobacter cloacae ArHS (4.8 Mb), along with three of its viral genotypes (n = 3; 42 kb, 49 kb and 50 kb), was assembled. These data were further used to analyse the population level evolutionary dynamics (taxonomical and functional) of reconstructed genotypes.

Phenotypic differentiation of gastrointestinal microbes is reflected in their encoded metabolic repertoires

BACKGROUND

The human gastrointestinal tract harbors a diverse microbial community, in which metabolic phenotypes play important roles for the human host. Recent developments in meta-omics attempt to unravel metabolic roles of microbes by linking genotypic and phenotypic characteristics. This connection, however, still remains poorly understood with respect to its evolutionary and ecological context.

RESULTS

We generated automatically refined draft genome-scale metabolic models of 301 representative intestinal microbes in silico. We applied a combination of unsupervised machine-learning and systems biology techniques to study individual and global differences in genomic content and inferred metabolic capabilities. Based on the global metabolic differences, we found that energy metabolism and membrane synthesis play important roles in delineating different taxonomic groups. Furthermore, we found an exponential relationship between phylogeny and the reaction composition, meaning that closely related microbes of the same genus can exhibit pronounced differences with respect to their metabolic capabilities while at the family level only marginal metabolic differences can be observed. This finding was further substantiated by the metabolic divergence within different genera. In particular, we could distinguish three sub-type clusters based on membrane and energy metabolism within the Lactobacilli as well as two clusters within the Bifidobacteria and Bacteroides.

CONCLUSIONS

We demonstrate that phenotypic differentiation within closely related species could be explained by their metabolic repertoire rather than their phylogenetic relationships. These results have important implications in our understanding of the ecological and evolutionary complexity of the human gastrointestinal microbiome.

Comparative Genomics of Candidate Phylum TM6 Suggests That Parasitism Is Widespread and Ancestral in This Lineage

Candidate phylum TM6 is a major bacterial lineage recognized through culture-independent rRNA surveys to be low abundance members in a wide range of habitats; however, they are poorly characterized due to a lack of pure culture representatives. Two recent genomic studies of TM6 bacteria revealed small genomes and limited gene repertoire, consistent with known or inferred dependence on eukaryotic hosts for their metabolic needs. Here, we obtained additional near-complete genomes of TM6 populations from agricultural soil and upflow anaerobic sludge blanket reactor metagenomes which, together with the two publicly available TM6 genomes, represent seven distinct family level lineages in the TM6 phylum. Genome-based phylogenetic analysis confirms that TM6 is an independent phylum level lineage in the bacterial domain, possibly affiliated with the Patescibacteria superphylum. All seven genomes are small (1.0–1.5 Mb) and lack complete biosynthetic pathways for various essential cellular building blocks including amino acids, lipids, and nucleotides. These and other features identified in the TM6 genomes such as a degenerated cell envelope, ATP/ADP translocases for parasitizing host ATP pools, and protein motifs to facilitate eukaryotic host interactions indicate that parasitism is widespread in this phylum. Phylogenetic analysis of ATP/ADP translocase genes suggests that the ancestral TM6 lineage was also parasitic. We propose the name Dependentiae (phyl. nov.) to reflect dependence of TM6 bacteria on host organisms.

Complete nitrification by a single microorganism

Nitrification is a two-step process where ammonia is considered to first be oxidized to nitrite by ammonia-oxidizing bacteria (AOB) and/or archaea (AOA), and subsequently to nitrate by nitrite-oxidizing bacteria (NOB). Described by Winogradsky already in 18901, this division of labour between the two functional groups is a generally accepted characteristic of the biogeochemical nitrogen cycle2. Complete oxidation of ammonia to nitrate in one organism (complete ammonia oxidation; comammox) is energetically feasible and it was postulated that this process could occur under conditions selecting for species with lower growth-rates but higher growth-yields than canonical ammonia-oxidizing microorganisms3. Still, organisms catalysing this process have not yet been discovered. Here, we report the enrichment and initial characterization of two Nitrospira species that encode all enzymes necessary for ammonia oxidation via nitrite to nitrate in their genomes, and indeed completely oxidize ammonium to nitrate to conserve energy. Their ammonia monooxygenase (AMO) enzymes are phylogenetically distinct from currently identified AMOs, rendering recent acquisition by horizontal gene transfer from known ammonia-oxidizing microorganisms unlikely. We also found highly similar amoA sequences (encoding the AMO subunit A) in public sequence databases, which were apparently misclassified as methane monooxygenases. This recognition of a novel amoA sequence group will lead to an improved understanding on the environmental abundance and distribution of ammonia-oxidizing microorganisms. Furthermore, the discovery of the long-sought-after comammox process will change our perception of the nitrogen cycle.

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.

Single-cell technologies in environmental omics

Environmental studies are primarily done by culturing isolated microorganisms or by amplifying and sequencing conserved genes. Difficulties understanding the complexity of large numbers of various microorganisms in an environment led to the development of techniques to enrich specific microorganisms for upstream analysis, ultimately leading to single-cell isolation and analyses. We discuss the significance of single-cell technologies in omics studies with focus on metagenomics and metatranscriptomics. We propose that by reducing sample heterogeneity using single-cell genomics, metaomic studies can be simplified.

Spatially extensive microbial biogeography of the Indian Ocean provides insights into the unique community structure of a pristine coral atoll

Microorganisms act both as drivers and indicators of perturbations in the marine environment. In an effort to establish baselines to predict the response of marine habitats to environmental change, here we report a broad survey of microbial diversity across the Indian Ocean, including the first microbial samples collected in the pristine lagoon of Salomon Islands, Chagos Archipelago. This was the first large-scale ecogenomic survey aboard a private yacht employing a ‘citizen oceanography’ approach and tools and protocols easily adapted to ocean going sailboats. Our data highlighted biogeographic patterns in microbial community composition across the Indian Ocean. Samples from within the Salomon Islands lagoon contained a community which was different even from adjacent samples despite constant water exchange, driven by the dominance of the photosynthetic cyanobacterium Synechococcus. In the lagoon, Synechococcus was also responsible for driving shifts in the metatranscriptional profiles. Enrichment of transcripts related to photosynthesis and nutrient cycling indicated bottom-up controls of community structure. However a five-fold increase in viral transcripts within the lagoon during the day, suggested a concomitant top-down control by bacteriophages. Indeed, genome recruitment against Synechococcus reference genomes suggested a role of viruses in providing the ecological filter for determining the β-diversity patterns in this system.

Anvi'o: an advanced analysis and visualization platform for 'omics data

Advances in high-throughput sequencing and ‘omics technologies are revolutionizing studies of naturally occurring microbial communities. Comprehensive investigations of microbial lifestyles require the ability to interactively organize and visualize genetic information and to incorporate subtle differences that enable greater resolution of complex data. Here we introduce anvi’o, an advanced analysis and visualization platform that offers automated and human-guided characterization of microbial genomes in metagenomic assemblies, with interactive interfaces that can link ‘omics data from multiple sources into a single, intuitive display. Its extensible visualization approach distills multiple dimensions of information about each contig, offering a dynamic and unified work environment for data exploration, manipulation, and reporting. Using anvi’o, we re-analyzed publicly available datasets and explored temporal genomic changes within naturally occurring microbial populations through de novo characterization of single nucleotide variations, and linked cultivar and single-cell genomes with metagenomic and metatranscriptomic data. Anvi’o is an open-source platform that empowers researchers without extensive bioinformatics skills to perform and communicate in-depth analyses on large ‘omics datasets.

Ocean plankton. Structure and function of the global ocean microbiome

Microbes are dominant drivers of biogeochemical processes, yet drawing a global picture of functional diversity, microbial community structure, and their ecological determinants remains a grand challenge. We analyzed 7.2 terabases of metagenomic data from 243 Tara Oceans samples from 68 locations in epipelagic and mesopelagic waters across the globe to generate an ocean microbial reference gene catalog with >40 million nonredundant, mostly novel sequences from viruses, prokaryotes, and picoeukaryotes. Using 139 prokaryote-enriched samples, containing >35,000 species, we show vertical stratification with epipelagic community composition mostly driven by temperature rather than other environmental factors or geography. We identify ocean microbial core functionality and reveal that >73% of its abundance is shared with the human gut microbiome despite the physicochemical differences between these two ecosystems.

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.

Needles in the blue sea: sub-species specificity in targeted protein biomarker analyses within the vast oceanic microbial metaproteome

Proteomics has great potential for studies of marine microbial biogeochemistry, yet high microbial diversity in many locales presents us with unique challenges. We addressed this challenge with a targeted metaproteomics workflow for NtcA and P-II, two nitrogen regulatory proteins, and demonstrated its application for cyanobacterial taxa within microbial samples from the Central Pacific Ocean. Using METATRYP, an open-source Python toolkit, we examined the number of shared (redundant) tryptic peptides in representative marine microbes, with the number of tryptic peptides shared between different species typically being 1% or less. The related cyanobacteria Prochlorococcus and Synechococcus shared an average of 4.8 ± 1.9% of their tryptic peptides, while shared intraspecies peptides were higher, 13 ± 15% shared peptides between 12 Prochlorococcus genomes. An NtcA peptide was found to target multiple cyanobacteria species, whereas a P-II peptide showed specificity to the high-light Prochlorococcus ecotype. Distributions of NtcA and P-II in the Central Pacific Ocean were similar except at the Equator likely due to differential nitrogen stress responses between Prochlorococcus and Synechococcus. The number of unique tryptic peptides coded for within three combined oceanic microbial metagenomes was estimated to be ∼4 × 10(7) , 1000-fold larger than an individual microbial proteome and 27-fold larger than the human proteome, yet still 20 orders of magnitude lower than the peptide diversity possible in all protein space, implying that peptide mapping algorithms should be able to withstand the added level of complexity in metaproteomic samples.

Evolutionary ecology of the marine Roseobacter clade

SUMMARY

Members of the Roseobacter clade are equipped with a tremendous diversity of metabolic capabilities, which in part explains their success in so many different marine habitats. Ideas on how this diversity evolved and is maintained are reviewed, focusing on recent evolutionary studies exploring the timing and mechanisms of Roseobacter ecological diversification.

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.

Widespread Distribution and Expression of Gamma A (UMB), an Uncultured, Diazotrophic, γ-Proteobacterial nifH Phylotype

Marine dinitrogen (N2) fixation studies have focused nearly exclusively on cyanobacterial diazotrophs; however γ-proteobacteria are an abundant component of the marine community and have been largely overlooked until recently. Here we present a phylogenetic analysis of all nifH γ-proteobacterial sequences available in public databases and qPCR data of a γ-proteobacterial phylotype, Gamma A (UMB), obtained during several research cruises. Our analysis revealed a complex diversity of diazotrophic γ-proteobacteria. One phylotype in particular, Gamma A, was described in several traditional and quantitative PCR studies. Though several γ-proteobacterial nifH sequences have been described as laboratory contaminants, Gamma A is part of a large cluster of sequences isolated from marine environments and distantly related to the clade of contaminants. Using a TaqMan probe and primer set, Gamma A nifH DNA abundance and expression were analyzed in nearly 1000 samples collected during 15 cruises to the Atlantic and Pacific Oceans. The data showed that Gamma A is an active, cosmopolitan diazotroph found throughout oxygenated, oligotrophic waters reaching maximum abundances of 8 x 104 nifH DNA copies l-1 and 5 x 105 nifH transcript copies l-1. Gamma A nifH transcript abundances were on average 3 fold higher than nifH DNA abundances. The widespread distribution and activity of Gamma A indicate that it has potential to be a globally important N2 fixing organism.

Capturing prokaryotic dark matter genomes

Prokaryotes are the most diverse and abundant cellular life forms on Earth. Most of them, identified by indirect molecular approaches, belong to microbial dark matter. The advent of metagenomic and single-cell genomic approaches has highlighted the metabolic capabilities of numerous members of this dark matter through genome reconstruction. Thus, linking functions back to the species has revolutionized our understanding of how ecosystem function is sustained by the microbial world. This review will present discoveries acquired through the illumination of prokaryotic dark matter genomes by these innovative approaches.

Zonation of bacterioplankton communities along aging upwelled water in the northern Benguela upwelling

Upwelling areas are shaped by enhanced primary production in surface waters, accompanied by a well-investigated planktonic succession. Although bacteria play an important role in biogeochemical cycles of upwelling systems, little is known about bacterial community composition and its development during upwelling events. The aim of this study was to investigate the succession of bacterial assemblages in aging upwelled water of the Benguela upwelling from coastal to offshore sites. Water from the upper mixed layer at 12 stations was sampled along two transects from the origin of the upwelling to a distance of 220 km. 16S rRNA gene amplicon sequencing was then used in a bacterial diversity analysis and major bacterial taxa were quantified by catalyzed reporter deposition-fluorescence in situ hybridization. Additionally, bacterial cell numbers and bacterial production were assessed. Community statistical analysis revealed a reproducible zonation along the two transects, with four clusters of significantly different microbial assemblages. Clustering was mainly driven by phytoplankton composition and abundance. Similar to the temporal succession that occurs during phytoplankton blooms in temperate coastal waters, operational taxonomic units (OTUs) affiliated with Bacteroidetes and Gammaproteobacteria were dominant during algal blooming whereas "Pelagibacterales" were highly abundant in regions with low algal abundance. The most dominant heterotrophic OTU (9% of all reads) was affiliated with "Pelagibacterales" and showed a strong negative correlation with phytoplankton. By contrast, the second most abundant heterotrophic OTU (6% of all reads) was affiliated with the phylum Verrucomicrobia and correlated positively with phytoplankton. Together with the close relation of bacterial production and phytoplankton abundance, our results showed that bacterial community dynamics is strongly driven by the development and composition of the phytoplankton community.

Characterization of an Unconventional Rhodopsin from the Freshwater Actinobacterium Rhodoluna lacicola

Rhodopsin-encoding microorganisms are common in many environments. However, knowing that rhodopsin genes are present provides little insight into how the host cells utilize light. The genome of the freshwater actinobacterium Rhodoluna lacicola encodes a rhodopsin of the uncharacterized actinorhodopsin family. We hypothesized that actinorhodopsin was a light-activated proton pump and confirmed this by heterologously expressing R. lacicola actinorhodopsin in retinal-producing Escherichia coli. However, cultures of R. lacicola did not pump protons, even though actinorhodopsin mRNA and protein were both detected. Proton pumping in R. lacicola was induced by providing exogenous retinal, suggesting that the cells lacked the retinal cofactor. We used high-performance liquid chromatography (HPLC) and oxidation of accessory pigments to confirm that R. lacicola does not synthesize retinal. These results suggest that in some organisms, the actinorhodopsin gene is constitutively expressed, but rhodopsin-based light capture may require cofactors obtained from the environment.

IMPORTANCE

Up to 70% of microbial genomes in some environments are predicted to encode rhodopsins. Because most microbial rhodopsins are light-activated proton pumps, the prevalence of this gene suggests that in some environments, most microorganisms respond to or utilize light energy. Actinorhodopsins were discovered in an analysis of freshwater metagenomic data and subsequently identified in freshwater actinobacterial cultures. We hypothesized that actinorhodopsin from the freshwater actinobacterium Rhodoluna lacicola was a light-activated proton pump and confirmed this by expressing actinorhodopsin in retinal-producing Escherichia coli. Proton pumping in R. lacicola was induced only after both light and retinal were provided, suggesting that the cells lacked the retinal cofactor. These results indicate that photoheterotrophy in this organism and others may require cofactors obtained from the environment.

How do divergent ecological strategies emerge among marine bacterioplankton lineages?

Heterotrophic bacteria in pelagic marine environments are frequently categorized into two canonical ecological groups: patch-associated and free-living. This framework provides a conceptual basis for understanding bacterial utilization of oceanic organic matter. Some patch-associated bacteria are ecologically linked with eukaryotic phytoplankton, and this observation fits with predicted coincidence of their genome expansion with marine phytoplankton diversification. By contrast, free-living bacteria in today's oceans typically live singly with streamlined metabolic and regulatory functions that allow them to grow in nutrient-poor seawater. Recent analyses of marine Alphaproteobacteria suggest that some free-living bacterioplankton lineages evolved from patch-associated ancestors up to several hundred million years ago. While evolutionary analyses agree with the hypothesis that natural selection has maintained these distinct ecological strategies and genomic traits in present-day populations, they do not rule out a major role for genetic drift in driving ancient ecological switches. These two evolutionary forces may have acted on ocean bacteria at different geological time scales and under different geochemical constraints, with possible implications for future adaptations to a changing ocean. New evolutionary models and genomic data are leading to a more comprehensive understanding of marine bacterioplankton evolutionary history.

Consequences of increased terrestrial dissolved organic matter and temperature on bacterioplankton community composition during a Baltic Sea mesocosm experiment

Predicted increases in runoff of terrestrial dissolved organic matter (DOM) and sea surface temperatures implicate substantial changes in energy fluxes of coastal marine ecosystems. Despite marine bacteria being critical drivers of marine carbon cycling, knowledge of compositional responses within bacterioplankton communities to such disturbances is strongly limited. Using 16S rRNA gene pyrosequencing, we examined bacterioplankton population dynamics in Baltic Sea mesocosms with treatments combining terrestrial DOM enrichment and increased temperature. Among the 200 most abundant taxa, 62 % either increased or decreased in relative abundance under changed environmental conditions. For example, SAR11 and SAR86 populations proliferated in combined increased terrestrial DOM/temperature mesocosms, while the hgcI and CL500-29 clades (Actinobacteria) decreased in the same mesocosms. Bacteroidetes increased in both control mesocosms and in the combined increased terrestrial DOM/temperature mesocosms. These results indicate considerable and differential responses among distinct bacterial populations to combined climate change effects, emphasizing the potential of such effects to induce shifts in ecosystem function and carbon cycling in the future Baltic Sea.

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.

Single cell genomics indicates horizontal gene transfer and viral infections in a deep subsurface Firmicutes population

A major fraction of Earth's prokaryotic biomass dwells in the deep subsurface, where cellular abundances per volume of sample are lower, metabolism is slower, and generation times are longer than those in surface terrestrial and marine environments. How these conditions impact biotic interactions and evolutionary processes is largely unknown. Here we employed single cell genomics to analyze cell-to-cell genome content variability and signatures of horizontal gene transfer (HGT) and viral infections in five cells of Candidatus Desulforudis audaxviator, which were collected from a 3 km-deep fracture water in the 2.9 Ga-old Witwatersrand Basin of South Africa. Between 0 and 32% of genes recovered from single cells were not present in the original, metagenomic assembly of Desulforudis, which was obtained from a neighboring subsurface fracture. We found a transposable prophage, a retron, multiple clustered regularly interspaced short palindromic repeats (CRISPRs) and restriction-modification systems, and an unusually high frequency of transposases in the analyzed single cell genomes. This indicates that recombination, HGT and viral infections are prevalent evolutionary events in the studied population of microorganisms inhabiting a highly stable deep subsurface environment.

Pan-genome dynamics of Pseudomonas gene complements enriched across hexachlorocyclohexane dumpsite

Background

Phylogenetic heterogeneity across Pseudomonas genus is complemented by its diverse genome architecture enriched by accessory genetic elements (plasmids, transposons, and integrons) conferring resistance across this genus. Here, we sequenced a stress tolerant genotype i.e. Pseudomonas sp. strain RL isolated from a hexachlorocyclohexane (HCH) contaminated pond (45 mg of total HCH g⁻¹ sediment) and further compared its gene repertoire with 17 reference ecotypes belonging to P. stutzeri, P. mendocina, P. aeruginosa, P. psychrotolerans and P. denitrificans, representing metabolically diverse ecosystems (i.e. marine, clinical, and soil/sludge). Metagenomic data from HCH contaminated pond sediment and similar HCH contaminated sites were further used to analyze the pan-genome dynamics of Pseudomonas genotypes enriched across increasing HCH gradient.

Results

Although strain RL demonstrated clear species demarcation (ANI ≤ 80.03%) from the rest of its phylogenetic relatives, it was found to be closest to P. stutzeri clade which was further complemented functionally. Comparative functional analysis elucidated strain specific enrichment of metabolic pathways like α-linoleic acid degradation and carbazole degradation in Pseudomonas sp. strain RL and P. stutzeri XLDN-R, respectively. Composition based methods (%codon bias and %G + C difference) further highlighted the significance of horizontal gene transfer (HGT) in evolution of nitrogen metabolism, two-component system (TCS) and methionine metabolism across the Pseudomonas genomes used in this study. An intact mobile class-I integron (3,552 bp) with a captured gene cassette encoding for dihydrofolate reductase (dhfra1) was detected in strain RL, distinctly demarcated from other integron harboring species (i.e. P. aeruginosa, P. stutzeri, and P. putida). Mobility of this integron was confirmed by its association with Tnp21-like transposon (95% identity) suggesting stress specific mobilization across HCH contaminated sites. Metagenomics data from pond sediment and recently surveyed HCH adulterated soils revealed the in situ enrichment of integron associated transposase gene (TnpA6100) across increasing HCH contamination (0.7 to 450 mg HCH g⁻¹ of soil).

Conclusions

Unlocking the potential of comparative genomics supplemented with metagenomics, we have attempted to resolve the environment and strain specific demarcations across 18 Pseudomonas gene complements. Pan-genome analyses of these strains indicate at astoundingly diverse metabolic strategies and provide genetic basis for the cosmopolitan existence of this taxon.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1488-2) contains supplementary material, which is available to authorized users.

A taxonomic framework for emerging groups of ecologically important marine gammaproteobacteria based on the reconstruction of evolutionary relationships using genome-scale data

In recent years a large number of isolates were obtained from saline environments that are phylogenetically related to distinct clades of oligotrophic marine gammaproteobacteria, which were originally identified in seawater samples using cultivation independent methods and are characterized by high seasonal abundances in coastal environments. To date a sound taxonomic framework for the classification of these ecologically important isolates and related species in accordance with their evolutionary relationships is missing. In this study we demonstrate that a reliable allocation of members of the oligotrophic marine gammaproteobacteria (OMG) group and related species to higher taxonomic ranks is possible by phylogenetic analyses of whole proteomes but also of the RNA polymerase beta subunit, whereas phylogenetic reconstructions based on 16S rRNA genes alone resulted in unstable tree topologies with only insignificant bootstrap support. The identified clades could be correlated with distinct phenotypic traits illustrating an adaptation to common environmental factors in their evolutionary history. Genome wide gene-content analyses revealed the existence of two distinct ecological guilds within the analyzed lineage of marine gammaproteobacteria which can be distinguished by their trophic strategies. Based on our results a novel order within the class Gammaproteobacteria is proposed, which is designated Cellvibrionales ord. nov. and comprises the five novel families Cellvibrionaceae fam. nov., Halieaceae fam. nov., Microbulbiferaceae fam. nov., Porticoccaceae fam. nov., and Spongiibacteraceae fam. nov.

Trichodesmium genome maintains abundant, widespread noncoding DNA in situ, despite oligotrophic lifestyle

Understanding the evolution of the free-living, cyanobacterial, diazotroph Trichodesmium is of great importance because of its critical role in oceanic biogeochemistry and primary production. Unlike the other >150 available genomes of free-living cyanobacteria, only 63.8% of the Trichodesmium erythraeum (strain IMS101) genome is predicted to encode protein, which is 20-25% less than the average for other cyanobacteria and nonpathogenic, free-living bacteria. We use distinctive isolates and metagenomic data to show that low coding density observed in IMS101 is a common feature of the Trichodesmium genus, both in culture and in situ. Transcriptome analysis indicates that 86% of the noncoding space is expressed, although the function of these transcripts is unclear. The density of noncoding, possible regulatory elements predicted in Trichodesmium, when normalized per intergenic kilobase, was comparable and twofold higher than that found in the gene-dense genomes of the sympatric cyanobacterial genera Synechococcus and Prochlorococcus, respectively. Conserved Trichodesmium noncoding RNA secondary structures were predicted between most culture and metagenomic sequences, lending support to the structural conservation. Conservation of these intergenic regions in spatiotemporally separated Trichodesmium populations suggests possible genus-wide selection for their maintenance. These large intergenic spacers may have developed during intervals of strong genetic drift caused by periodic blooms of a subset of genotypes, which may have reduced effective population size. Our data suggest that transposition of selfish DNA, low effective population size, and high-fidelity replication allowed the unusual "inflation" of noncoding sequence observed in Trichodesmium despite its oligotrophic lifestyle.

The use of evolutionary approaches to understand single cell genomes

The vast majority of environmental bacteria and archaea remain uncultivated, yet their genome sequences are rapidly becoming available through single cell sequencing technologies. Reconstructing metabolism is one common way to make use of genome sequences of ecologically important bacteria, but molecular evolutionary analysis is another approach that, while currently underused, can reveal important insights into the function of these uncultivated microbes in nature. Because genome sequences from single cells are often incomplete, metabolic reconstruction based on genome content can be compromised. However, this problem does not necessarily impede the use of phylogenomic and population genomic approaches that are based on patterns of polymorphisms and substitutions at nucleotide and amino acid sites. These approaches explore how various evolutionary forces act to assemble genetic diversity within and between lineages. In this mini-review, I present examples illustrating the benefits of analyzing single cell genomes using evolutionary approaches.

Microbial community successional patterns in beach sands impacted by the Deepwater Horizon oil spill

Although petroleum hydrocarbons discharged from the Deepwater Horizon (DWH) blowout were shown to have a pronounced impact on indigenous microbial communities in the Gulf of Mexico, effects on nearshore or coastal ecosystems remain understudied. This study investigated the successional patterns of functional and taxonomic diversity for over 1 year after the DWH oil was deposited on Pensacola Beach sands (FL, USA), using metagenomic and 16S rRNA gene amplicon techniques. Gamma- and Alphaproteobacteria were enriched in oiled sediments, in corroboration of previous studies. In contrast to previous studies, we observed an increase in the functional diversity of the community in response to oil contamination and a functional transition from generalist populations within 4 months after oil came ashore to specialists a year later, when oil was undetectable. At the latter time point, a typical beach community had reestablished that showed little to no evidence of oil hydrocarbon degradation potential, was enriched in archaeal taxa known to be sensitive to xenobiotics, but differed significantly from the community before the oil spill. Further, a clear succession pattern was observed, where early responders to oil contamination, likely degrading aliphatic hydrocarbons, were replaced after 3 months by populations capable of aromatic hydrocarbon decomposition. Collectively, our results advance the understanding of how natural benthic microbial communities respond to crude oil perturbation, supporting the specialization-disturbance hypothesis; that is, the expectation that disturbance favors generalists, while providing (microbial) indicator species and genes for the chemical evolution of oil hydrocarbons during degradation and weathering.

Illumina sequencing-based analysis of free-living bacterial community dynamics during an Akashiwo sanguine bloom in Xiamen sea, China

Although phytoplankton are the major source of marine dissolved organic matter (DOM), their blooms are a global problem that can greatly affect marine ecological systems, especially free-living bacteria, which are the primary DOM degraders. In this study, we analyzed free-living bacterial communities from Xiamen sea during an Akashiwo sanguine bloom using Illumina MiSeq sequencing of 16S rRNA gene amplicons. The bloom was probably stimulated by low salinity and ended after abatement of eutrophication pollution. A total of 658,446 sequence reads and 11,807 OTUs were obtained in both bloom and control samples with Alpha-proteobacteria and Gamma-proteobacteria being the predominant classes detected. The bloom decreased bacterial diversity, increased species evenness, and significantly changed the bacterial community structure. Bacterial communities within the bloom were more homogeneous than those within the control area. The bacteria stimulated by this bloom included the SAR86 and SAR116 clades and the AEGEAN-169 marine group, but a few were suppressed. In addition, many bacteria known to be associated with phytoplankton were detected only in the bloom samples. This study revealed the great influence of an A. sanguinea bloom on free-living bacterial communities, and provided new insights into the relationship between bacteria and A. sanguinea in marine ecosystems.

Disentangling seasonal bacterioplankton population dynamics by high-frequency sampling

Multiyear comparisons of bacterioplankton succession reveal that environmental conditions drive community shifts with repeatable patterns between years. However, corresponding insight into bacterioplankton dynamics at a temporal resolution relevant for detailed examination of variation and characteristics of specific populations within years is essentially lacking. During 1 year, we collected 46 samples in the Baltic Sea for assessing bacterial community composition by 16S rRNA gene pyrosequencing (nearly twice weekly during productive season). Beta-diversity analysis showed distinct clustering of samples, attributable to seemingly synchronous temporal transitions among populations (populations defined by 97% 16S rRNA gene sequence identity). A wide spectrum of bacterioplankton dynamics was evident, where divergent temporal patterns resulted both from pronounced differences in relative abundance and presence/absence of populations. Rates of change in relative abundance calculated for individual populations ranged from 0.23 to 1.79 day(-1) . Populations that were persistently dominant, transiently abundant or generally rare were found in several major bacterial groups, implying evolution has favoured a similar variety of life strategies within these groups. These findings suggest that high temporal resolution sampling allows constraining the timescales and frequencies at which distinct populations transition between being abundant or rare, thus potentially providing clues about physical, chemical or biological forcing on bacterioplankton community structure.

First genomic insights into members of a candidate bacterial phylum responsible for wastewater bulking

Filamentous cells belonging to the candidate bacterial phylum KSB3 were previously identified as the causative agent of fatal filament overgrowth (bulking) in a high-rate industrial anaerobic wastewater treatment bioreactor. Here, we obtained near complete genomes from two KSB3 populations in the bioreactor, including the dominant bulking filament, using differential coverage binning of metagenomic data. Fluorescence in situ hybridization with 16S rRNA-targeted probes specific for the two populations confirmed that both are filamentous organisms. Genome-based metabolic reconstruction and microscopic observation of the KSB3 filaments in the presence of sugar gradients indicate that both filament types are Gram-negative, strictly anaerobic fermenters capable of non-flagellar based gliding motility, and have a strikingly large number of sensory and response regulator genes. We propose that the KSB3 filaments are highly sensitive to their surroundings and that cellular processes, including those causing bulking, are controlled by external stimuli. The obtained genomes lay the foundation for a more detailed understanding of environmental cues used by KSB3 filaments, which may lead to more robust treatment options to prevent bulking.

Genomic and proteomic characterization of "Candidatus Nitrosopelagicus brevis": an ammonia-oxidizing archaeon from the open ocean

Thaumarchaeota are among the most abundant microbial cells in the ocean, but difficulty in cultivating marine Thaumarchaeota has hindered investigation into the physiological and evolutionary basis of their success. We report here a closed genome assembled from a highly enriched culture of the ammonia-oxidizing pelagic thaumarchaeon CN25, originating from the open ocean. The CN25 genome exhibits strong evidence of genome streamlining, including a 1.23-Mbp genome, a high coding density, and a low number of paralogous genes. Proteomic analysis recovered nearly 70% of the predicted proteins encoded by the genome, demonstrating that a high fraction of the genome is translated. In contrast to other minimal marine microbes that acquire, rather than synthesize, cofactors, CN25 encodes and expresses near-complete biosynthetic pathways for multiple vitamins. Metagenomic fragment recruitment indicated the presence of DNA sequences >90% identical to the CN25 genome throughout the oligotrophic ocean. We propose the provisional name "Candidatus Nitrosopelagicus brevis" str. CN25 for this minimalist marine thaumarchaeon and suggest it as a potential model system for understanding archaeal adaptation to the open ocean.

Statistical research on the bioactivity of new marine natural products discovered during the 28 years from 1985 to 2012

Every year, hundreds of new compounds are discovered from the metabolites of marine organisms. Finding new and useful compounds is one of the crucial drivers for this field of research. Here we describe the statistics of bioactive compounds discovered from marine organisms from 1985 to 2012. This work is based on our database, which contains information on more than 15,000 chemical substances including 4196 bioactive marine natural products. We performed a comprehensive statistical analysis to understand the characteristics of the novel bioactive compounds and detail temporal trends, chemical structures, species distribution, and research progress. We hope this meta-analysis will provide useful information for research into the bioactivity of marine natural products and drug development.

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.