Comparative genomic analysis of phylogenetically closely related Hydrogenobaculum sp. isolates from Yellowstone National Park

Genome sequence and characterisation of a freshwater photoarsenotroph, Cereibacter azotoformans strain ORIO, isolated from sediments capable of cyclic light-dark arsenic oxidation and reduction

A freshwater photosynthetic arsenite-oxidizing bacterium, Cereibacter azotoformans strain ORIO, was isolated from Owens River, CA, USA. The waters from Owens River are elevated in arsenic and serve as the headwaters to the Los Angeles Aqueduct. The complete genome sequence of strain ORIO is 4.8 Mb genome (68% G + C content) and comprises two chromosomes and six plasmids. Taxonomic analysis placed ORIO within the Cereibacter genus (formerly Rhodobacter). The ORIO genome contains arxB2 AB1 CD (encoding an arsenite oxidase), arxXSR (regulators) and several ars arsenic resistance genes all co-localised on a 136 kb plasmid, named pORIO3. Phylogenetic analysis of ArxA, the molybdenum-containing arsenite oxidase catalytic subunit, demonstrated photoarsenotrophy is likely to occur within members of the Alphaproteobacteria. ORIO is a mixotroph, oxidises arsenite to arsenate (As(V)) photoheterotrophically, and expresses arxA in cultures grown with arsenite. Further ecophysiology studies with Owens River sediment demonstrated the interconversion of arsenite and As(V) was dependent on light-dark cycling. arxA and arrA (As(V) respiratory reductase) genes were detected in the light-dark cycled sediment metagenomes suggesting syntrophic interactions among arsenotrophs. This work establishes C. azotoformans str. ORIO as a new model organism for studying photoarsenotrophy and light-dark arsenic biogeochemical cycling.

Effects of hydrogeochemistry on the microbial ecology of terrestrial hot springs

Temperature, pH, and hydrochemistry of terrestrial hot springs play a critical role in shaping thermal microbial communities. However, the interactions of biotic and abiotic factors at this terrestrial-aquatic interface are still not well understood on a global scale, and the question of how underground events influence microbial communities remains open. To answer this, 11 new samples obtained from the El Tatio geothermal field were analyzed by 16S rRNA amplicon sequencing (V4 region), along with 191 samples from previous publications obtained from the Taupo Volcanic Zone, the Yellowstone Plateau Volcanic Field, and the Eastern Tibetan Plateau, with their temperature, pH, and major ion concentration. Microbial alpha diversity was lower in acid-sulfate waters, and no significant correlations were found with temperature. However, moderate correlations were observed between chemical parameters such as pH (mostly constrained to temperatures below 70°C), SO4 2- and abundances of members of the phyla Armatimonadota, Deinococcota, Chloroflexota, Campilobacterota, and Thermoplasmatota. pH and SO4 2- gradients were explained by phase separation of sulfur-rich hydrothermal fluids and oxidation of reduced sulfur in the steam phase, which were identified as key processes shaping these communities. Ordination and permutational analysis of variance showed that temperature, pH, and major element hydrochemistry explain only 24% of the microbial community structure. Therefore, most of the variance remained unexplained, suggesting that other environmental or biotic factors are also involved and highlighting the environmental complexity of the ecosystem and its great potential to test niche theory ecological associated questions. IMPORTANCE This is the first approach to investigate whether geothermal processes could have an influence on the ecology of thermal microbial communities on a global scale. In addition to temperature and pH, microbial communities are structured by sulfate concentrations, which depends on the tectono-magmatic settings (such as the depth of magmatic chambers) and the local settings (such as the availability of a confining layer separating NaCl waters from steam after phase separation) and the possibility of mixing with more diluted fluids. Comparison of microbial communities from different geothermal areas by homogeneous sequence processing showed that no significant geographic distance decay was detected on the microbial communities according to Bray-Curtis, Jaccard, unweighted, and weighted Unifrac similarity/dissimilarity indices. Instead, an ancient potential divergence in the same taxonomic groups is suggested between globally distant thermal zones.

Seasonal hydrologic and geologic forcing drive hot spring geochemistry and microbial biodiversity

Hot springs integrate hydrologic and geologic processes that vary over short- and long-term time scales. However, the influence of temporal hydrologic and geologic change on hot spring biodiversity is unknown. Here, we coordinated near-weekly, cross-seasonal (~140 days) geochemical and microbial community analyses of three widely studied hot springs with local precipitation data in Yellowstone National Park. One spring ('HFS') exhibited statistically significant, coupled microbial and geochemical variation across seasons that was associated with recent precipitation patterns. Two other spring communities, 'CP' and 'DS', exhibited minimal to no variation across seasons. Variability in the seasonal response of springs is attributed to differences in the timing and extent of aquifer recharge with oxidized near-surface water from precipitation. This influx of oxidized water is associated with changes in community composition, and in particular, the abundances of aerobic sulfide-/sulfur-oxidizers that can acidify waters. During sampling, a new spring formed after a period of heavy precipitation and its successional dynamics were also influenced by surface water recharge. Collectively, these results indicate that changes in short-term hydrology associated with precipitation can impact hot spring geochemistry and microbial biodiversity. These results point to potential susceptibility of certain hot springs and their biodiversity to sustained, longer-term hydrologic changes.

Diversity and Distribution of a Novel Genus of Hyperthermophilic Aquificae Viruses Encoding a Proof-Reading Family-A DNA Polymerase

Despite the high abundance of Aquificae in many geothermal systems, these bacteria are difficult to culture and no viruses infecting members of this phylum have been isolated. Here, we describe the complete, circular dsDNA Uncultivated Virus Genome (UViG) of Thermocrinis Octopus Spring virus (TOSV), derived from metagenomic data, along with eight related UViGs representing three additional viral species. Despite low overall similarity among viruses from different hot springs, the genomes shared a high degree of synteny, and encoded numerous genes for nucleotide metabolism, including a PolA-type DNA polymerase polyprotein with likely accessory functions, a DNA Pol III sliding clamp, a thymidylate kinase, a DNA gyrase, a helicase, and a DNA methylase. Also present were conserved genes predicted to code for phage capsid, large and small subunits of terminase, portal protein, holin, and lytic transglycosylase, all consistent with a distant relatedness to cultivated Caudovirales. These viruses are predicted to infect Aquificae, as multiple CRISPR spacers matching the viral genomes were identified within the genomes and metagenomic contigs from these bacteria. Based on the predicted atypical bi-directional replication strategy, low sequence similarity to known viral genomes, and unique position in gene-sharing networks, we propose a new putative genus, "Pyrovirus," in the order Caudovirales.

Phylogenomic analysis of novel Diaforarchaea is consistent with sulfite but not sulfate reduction in volcanic environments on early Earth

The origin(s) of dissimilatory sulfate and/or (bi)sulfite reducing organisms (SRO) remains enigmatic despite their importance in global carbon and sulfur cycling since at least 3.4 Ga. Here, we describe novel, deep-branching archaeal SRO populations distantly related to other Diaforarchaea from two moderately acidic thermal springs. Dissimilatory (bi)sulfite reductase homologs, DsrABC, encoded in metagenome assembled genomes (MAGs) from spring sediments comprise one of the earliest evolving Dsr lineages. DsrA homologs were expressed in situ under moderately acidic conditions. MAGs lacked genes encoding proteins that activate sulfate prior to (bi)sulfite reduction. This is consistent with sulfide production in enrichment cultures provided sulfite but not sulfate. We suggest input of volcanic sulfur dioxide to anoxic spring-water yields (bi)sulfite and moderately acidic conditions that favor its stability and bioavailability. The presence of similar volcanic springs at the time SRO are thought to have originated (>3.4 Ga) may have supplied (bi)sulfite that supported ancestral SRO. These observations coincide with the lack of inferred SO42- reduction capacity in nearly all organisms with early-branching DsrAB and which are near universally found in hydrothermal environments.

Biogeochemical probing of microbial communities in a basalt-hosted hot spring at Kverkfjöll volcano, Iceland

We investigated bacterial and archaeal communities along an ice-fed surficial hot spring at Kverkfjöll volcano-a partially ice-covered basaltic volcano at Vatnajökull glacier, Iceland, using biomolecular (16S rRNA, apsA, mcrA, amoA, nifH genes) and stable isotope techniques. The hot spring environment is characterized by high temperatures and low dissolved oxygen concentrations at the source (68°C and <1 mg/L (±0.1%)) changing to lower temperatures and higher dissolved oxygen downstream (34.7°C and 5.9 mg/L), with sulfate the dominant anion (225 mg/L at the source). Sediments are comprised of detrital basalt, low-temperature alteration phases and pyrite, with <0.4 wt. % total organic carbon (TOC). 16S rRNA gene profiles reveal that organisms affiliated with Hydrogenobaculum (54%-87% bacterial population) and Thermoproteales (35%-63% archaeal population) dominate the micro-oxic hot spring source, while sulfur-oxidizing archaea (Sulfolobales, 57%-82%), and putative sulfur-oxidizing and heterotrophic bacterial groups dominate oxic downstream environments. The δ¹³ C_org (‰ V-PDB) values for sediment TOC and microbial biomass range from -9.4‰ at the spring's source decreasing to -12.6‰ downstream. A reverse effect isotope fractionation of ~3‰ between sediment sulfide (δ³⁴ S ~0‰) and dissolved water sulfate (δ³⁴ S +3.2‰), and δ¹⁸ O values of ~ -5.3‰ suggest pyrite forms abiogenically from volcanic sulfide, followed by abiogenic and microbial oxidation. These environments represent an unexplored surficial geothermal environment analogous to transient volcanogenic habitats during putative "snowball Earth" scenarios and volcano-ice geothermal environments on Mars.

MEBS, a software platform to evaluate large (meta)genomic collections according to their metabolic machinery: unraveling the sulfur cycle

The increasing number of metagenomic and genomic sequences has dramatically improved our understanding of microbial diversity, yet our ability to infer metabolic capabilities in such datasets remains challenging. We describe the Multigenomic Entropy Based Score pipeline (MEBS), a software platform designed to evaluate, compare, and infer complex metabolic pathways in large “omic” datasets, including entire biogeochemical cycles. MEBS is open source and available through https://github.com/eead-csic-compbio/metagenome_Pfam_score. To demonstrate its use, we modeled the sulfur cycle by exhaustively curating the molecular and ecological elements involved (compounds, genes, metabolic pathways, and microbial taxa). This information was reduced to a collection of 112 characteristic Pfam protein domains and a list of complete-sequenced sulfur genomes. Using the mathematical framework of relative entropy (H^΄), we quantitatively measured the enrichment of these domains among sulfur genomes. The entropy of each domain was used both to build up a final score that indicates whether a (meta)genomic sample contains the metabolic machinery of interest and to propose marker domains in metagenomic sequences such as DsrC (PF04358). MEBS was benchmarked with a dataset of 2107 non-redundant microbial genomes from RefSeq and 935 metagenomes from MG-RAST. Its performance, reproducibility, and robustness were evaluated using several approaches, including random sampling, linear regression models, receiver operator characteristic plots, and the area under the curve metric (AUC). Our results support the broad applicability of this algorithm to accurately classify (AUC = 0.985) hard-to-culture genomes (e.g., Candidatus Desulforudis audaxviator), previously characterized ones, and metagenomic environments such as hydrothermal vents, or deep-sea sediment. Our benchmark indicates that an entropy-based score can capture the metabolic machinery of interest and can be used to efficiently classify large genomic and metagenomic datasets, including uncultivated/unexplored taxa.

Microbial community dynamics in Inferno Crater Lake, a thermally fluctuating geothermal spring

Understanding how microbial communities respond and adjust to ecosystem perturbation is often difficult to interpret due to multiple and often simultaneous variations in observed conditions. In this research, we investigated the microbial community dynamics of Inferno Crater Lake, an acidic geothermal spring in New Zealand with a unique thermal cycle that varies between 30 and 80 °C over a period of 40-60 days. Using a combination of next-generation sequencing, geochemical analysis and quantitative PCR we found that the microbial community composition was predominantly chemolithotrophic and strongly associated with the thermal cycle. At temperatures >65 °C, the microbial community was dominated almost exclusively by sulphur-oxidising archaea (Sulfolobus-like spp.). By contrast, at mesophilic temperatures the community structure was more mixed, comprising both archaea and bacteria but dominated primarily by chemolithotrophic sulphur and hydrogen oxidisers. Multivariate analysis of physicochemical data confirmed that temperature was the only significant variable associated with community turnover. This research contributes to our understanding of microbial community dynamics in variable environments, using a naturally alternating system as a model and extends our limited knowledge of acidophile ecology in geothermal habitats.

Microbes in mercury-enriched geothermal springs in western North America

Because geothermal environments contain mercury (Hg) from natural sources, microorganisms that evolved in these systems have likely adapted to this element. Knowledge of the interactions between microorganisms and Hg in geothermal systems may assist in understanding the long-term evolution of microbial adaptation to Hg with relevance to other environments where Hg is introduced from anthropogenic sources. A number of microbiological studies with supporting geochemistry have been conducted in geothermal systems across western North America. Approximately 1 in 5 study sites include measurements of Hg. Of all prokaryotic taxa reported across sites with microbiological and accompanying physicochemical data, 42% have been detected at sites in which Hg was measured. Genes specifying Hg reduction and detoxification by microorganisms were detected in a number of hot springs across the region. Archaeal-like sequences, representing two crenarchaeal orders and one order each of the Euryarchaeota and Thaumarchaeota, dominated in metagenomes' MerA (the mercuric reductase protein) inventories, while bacterial homologs were mostly found in one deeply sequenced metagenome. MerA homologs were more frequently found in metagenomes of microbial communities in acidic springs than in circumneutral or high pH geothermal systems, possibly reflecting higher bioavailability of Hg under acidic conditions. MerA homologs were found in hot springs prokaryotic isolates affiliated with Bacteria and Archaea taxa. Acidic sites with high Hg concentrations contain more of Archaea than Bacteria taxa, while the reverse appears to be the case in circumneutral and high pH sites with high Hg concentrations. However, MerA was detected in only a small fraction of the Archaea and Bacteria taxa inhabiting sites containing Hg. Nevertheless, the presence of MerA homologs and their distribution patterns in systems, in which Hg has yet to be measured, demonstrates the potential for detoxification by Hg reduction in these geothermal systems, particularly the low pH springs that are dominated by Archaea.

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⁻¹, and in situ microelectrode measurements showed that the oxidation of Fe(II) is limited by the diffusion of O₂ into the Fe(III)-oxide mat. The formation of microterracettes also implicated O₂ 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.

Novel viral genomes identified from six metagenomes reveal wide distribution of archaeal viruses and high viral diversity in terrestrial hot springs

Limited by culture-dependent methods the number of viruses identified from thermophilic Archaea and Bacteria is still very small. In this study we retrieved viral sequences from six hot spring metagenomes isolated worldwide, revealing a wide distribution of four archaeal viral families, Ampullaviridae, Bicaudaviridae, Lipothrixviridae and Rudiviridae. Importantly, we identified 10 complete or near complete viral genomes allowing, for the first time, an assessment of genome conservation and evolution of the Ampullaviridae family as well as Sulfolobus Monocaudavirus 1 (SMV1)-related viruses. Among the novel genomes, one belongs to a putative thermophilic virus infecting the bacterium Hydrogenobaculum, for which no virus has been reported in the literature. Moreover, a high viral diversity was observed in the metagenomes, especially among the Lipothrixviridae, as indicated by the large number of unique contigs and the lack of a completely assembled genome for this family. This is further supported by the large number of novel genes in the complete and partial genomes showing no sequence similarities to public databases. CRISPR analysis revealed hundreds of novel CRISPR loci and thousands of novel CRISPR spacers from each metagenome, reinforcing the notion of high viral diversity in the thermal environment.

Isolation of diverse members of the Aquificales from geothermal springs in Tengchong, China

The order Aquificales (phylum Aquificae) consists of thermophilic and hyperthermophilic bacteria that are prominent in many geothermal systems, including those in Tengchong, Yunnan Province, China. However, Aquificales have not previously been isolated from Tengchong. We isolated five strains of Aquificales from diverse springs (temperature 45.2–83.3°C and pH 2.6–9.1) in the Rehai Geothermal Field from sites in which Aquificales were abundant. Phylogenetic analysis showed that four of the strains belong to the genera Hydrogenobacter, Hydrogenobaculum, and Sulfurihydrogenibium, including strains distant enough to likely justify new species of Hydrogenobacter and Hydrogenobaculum. The additional strain may represent a new genus in the Hydrogenothermaceae. All strains were capable of aerobic respiration under microaerophilic conditions; however, they had variable capacity for chemolithotrophic oxidation of hydrogen and sulfur compounds and nitrate reduction.

Uncovering rare NADH-preferring ketol-acid reductoisomerases

All members of the ketol-acid reductoisomerase (KARI) enzyme family characterized to date have been shown to prefer the nicotinamide adenine dinucleotide phosphate hydride (NADPH) cofactor to nicotinamide adenine dinucleotide hydride (NADH). However, KARIs with the reversed cofactor preference are desirable for industrial applications, including anaerobic fermentation to produce branched-chain amino acids. By applying insights gained from structural and engineering studies of this enzyme family to a comprehensive multiple sequence alignment of KARIs, we identified putative NADH-utilizing KARIs and characterized eight whose catalytic efficiencies using NADH were equal to or greater than NADPH. These are the first naturally NADH-preferring KARIs reported and demonstrate that this property has evolved independently multiple times, using strategies unlike those used previously in the laboratory to engineer a KARI cofactor switch.

Strain-level genomic variation in natural populations of Lebetimonas from an erupting deep-sea volcano

Chemolithoautotrophic Epsilonproteobacteria are ubiquitous in sulfidic, oxygen-poor habitats, including hydrothermal vents, marine oxygen minimum zones, marine sediments and sulfidic caves and have a significant role in cycling carbon, hydrogen, nitrogen and sulfur in these environments. The isolation of diverse strains of Epsilonproteobacteria and the sequencing of their genomes have revealed that this group has the metabolic potential to occupy a wide range of niches, particularly at dynamic deep-sea hydrothermal vents. We expand on this body of work by examining the population genomics of six strains of Lebetimonas, a vent-endemic, thermophilic, hydrogen-oxidizing Epsilonproteobacterium, from a single seamount in the Mariana Arc. Using Lebetimonas as a model for anaerobic, moderately thermophilic organisms in the warm, anoxic subseafloor environment, we show that genomic content is highly conserved and that recombination is limited between closely related strains. The Lebetimonas genomes are shaped by mobile genetic elements and gene loss as well as the acquisition of novel functional genes by horizontal gene transfer, which provide the potential for adaptation and microbial speciation in the deep sea. In addition, these Lebetimonas genomes contain two operons of nitrogenase genes with different evolutionary origins. Lebetimonas expressed nifH during growth with nitrogen gas as the sole nitrogen source, thus providing the first evidence of nitrogen fixation in any Epsilonproteobacteria from deep-sea hydrothermal vents. In this study, we provide a comparative overview of the genomic potential within the Nautiliaceae as well as among more distantly related hydrothermal vent Epsilonproteobacteria to broaden our understanding of microbial adaptation and diversity in the deep sea.

New insights into the evolutionary history of biological nitrogen fixation

Nitrogenase, which catalyzes the ATP-dependent reduction of dinitrogen (N₂) to ammonia (NH₃), accounts for roughly half of the bioavailable nitrogen supporting extant life. The fundamental requirement for fixed forms of nitrogen for life on Earth, both at present and in the past, has led to broad and significant interest in the origin and evolution of biological N₂ fixation. One key question is whether the limited availability of fixed nitrogen was a factor in life's origin or whether there were ample sources of fixed nitrogen produced by abiotic processes or delivered through the weathering of bolide impact materials to support this early life. If the latter, the key questions become what were the characteristics of the environment that precipitated the evolution of this oxygen sensitive process, when did this occur, and how was its subsequent evolutionary history impacted by the advent of oxygenic photosynthesis and the rise of oxygen in the Earth's biosphere. Since the availability of fixed sources of nitrogen capable of supporting early life is difficult to glean from the geologic record, there are limited means to get direct insights into these questions. Indirect insights, however, can be gained through phylogenetic studies of nitrogenase structural gene products and additional gene products involved in the biosynthesis of the complex metal-containing prosthetic groups associated with this enzyme complex. Insights gained from such studies, as reviewed herein, challenge traditional models for the evolution of biological nitrogen fixation and provide the basis for the development of new conceptual models that explain the stepwise evolution of this highly complex life sustaining process.

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.

Evidence for the widespread distribution of CRISPR-Cas system in the Phylum Cyanobacteria

Members of the phylum Cyanobacteria inhabit ecologically diverse environments. However, the CRISPR-Cas (clustered regularly interspaced short palindromic repeats, CRISPR associated genes), an extremely adaptable defense system, has not been surveyed in this phylum. We analyzed 126 cyanobacterial genomes and, surprisingly, found CRISPR-Cas in the majority except the marine subclade (Synechococcus and Prochlorococcus), in which cyanophages are a known force shaping their evolution. Multiple observations of CRISPR loci in the absence of cas1/cas2 genes may represent an early stage of losing a CRISPR-Cas locus. Our findings reveal the widespread distribution of their role in the phylum Cyanobacteria and provide a first step to systematically understanding CRISPR-Cas systems in cyanobacteria.