Zackary J. Jay

Metagenomes from high-temperature chemotrophic systems reveal geochemical controls on microbial community structure and function.

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.

In silico approaches to study mass and energy flows in microbial consortia: a syntrophic case study

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 communitys 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 YNP metagenome project: environmental parameters responsible for microbial distribution in the Yellowstone geothermal ecosystem.

The Yellowstone geothermal complex contains over 10,000 diverse geothermal features that host numerous phylogenetically deeply rooted and poorly understood archaea, bacteria, and viruses. Microbial communities in high-temperature environments are generally less diverse than soil, marine, sediment, or lake habitats and therefore offer a tremendous opportunity for studying the structure and function of different model microbial communities using environmental metagenomics. One of the broader goals of this study was to establish linkages among microbial distribution, metabolic potential, and environmental variables. Twenty geochemically distinct geothermal ecosystems representing a broad spectrum of Yellowstone hot-spring environments were used for metagenomic and geochemical analysis and included approximately equal numbers of: (1) phototrophic mats, (2) “filamentous streamer” communities, and (3) archaeal-dominated sediments. The metagenomes were analyzed using a suite of complementary and integrative bioinformatic tools, including phylogenetic and functional analysis of both individual sequence reads and assemblies of predominant phylotypes. This volume identifies major environmental determinants of a large number of thermophilic microbial lineages, many of which have not been fully described in the literature nor previously cultivated to enable functional and genomic analyses. Moreover, protein family abundance comparisons and in-depth analyses of specific genes and metabolic pathways relevant to these hot-spring environments reveal hallmark signatures of metabolic capabilities that parallel the distribution of phylotypes across specific types of geochemical environments.

Metagenome Sequence Analysis of Filamentous Microbial Communities Obtained from Geochemically Distinct Geothermal Channels Reveals Specialization of Three Aquificales Lineages

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.

Community Structure and Function of High-Temperature Chlorophototrophic Microbial Mats Inhabiting Diverse Geothermal Environments

Six phototrophic microbial mat communities from different geothermal springs (YNP) were studied using metagenome sequencing and geochemical analyses. The primary goals of this work were to determine differences in community composition of high-temperature phototrophic mats distributed across the Yellowstone geothermal ecosystem, and to identify metabolic attributes of predominant organisms present in these communities that may correlate with environmental attributes important in niche differentiation. Random shotgun metagenome sequences from six phototrophic communities (average ∼53 Mbp/site) were subjected to multiple taxonomic, phylogenetic, and functional analyses. All methods, including G + C content distribution, MEGAN analyses, and oligonucleotide frequency-based clustering, provided strong support for the dominant community members present in each site. Cyanobacteria were only observed in non-sulfidic sites; de novo assemblies were obtained for Synechococcus-like populations at Chocolate Pots (CP_7) and Fischerella-like populations at White Creek (WC_6). Chloroflexi-like sequences (esp. Roseiflexus and/or Chloroflexus spp.) were observed in all six samples and contained genes involved in bacteriochlorophyll biosynthesis and the 3-hydroxypropionate carbon fixation pathway. Other major sequence assemblies were obtained for a Chlorobiales population from CP_7 (proposed family Thermochlorobacteriaceae), and an anoxygenic, sulfur-oxidizing Thermochromatium-like (Gamma-proteobacteria) population from Bath Lake Vista Annex (BLVA_20). Additional sequence coverage is necessary to establish more complete assemblies of other novel bacteria in these sites (e.g., Bacteroidetes and Firmicutes); however, current assemblies suggested that several of these organisms play important roles in heterotrophic and fermentative metabolisms. Definitive linkages were established between several of the dominant phylotypes present in these habitats and important functional processes such as photosynthesis, carbon fixation, sulfur oxidation, and fermentation.

Phylogenetic and functional analysis of metagenome sequence from high-temperature archaeal habitats demonstrate linkages between metabolic potential and geochemistry

Geothermal habitats in Yellowstone National Park (YNP) provide an unparalleled opportunity to understand the environmental factors that control the distribution of archaea in thermal habitats. Here we describe, analyze, and synthesize metagenomic and geochemical data collected from seven high-temperature sites that contain microbial communities dominated by archaea relative to bacteria. The specific objectives of the study were to use metagenome sequencing to determine the structure and functional capacity of thermophilic archaeal-dominated microbial communities across a pH range from 2.5 to 6.4 and to discuss specific examples where the metabolic potential correlated with measured environmental parameters and geochemical processes occurring in situ. Random shotgun metagenome sequence (∼40–45 Mb Sanger sequencing per site) was obtained from environmental DNA extracted from high-temperature sediments and/or microbial mats and subjected to numerous phylogenetic and functional analyses. Analysis of individual sequences (e.g., MEGAN and G + C content) and assemblies from each habitat type revealed the presence of dominant archaeal populations in all environments, 10 of whose genomes were largely reconstructed from the sequence data. Analysis of protein family occurrence, particularly of those involved in energy conservation, electron transport, and autotrophic metabolism, revealed significant differences in metabolic strategies across sites consistent with differences in major geochemical attributes (e.g., sulfide, oxygen, pH). These observations provide an ecological basis for understanding the distribution of indigenous archaeal lineages across high-temperature systems of YNP.

Niche specialization of novel Thaumarchaeota to oxic and hypoxic acidic geothermal springs of Yellowstone National Park

Novel lineages of the phylum Thaumarchaeota are endemic to thermal habitats, and may exhibit physiological capabilities that are not yet observed in members of this phylum. The primary goals of this study were to conduct detailed phylogenetic and functional analyses of metagenome sequence assemblies of two different thaumarchaeal populations found in high-temperature (65–72 °C), acidic (pH∼3) iron oxide and sulfur sediment environments of Yellowstone National Park (YNP). Metabolic reconstruction was coupled with detailed geochemical measurements of each geothermal habitat and reverse-transcriptase PCR to confirm the in situ activity of these populations. Phylogenetic analyses of ribosomal and housekeeping proteins place these archaea near the root of the thaumarchaeal branch. Metabolic reconstruction suggests that these populations are chemoorganotrophic and couple growth with the reduction of oxygen or nitrate in iron oxide habitats, or sulfur in hypoxic sulfur sediments. The iron oxide population has the potential for growth via the oxidation of sulfide to sulfate using a novel reverse sulfate reduction pathway. Possible carbon sources include aromatic compounds (for example, 4-hydroxyphenylacetate), complex carbohydrates (for example, starch), oligopeptides and amino acids. Both populations contain a type III ribulose bisphosphate carboxylase/oxygenase used for carbon dioxide fixation or adenosine monophosphate salvage. No evidence for the oxidation of ammonia was obtained from de novo sequence assemblies. Our results show that thermoacidophilic Thaumarchaeota from oxic iron mats and hypoxic sulfur sediments exhibit different respiratory machinery depending on the presence of oxygen versus sulfide, represent deeply rooted lineages within the phylum Thaumarchaeota and are endemic to numerous sites in YNP.

Microbial community structure and sulfur biogeochemistry in mildly‐acidic sulfidic geothermal springs in Yellowstone National Park

Geothermal and hydrothermal waters often contain high concentrations of dissolved sulfide, which reacts with oxygen (abiotically or biotically) to yield elemental sulfur and other sulfur species that may support microbial metabolism. The primary goal of this study was to elucidate predominant biogeochemical processes important in sulfur biogeochemistry by identifying predominant sulfur species and describing microbial community structure within high-temperature, hypoxic, sulfur sediments ranging in pH from 4.2 to 6.1. Detailed analysis of aqueous species and solid phases present in hypoxic sulfur sediments revealed unique habitats containing high concentrations of dissolved sulfide, thiosulfate, and arsenite, as well as rhombohedral and spherical elemental sulfur and/or sulfide phases such as orpiment, stibnite, and pyrite, as well as alunite and quartz. Results from 16S rRNA gene sequencing show that these sediments are dominated by Crenarchaeota of the orders Desulfurococcales and Thermoproteales. Numerous cultivated representatives of these lineages, as well as the Thermoproteales strain (WP30) isolated in this study, require complex sources of carbon and respire elemental sulfur. We describe a new archaeal isolate (strain WP30) belonging to the order Thermoproteales (phylum Crenarchaeota, 98% identity to Pyrobaculum/Thermoproteus spp. 16S rRNA genes), which was obtained from sulfur sediments using in situ geochemical composition to design cultivation medium. This isolate produces sulfide during growth, which further promotes the formation of sulfide phases including orpiment, stibnite, or pyrite, depending on solution conditions. Geochemical, molecular, and physiological data were integrated to suggest primary factors controlling microbial community structure and function in high-temperature sulfur sediments.

Assembly and Succession of Iron Oxide Microbial Mat Communities in Acidic Geothermal Springs

Biomineralized ferric oxide microbial mats are ubiquitous features on Earth, are common in hot springs of Yellowstone National Park (YNP, WY, USA), and form due to direct interaction between microbial and physicochemical processes. The overall goal of this study was to determine the contribution of different community members to the assembly and succession of acidic high-temperature Fe(III)-oxide mat ecosystems. Spatial and temporal changes in Fe(III)-oxide accretion and the abundance of relevant community members were monitored over 70 days using sterile glass microscope slides incubated in the outflow channels of two acidic geothermal springs (pH = 3–3.5; temperature = 68–75°C) in YNP. Hydrogenobaculum spp. were the most abundant taxon identified during early successional stages (4–40 days), and have been shown to oxidize arsenite, sulfide, and hydrogen coupled to oxygen reduction. Iron-oxidizing populations of Metallosphaera yellowstonensis were detected within 4 days, and reached steady-state levels within 14–30 days, corresponding to visible Fe(III)-oxide accretion. Heterotrophic archaea colonized near 30 days, and emerged as the dominant functional guild after 70 days and in mature Fe(III)-oxide mats (1–2 cm thick). First-order rate constants of Fe(III)-oxide accretion ranged from 0.046 to 0.05 day−1, and in situ microelectrode measurements showed that the oxidation of Fe(II) is limited by the diffusion of O2 into the Fe(III)-oxide mat. The formation of microterracettes also implicated O2 as a major variable controlling microbial growth and subsequent mat morphology. The assembly and succession of Fe(III)-oxide mat communities follows a repeatable pattern of colonization by lithoautotrophic organisms, and the subsequent growth of diverse organoheterotrophs. The unique geochemical signatures and micromorphology of extant biomineralized Fe(III)-oxide mats are also useful for understanding other Fe(II)-oxidizing systems.

Predominant Acidilobus-Like Populations from Geothermal Environments in Yellowstone National Park Exhibit Similar Metabolic Potential in Different Hypoxic Microbial Communities

ABSTRACT High-temperature (>70°C) ecosystems in Yellowstone National Park (YNP) provide an unparalleled opportunity to study chemotrophic archaea and their role in microbial community structure and function under highly constrained geochemical conditions. Acidilobus spp. (order Desulfurococcales) comprise one of the dominant phylotypes in hypoxic geothermal sulfur sediment and Fe(III)-oxide environments along with members of the Thermoproteales and Sulfolobales. Consequently, the primary goals of the current study were to analyze and compare replicate de novo sequence assemblies of Acidilobus-like populations from four different mildly acidic (pH 3.3 to 6.1) high-temperature (72°C to 82°C) environments and to identify metabolic pathways and/or protein-encoding genes that provide a detailed foundation of the potential functional role of these populations in situ. De novo assemblies of the highly similar Acidilobus-like populations (>99% 16S rRNA gene identity) represent near-complete consensus genomes based on an inventory of single-copy genes, deduced metabolic potential, and assembly statistics generated across sites. Functional analysis of coding sequences and confirmation of gene transcription by Acidilobus-like populations provide evidence that they are primarily chemoorganoheterotrophs, generating acetyl coenzyme A (acetyl-CoA) via the degradation of carbohydrates, lipids, and proteins, and auxotrophic with respect to several external vitamins, cofactors, and metabolites. No obvious pathways or protein-encoding genes responsible for the dissimilatory reduction of sulfur were identified. The presence of a formate dehydrogenase (Fdh) and other protein-encoding genes involved in mixed-acid fermentation supports the hypothesis that Acidilobus spp. function as degraders of complex organic constituents in high-temperature, mildly acidic, hypoxic geothermal systems.