The Yellowstone caldera contains the most numerous and diverse geothermal systems on Earth, yielding an extensive array of unique high-temperature environments that host a variety of deeply-rooted and understudied Archaea, Bacteria and Eukarya. The combination of extreme temperature and chemical conditions encountered in geothermal environments often results in considerably less microbial diversity than other terrestrial habitats and offers a tremendous opportunity for studying the structure and function of indigenous microbial communities and for establishing linkages between putative metabolisms and element cycling. Metagenome sequence (14–15,000 Sanger reads per site) was obtained for five high-temperature (>65°C) chemotrophic microbial communities sampled from geothermal springs (or pools) in Yellowstone National Park (YNP) that exhibit a wide range in geochemistry including pH, dissolved sulfide, dissolved oxygen and ferrous iron. Metagenome data revealed significant differences in the predominant phyla associated with each of these geochemical environments. Novel members of the Sulfolobales are dominant in low pH environments, while other Crenarchaeota including distantly-related Thermoproteales and Desulfurococcales populations dominate in suboxic sulfidic sediments. Several novel archaeal groups are well represented in an acidic (pH 3) Fe-oxyhydroxide mat, where a higher O2 influx is accompanied with an increase in archaeal diversity. The presence or absence of genes and pathways important in S oxidation-reduction, H2-oxidation, and aerobic respiration (terminal oxidation) provide insight regarding the metabolic strategies of indigenous organisms present in geothermal systems. Multiple-pathway and protein-specific functional analysis of metagenome sequence data corroborated results from phylogenetic analyses and clearly demonstrate major differences in metabolic potential across sites. The distribution of functional genes involved in electron transport is consistent with the hypothesis that geochemical parameters (e.g., pH, sulfide, Fe, O2) control microbial community structure and function in YNP geothermal springs.

The function of cells in their native habitat often cannot be reliably predicted from genomic data or from physiology studies of isolates. Traditional experimental approaches to study the function of taxonomically and metabolically diverse microbiomes are limited by their destructive nature, low spatial resolution, or low throughput. Recently developed technologies can offer new insights into cellular function in natural and human-made systems and how microorganisms interact with and shape the environments that they inhabit. In this Review, we provide an overview of these next-generation physiology approaches and discuss how the non-destructive analysis of cellular phenotypes, in combination with the separation of the target cells for downstream analyses, provide powerful new, complementary ways to study microbiome function. We anticipate that the widespread application of next-generation physiology approaches will transform the field of microbial ecology and dramatically improve our understanding of how microorganisms function in their native environment. ToC blurbIn this Review, Hatzenpichler et al. introduce next-generation physiology, which is a suite of new techniques that enable to investigate the phenotypes of individual cells in a non-destructive manner. Next-generation physiology complements genomics and culturing and provides new insights into microbiome function.

BackgroundThree methods were developed for the application of stoichiometry-based network analysis approaches including elementary mode analysis to the study of mass and energy flows in microbial communities. Each has distinct advantages and disadvantages suitable for analyzing systems with different degrees of complexity and a priori knowledge. These approaches were tested and compared using data from the thermophilic, phototrophic mat communities from Octopus and Mushroom Springs in Yellowstone National Park (USA). The models were based on three distinct microbial guilds: oxygenic phototrophs, filamentous anoxygenic phototrophs, and sulfate-reducing bacteria. Two phases, day and night, were modeled to account for differences in the sources of mass and energy and the routes available for their exchange.ResultsThe in silico models were used to explore fundamental questions in ecology including the prediction of and explanation for measured relative abundances of primary producers in the mat, theoretical tradeoffs between overall productivity and the generation of toxic by-products, and the relative robustness of various guild interactions.ConclusionThe three modeling approaches represent a flexible toolbox for creating cellular metabolic networks to study microbial communities on scales ranging from cells to ecosystems. A comparison of the three methods highlights considerations for selecting the one most appropriate for a given microbial system. For instance, communities represented only by metagenomic data can be modeled using the pooled method which analyzes a community's total metabolic potential without attempting to partition enzymes to different organisms. Systems with extensive a priori information on microbial guilds can be represented using the compartmentalized technique, employing distinct control volumes to separate guild-appropriate enzymes and metabolites. If the complexity of a compartmentalized network creates an unacceptable computational burden, the nested analysis approach permits greater scalability at the cost of more user intervention through multiple rounds of pathway analysis.

The biological reduction of dinitrogen (N) to ammonia (NH) by nitrogenase is an energetically demanding reaction that requires low-potential electrons and ATP; however, pathways used to deliver the electrons from central metabolism to the reductants of nitrogenase, ferredoxin or flavodoxin, remain unknown for many diazotrophic microbes. The FixABCX protein complex has been proposed to reduce flavodoxin or ferredoxin using NADH as the electron donor in a process known as electron bifurcation. Herein, the FixABCX complex from Azotobacter vinelandii was purified and demonstrated to catalyze an electron bifurcation reaction: oxidation of NADH (E = -320 mV) coupled to reduction of flavodoxin semiquinone (E = -460 mV) and reduction of coenzyme Q (E = 10 mV). Knocking out fix genes rendered Δrnf A. vinelandii cells unable to fix dinitrogen, confirming that the FixABCX system provides another route for delivery of electrons to nitrogenase. Characterization of the purified FixABCX complex revealed the presence of flavin and iron-sulfur cofactors confirmed by native mass spectrometry, electron paramagnetic resonance spectroscopy, and transient absorption spectroscopy. Transient absorption spectroscopy further established the presence of a short-lived flavin semiquinone radical, suggesting that a thermodynamically unstable flavin semiquinone may participate as an intermediate in the transfer of an electron to flavodoxin. A structural model of FixABCX, generated using chemical cross-linking in conjunction with homology modeling, revealed plausible electron transfer pathways to both high- and low-potential acceptors. Overall, this study informs a mechanism for electron bifurcation, offering insight into a unique method for delivery of low-potential electrons required for energy-intensive biochemical conversions.

The Aquificales are thermophilic microorganisms that inhabit hydrothermal systems worldwide and are considered one of the earliest lineages of the domain Bacteria. We analyzed metagenome sequence obtained from six thermal “filamentous streamer” communities (∼40 Mbp per site), which targeted three different groups of Aquificales found in Yellowstone National Park (YNP). Unassembled metagenome sequence and PCR-amplified 16S rRNA gene libraries revealed that acidic, sulfidic sites were dominated by Hydrogenobaculum (Aquificaceae) populations, whereas the circum-neutral pH (6.5–7.8) sites containing dissolved sulfide were dominated by Sulfurihydrogenibium spp. (Hydrogenothermaceae). Thermocrinis (Aquificaceae) populations were found primarily in the circum-neutral sites with undetectable sulfide, and to a lesser extent in one sulfidic system at pH 8. Phylogenetic analysis of assembled sequence containing 16S rRNA genes as well as conserved protein-encoding genes revealed that the composition and function of these communities varied across geochemical conditions. Each Aquificales lineage contained genes for CO2 fixation by the reverse-TCA cycle, but only the Sulfurihydrogenibium populations perform citrate cleavage using ATP citrate lyase (Acl). The Aquificaceae populations use an alternative pathway catalyzed by two separate enzymes, citryl-CoA synthetase (Ccs), and citryl-CoA lyase (Ccl). All three Aquificales lineages contained evidence of aerobic respiration, albeit due to completely different types of heme Cu oxidases (subunit I) involved in oxygen reduction. The distribution of Aquificales populations and differences among functional genes involved in energy generation and electron transport is consistent with the hypothesis that geochemical parameters (e.g., pH, sulfide, H2, O2) have resulted in niche specialization among members of the Aquificales.