Complete genome of the thermophilic purple sulfur Bacterium Thermochromatium tepidum compared to Allochromatium vinosum and other Chromatiaceae.

The complete genome sequence of the thermophilic purple sulfur bacterium Thermochromatium tepidum strain MCT (DSM 3771T) is described and contrasted with that of its mesophilic relative Allochromatium vinosum strain D (DSM 180T) and other Chromatiaceae. The Tch. tepidum genome is a single circular chromosome of 2,958,290 base pairs with no plasmids and is substantially smaller than the genome of Alc. vinosum. The Tch. tepidum genome encodes two forms of RuBisCO and contains nifHDK and several other genes encoding a molybdenum nitrogenase but lacks a gene encoding a protein that assembles the Fe-S cluster required to form a functional nitrogenase molybdenum-iron cofactor, leaving the phototroph phenotypically Nif-. Tch. tepidum contains genes necessary for oxidizing sulfide to sulfate as photosynthetic electron donor but is genetically unequipped to either oxidize thiosulfate as an electron donor or carry out assimilative sulfate reduction, both of which are physiological hallmarks of Alc. vinosum. Also unlike Alc. vinosum, Tch. tepidum is obligately phototrophic and unable to grow chemotrophically in darkness by respiration. Several genes present in the Alc. vinosum genome that are absent from the genome of Tch. tepidum likely contribute to the major physiological differences observed between these related purple sulfur bacteria that inhabit distinct ecological niches.

tRNA functional signatures classify plastids as late-branching cyanobacteria.

BACKGROUND: Eukaryotes acquired the trait of oxygenic photosynthesis through endosymbiosis of the cyanobacterial progenitor of plastid organelles. Despite recent advances in the phylogenomics of Cyanobacteria, the phylogenetic root of plastids remains controversial. Although a single origin of plastids by endosymbiosis is broadly supported, recent phylogenomic studies are contradictory on whether plastids branch early or late within Cyanobacteria. One underlying cause may be poor fit of evolutionary models to complex phylogenomic data. RESULTS: Using Posterior Predictive Analysis, we show that recently applied evolutionary models poorly fit three phylogenomic datasets curated from cyanobacteria and plastid genomes because of heterogeneities in both substitution processes across sites and of compositions across lineages. To circumvent these sources of bias, we developed CYANO-MLP, a machine learning algorithm that consistently and accurately phylogenetically classifies ("phyloclassifies") cyanobacterial genomes to their clade of origin based on bioinformatically predicted function-informative features in tRNA gene complements. Classification of cyanobacterial genomes with CYANO-MLP is accurate and robust to deletion of clades, unbalanced sampling, and compositional heterogeneity in input tRNA data. CYANO-MLP consistently classifies plastid genomes into a late-branching cyanobacterial sub-clade containing single-cell, starch-producing, nitrogen-fixing ecotypes, consistent with metabolic and gene transfer data. CONCLUSIONS: Phylogenomic data of cyanobacteria and plastids exhibit both site-process heterogeneities and compositional heterogeneities across lineages. These aspects of the data require careful modeling to avoid bias in phylogenomic estimation. Furthermore, we show that amino acid recoding strategies may be insufficient to mitigate bias from compositional heterogeneities. However, the combination of our novel tRNA-specific strategy with machine learning in CYANO-MLP appears robust to these sources of bias with high accuracy in phyloclassification of cyanobacterial genomes. CYANO-MLP consistently classifies plastids as late-branching Cyanobacteria, consistent with independent evidence from signature-based approaches and some previous phylogenetic studies.

Soil microbial community composition in tallgrass prairie restorations converge with remnants across a 27-year chronosequence.

Restoration and management of natural ecosystems is a critical strategy in mitigating global biodiversity loss. This is exemplified in the American Midwest by efforts aimed at reclaiming historical grasslands lost to high-yield agriculture. While restorations traditionally take the form of plant reintroduction and management, advances in microbial analyses suggest that soil communities could be indicators restoration success. However, current understanding of key microbial taxa and functional activities in both natural and restored ecosystems is limited. Here, we investigated the impact of nearly 30 years of carefully managed restoration on soil microbial communities at the Nachusa Grasslands in northern Illinois, USA. We characterized bacterial and archaeal communities in a chronosequence of restored tallgrass prairies ranging from 1 to 27 years old across a growing season and compared them to communities in pre-restoration agricultural fields and remnant prairies. Results indicate that older restorations harboured communities statistically distinct from newer restorations. These communities converged toward those in local prairie remnants, suggesting that plant-focussed restoration has yielded soil bacterial communities reflective of a successful restoration. Recovery of microbial clades within the Verrucomicrobia and Acidobacteria are an important feature of this convergence, and these groups could be targeted for future soil-focussed, bottom-up restoration studies.

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.

tRNA signatures reveal a polyphyletic origin of SAR11 strains among alphaproteobacteria.

Molecular phylogenetics and phylogenomics are subject to noise from horizontal gene transfer (HGT) and bias from convergence in macromolecular compositions. Extensive variation in size, structure and base composition of alphaproteobacterial genomes has complicated their phylogenomics, sparking controversy over the origins and closest relatives of the SAR11 strains. SAR11 are highly abundant, cosmopolitan aquatic Alphaproteobacteria with streamlined, A+T-biased genomes. A dominant view holds that SAR11 are monophyletic and related to both Rickettsiales and the ancestor of mitochondria. Other studies dispute this, finding evidence of a polyphyletic origin of SAR11 with most strains distantly related to Rickettsiales. Although careful evolutionary modeling can reduce bias and noise in phylogenomic inference, entirely different approaches may be useful to extract robust phylogenetic signals from genomes. Here we develop simple phyloclassifiers from bioinformatically derived tRNA Class-Informative Features (CIFs), features predicted to target tRNAs for specific interactions within the tRNA interaction network. Our tRNA CIF-based model robustly and accurately classifies alphaproteobacterial genomes into one of seven undisputed monophyletic orders or families, despite great variability in tRNA gene complement sizes and base compositions. Our model robustly rejects monophyly of SAR11, classifying all but one strain as Rhizobiales with strong statistical support. Yet remarkably, conventional phylogenetic analysis of tRNAs classifies all SAR11 strains identically as Rickettsiales. We attribute this discrepancy to convergence of SAR11 and Rickettsiales tRNA base compositions. Thus, tRNA CIFs appear more robust to compositional convergence than tRNA sequences generally. Our results suggest that tRNA-CIF-based phyloclassification is robust to HGT of components of the tRNA interaction network, such as aminoacyl-tRNA synthetases. We explain why tRNAs are especially advantageous for prediction of traits governing macromolecular interactions from genomic data, and why such traits may be advantageous in the search for robust signals to address difficult problems in classification and phylogeny.

Evolution: Spectral speciation.

Speciation is a complex process sparked by multitudes of environmental stressors and culminating in adaptive, and perhaps novel, phenotypic traits. A new study presents evidence supporting spectral niche-partitioning in a cyanobacterial clade specializing in far-red photosynthesis.

Restoration age and reintroduced bison may shape soil bacterial communities in restored tallgrass prairies.

Knowledge of how habitat restoration shapes soil microbial communities often is limited despite their critical roles in ecosystem function. Soil community diversity and composition change after restoration, but the trajectory of these successional changes may be influenced by disturbances imposed for habitat management. We studied soil bacterial communities in a restored tallgrass prairie chronosequence for >6 years to document how diversity and composition changed with age, management through fire, and grazing by reintroduced bison, and in comparison to pre-restoration agricultural fields and remnant prairies. Soil C:N increased with restoration age and bison, and soil pH first increased and then declined with age, although bison weakened this pattern. Bacterial richness and diversity followed a similar hump-shaped pattern as soil pH, such that the oldest restorations approached the low diversity of remnant prairies. ß-diversity patterns indicated that composition in older restorations with bison resembled bison-free sites, but over time they became more distinct. In contrast, younger restorations with bison maintained unique compositions throughout the study, suggesting bison disturbances may cause a different successional trajectory. We used a novel random forest approach to identify taxa that indicate these differences, finding that they were frequently associated with bacteria that respond to grazing in other grasslands.

Taxonomic and functional restoration of tallgrass prairie soil microbial communities in comparison to remnant and agricultural soils.

Restoring ecosystems requires the re-establishment of diverse soil microbial communities that drive critical ecosystem functions. In grasslands, restoration and management require the application of disturbances like fire and grazing. Disturbances can shape microbial taxonomic composition and potentially functional composition as well. We characterized taxonomic and functional gene composition of soil communities using whole genome shotgun metagenomic sequencing to determine how restored soil communities differed from pre-restoration agricultural soils and original remnant soils, how management affects soil microbes, and whether restoration and management affect the number of microbial genes associated with carbohydrate degradation. We found distinct differences in both taxonomic and functional diversity and composition among restored, remnant, and agricultural soils. Remnant soils had low taxonomic and functional richness and diversity, as well as distinct composition, indicating that restoration of agricultural soils does not re-create soil microbial communities that match remnants. Prescribed fire management increased functional diversity, which also was higher in more recently planted restorations. Finally, restored and post-fire soils included high abundances of genes encoding cellulose-degrading enzymes, so restorations and their ongoing management can potentially support functions important in carbon cycling.

Discovery of Chlorophyll d: Isolation and Characterization of a Far-Red Cyanobacterium from the Original Site of Manning and Strain (1943) at Moss Beach, California.

We have isolated a chlorophyll-d-containing cyanobacterium from the intertidal field site at Moss Beach, on the coast of Central California, USA, where Manning and Strain (1943) originally discovered this far-red chlorophyll. Here, we present the cyanobacterium's environmental description, culturing procedure, pigment composition, ultrastructure, and full genome sequence. Among cultures of far-red cyanobacteria obtained from red algae from the same site, this strain was an epiphyte on a brown macroalgae. Its Qyin vivo absorbance peak is centered at 704-705 nm, the shortest wavelength observed thus far among the various known Acaryochloris strains. Its Chl a/Chl d ratio was 0.01, with Chl d accounting for 99% of the total Chl d and Chl a mass. TEM imagery indicates the absence of phycobilisomes, corroborated by both pigment spectra and genome analysis. The Moss Beach strain codes for only a single set of genes for producing allophycocyanin. Genomic sequencing yielded a 7.25 Mbp circular chromosome and 10 circular plasmids ranging from 16 kbp to 394 kbp. We have determined that this strain shares high similarity with strain S15, an epiphyte of red algae, while its distinct gene complement and ecological niche suggest that this strain could be the closest known relative to the original Chl d source of Manning and Strain (1943). The Moss Beach strain is designated Acaryochloris sp. (marina) strain Moss Beach.

Characterization of photosystem I antenna proteins in the prasinophyte Ostreococcus tauri.

Prasinophyceae are a broad class of early-branching eukaryotic green algae. These picophytoplankton are found ubiquitously throughout the ocean and contribute considerably to global carbon-fixation. Ostreococcus tauri, as the first sequenced prasinophyte, is a model species for studying the functional evolution of light-harvesting systems in photosynthetic eukaryotes. In this study we isolated and characterized O. tauri pigment-protein complexes. Two photosystem I (PSI) fractions were obtained by sucrose density gradient centrifugation in addition to free light-harvesting complex (LHC) fraction and photosystem II (PSII) core fractions. The smaller PSI fraction contains the PSI core proteins, LHCI, which are conserved in all green plants, Lhcp1, a prasinophyte-specific LHC protein, and the minor, monomeric LHCII proteins CP26 and CP29. The larger PSI fraction contained the same antenna proteins as the smaller, with the addition of Lhca6 and Lhcp2, and a 30% larger absorption cross-section. When O. tauri was grown under high-light conditions, only the smaller PSI fraction was present. The two PSI preparations were also found to be devoid of the far-red chlorophyll fluorescence (715-730 nm), a signature of PSI in oxygenic phototrophs. These unique features of O. tauri PSI may reflect primitive light-harvesting systems in green plants and their adaptation to marine ecosystems. Possible implications for the evolution of the LHC-superfamily in photosynthetic eukaryotes are discussed.

Niche adaptation and genome expansion in the chlorophyll d-producing cyanobacterium Acaryochloris marina.

Acaryochloris marina is a unique cyanobacterium that is able to produce chlorophyll d as its primary photosynthetic pigment and thus efficiently use far-red light for photosynthesis. Acaryochloris species have been isolated from marine environments in association with other oxygenic phototrophs, which may have driven the niche-filling introduction of chlorophyll d. To investigate these unique adaptations, we have sequenced the complete genome of A. marina. The DNA content of A. marina is composed of 8.3 million base pairs, which is among the largest bacterial genomes sequenced thus far. This large array of genomic data is distributed into nine single-copy plasmids that code for >25% of the putative ORFs. Heavy duplication of genes related to DNA repair and recombination (primarily recA) and transposable elements could account for genetic mobility and genome expansion. We discuss points of interest for the biosynthesis of the unusual pigments chlorophyll d and alpha-carotene and genes responsible for previously studied phycobilin aggregates. Our analysis also reveals that A. marina carries a unique complement of genes for these phycobiliproteins in relation to those coding for antenna proteins related to those in Prochlorococcus species. The global replacement of major photosynthetic pigments appears to have incurred only minimal specializations in reaction center proteins to accommodate these alternate pigments. These features clearly show that the genus Acaryochloris is a fitting candidate for understanding genome expansion, gene acquisition, ecological adaptation, and photosystem modification in the cyanobacteria.

Metabolic flexibility revealed in the genome of the cyst-forming alpha-1 proteobacterium Rhodospirillum centenum.

BACKGROUND: Rhodospirillum centenum is a photosynthetic non-sulfur purple bacterium that favors growth in an anoxygenic, photosynthetic N2-fixing environment. It is emerging as a genetically amenable model organism for molecular genetic analysis of cyst formation, photosynthesis, phototaxis, and cellular development. Here, we present an analysis of the genome of this bacterium. RESULTS: R. centenum contains a singular circular chromosome of 4,355,548 base pairs in size harboring 4,105 genes. It has an intact Calvin cycle with two forms of Rubisco, as well as a gene encoding phosphoenolpyruvate carboxylase (PEPC) for mixotrophic CO2 fixation. This dual carbon-fixation system may be required for regulating internal carbon flux to facilitate bacterial nitrogen assimilation. Enzymatic reactions associated with arsenate and mercuric detoxification are rare or unique compared to other purple bacteria. Among numerous newly identified signal transduction proteins, of particular interest is a putative bacteriophytochrome that is phylogenetically distinct from a previously characterized R. centenum phytochrome, Ppr. Genes encoding proteins involved in chemotaxis as well as a sophisticated dual flagellar system have also been mapped. CONCLUSIONS: Remarkable metabolic versatility and a superior capability for photoautotrophic carbon assimilation is evident in R. centenum.

Integrating Markov clustering and molecular phylogenetics to reconstruct the cyanobacterial species tree from conserved protein families.

Attempts to classify living organisms by their physical characteristics are as old as biology itself. The advent of protein and DNA sequencing--most notably the use of 16S ribosomal RNA--defined a new level of classification that now forms our basic understanding of the history of life on earth. High-throughput sequencing currently provides DNA sequences at an unprecedented rate, not only providing a wealth of information but also posing considerable analytical challenges. Here we present comparative genomics-based methods useful for automating evolutionary analysis between any number of species. As a practical example, we applied our method to the well-studied cyanobacterial lineage. The 24 cyanobacterial genomes compared here occupy a wide variety of environmental niches and play major roles in global carbon and nitrogen cycles. By integrating phylogenetic data inferred for upward of 1,000 protein-coding genes common to all or most cyanobacteria, we have reconstructed an evolutionary history of the phylum, establishing a framework for resolving key issues regarding the evolution of their metabolic and phenotypic diversity. Greater resolution on individual branches can be attained by telescoping inward to the larger set of conserved proteins between fewer taxa. The construction of all individual protein phylogenies allows for quantitative tree scoring, providing insight into the evolutionary history of each protein family as well as probing the limits of phylogenetic resolution. The tools incorporated here are fast, computationally tractable, and easily extendable to other phyla and provide a scaleable framework for contrasting and integrating the information present in thousands of protein-coding genes within related genomes.

The genome of Heliobacterium modesticaldum, a phototrophic representative of the Firmicutes containing the simplest photosynthetic apparatus.

Despite the fact that heliobacteria are the only phototrophic representatives of the bacterial phylum Firmicutes, genomic analyses of these organisms have yet to be reported. Here we describe the complete sequence and analysis of the genome of Heliobacterium modesticaldum, a thermophilic species belonging to this unique group of phototrophs. The genome is a single 3.1-Mb circular chromosome containing 3,138 open reading frames. As suspected from physiological studies of heliobacteria that have failed to show photoautotrophic growth, genes encoding enzymes for known autotrophic pathways in other phototrophic organisms, including ribulose bisphosphate carboxylase (Calvin cycle), citrate lyase (reverse citric acid cycle), and malyl coenzyme A lyase (3-hydroxypropionate pathway), are not present in the H. modesticaldum genome. Thus, heliobacteria appear to be the only known anaerobic anoxygenic phototrophs that are not capable of autotrophy. Although for some cellular activities, such as nitrogen fixation, there is a full complement of genes in H. modesticaldum, other processes, including carbon metabolism and endosporulation, are more genetically streamlined than they are in most other low-G+C gram-positive bacteria. Moreover, several genes encoding photosynthetic functions in phototrophic purple bacteria are not present in the heliobacteria. In contrast to the nutritional flexibility of many anoxygenic phototrophs, the complete genome sequence of H. modesticaldum reveals an organism with a notable degree of metabolic specialization and genomic reduction.

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.

Genome Sequence of Rhodoferax antarcticus ANT.BRT; A Psychrophilic Purple Nonsulfur Bacterium from an Antarctic Microbial Mat.

Rhodoferax antarcticus is an Antarctic purple nonsulfur bacterium and the only characterized anoxygenic phototroph that grows best below 20 °C. We present here a high-quality draft genome of Rfx. antarcticus strain ANT.BRT, isolated from an Antarctic microbial mat. The circular chromosome (3.8 Mbp) of Rfx. antarcticus has a 59.1% guanine + cytosine (GC) content and contains 4036 open reading frames. In addition, the bacterium contains a sizable plasmid (198.6 kbp, 48.4% GC with 226 open reading frames) that comprises about 5% of the total genetic content. Surprisingly, genes encoding light-harvesting complexes 1 and 3 (LH1 and LH3), but not light-harvesting complex 2 (LH2), were identified in the photosynthesis gene cluster of the Rfx. antarcticus genome, a feature that is unique among purple phototrophs. Consistent with physiological studies that showed a strong capacity for nitrogen fixation in Rfx. antarcticus, a nitrogen fixation gene cluster encoding a molybdenum-type nitrogenase was present, but no alternative nitrogenases were identified despite the cold-active phenotype of this phototroph. Genes encoding two forms of ribulose 1,5-bisphosphate carboxylase/oxygenase were present in the Rfx. antarcticus genome, a feature that likely provides autotrophic flexibility under varying environmental conditions. Lastly, genes for assembly of both type IV pili and flagella are present, with the latter showing an unusual degree of clustering. This report represents the first genomic analysis of a psychrophilic anoxygenic phototroph and provides a glimpse of the genetic basis for maintaining a phototrophic lifestyle in a permanently cold, yet highly variable, environment.

Correction: The First Complete Plastid Genome from Joinvilleaceae (J. ascendens; Poales) Shows Unique and Unpredicted Rearrangements.

[This corrects the article DOI: 10.1371/journal.pone.0163218.].

The First Complete Plastid Genome from Joinvilleaceae (J. ascendens; Poales) Shows Unique and Unpredicted Rearrangements.

Joinvilleaceae is a family of tropical grass-like monocots that comprises only the genus Joinvillea. Previous studies have placed Joinvilleaceae in close phylogenetic proximity to the well-studied grass family. A full plastome sequence was determined and characterized for J. ascendens. The plastome was sequenced with next generation methods, fully assembled de novo and annotated. The assembly revealed two novel inversions specific to the Joinvilleaceae lineage and at least one novel plastid inversion in the Joinvilleaceae-Poaceae lineage. Two previously documented inversions in the Joinvilleaceae-Poaceae lineage and one previously documented inversion in the Poaceae lineage were also verified. Inversion events were identified visually and verified computationally by simulation mutations. Additionally, the loss and subsequent degradation of the accD gene in order Poales was explored extensively in Poaceae and J. ascendens. The two novel inversions along with changes in gene composition between families better delimited lineages in the Poales. The presence of large inversions and subsequent reversals in this small family suggested a high potential for large-scale rearrangements to occur in plastid genomes.

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

OP9 is a yet-uncultivated bacterial lineage found in geothermal systems, petroleum reservoirs, anaerobic digesters and wastewater treatment facilities. Here we use single-cell and metagenome sequencing to obtain two distinct, nearly complete OP9 genomes, one constructed from single cells sorted from hot spring sediments and the other derived from binned metagenomic contigs from an in situ-enriched cellulolytic, thermophilic community. Phylogenomic analyses support the designation of OP9 as a candidate phylum for which we propose the name 'Atribacteria'. Although a plurality of predicted proteins is most similar to those from Firmicutes, the presence of key genes suggests a diderm cell envelope. Metabolic reconstruction from the core genome suggests an anaerobic lifestyle based on sugar fermentation by Embden-Meyerhof glycolysis with production of hydrogen, acetate and ethanol. Putative glycohydrolases and an endoglucanase may enable catabolism of (hemi)cellulose in thermal environments. This study lays a foundation for understanding the physiology and ecological role of the 'Atribacteria'.

Coordinating environmental genomics and geochemistry reveals metabolic transitions in a hot spring ecosystem.

We have constructed a conceptual model of biogeochemical cycles and metabolic and microbial community shifts within a hot spring ecosystem via coordinated analysis of the "Bison Pool" (BP) Environmental Genome and a complementary contextual geochemical dataset of ~75 geochemical parameters. 2,321 16S rRNA clones and 470 megabases of environmental sequence data were produced from biofilms at five sites along the outflow of BP, an alkaline hot spring in Sentinel Meadow (Lower Geyser Basin) of Yellowstone National Park. This channel acts as a >22 m gradient of decreasing temperature, increasing dissolved oxygen, and changing availability of biologically important chemical species, such as those containing nitrogen and sulfur. Microbial life at BP transitions from a 92 °C chemotrophic streamer biofilm community in the BP source pool to a 56 °C phototrophic mat community. We improved automated annotation of the BP environmental genomes using BLAST-based Markov clustering. We have also assigned environmental genome sequences to individual microbial community members by complementing traditional homology-based assignment with nucleotide word-usage algorithms, allowing more than 70% of all reads to be assigned to source organisms. This assignment yields high genome coverage in dominant community members, facilitating reconstruction of nearly complete metabolic profiles and in-depth analysis of the relation between geochemical and metabolic changes along the outflow. We show that changes in environmental conditions and energy availability are associated with dramatic shifts in microbial communities and metabolic function. We have also identified an organism constituting a novel phylum in a metabolic "transition" community, located physically between the chemotroph- and phototroph-dominated sites. The complementary analysis of biogeochemical and environmental genomic data from BP has allowed us to build ecosystem-based conceptual models for this hot spring, reconstructing whole metabolic networks in order to illuminate community roles in shaping and responding to geochemical variability.