Mark A. Kozubal

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.

Linking microbial oxidation of arsenic with detection and phylogenetic analysis of arsenite oxidase genes in diverse geothermal environments.

The identification and characterization of genes involved in the microbial oxidation of arsenite will contribute to our understanding of factors controlling As cycling in natural systems. Towards this goal, we recently characterized the widespread occurrence of aerobic arsenite oxidase genes (aroA-like) from pure-culture bacterial isolates, soils, sediments and geothermal mats, but were unable to detect these genes in all geothermal systems where we have observed microbial arsenite oxidation. Consequently, the objectives of the current study were to measure arsenite-oxidation rates in geochemically diverse thermal habitats in Yellowstone National Park (YNP) ranging in pH from 2.6 to 8, and to identify corresponding 16S rRNA and aroA genotypes associated with these arsenite-oxidizing environments. Geochemical analyses, including measurement of arsenite-oxidation rates within geothermal outflow channels, were combined with 16S rRNA gene and aroA functional gene analysis using newly designed primers to capture previously undescribed aroA-like arsenite oxidase gene diversity. The majority of bacterial 16S rRNA gene sequences found in acidic (pH 2.6-3.6) Fe-oxyhydroxide microbial mats were closely related to Hydrogenobaculum spp. (members of the bacterial order Aquificales), while the predominant sequences from near-neutral (pH 6.2-8) springs were affiliated with other Aquificales including Sulfurihydrogenibium spp., Thermocrinis spp. and Hydrogenobacter spp., as well as members of the Deinococci, Thermodesulfobacteria and beta-Proteobacteria. Modified primers designed around previously characterized and newly identified aroA-like genes successfully amplified new lineages of aroA-like genes associated with members of the Aquificales across all geothermal systems examined. The expression of Aquificales aroA-like genes was also confirmed in situ, and the resultant cDNA sequences were consistent with aroA genotypes identified in the same environments. The aroA sequences identified in the current study expand the phylogenetic distribution of known Mo-pterin arsenite oxidase genes, and suggest the importance of three prominent genera of the order Aquificales in arsenite oxidation across geochemically distinct geothermal habitats ranging in pH from 2.6 to 8.

Geoarchaeota: a new candidate phylum in the Archaea from high-temperature acidic iron mats in Yellowstone National Park

Geothermal systems in Yellowstone National Park (YNP) provide an outstanding opportunity to understand the origin and evolution of metabolic processes necessary for life in extreme environments including low pH, high temperature, low oxygen and elevated concentrations of reduced iron. Previous phylogenetic studies of acidic ferric iron mats from YNP have revealed considerable diversity of uncultivated and undescribed archaea. The goal of this study was to obtain replicate de novo genome assemblies for a dominant archaeal population inhabiting acidic iron-oxide mats in YNP. Detailed analysis of conserved ribosomal and informational processing genes indicates that the replicate assemblies represent a new candidate phylum within the domain Archaea referred to here as ‘Geoarchaeota’ or ‘novel archaeal group 1 (NAG1)’. The NAG1 organisms contain pathways necessary for the catabolism of peptides and complex carbohydrates as well as a bacterial-like Form I carbon monoxide dehydrogenase complex likely used for energy conservation. Moreover, this novel population contains genes involved in the metabolism of oxygen including a Type A heme copper oxidase, a bd-type terminal oxidase and a putative oxygen-sensing protoglobin. NAG1 has a variety of unique bacterial-like cofactor biosynthesis and transport genes and a Type3-like CRISPR system. Discovery of NAG1 is critical to our understanding of microbial community structure and function in extant thermophilic iron-oxide mats of YNP, and will provide insight regarding the evolution of Archaea in early Earth environments that may have important analogs active in YNP today.

Microbial Iron Cycling in Acidic Geothermal Springs of Yellowstone National Park: Integrating Molecular Surveys, Geochemical Processes, and Isolation of Novel Fe-Active Microorganisms

Geochemical, molecular, and physiological analyses of microbial isolates were combined to study the geomicrobiology of acidic iron oxide mats in Yellowstone National Park. Nineteen sampling locations from 11 geothermal springs were studied ranging in temperature from 53 to 88°C and pH 2.4 to 3.6. All iron oxide mats exhibited high diversity of crenarchaeal sequences from the Sulfolobales, Thermoproteales, and Desulfurococcales. The predominant Sulfolobales sequences were highly similar to Metallosphaera yellowstonensis str. MK1, previously isolated from one of these sites. Other groups of archaea were consistently associated with different types of iron oxide mats, including undescribed members of the phyla Thaumarchaeota and Euryarchaeota. Bacterial sequences were dominated by relatives of Hydrogenobaculum spp. above 65–70°C, but increased in diversity below 60°C. Cultivation of relevant iron-oxidizing and iron-reducing microbial isolates included Sulfolobus str. MK3, Sulfobacillus str. MK2, Acidicaldus str. MK6, and a new candidate genus in the Sulfolobales referred to as Sulfolobales str. MK5. Strains MK3 and MK5 are capable of oxidizing ferrous iron autotrophically, while strain MK2 oxidizes iron mixotrophically. Similar rates of iron oxidation were measured for M. yellowstonensis str. MK1 and Sulfolobales str. MK5. Biomineralized phases of ferric iron varied among cultures and field sites, and included ferric oxyhydroxides, K-jarosite, goethite, hematite, and scorodite depending on geochemical conditions. Strains MK5 and MK6 are capable of reducing ferric iron under anaerobic conditions with complex carbon sources. The combination of geochemical and molecular data as well as physiological observations of isolates suggests that the community structure of acidic Fe mats is linked with Fe cycling across temperatures ranging from 53 to 88°C.

Isolation and Distribution of a Novel Iron-Oxidizing Crenarchaeon from Acidic Geothermal Springs in Yellowstone National Park

ABSTRACT Novel thermophilic crenarchaea have been observed in Fe(III) oxide microbial mats of Yellowstone National Park (YNP); however, no definitive work has identified specific microorganisms responsible for the oxidation of Fe(II). The objectives of the current study were to isolate and characterize an Fe(II)-oxidizing member of the Sulfolobales observed in previous 16S rRNA gene surveys and to determine the abundance and distribution of close relatives of this organism in acidic geothermal springs containing high concentrations of dissolved Fe(II). Here we report the isolation and characterization of the novel, Fe(II)-oxidizing, thermophilic, acidophilic organism Metallosphaera sp. strain MK1 obtained from a well-characterized acid-sulfate-chloride geothermal spring in Norris Geyser Basin, YNP. Full-length 16S rRNA gene sequence analysis revealed that strain MK1 exhibits only 94.9 to 96.1% sequence similarity to other known Metallosphaera spp. and less than 89.1% similarity to known Sulfolobus spp. Strain MK1 is a facultative chemolithoautotroph with an optimum pH range of 2.0 to 3.0 and an optimum temperature range of 65 to 75°C. Strain MK1 grows optimally on pyrite or Fe(II) sorbed onto ferrihydrite, exhibiting doubling times between 10 and 11 h under aerobic conditions (65°C). The distribution and relative abundance of MK1-like 16S rRNA gene sequences in 14 acidic geothermal springs containing Fe(III) oxide microbial mats were evaluated. Highly related MK1-like 16S rRNA gene sequences (>99% sequence similarity) were consistently observed in Fe(III) oxide mats at temperatures ranging from 55 to 80°C. Quantitative PCR using Metallosphaera-specific primers confirmed that organisms highly similar to strain MK1 comprised up to 40% of the total archaeal community at selected sites. The broad distribution of highly related MK1-like 16S rRNA gene sequences in acidic Fe(III) oxide microbial mats is consistent with the observed characteristics and growth optima of Metallosphaera-like strain MK1 and emphasizes the importance of this newly described taxon in Fe(II) chemolithotrophy in acidic high-temperature environments of YNP.

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.

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.

Terminal Oxidase Diversity and Function in “Metallosphaera yellowstonensis”: Gene Expression and Protein Modeling Suggest Mechanisms of Fe(II) Oxidation in the Sulfolobales

ABSTRACT “Metallosphaera yellowstonensis” is a thermoacidophilic archaeon isolated from Yellowstone National Park that is capable of autotrophic growth using Fe(II), elemental S, or pyrite as electron donors. Analysis of the draft genome sequence from M. yellowstonensis strain MK1 revealed seven different copies of heme copper oxidases (subunit I) in a total of five different terminal oxidase complexes, including doxBCEF, foxABCDEFGHIJ, soxABC, and the soxM supercomplex, as well as a novel hypothetical two-protein doxB-like polyferredoxin complex. Other genes found in M. yellowstonensis with possible roles in S and or Fe cycling include a thiosulfate oxidase (tqoAB), a sulfite oxidase (som), a cbsA cytochrome b 558/566, several small blue copper proteins, and a novel gene sequence coding for a putative multicopper oxidase (Mco). Results from gene expression studies, including reverse transcriptase (RT) quantitative PCR (qPCR) of cultures grown autotrophically on either Fe(II), pyrite, or elemental S showed that the fox gene cluster and mco are highly expressed under conditions where Fe(II) is an electron donor. Metagenome sequence and gene expression studies of Fe-oxide mats confirmed the importance of fox genes (e.g., foxA and foxC) and mco under Fe(II)-oxidizing conditions. Protein modeling of FoxC suggests a novel lysine-lysine or lysine-arginine heme B binding domain, indicating that it is likely the cytochrome component of a heterodimer complex with foxG as a ferredoxin subunit. Analysis of mco shows that it encodes a novel multicopper blue protein with two plastocyanin type I copper domains that may play a role in the transfer of electrons within the Fox protein complex. An understanding of metabolic pathways involved in aerobic iron and sulfur oxidation in Sulfolobales has broad implications for understanding the evolution and niche diversification of these thermophiles as well as practical applications in fields such as bioleaching of trace metals from pyritic ores.

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.