Single-cell and metagenomic analyses indicate a fermentative and saccharolytic lifestyle for members of the OP9 lineage

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.

Metabolic and microbial community dynamics during the anaerobic digestion of maize silage in a two-phase process

Two-phasic anaerobic digestion processes (hydrolysis/acidogenesis separated from acetogenesis/methanogenesis) can be used for biogas production on demand or a combined chemicals/bioenergy production. For an effective process control, detailed knowledge about the microbial catalysts and their correlation to process conditions is crucial. In this study, maize silage was digested in a two-phase process and interrelationships between process parameters and microbial communities were revealed. In the first-phase reactor, alternating metabolic periods were observed which emerged independently from the feeding frequency. During the L-period, up to 11.8 g L(-1) lactic acid was produced which significantly correlated to lactic acid bacteria of the genus Lactobacillus as the most abundant community members. During the alternating G-period, the production of volatile fatty acids (up to 5.3, 4.0 and 3.1 g L(-1) for propionic, n-butyric and n-caproic acid, respectively) dominated accompanied by a high gas production containing up to 28 % hydrogen. The relative abundance of various Clostridiales increased during this metabolic period. In the second-phase reactor, the metabolic fluctuations of the first phase were smoothed out resulting in a stable biogas production as well as stable bacterial and methanogenic communities. However, the biogas composition followed the metabolic dynamics of the first phase: the hydrogen content increased during the L-period whereas highest CH4/CO2 ratios (up to 2.8) were reached during the G-period. Aceticlastic Methanosaeta as well as hydrogenotrophic Methanoculleus and Methanobacteriaceae were identified as dominant methanogens. Consequently, a directed control of the first-phase stabilizing desired metabolic states can lead to an enhanced productivity regarding chemicals and bioenergy.

Monodisperse Picoliter Droplets for Low-Bias and Contamination-Free Reactions in Single-Cell Whole Genome Amplification

Whole genome amplification (WGA) is essential for obtaining genome sequences from single bacterial cells because the quantity of template DNA contained in a single cell is very low. Multiple displacement amplification (MDA), using Phi29 DNA polymerase and random primers, is the most widely used method for single-cell WGA. However, single-cell MDA usually results in uneven genome coverage because of amplification bias, background amplification of contaminating DNA, and formation of chimeras by linking of non-contiguous chromosomal regions. Here, we present a novel MDA method, termed droplet MDA, that minimizes amplification bias and amplification of contaminants by using picoliter-sized droplets for compartmentalized WGA reactions. Extracted DNA fragments from a lysed cell in MDA mixture are divided into 10⁵ droplets (67 pL) within minutes via flow through simple microfluidic channels. Compartmentalized genome fragments can be individually amplified in these droplets without the risk of encounter with reagent-borne or environmental contaminants. Following quality assessment of WGA products from single Escherichia coli cells, we showed that droplet MDA minimized unexpected amplification and improved the percentage of genome recovery from 59% to 89%. Our results demonstrate that microfluidic-generated droplets show potential as an efficient tool for effective amplification of low-input DNA for single-cell genomics and greatly reduce the cost and labor investment required for determination of nearly complete genome sequences of uncultured bacteria from environmental samples.

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.

Abundant Atribacteria in deep marine sediment from the Adélie Basin, Antarctica

Bacteria belonging to the newly classified candidate phylum “Atribacteria” (formerly referred to as “OP9” and “JS1”) are common in anoxic methane-rich sediments. However, the metabolic functions and biogeochemical role of these microorganisms in the subsurface remains unrealized due to the lack of pure culture representatives. In this study of deep sediment from Antarctica’s Adélie Basin, collected during Expedition 318 of the Integrated Ocean Drilling Program (IODP), Atribacteria-related sequences of the 16S rRNA gene were abundant (up to 51% of the sequences) and steadily increased in relative abundance with depth throughout the methane-rich zones. To better understand the metabolic potential of Atribacteria within this environment, and to compare with phylogenetically distinct Atribacteria from non-deep-sea environments, individual cells were sorted for single cell genomics from sediment collected from 97.41 m below the seafloor from IODP Hole U1357C. As observed for non-marine Atribacteria, a partial single cell genome suggests a heterotrophic metabolism, with Atribacteria potentially producing fermentation products such as acetate, ethanol, and CO₂. These products may in turn support methanogens within the sediment microbial community and explain the frequent occurrence of Atribacteria in anoxic methane-rich sediments. This first report of a single cell genome from deep sediment broadens the known diversity within the Atribacteria phylum and highlights the potential role of Atribacteria in carbon cycling in deep sediment.

DEEP BIOSPHERE. Exploring deep microbial life in coal-bearing sediment down to ~2.5 km below the ocean floor

Microbial life inhabits deeply buried marine sediments, but the extent of this vast ecosystem remains poorly constrained. Here we provide evidence for the existence of microbial communities in ~40° to 60°C sediment associated with lignite coal beds at ~1.5 to 2.5 km below the seafloor in the Pacific Ocean off Japan. Microbial methanogenesis was indicated by the isotopic compositions of methane and carbon dioxide, biomarkers, cultivation data, and gas compositions. Concentrations of indigenous microbial cells below 1.5 km ranged from <10 to ~10(4) cells cm(-3). Peak concentrations occurred in lignite layers, where communities differed markedly from shallower subseafloor communities and instead resembled organotrophic communities in forest soils. This suggests that terrigenous sediments retain indigenous community members tens of millions of years after burial in the seabed.

Uncultivated thermophiles: current status and spotlight on 'Aigarchaeota'

Meta-analysis of cultivation-independent sequence data shows that geothermal systems host an abundance of novel organisms, representing a vast unexplored phylogenetic and functional diversity among yet-uncultivated thermophiles. A number of thermophiles have recently been interrogated using metagenomic and/or single-cell genomic approaches, including members of taxonomic groups that inhabit both thermal and non-thermal environments, such as 'Acetothermia' (OP1) and 'Atribacteria' (OP9/JS1), as well as the exclusively thermophilic lineages 'Korarchaeota', 'Calescamantes' (EM19), 'Fervidibacteria' (OctSpA1-106), and 'Aigarchaeota' (HWCG-I). The 'Aigarchaeota', a sister lineage to the Thaumarchaeota, likely includes both hyperthermophiles and moderate thermophiles. They inhabit terrestrial, marine, and subsurface thermal environments and comprise at least nine genus-level lineages, several of which are globally distributed.

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.

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.

Manganese and iron as structuring parameters of microbial communities in Arctic marine sediments from the Baffin Bay

The Arctic Baffin Bay between Canada and Greenland is sea ice-covered during the majority of the year, restricting primary production to the summer months. Sediments receive low amounts of mostly terrestrial- and less marine-derived organic matter. To study microbial communities constrained by physicochemical conditions changing with distance from land and ocean depth, we applied high-throughput 16S rRNA gene sequencing and compared sequence diversity with biogeochemical parameters in 40 different sediment samples. Samples originated from seven cores down to 470 cm below seafloor along a shelf-to-basin transect. Bacterial diversity decreased faster with depth in basin than in shelf sediments, suggesting higher organic matter content sustained diversity into greater depths. All samples were dominated by Betaproteobacteria (mostly order Burkholderiales), which were especially abundant in basin sediments with low organic carbon and high Mn and Fe pore water concentrations. Strong statistical correlations between concentrations of reduced Mn and/or Fe and the relative abundances of Betaproteobacteria suggest that this group is involved in metal reduction in Baffin Bay sediments. Dehalococcoidia (phylum Chloroflexi) were abundant in all samples, especially in shelf sediments with high organic content. This study indicates that Mn and/or Fe play important roles structuring microbial communities in Arctic sediments poor in organic matter.

Distinct microbial populations are tightly linked to the profile of dissolved iron in the methanic sediments of the Helgoland mud area, North Sea

Iron reduction in subseafloor sulfate-depleted and methane-rich marine sediments is currently a subject of interest in subsurface geomicrobiology. While iron reduction and microorganisms involved have been well studied in marine surface sediments, little is known about microorganisms responsible for iron reduction in deep methanic sediments. Here, we used quantitative PCR-based 16S rRNA gene copy numbers and pyrosequencing-based relative abundances of bacteria and archaea to investigate covariance between distinct microbial populations and specific geochemical profiles in the top 5 m of sediment cores from the Helgoland mud area, North Sea. We found that gene copy numbers of bacteria and archaea were specifically higher around the peak of dissolved iron in the methanic zone (250–350 cm). The higher copy numbers at these depths were also reflected by the relative sequence abundances of members of the candidate division JS1, methanogenic and Methanohalobium/ANME-3 related archaea. The distribution of these populations was strongly correlated to the profile of pore-water Fe²⁺ while that of Desulfobacteraceae corresponded to the pore-water sulfate profile. Furthermore, specific JS1 populations also strongly co-varied with the distribution of Methanosaetaceae in the methanic zone. Our data suggest that the interplay among JS1 bacteria, methanogenic archaea and Methanohalobium/ANME-3-related archaea may be important for iron reduction and methane cycling in deep methanic sediments of the Helgoland mud area and perhaps in other methane-rich depositional environments.

Thermophiles in the genomic era: Biodiversity, science, and applications

Thermophiles and hyperthermophiles are present in various regions of the Earth, including volcanic environments, hot springs, mud pots, fumaroles, geysers, coastal thermal springs, and even deep-sea hydrothermal vents. They are also found in man-made environments, such as heated compost facilities, reactors, and spray dryers. Thermophiles, hyperthermophiles, and their bioproducts facilitate various industrial, agricultural, and medicinal applications and offer potential solutions to environmental damages and the demand for biofuels. Intensified efforts to sequence the entire genome of hyperthermophiles and thermophiles are increasing rapidly, as evidenced by the fact that over 120 complete genome sequences of the hyperthermophiles Aquificae, Thermotogae, Crenarchaeota, and Euryarchaeota are now available. In this review, we summarise the major current applications of thermophiles and thermozymes. In addition, emphasis is placed on recent progress in understanding the biodiversity, genomes, transcriptomes, metagenomes, and single-cell sequencing of thermophiles in the genomic era.

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.

Low frequency of endospore-specific genes in subseafloor sedimentary metagenomes

Spore formation is considered to be one of the microbial strategies for long-term survival in subseafloor sedimentary habitats. However, our knowledge of the genetic and physiological characteristics of subseafloor microbes is limited. Here, we studied the distribution and frequency of genes that are related to endospore formation in 10 subseafloor sedimentary metagenomes from Site C9001 off Japan and Site 1229 off Peru. None or very low frequencies of endospore-specific genes (e.g. dpaA, dpaB, sspA, spo0A, spoIIGA, spoIIM, spoIIIAB, spoIVA, spoIVB, yabP, yunB, spoVM) were observed in the subseafloor metagenomes. Based on the number of universally conserved single copy genes, the frequency ratio of putative endospore-formers was estimated to be < 10%, which is consistent with the frequency of Clostridia-derived genomes (2-4%) but is lower than previous estimates based on the concentration of dipicolinic acid. Conceivable explanations for this discrepancy are as follows: the efficiency of lysis and DNA extraction of subseafloor endospore cells may have been lower than those of vegetative cells, conversion factor of dipicolinic acid content per cell may differ, and/or sporulation-related genes and other functional strategies for long-term survival in the deep subseafloor biosphere are evolutionarily distinct from known spore-forming gene repertoires.

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.

Metagenomic profiling of known and unknown microbes with microbeGPS

Microbial community profiling identifies and quantifies organisms in metagenomic sequencing data using either reference based or unsupervised approaches. However, current reference based profiling methods only report the presence and abundance of single reference genomes that are available in databases. Since only a small fraction of environmental genomes is represented in genomic databases, these approaches entail the risk of false identifications and often suggest a higher precision than justified by the data. Therefore, we developed MicrobeGPS, a novel metagenomic profiling approach that overcomes these limitations. MicrobeGPS is the first method that identifies microbiota in the sample and estimates their genomic distances to known reference genomes. With this strategy, MicrobeGPS identifies organisms down to the strain level and highlights possibly inaccurate identifications when the correct reference genome is missing. We demonstrate on three metagenomic datasets with different origin that our approach successfully avoids misleading interpretation of results and additionally provides more accurate results than current profiling methods. Our results indicate that MicrobeGPS can enable reference based taxonomic profiling of complex and less characterized microbial communities. MicrobeGPS is open source and available from https://sourceforge.net/projects/microbegps/ as source code and binary distribution for Windows and Linux operating systems.

Reconstructing each cell's genome within complex microbial communities-dream or reality?

As the vast majority of microorganisms have yet to be cultivated in a laboratory setting, access to their genetic makeup has largely been limited to cultivation-independent methods. These methods, namely metagenomics and more recently single-cell genomics, have become cornerstones for microbial ecology and environmental microbiology. One ultimate goal is the recovery of genome sequences from each cell within an environment to move toward a better understanding of community metabolic potential and to provide substrate for experimental work. As single-cell sequencing has the ability to decipher all sequence information contained in an individual cell, this method holds great promise in tackling such challenge. Methodological limitations and inherent biases however do exist, which will be discussed here based on environmental and benchmark data, to assess how far we are from reaching this goal.

Microbial population dynamics during long-term sludge adaptation of thermophilic and mesophilic sequencing batch digesters treating sewage fine sieved fraction at varying organic loading rates

BACKGROUND

In this research, the feasibility of, and population dynamics in, one-step anaerobic sequencing batch reactor systems treating the fine sieved fraction (FSF) from raw municipal wastewater was studied under thermophilic (55 °C) and mesophilic (35 °C) conditions. FSF was sequestered from raw municipal wastewater, in the Netherlands, using a rotating belt filter (mesh size 350 micron). FSF is a heterogeneous substrate that mainly consists of fibres originating from toilet paper and thus contains a high cellulosic fraction (60-80 % of total solids content), regarded as an energy-rich material.

RESULTS

Results of the 656-day fed-batch operation clearly showed that thermophilic digestion was more stable, applying high organic loading rates (OLR) up to 22 kg COD/(m(3) day). In contrast, the mesophilic digester already failed applying an OLR of 5.5 kg COD/(m(3) day), indicated by a drop in pH and increase in volatile fatty acids (VFAs). The observed viscosity values of the mesophilic sludge were more than tenfold higher than the thermophilic sludge. 454-pyrosequencing of eight mesophilic and eight thermophilic biomass samples revealed that Bacteroides and aceticlastic methanogen Methanosaeta were the dominant genera in the mesophilic digester, whereas OP9 lineages, Clostridium and the hydrogenotrophic methanogen Methanothermobacter dominated the thermophilic one.

CONCLUSIONS

Our study suggests that applying thermophilic conditions for FSF digestion would result in a higher biogas production rate and/or a smaller required reactor volume, comparing to mesophilic conditions.

Mapping the inner workings of the microbiome: genomic- and metagenomic-based study of metabolism and metabolic interactions in the human microbiome

The human gut microbiome is a major contributor to human metabolism and health, yet the metabolic processes that are carried out by various community members, the way these members interact with each other and with the host, and the impact of such interactions on the overall metabolic machinery of the microbiome have not yet been mapped. Here, we discuss recent efforts to study the metabolic inner workings of this complex ecosystem. We will specifically highlight two interrelated lines of work, the first aiming to deconvolve the microbiome and to characterize the metabolic capacity of various microbiome species and the second aiming to utilize computational modeling to infer and study metabolic interactions between these species.

Improved assemblies using a source-agnostic pipeline for MetaGenomic Assembly by Merging (MeGAMerge) of contigs

Assembly of metagenomic samples is a very complex process, with algorithms designed to address sequencing platform-specific issues, (read length, data volume, and/or community complexity), while also faced with genomes that differ greatly in nucleotide compositional biases and in abundance. To address these issues, we have developed a post-assembly process: MetaGenomic Assembly by Merging (MeGAMerge). We compare this process to the performance of several assemblers, using both real, and in-silico generated samples of different community composition and complexity. MeGAMerge consistently outperforms individual assembly methods, producing larger contigs with an increased number of predicted genes, without replication of data. MeGAMerge contigs are supported by read mapping and contig alignment data, when using synthetically-derived and real metagenomic data, as well as by gene prediction analyses and similarity searches. MeGAMerge is a flexible method that generates improved metagenome assemblies, with the ability to accommodate upcoming sequencing platforms, as well as present and future assembly algorithms.

Recent advances in genomic DNA sequencing of microbial species from single cells

The vast majority of microbial species remain uncultivated and, until recently, about half of all known bacterial phyla were identified only from their 16S ribosomal RNA gene sequence. With the advent of single-cell sequencing, genomes of uncultivated species are rapidly filling in unsequenced branches of the microbial phylogenetic tree. The wealth of new insights gained from these previously inaccessible groups is providing a deeper understanding of their basic biology, taxonomy and evolution, as well as their diverse roles in environmental ecosystems and human health.

Illuminating microbial dark matter in meromictic Sakinaw Lake

Despite recent advances in metagenomic and single-cell genomic sequencing to investigate uncultivated microbial diversity and metabolic potential, fundamental questions related to population structure, interactions, and biogeochemical roles of candidate divisions remain. Numerous molecular surveys suggest that stratified ecosystems manifesting anoxic, sulfidic, and/or methane-rich conditions are enriched in these enigmatic microbes. Here we describe diversity, abundance, and cooccurrence patterns of uncultivated microbial communities inhabiting the permanently stratified waters of meromictic Sakinaw Lake, British Columbia, Canada, using 454 sequencing of the small-subunit rRNA gene with three-domain resolution. Operational taxonomic units (OTUs) were affiliated with 64 phyla, including more than 25 candidate divisions. Pronounced trends in community structure were observed for all three domains with eukaryotic sequences vanishing almost completely below the mixolimnion, followed by a rapid and sustained increase in methanogen-affiliated (∼10%) and unassigned (∼60%) archaeal sequences as well as bacterial OTUs affiliated with Chloroflexi (∼22%) and candidate divisions (∼28%). Network analysis revealed highly correlated, depth-dependent cooccurrence patterns between Chloroflexi, candidate divisions WWE1, OP9/JS1, OP8, and OD1, methanogens, and unassigned archaeal OTUs indicating niche partitioning and putative syntrophic growth modes. Indeed, pathway reconstruction using recently published Sakinaw Lake single-cell genomes affiliated with OP9/JS1 and OP8 revealed complete coverage of the Wood-Ljungdahl pathway with potential to drive syntrophic acetate oxidation to hydrogen and carbon dioxide under methanogenic conditions. Taken together, these observations point to previously unrecognized syntrophic networks in meromictic lake ecosystems with the potential to inform design and operation of anaerobic methanogenic bioreactors.

A quantitative comparison of single-cell whole genome amplification methods

Single-cell sequencing is emerging as an important tool for studies of genomic heterogeneity. Whole genome amplification (WGA) is a key step in single-cell sequencing workflows and a multitude of methods have been introduced. Here, we compare three state-of-the-art methods on both bulk and single-cell samples of E. coli DNA: Multiple Displacement Amplification (MDA), Multiple Annealing and Looping Based Amplification Cycles (MALBAC), and the PicoPLEX single-cell WGA kit (NEB-WGA). We considered the effects of reaction gain on coverage uniformity, error rates and the level of background contamination. We compared the suitability of the different WGA methods for the detection of copy-number variations, for the detection of single-nucleotide polymorphisms and for de-novo genome assembly. No single method performed best across all criteria and significant differences in characteristics were observed; the choice of which amplifier to use will depend strongly on the details of the type of question being asked in any given experiment.

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.

Temperature and injection water source influence microbial community structure in four Alaskan North Slope hydrocarbon reservoirs

A fundamental knowledge of microbial community structure in petroleum reservoirs can improve predictive modeling of these environments. We used hydrocarbon profiles, stable isotopes, and high-density DNA microarray analysis to characterize microbial communities in produced water from four Alaskan North Slope hydrocarbon reservoirs. Produced fluids from Schrader Bluff (24–27°C), Kuparuk (47–70°C), Sag River (80°C), and Ivishak (80–83°C) reservoirs were collected, with paired soured/non-soured wells sampled from Kuparuk and Ivishak. Chemical and stable isotope data suggested Schrader Bluff had substantial biogenic methane, whereas methane was mostly thermogenic in deeper reservoirs. Acetoclastic methanogens (Methanosaeta) were most prominent in Schrader Bluff samples, and the combined δD and δ¹³C values of methane also indicated acetoclastic methanogenesis could be a primary route for biogenic methane. Conversely, hydrogenotrophic methanogens (e.g., Methanobacteriaceae) and sulfide-producing Archaeoglobus and Thermococcus were more prominent in Kuparuk samples. Sulfide-producing microbes were detected in all reservoirs, uncoupled from souring status (e.g., the non-soured Kuparuk samples had higher relative abundances of many sulfate-reducers compared to the soured sample, suggesting sulfate-reducers may be living fermentatively/syntrophically when sulfate is limited). Sulfate abundance via long-term seawater injection resulted in greater relative abundances of Desulfonauticus, Desulfomicrobium, and Desulfuromonas in the soured Ivishak well compared to the non-soured well. In the non-soured Ivishak sample, several taxa affiliated with Thermoanaerobacter and Halomonas predominated. Archaea were not detected in the deepest reservoirs. Functional group taxa differed in relative abundance among reservoirs, likely reflecting differing thermal and/or geochemical influences.

Patterns of Microbially Driven Carbon Cycling in the Ocean: Links between Extracellular Enzymes and Microbial Communities

Heterotrophic microbial communities play a central role in the marine carbon cycle. They are active in nearly all known environments, from the surface to the deep ocean, in the sediments, and from the equator to the Poles. In order to process complex organic matter, these communities produce extracellular enzymes of the correct structural specificity to hydrolyze substrates to sizes sufficiently small for uptake. Extracellular enzymatic hydrolysis thus initiates heterotrophic carbon cycling. Our knowledge of the enzymatic capabilities of microbial communities in the ocean is still underdeveloped. Recent studies, however, suggest that there may be large-scale patterns of enzymatic function in the ocean, patterns of community function that may be connected to emerging patterns of microbial community composition. Here I review some of these large-scale contrasts in microbial enzyme activities, between high-latitude and temperate surface ocean waters, contrasts between inshore and offshore waters, changes with depth gradients in the ocean, and contrasts between the water column and underlying sediments. These contrasting patterns are set in the context of recent studies of microbial communities and patterns of microbial biogeography. Focusing on microbial community function as well as composition and potential should yield clearer understanding of the factors driving carbon cycling in the ocean.

Obtaining genomes from uncultivated environmental microorganisms using FACS-based single-cell genomics

Single-cell genomics is a powerful tool for exploring the genetic makeup of environmental microorganisms, the vast majority of which are difficult, if not impossible, to cultivate with current approaches. Here we present a comprehensive protocol for obtaining genomes from uncultivated environmental microbes via high-throughput single-cell isolation by FACS. The protocol encompasses the preservation and pretreatment of differing environmental samples, followed by the physical separation, lysis, whole-genome amplification and 16S rRNA-based identification of individual bacterial and archaeal cells. The described procedure can be performed with standard molecular biology equipment and a FACS machine. It takes <12 h of bench time over a 4-d time period, and it generates up to 1 μg of genomic DNA from an individual microbial cell, which is suitable for downstream applications such as PCR amplification and shotgun sequencing. The completeness of the recovered genomes varies, with an average of ∼50%.

Current challenges in the bioinformatics of single cell genomics

Single cell genomics is a rapidly growing field with many new techniques emerging in the past few years. However, few bioinformatics tools specific for single cell genomics analysis are available. Single cell DNA/RNA sequencing data usually have low genome coverage and high amplification bias, which makes bioinformatics analysis challenging. Many current bioinformatics tools developed for bulk cell sequencing do not work well with single cell sequencing data. Here, we summarize current challenges in the bioinformatics analysis of single cell genomic DNA sequencing and single cell transcriptomes. These challenges include calling copy number variations, identifying mutated genes in tumor samples, reconstructing cell lineages, recovering low abundant transcripts, and improving the accuracy of quantitative analysis of transcripts. Development in single cell genomics bioinformatics analysis will promote the application of this technology to basic biology and medical research.

Kallotenue papyrolyticum gen. nov., sp. nov., a cellulolytic and filamentous thermophile that represents a novel lineage (Kallotenuales ord. nov., Kallotenuaceae fam. nov.) within the class Chloroflexia

Several closely related, thermophilic and cellulolytic bacterial strains, designated JKG1T, JKG2, JKG3, JKG4 and JKG5, were isolated from a cellulolytic enrichment (corn stover) incubated in the water column of Great Boiling Spring, NV. Strain JKG1T had cells of diameter 0.7–0.9 µm and length ~2.0 µm that formed non-branched, multicellular filaments reaching >300 µm. Spores were not formed and dense liquid cultures were red. The temperature range for growth was 45–65 °C, with an optimum of 55 °C. The pH range for growth was pH 5.6–9.0, with an optimum of pH 7.5. JKG1T grew as an aerobic heterotroph, utilizing glucose, sucrose, xylose, arabinose, cellobiose, CM-cellulose, filter paper, microcrystalline cellulose, xylan, starch, Casamino acids, tryptone, peptone, yeast extract, acetate, citrate, lactate, pyruvate and glycerol as sole carbon sources, and was not observed to photosynthesize. The cells stained Gram-negative. Phylogenetic analysis using 16S rRNA gene sequences placed the new isolates in the class Chloroflexia , but distant from other cultivated members, with the highest sequence identity of 82.5 % to Roseiflexus castenholzii . The major quinone was menaquinone-9; no ubiquinones were detected. The major cellular fatty acids (>5 %) were C18 : 0, anteiso-C17 : 0, iso-C18 : 0, iso-C17 : 0, C16 : 0, iso-C16 : 0 and C17 : 0. The peptidoglycan amino acids were alanine, ornithine, glutamic acid, serine and asparagine. Whole-cell sugars included mannose, rhamnose, glucose, galactose, ribose, arabinose and xylose. Morphological, phylogenetic and chemotaxonomic results suggest that JKG1T is representative of a new lineage within the class Chloroflexia , which we propose to designate Kallotenue papyrolyticum gen. nov., sp. nov., Kallotenuaceae fam. nov., Kallotenuales ord. nov. The type strain of Kallotenue papyrolyticum gen. nov., sp. nov. is JKG1T ( = DSM 26889T = JCM 19132T).

Testing the reproducibility of multiple displacement amplification on genomes of clonal endosymbiont populations

The multiple displacement amplification method has revolutionized genomic studies of uncultured bacteria, where the extraction of pure DNA in sufficient quantity for next-generation sequencing is challenging. However, the method is problematic in that it amplifies the target DNA unevenly, induces the formation of chimeric reads and also amplifies contaminating DNA. Here, we have tested the reproducibility of the multiple displacement amplification method using serial dilutions of extracted genomic DNA and intact cells from the cultured endosymbiont Bartonella australis. The amplified DNA was sequenced with the Illumina sequencing technology, and the results were compared to sequence data obtained from unamplified DNA in this study as well as from a previously published genome project. We show that artifacts such as the extent of the amplification bias, the percentage of chimeric reads and the relative fraction of contaminating DNA increase dramatically for the smallest amounts of template DNA. The pattern of read coverage was reproducibly obtained for samples with higher amounts of template DNA, suggesting that the bias is non-random and genome-specific. A re-analysis of previously published sequence data obtained after amplification from clonal endosymbiont populations confirmed these predictions. We conclude that many of the artifacts associated with the use of the multiple displacement amplification method can be alleviated or much reduced by using multiple cells as the template for the amplification. These findings should be particularly useful for researchers studying the genomes of endosymbionts and other uncultured bacteria, for which a small clonal population of cells can be isolated.

Is metagenomics resolving identification of functions in microbial communities?

We are coming up on the tenth anniversary of the broad use of the method involving whole metagenome shotgun sequencing, referred to as metagenomics. The application of this approach has definitely revolutionized microbiology and the related fields, including the realization of the importance of the human microbiome. As such, metagenomics has already provided a novel outlook on the complexity and dynamics of microbial communities that are an important part of the biosphere of the planet. Accumulation of massive amounts of sequence data also caused a surge in the development of bioinformatics tools specially designed to provide pipelines for data analysis and visualization. However, a critical outlook into the field is required to appreciate what could be and what has currently been gained from the massive sequence databases that are being generated with ever-increasing speed.

Wide distribution of autochthonous branched glycerol dialkyl glycerol tetraethers (bGDGTs) in U.S. Great Basin hot springs

Branched glycerol dialkyl glycerol tetraethers (bGDGTs) are membrane-spanning lipids that likely stabilize membranes of some bacteria. Although bGDGTs have been reported previously in certain geothermal environments, it has been suggested that they may derive from surrounding soils since bGDGTs are known to be produced by soil bacteria. To test the hypothesis that bGDGTs can be produced by thermophiles in geothermal environments, we examined the distribution and abundance of bGDGTs, along with extensive geochemical data, in 40 sediment and mat samples collected from geothermal systems in the U.S. Great Basin (temperature: 31–95°C; pH: 6.8–10.7). bGDGTs were found in 38 out of 40 samples at concentrations up to 824 ng/g sample dry mass and comprised up to 99.5% of total GDGTs (branched plus isoprenoidal). The wide distribution of bGDGTs in hot springs, strong correlation between core and polar lipid abundances, distinctness of bGDGT profiles compared to nearby soils, and higher concentration of bGDGTs in hot springs compared to nearby soils provided evidence of in situ production, particularly for the minimally methylated bGDGTs I, Ib, and Ic. Polar bGDGTs were found almost exclusively in samples ≤70°C and the absolute abundance of polar bGDGTs correlated negatively with properties of chemically reduced, high temperature spring sources (temperature, H₂S/HS⁻) and positively with properties of oxygenated, low temperature sites (O₂, NO⁻₃). Two-way cluster analysis and nonmetric multidimensional scaling based on relative abundance of polar bGDGTs supported these relationships and showed a negative relationship between the degree of methylation and temperature, suggesting a higher abundance for minimally methylated bGDGTs at high temperature. This study presents evidence of the widespread production of bGDGTs in mats and sediments of natural geothermal springs in the U.S. Great Basin, especially in oxygenated, low-temperature sites (≤70°C).

In situ production of branched glycerol dialkyl glycerol tetraethers in a great basin hot spring (USA)

Branched glycerol dialkyl glycerol tetraethers (bGDGTs) are predominantly found in soils and peat bogs. In this study, we analyzed core (C)-bGDGTs after hydrolysis of polar fractions using liquid chromatography-atmospheric pressure chemical ionization-mass spectrometry and analyzed intact P-bGDGTs using total lipid extract (TLE) without hydrolysis by liquid chromatography-electrospray ionization-multiple stage mass spectrometry. Our results show multiple lines of evidence for the production of bGDGTs in sediments and cellulolytic enrichments in a hot spring (62–86°C) in the Great Basin (USA). First, in situ cellulolytic enrichment led to an increase in the relative abundance of hydrolysis-derived P-bGDGTs over their C-bGDGT counterparts. Second, the hydrolysis-derived P- and C-bGDGT profiles in the hot spring were different from those of the surrounding soil samples; in particular, a monoglycosidic bGDGT Ib containing 13,16-dimethyloctacosane and one cyclopentane moiety was detected in the TLE but it was undetectable in surrounding soil samples even after sample enrichments. Third, previously published 16S rRNA gene pyrotag analysis from the same lignocellulose samples demonstrated the enrichment of thermophiles, rather than mesophiles, and total bGDGT abundance in cellulolytic enrichments correlated with the relative abundance of 16S rRNA gene pyrotags from thermophilic bacteria in the phyla Bacteroidetes, Dictyoglomi, EM3, and OP9 (“Atribacteria”). These observations conclusively demonstrate the production of bGDGTs in this hot spring; however, the identity of organisms that produce bGDGTs in the geothermal environment remains unclear.

Optofluidic cell selection from complex microbial communities for single-genome analysis

Genetic analysis of single cells is emerging as a powerful approach for studies of heterogeneous cell populations. Indeed, the notion of homogeneous cell populations is receding as approaches to resolve genetic and phenotypic variation between single cells are applied throughout the life sciences. A key step in single-cell genomic analysis today is the physical isolation of individual cells from heterogeneous populations, particularly microbial populations, which often exhibit high diversity. Here, we detail the construction and use of instrumentation for optical trapping inside microfluidic devices to select individual cells for analysis by methods including nucleic acid sequencing. This approach has unique advantages for analyses of rare community members, cells with irregular morphologies, small quantity samples, and studies that employ advanced optical microscopy.