Relationship between Microorganisms Inhabiting Alkaline Siliceous Hot Spring Mat Communities and Overflowing Water.

The compositions of Octopus Spring and Mushroom Spring (Yellowstone National Park, Wyoming, USA) microbial mats have been thoroughly studied, but the compositions of the effluent waters that flow above the mats have not. In this study, cells in the mats and overflowing waters of both springs were investigated at multiple sites where Synechococcus spp. are the dominant cyanobacteria (ca. 72°C to ca. 50°C), and on several dates. In addition to microscopic analyses of stained and autofluorescent cells, 16S rRNA gene sequencing was used to characterize the major taxa present and a protein-encoding gene (psaA) was sequenced and analyzed by ecotype simulation to predict species of Synechococcus The mats of both springs were similar in terms of the downstream distribution of predominant taxa detected previously. However, waters above these mats were predominated by taxa that reside in upstream mats or communities above the upper-temperature limit of the mat. A disturbance/recolonization study was performed at a site normally predominated by Synechococcus species adapted to low temperatures. After removing indigenous Synechococcus cells, Synechococcus species adapted to higher temperatures, which were predominant in the water overflowing this site, colonized the newly forming mat. Differences in recolonization under reduced and UV-screened irradiance suggested that, in addition to physical transport, environmental conditions likely select for species that are better adapted to these different conditions and can influence mat recovery. A transport model was developed and used to predict that, in Mushroom Spring, erosion predominates in the narrower and deeper upstream effluents and deposition predominates over erosion in wider and shallower downstream effluents.IMPORTANCE In flowing aquatic systems, cell erosion and deposition are important to the dispersal of cells from one location to another. Very little is known about microbial dispersal and the physical processes that underlie it. This study demonstrates its importance to colonization of downstream surfaces and especially to the recolonization and functioning of disturbed sites. Ecological systems in flowing environments are often, roughly speaking, pseudosteady, in that nutrients enter the system and by-products leave at relatively steady rates. Over time, material inputs and outputs must balance. Measurements of input fluxes (e.g., growth rates and proxies, such as photosynthesis rates) are frequent. However, erosion and deposition of cells are seldom measured and ecological significance is sometimes neglected. The importance of these parameters is immediately evident in any attempt to construct a model of long-time community behavior, as spatial ecological structure is significantly impacted and can be dominated by migration of organisms, even in small numbers.

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.

Geographic and Ecological Diversity of Green Sulfur Bacteria in Hot Spring Mat Communities.

Three strains of thermophilic green sulfur bacteria (GSB) are known; all are from microbial mats in hot springs in Rotorua, New Zealand (NZ) and belong to the species Chlorobaculum tepidum. Here, we describe diverse populations of GSB inhabiting Travel Lodge Spring (TLS) (NZ) and hot springs ranging from 36.1 °C to 51.1 °C in the Republic of the Philippines (PHL) and Yellowstone National Park (YNP), Wyoming, USA. Using targeted amplification and restriction fragment length polymorphism analysis, GSB 16S rRNA sequences were detected in mats in TLS, one PHL site, and three regions of YNP. GSB enrichments from YNP and PHL mats contained small, green, nonmotile rods possessing chlorosomes, chlorobactene, and bacteriochlorophyll c. Partial 16S rRNA gene sequences from YNP, NZ, and PHL mats and enrichments from YNP and PHL samples formed distinct phylogenetic clades, suggesting geographic isolation, and were associated with samples differing in temperature and pH, suggesting adaptations to these parameters. Sequences from enrichments and corresponding mats formed clades that were sometimes distinct, increasing the diversity detected. Sequence differences, monophyly, distribution patterns, and evolutionary simulation modeling support our discovery of at least four new putative moderately thermophilic Chlorobaculum species that grew rapidly at 40 °C to 44 °C.

Hyperactive nanobacteria with host-dependent traits pervade Omnitrophota.

Candidate bacterial phylum Omnitrophota has not been isolated and is poorly understood. We analysed 72 newly sequenced and 349 existing Omnitrophota genomes representing 6 classes and 276 species, along with Earth Microbiome Project data to evaluate habitat, metabolic traits and lifestyles. We applied fluorescence-activated cell sorting and differential size filtration, and showed that most Omnitrophota are ultra-small (~0.2 µm) cells that are found in water, sediments and soils. Omnitrophota genomes in 6 classes are reduced, but maintain major biosynthetic and energy conservation pathways, including acetogenesis (with or without the Wood-Ljungdahl pathway) and diverse respirations. At least 64% of Omnitrophota genomes encode gene clusters typical of bacterial symbionts, suggesting host-associated lifestyles. We repurposed quantitative stable-isotope probing data from soils dominated by andesite, basalt or granite weathering and identified 3 families with high isotope uptake consistent with obligate bacterial predators. We propose that most Omnitrophota inhabit various ecosystems as predators or parasites.

Fine-scale distribution patterns of Synechococcus ecological diversity in microbial mats of Mushroom Spring, Yellowstone National Park.

Past analyses of sequence diversity in high-resolution protein-encoding genes have identified putative ecological species of unicellular cyanobacteria in the genus Synechococcus, which are specialized to 60°C but not 65°C in Mushroom Spring microbial mats. Because these studies were limited to only two habitats, we studied the distribution of Synechococcus sequence variants at 1°C intervals along the effluent flow channel and at 80-µm vertical-depth intervals throughout the upper photic layer of the microbial mat. Diversity at the psaA locus, which encodes a photosynthetic reaction center protein (PsaA), was sampled by PCR amplification, cloning, and sequencing methods at 60, 63, and 65°C sites. The evolutionary simulation programs Ecotype Simulation and AdaptML were used to identify putative ecologically distinct populations (ecotypes). Ecotype Simulation predicted a higher number of putative ecotypes in cases where habitat variation was limited, while AdaptML predicted a higher number of ecologically distinct phylogenetic clades in cases where habitat variation was high. Denaturing gradient gel electrophoresis was used to track the distribution of dominant sequence variants of ecotype populations relative to temperature variation and to O2, pH, and spectral irradiance variation, as measured using microsensors. Different distributions along effluent channel flow and vertical gradients, where temperature, light, and O2 concentrations are known to vary, confirmed the ecological distinctness of putative ecotypes.

Synthase-Selective Exploration of a Tunicate Microbiome by Activity-Guided Single-Cell Genomics.

While thousands of environmental metagenomes have been mined for the presence of novel biosynthetic gene clusters, such computational predictions do not provide evidence of their in vivo biosynthetic functionality. Using fluorescent in situ enzyme assay targeting carrier proteins common to polyketide (PKS) and nonribosomal peptide synthetases (NRPS), we applied fluorescence-activated cell sorting to tunicate microbiome to enrich for microbes with active secondary metabolic capabilities. Single-cell genomics uncovered the genetic basis for a wide biosynthetic diversity in the enzyme-active cells and revealed a member of marine Oceanospirillales harboring a novel NRPS gene cluster with high similarity to phylogenetically distant marine and terrestrial bacteria. Interestingly, this synthase belongs to a larger class of siderophore biosynthetic gene clusters commonly associated with pestilence and disease. This demonstrates activity-guided single-cell genomics as a tool to guide novel biosynthetic discovery.

Evolutionary stasis of a deep subsurface microbial lineage.

Sulfate-reducing bacteria Candidatus Desulforudis audaxviator (CDA) were originally discovered in deep fracture fluids accessed via South African gold mines and have since been found in geographically widespread deep subsurface locations. In order to constrain models for subsurface microbial evolution, we compared CDA genomes from Africa, North America and Eurasia using single cell genomics. Unexpectedly, 126 partial single amplified genomes from the three continents, a complete genome from of an isolate from Eurasia, and metagenome-assembled genomes from Africa and Eurasia shared >99.2% average nucleotide identity, low frequency of SNP's, and near-perfectly conserved prophages and CRISPRs. Our analyses reject sample cross-contamination, recent natural dispersal, and unusually strong purifying selection as likely explanations for these unexpected results. We therefore conclude that the analyzed CDA populations underwent only minimal evolution since their physical separation, potentially as far back as the breakup of Pangea between 165 and 55 Ma ago. High-fidelity DNA replication and repair mechanisms are the most plausible explanation for the highly conserved genome of CDA. CDA presents a stark contrast to the current model organisms in microbial evolutionary studies, which often develop adaptive traits over far shorter periods of time.

Biogeography of American Northwest Hot Spring A/B'-Lineage Synechococcus Populations.

Previous analyses have shown how diversity among unicellular cyanobacteria inhabiting island-like hot springs is structured relative to physical separation and physiochemical differences among springs, especially at local to regional scales. However, these studies have been limited by the low resolution provided by the molecular markers surveyed. We analyzed large datasets obtained by high-throughput sequencing of a segment of the photosynthesis gene psaA from samples collected in hot springs from geothermal basins in Yellowstone National Park, Montana, and Oregon, all known from previous studies to contain populations of A/B'-lineage Synechococcus. The fraction of identical sequences was greater among springs separated by <50 km than among springs separated by >50 km, and springs separated by >800 km shared sequence variants only rarely. Phylogenetic analyses provided evidence for endemic lineages that could be related to geographic isolation and/or geochemical differences on regional scales. Ecotype Simulation 2 was used to predict putative ecotypes (ecologically distinct populations), and their membership, and canonical correspondence analysis was used to examine the geographical and geochemical bases for variation in their distribution. Across the range of Oregon and Yellowstone, geographical separation explained the largest percentage of the differences in distribution of ecotypes (9.5% correlated to longitude; 9.4% to latitude), with geochemical differences explaining the largest percentage of the remaining differences in distribution (7.4-9.3% correlated to magnesium, sulfate, and sulfide). Among samples within the Greater Yellowstone Ecosystem, geochemical differences significantly explained the distribution of ecotypes (6.5-9.3% correlated to magnesium, boron, sulfate, silicon dioxide, chloride, and pH). Nevertheless, differences in the abundance and membership of ecotypes in Yellowstone springs with similar chemistry suggested that allopatry may be involved even at local scales. Synechococcus populations have diverged both by physical isolation and physiochemical differences, and populations on surprisingly local scales have been evolving independently.

Hiding in Plain Sight: The Globally Distributed Bacterial Candidate Phylum PAUC34f.

Bacterial candidate phylum PAUC34f was originally discovered in marine sponges and is widely considered to be composed of sponge symbionts. Here, we report 21 single amplified genomes (SAGs) of PAUC34f from a variety of environments, including the dark ocean, lake sediments, and a terrestrial aquifer. The diverse origins of the SAGs and the results of metagenome fragment recruitment suggest that some PAUC34f lineages represent relatively abundant, free-living cells in environments other than sponge microbiomes, including the deep ocean. Both phylogenetic and biogeographic patterns, as well as genome content analyses suggest that PAUC34f associations with hosts evolved independently multiple times, while free-living lineages of PAUC34f are distinct and relatively abundant in a wide range of environments.

Ancestral Absence of Electron Transport Chains in Patescibacteria and DPANN.

Recent discoveries suggest that the candidate superphyla Patescibacteria and DPANN constitute a large fraction of the phylogenetic diversity of Bacteria and Archaea. Their small genomes and limited coding potential have been hypothesized to be ancestral adaptations to obligate symbiotic lifestyles. To test this hypothesis, we performed cell-cell association, genomic, and phylogenetic analyses on 4,829 individual cells of Bacteria and Archaea from 46 globally distributed surface and subsurface field samples. This confirmed the ubiquity and abundance of Patescibacteria and DPANN in subsurface environments, the small size of their genomes and cells, and the divergence of their gene content from other Bacteria and Archaea. Our analyses suggest that most Patescibacteria and DPANN in the studied subsurface environments do not form specific physical associations with other microorganisms. These data also suggest that their unusual genomic features and prevalent auxotrophies may be a result of ancestral, minimal cellular energy transduction mechanisms that lack respiration, thus relying solely on fermentation for energy conservation.

Four Draft Single-Cell Genome Sequences of Novel, Nearly Identical Kiritimatiellaeota Strains Isolated from the Continental Deep Subsurface.

The recently proposed bacterial phylum Kiritimatiellaeota represents a globally distributed monophyletic clade distinct from other members of the Planctomycetes, Verrucomicrobia, and Chlamydiae (PVC) superphylum. Here, we present four phylogenetically distinct single-cell genome sequences from within the Kiritimatiellaeota lineage sampled from deep continental subsurface aquifer fluids of the Death Valley Regional Flow System in the United States.

Hydrogen-based metabolism as an ancestral trait in lineages sibling to the Cyanobacteria.

The evolution of aerobic respiration was likely linked to the origins of oxygenic Cyanobacteria. Close phylogenetic neighbors to Cyanobacteria, such as Margulisbacteria (RBX-1 and ZB3), Saganbacteria (WOR-1), Melainabacteria and Sericytochromatia, may constrain the metabolic platform in which aerobic respiration arose. Here, we analyze genomic sequences and predict that sediment-associated Margulisbacteria have a fermentation-based metabolism featuring a variety of hydrogenases, a streamlined nitrogenase, and electron bifurcating complexes involved in cycling of reducing equivalents. The genomes of ocean-associated Margulisbacteria encode an electron transport chain that may support aerobic growth. Some Saganbacteria genomes encode various hydrogenases, and others may be able to use O2 under certain conditions via a putative novel type of heme copper O2 reductase. Similarly, Melainabacteria have diverse energy metabolisms and are capable of fermentation and aerobic or anaerobic respiration. The ancestor of all these groups may have been an anaerobe in which fermentation and H2 metabolism were central metabolic features. The ability to use O2 as a terminal electron acceptor must have been subsequently acquired by these lineages.

Author Correction: Hydrogen-based metabolism as an ancestral trait in lineages sibling to the Cyanobacteria.

The original version of this Article contained errors in Fig. 4. In panel a, the labels 'F420-reducing NiFe hydrogenase (group 3a)' and 'Group 2 NiFe hydrogenase' were misplaced. These errors have been corrected in both the PDF and HTML versions of the Article.

Genomic Comparison of Two Family-Level Groups of the Uncultivated NAG1 Archaeal Lineage from Chemically and Geographically Disparate Hot Springs.

Recent progress based on single-cell genomics and metagenomic investigations of archaea in a variety of extreme environments has led to significant advances in our understanding of the diversity, evolution, and metabolic potential of archaea, yet the vast majority of archaeal diversity remains undersampled. In this work, we coordinated single-cell genomics with metagenomics in order to construct a near-complete genome from a deeply branching uncultivated archaeal lineage sampled from Great Boiling Spring (GBS) in the U.S. Great Basin, Nevada. This taxon is distantly related (distinct families) to an archaeal genome, designated "Novel Archaeal Group 1" (NAG1), which was extracted from a metagenome recovered from an acidic iron spring in Yellowstone National Park (YNP). We compared the metabolic predictions of the NAG1 lineage to better understand how these archaea could inhabit such chemically distinct environments. Similar to the NAG1 population previously studied in YNP, the NAG1 population from GBS is predicted to utilize proteins as a primary carbon source, ferment simple carbon sources, and use oxygen as a terminal electron acceptor under oxic conditions. However, GBS NAG1 populations contained distinct genes involved in central carbon metabolism and electron transfer, including nitrite reductase, which could confer the ability to reduce nitrite under anaerobic conditions. Despite inhabiting chemically distinct environments with large variations in pH, GBS NAG1 populations shared many core genomic and metabolic features with the archaeon identified from YNP, yet were able to carve out a distinct niche at GBS.

Improved genome recovery and integrated cell-size analyses of individual uncultured microbial cells and viral particles.

Microbial single-cell genomics can be used to provide insights into the metabolic potential, interactions, and evolution of uncultured microorganisms. Here we present WGA-X, a method based on multiple displacement amplification of DNA that utilizes a thermostable mutant of the phi29 polymerase. WGA-X enhances genome recovery from individual microbial cells and viral particles while maintaining ease of use and scalability. The greatest improvements are observed when amplifying high G+C content templates, such as those belonging to the predominant bacteria in agricultural soils. By integrating WGA-X with calibrated index-cell sorting and high-throughput genomic sequencing, we are able to analyze genomic sequences and cell sizes of hundreds of individual, uncultured bacteria, archaea, protists, and viral particles, obtained directly from marine and soil samples, in a single experiment. This approach may find diverse applications in microbiology and in biomedical and forensic studies of humans and other multicellular organisms.Single-cell genomics can be used to study uncultured microorganisms. Here, Stepanauskas et al. present a method combining improved multiple displacement amplification and FACS, to obtain genomic sequences and cell size information from uncultivated microbial cells and viral particles in environmental samples.

Author Correction: Improved genome recovery and integrated cell-size analyses of individual uncultured microbial cells and viral particles.

The original version of this Article contained errors in the units of concentration of three reagents listed in the Methods. These errors have all been corrected in both the PDF and HTML versions of the Article.

Rokubacteria: Genomic Giants among the Uncultured Bacterial Phyla.

Recent advances in single-cell genomic and metagenomic techniques have facilitated the discovery of numerous previously unknown, deep branches of the tree of life that lack cultured representatives. Many of these candidate phyla are composed of microorganisms with minimalistic, streamlined genomes lacking some core metabolic pathways, which may contribute to their resistance to growth in pure culture. Here we analyzed single-cell genomes and metagenome bins to show that the "Candidate phylum Rokubacteria," formerly known as SPAM, represents an interesting exception, by having large genomes (6-8 Mbps), high GC content (66-71%), and the potential for a versatile, mixotrophic metabolism. We also observed an unusually high genomic heterogeneity among individual Rokubacteria cells in the studied samples. These features may have contributed to the limited recovery of sequences of this candidate phylum in prior cultivation and metagenomic studies. Our analyses suggest that Rokubacteria are distributed globally in diverse terrestrial ecosystems, including soils, the rhizosphere, volcanic mud, oil wells, aquifers, and the deep subsurface, with no reports from marine environments to date.

Minimum information about a single amplified genome (MISAG) and a metagenome-assembled genome (MIMAG) of bacteria and archaea.

We present two standards developed by the Genomic Standards Consortium (GSC) for reporting bacterial and archaeal genome sequences. Both are extensions of the Minimum Information about Any (x) Sequence (MIxS). The standards are the Minimum Information about a Single Amplified Genome (MISAG) and the Minimum Information about a Metagenome-Assembled Genome (MIMAG), including, but not limited to, assembly quality, and estimates of genome completeness and contamination. These standards can be used in combination with other GSC checklists, including the Minimum Information about a Genome Sequence (MIGS), Minimum Information about a Metagenomic Sequence (MIMS), and Minimum Information about a Marker Gene Sequence (MIMARKS). Community-wide adoption of MISAG and MIMAG will facilitate more robust comparative genomic analyses of bacterial and archaeal diversity.

The molecular dimension of microbial species: 2. Synechococcus strains representative of putative ecotypes inhabiting different depths in the Mushroom Spring microbial mat exhibit different adaptive and acclimative responses to light.

Closely related strains of thermophilic Synechococcus were cultivated from the microbial mats found in the effluent channels of Mushroom Spring, Yellowstone National Park (YNP). These strains have identical or nearly identical 16S rRNA sequences but are representative of separate, predicted putative ecotype (PE) populations, which were identified by using the more highly resolving psaA locus and which predominate at different vertical positions within the 1-mm-thick upper-green layer of the mat. Pyrosequencing confirmed that each strain contained a single, predominant psaA genotype. Strains differed in growth rate as a function of irradiance. A strain with a psaA genotype corresponding to a predicted PE that predominates near the mat surface grew fastest at high irradiances, whereas strains with psaA genotypes representative of predominant subsurface populations grew faster at low irradiance and exhibited greater sensitivity to abrupt shifts to high light. The high-light-adapted and low-light-adapted strains also exhibited differences in pigment content and the composition of the photosynthetic apparatus (photosystem ratio) when grown under different light intensities. Cells representative of the different strains had similar morphologies under low-light conditions, but under high-light conditions, cells of low-light-adapted strains became elongated and formed short chains of cells. Collectively, the results presented here are consistent with the hypothesis that closely related, but distinct, ecological species of Synechococcus occupy different light niches in the Mushroom Spring microbial mat and acclimate differently to changing light environments.

The molecular dimension of microbial species: 1. Ecological distinctions among, and homogeneity within, putative ecotypes of Synechococcus inhabiting the cyanobacterial mat of Mushroom Spring, Yellowstone National Park.

Based on the Stable Ecotype Model, evolution leads to the divergence of ecologically distinct populations (e.g., with different niches and/or behaviors) of ecologically interchangeable membership. In this study, pyrosequencing was used to provide deep sequence coverage of Synechococcus psaA genes and transcripts over a large number of habitat types in the Mushroom Spring microbial mat. Putative ecological species [putative ecotypes (PEs)], which were predicted by an evolutionary simulation based on the Stable Ecotype Model (Ecotype Simulation), exhibited distinct distributions relative to temperature-defined positions in the effluent channel and vertical position in the upper 1 mm-thick mat layer. Importantly, in most cases variants predicted to belong to the same PE formed unique clusters relative to temperature and depth in the mat in canonical correspondence analysis, supporting the hypothesis that while the PEs are ecologically distinct, the members of each ecotype are ecologically homogeneous. PEs responded differently to experimental perturbations of temperature and light, but the genetic variation within each PE was maintained as the relative abundances of PEs changed, further indicating that each population responded as a set of ecologically interchangeable individuals. Compared to PEs that predominate deeper within the mat photic zone, the timing of transcript abundances for selected genes differed for PEs that predominate in microenvironments closer to upper surface of the mat with spatiotemporal differences in light and O2 concentration. All of these findings are consistent with the hypotheses that Synechococcus species in hot spring mats are sets of ecologically interchangeable individuals that are differently adapted, that these adaptations control their distributions, and that the resulting distributions constrain the activities of the species in space and time.