Pyrobaculum yellowstonensis Strain WP30 Respires on Elemental Sulfur and/or Arsenate in Circumneutral Sulfidic Geothermal Sediments of Yellowstone National Park

Tectonic and geological setting influence hot spring microbiology

Hydrothermal systems form at divergent and convergent boundaries of lithospheric plates and within plates due to weakened crust and mantle plumes, playing host to diverse microbial ecosystems. Little is known of how differences in tectonic setting influence the geochemical and microbial compositions of these hydrothermal ecosystems. Here, coordinated geochemical and microbial community analyses were conducted on 87 high-temperature (>65°C) water and sediment samples from hot springs in Yellowstone National Park, Wyoming, USA (n = 41; mantle plume setting), Iceland (n = 41, divergent boundary), and Japan (n = 5; convergent boundary). Region-specific variation in geochemistry and sediment-associated 16S rRNA gene amplicon sequence variant (ASV) composition was observed, with 16S rRNA gene assemblages being nearly completely distinguished by region and pH being the most explanatory parameter within regions. Several low abundance ASVs exhibited cosmopolitan distributions across regions, while most high-abundance ASVs were only identified in specific regions. The presence of some cosmopolitan ASVs across regions argues against dispersal limitation primarily shaping the distribution of taxa among regions. Rather, the results point to local tectonic and geologic characteristics shaping the geochemistry of continental hydrothermal systems that then select for distinct microbial assemblages. These results provide new insights into the co-evolution of hydrothermal systems and their microbial communities.

Evolution of the cytochrome bd oxygen reductase superfamily and the function of CydAA' in Archaea

Cytochrome bd-type oxygen reductases (cytbd) belong to one of three enzyme superfamilies that catalyze oxygen reduction to water. They are widely distributed in Bacteria and Archaea, but the full extent of their biochemical diversity is unknown. Here we used phylogenomics to identify three families and several subfamilies within the cytbd superfamily. The core architecture shared by all members of the superfamily consists of four transmembrane helices that bind two active site hemes, which are responsible for oxygen reduction. While previously characterized cytochrome bd-type oxygen reductases use quinol as an electron donor to reduce oxygen, sequence analysis shows that only one of the identified families has a conserved quinol binding site. The other families are missing this feature, suggesting that they use an alternative electron donor. Multiple gene duplication events were identified within the superfamily, resulting in significant evolutionary and structural diversity. The CydAA' cytbd, found exclusively in Archaea, is formed by the co-association of two superfamily paralogs. We heterologously expressed CydAA' from Caldivirga maquilingensis and demonstrated that it performs oxygen reduction with quinol as an electron donor. Strikingly, CydAA' is the first isoform of cytbd containing only b-type hemes shown to be active when isolated from membranes, demonstrating that oxygen reductase activity in this superfamily is not dependent on heme d.

Ecological Dichotomies Arise in Microbial Communities Due to Mixing of Deep Hydrothermal Waters and Atmospheric Gas in a Circumneutral Hot Spring

Little is known of how the confluence of subsurface and surface processes influences the assembly and habitability of hydrothermal ecosystems. To address this knowledge gap, the geochemical and microbial composition of a high-temperature, circumneutral hot spring in Yellowstone National Park was examined to identify the sources of solutes and their effect on the ecology of microbial inhabitants. Metagenomic analysis showed that populations comprising planktonic and sediment communities are archaeal dominated, are dependent on chemical energy (chemosynthetic), share little overlap in their taxonomic composition, and are differentiated by their inferred use of/tolerance to oxygen and mode of carbon metabolism. The planktonic community is dominated by putative aerobic/aerotolerant autotrophs, while the taxonomic composition of the sediment community is more evenly distributed and comprised of anaerobic heterotrophs. These observations are interpreted to reflect sourcing of the spring by anoxic, organic carbon-limited subsurface hydrothermal fluids and ingassing of atmospheric oxygen that selects for aerobic/aerotolerant organisms that have autotrophic capabilities in the water column. Autotrophy and consumption of oxygen by the planktonic community may influence the assembly of the anaerobic and heterotrophic sediment community. Support for this inference comes from higher estimated rates of genome replication in planktonic populations than sediment populations, indicating faster growth in planktonic populations. Collectively, these observations provide new insight into how mixing of subsurface waters and atmospheric oxygen create dichotomy in the ecology of hot spring communities and suggest that planktonic and sediment communities may have been less differentiated taxonomically and functionally prior to the rise of oxygen at ∼2.4 billion years ago (Gya). IMPORTANCE Understanding the source and availability of energy capable of supporting life in hydrothermal environments is central to predicting the ecology of microbial life on early Earth when volcanic activity was more widespread. Little is known of the substrates supporting microbial life in circumneutral to alkaline springs, despite their relevance to early Earth habitats. Using metagenomic and informatics approaches, water column and sediment habitats in a representative circumneutral hot spring in Yellowstone were shown to be dichotomous, with the former largely hosting aerobic/aerotolerant autotrophs and the latter primarily hosting anaerobic heterotrophs. This dichotomy is attributed to influx of atmospheric oxygen into anoxic deep hydrothermal spring waters. These results indicate that the ecology of microorganisms in circumneutral alkaline springs sourced by deep hydrothermal fluids was different prior to the rise of atmospheric oxygen ∼2.4 Gya, with planktonic and sediment communities likely to be less differentiated than contemporary circumneutral hot springs.

Deciphering Symbiotic Interactions of "Candidatus Aenigmarchaeota" with Inferred Horizontal Gene Transfers and Co-occurrence Networks

"Candidatus Aenigmarchaeota" ("Ca. Aenigmarchaeota") represents one of the earliest proposed evolutionary branches within the Diapherotrites, Parvarchaeota, Aenigmarchaeota, Nanoarchaeota, and Nanohaloarchaeota (DPANN) superphylum. However, their ecological roles and potential host-symbiont interactions are still poorly understood. Here, eight metagenome-assembled genomes (MAGs) were reconstructed from hot spring ecosystems, and further in-depth comparative and evolutionary genomic analyses were conducted on these MAGs and other genomes downloaded from public databases. Although with limited metabolic capacities, we reported that "Ca. Aenigmarchaeota" in thermal environments harbor more genes related to carbohydrate metabolism than "Ca. Aenigmarchaeota" in nonthermal environments. Evolutionary analyses suggested that members from the Thaumarchaeota, Aigarchaeota, Crenarchaeota, and Korarchaeota (TACK) superphylum and Euryarchaeota contribute substantially to the niche expansion of "Ca. Aenigmarchaeota" via horizontal gene transfer (HGT), especially genes related to virus defense and stress responses. Based on co-occurrence network results and recent genetic exchanges among community members, we conjectured that "Ca. Aenigmarchaeota" may be symbionts associated with one MAG affiliated with the genus Pyrobaculum, though host specificity might be wide and variable across different "Ca. Aenigmarchaeota" organisms. This study provides significant insight into possible DPANN-host interactions and ecological roles of "Ca. Aenigmarchaeota." IMPORTANCE Recent advances in sequencing technology promoted the blowout discovery of super tiny microbes in the Diapherotrites, Parvarchaeota, Aenigmarchaeota, Nanoarchaeota, and Nanohaloarchaeota (DPANN) superphylum. However, the unculturable properties of the majority of microbes impeded our investigation of their behavior and symbiotic lifestyle in the corresponding community. By integrating horizontal gene transfer (HGT) detection and co-occurrence network analysis on "Candidatus Aenigmarchaeota" ("Ca. Aenigmarchaeota"), we made one of the first attempts to infer their putative interaction partners and further decipher the potential functional and genetic interactions between the symbionts. We revealed that HGTs contributed by members from the Thaumarchaeota, Aigarchaeota, Crenarchaeota, and Korarchaeota (TACK) superphylum and Euryarchaeota conferred "Ca. Aenigmarchaeota" with the ability to survive under different environmental stresses, such as virus infection, high temperature, and oxidative stress. This study demonstrates that the interaction partners might be inferable by applying informatics analyses on metagenomic sequencing data.

Novel hyperthermophilic crenarchaeon Infirmifilum lucidum gen. nov. sp. nov., reclassification of Thermofilum uzonense as Infirmifilum uzonense comb. nov. and assignment of the family Thermofilaceae to the order Thermofilales ord. nov

A novel hyperthermophilic crenarchaeon, strain 3507LT^T, was isolated from a terrestrial hot spring near Tinguiririca volcano, Chile. Cells were non-motile thin, slightly curved filamentous rods. It grew at 73-93 °C and pH range of 5 to 7.5 with an optimum at 85 °C and pH 6.0-6.7. The presence of culture broth filtrate of another hyperthemophilic archaeon as well as yeast extract was obligatory for growth of the novel isolate. Strain 3507LT^T is an anaerobic chemoorganoheterotroph, fermenting monosaccharides, disaccharides and polysaccharides (lichenan, starch, xanthan gum, xyloglucan, alpha-cellulose and amorphous cellulose). No growth stimulation was detected when nitrate, thiosulfate, selenate or elemental sulfur were added as the electron acceptors. The complete genome of strain 3507LT^T consisted of a single circular chromosome with size of 1.63 Mbp. The DNA G+C content was 53.9%. According to the 16S rRNA gene sequence as well as conserved protein sequences phylogenetic analyses, strain 3507LT^T together with Thermofilum uzonense formed a separate cluster within a Thermofilaceae family (Thermoproteales/Thermoprotei/Crenarchaeota). Based on phenotypic characteristics, phylogeny as well as AAI comparisons, a novel genus and species Infirmifilum lucidum strain 3507LT^T (=VKM B-3376^T = KCTC 15797^T) gen. nov. sp. nov. was proposed. Its closest relative, Thermofilum uzonense strain 1807-2^T should be reclassified as Infirmifilum uzonense strain 1807-2^T comb. nov. Finally, based on phylogenomic and comparative genome analyses of representatives of Thermofilaceae family and other representatives of Thermoproteales order, a proposal of transfer of the family Thermofilaceae into a separate order Thermofilales ord. nov. was made.

Microbe-Mineral Interaction and Novel Proteins for Iron Oxide Mineral Reduction in the Hyperthermophilic Crenarchaeon Pyrodictium delaneyi

Understanding iron reduction in the hyperthermophilic crenarchaeon Pyrodictium delaneyi provides insight into the diversity of mechanisms used for this process and its potential impact in geothermal environments. The ability of P. delaneyi to reduce Fe(III) oxide minerals through direct contact potentially using a novel cytochrome respiratory complex and a membrane-bound molybdopterin respiratory complex sets iron reduction in this organism apart from previously described iron reduction processes.

Dissimilatory iron reduction by hyperthermophilic archaea occurs in many geothermal environments and generally relies on microbe-mineral interactions that transform various iron oxide minerals. In this study, the physiology of dissimilatory iron and nitrate reduction was examined in the hyperthermophilic crenarchaeon type strain Pyrodictium delaneyi Su06. Iron barrier experiments showed that P. delaneyi required direct contact with the Fe(III) oxide mineral ferrihydrite for reduction. The separate addition of an exogenous electron shuttle (anthraquinone-2,6-disulfonate), a metal chelator (nitrilotriacetic acid), and 75% spent cell-free supernatant did not stimulate growth with or without the barrier. Protein electrophoresis showed that the c-type cytochrome and general protein compositions of P. delaneyi changed when grown on ferrihydrite relative to nitrate. Differential proteomic analyses using tandem mass tagged protein fragments and mass spectrometry detected 660 proteins and differential production of 127 proteins. Among these, two putative membrane-bound molybdopterin-dependent oxidoreductase complexes increased in relative abundance 60- to 3,000-fold and 50- to 100-fold in cells grown on iron oxide. A putative 8-heme c-type cytochrome was 60-fold more abundant in iron-grown cells and was unique to the Pyrodictiaceae. There was also a >14,700-fold increase in a membrane transport protein in iron-grown cells. For flagellin proteins and a putative nitrate reductase, there were no changes in abundance, but a membrane nitric oxide reductase was more abundant on nitrate. These data help to elucidate the mechanisms by which hyperthermophilic crenarchaea generate energy and transfer electrons across the membrane to iron oxide minerals.

IMPORTANCE Understanding iron reduction in the hyperthermophilic crenarchaeon Pyrodictium delaneyi provides insight into the diversity of mechanisms used for this process and its potential impact in geothermal environments. The ability of P. delaneyi to reduce Fe(III) oxide minerals through direct contact potentially using a novel cytochrome respiratory complex and a membrane-bound molybdopterin respiratory complex sets iron reduction in this organism apart from previously described iron reduction processes. A model is presented where obligatory H₂ oxidation on the membrane coupled with electron transport and either Fe(III) oxide or nitrate reduction leads to the generation of a proton motive force and energy generation by oxidative phosphorylation. However, P. delaneyi cannot fix CO₂ and relies on organic compounds from its environment for biosynthesis.

Activity-based cell sorting reveals responses of uncultured archaea and bacteria to substrate amendment

Metagenomic studies have revolutionized our understanding of the metabolic potential of uncultured microorganisms in various ecosystems. However, many of these genomic predictions have yet to be experimentally tested, and the functional expression of genomic potential often remains unaddressed. In order to obtain a more thorough understanding of cell physiology, novel techniques capable of testing microbial metabolism under close to in situ conditions must be developed. Here, we provide a benchmark study to demonstrate that bioorthogonal non-canonical amino acid tagging (BONCAT) in combination with fluorescence-activated cell sorting (FACS) and 16S rRNA gene sequencing can be used to identify anabolically active members of a microbial community incubated in the presence of various growth substrates or under changing physicochemical conditions. We applied this approach to a hot spring sediment microbiome from Yellowstone National Park (Wyoming, USA) and identified several microbes that changed their activity levels in response to substrate addition, including uncultured members of the phyla Thaumarchaeota, Acidobacteria, and Fervidibacteria. Because shifts in activity in response to substrate amendment or headspace changes are indicative of microbial preferences for particular growth conditions, results from this and future BONCAT-FACS studies could inform the development of cultivation media to specifically enrich uncultured microbes. Most importantly, BONCAT-FACS is capable of providing information on the physiology of uncultured organisms at as close to in situ conditions as experimentally possible.

Arsenic-resistance mechanisms in bacterium Leclercia adecarboxylata strain As3-1: Biochemical and genomic analyses

Microbial arsenic transformation is important in As biogeochemical cycles in the environment. In this study, a new As-resistant bacterial strain Leclercia adecarboxylata As3-1 was isolated and its associated mechanisms in As resistance and detoxification were evaluated based on genome sequencing and gene annotations. After subjecting strain As3-1 to medium containing arsenate (AsV), AsV reduction occurred and an AsV-enhanced bacterial growth was observed. Strain As3-1 lacked arsenite (AsIII) oxidation ability and displayed lower AsIII resistance than AsV, probably due to its higher AsIII accumulation. Polymerase chain reaction and phylogenetic analysis showed that strain As3-1 harbored a typical AsV reductase gene (arsC) on the plasmids. Genome sequencing and gene annotations identified four operons phoUpstBACS, arsHRBC, arsCRDABC and ttrRSBCA, with 8 additional genes outside the operons that might have involved in As resistance and detoxification in strain As3-1. These included 5 arsC genes explaining why strain As3-1 tolerated high AsV concentrations. Besides ArsC, TtrB, TtrC and TtrA proteins could also be involved in AsV reduction and consequent energy acquisition for bacterial growth. Our data provided a new example of diverse As-regulating systems and AsV-enhanced growth without ArrA in bacteria. The information helps to understand the role of As in selecting microbial systems that can transform and utilize As.

Subsurface processes influence oxidant availability and chemoautotrophic hydrogen metabolism in Yellowstone hot springs

The geochemistry of hot springs and the availability of oxidants capable of supporting microbial metabolisms are influenced by subsurface processes including the separation of hydrothermal fluids into vapor and liquid phases. Here, we characterized the influence of geochemical variation and oxidant availability on the abundance, composition, and activity of hydrogen (H₂ )-dependent chemoautotrophs along the outflow channels of two-paired hot springs in Yellowstone National Park. The hydrothermal fluid at Roadside East (RSE; 82.4°C, pH 3.0) is acidic due to vapor-phase input while the fluid at Roadside West (RSW; 68.1°C, pH 7.0) is circumneutral due to liquid-phase input. Most chemotrophic communities exhibited net rates of H₂ oxidation, consistent with H₂ support of primary productivity, with one chemotrophic community exhibiting a net rate of H₂ production. Abundant H₂ -oxidizing chemoautotrophs were supported by reduction in oxygen, elemental sulfur, sulfate, and nitrate in RSW and oxygen and ferric iron in RSE; O₂ utilizing hydrogenotrophs increased in abundance down both outflow channels. Sequencing of 16S rRNA transcripts or genes from native sediments and dilution series incubations, respectively, suggests that members of the archaeal orders Sulfolobales, Desulfurococcales, and Thermoproteales are likely responsible for H₂ oxidation in RSE, whereas members of the bacterial order Thermoflexales and the archaeal order Thermoproteales are likely responsible for H₂ oxidation in RSW. These observations suggest that subsurface processes strongly influence spring chemistry and oxidant availability, which in turn select for unique assemblages of H₂ oxidizing microorganisms. Therefore, these data point to the role of oxidant availability in shaping the ecology and evolution of hydrogenotrophic organisms.

Archaeal S-Layers: Overview and Current State of the Art

In contrast to bacteria, all archaea possess cell walls lacking peptidoglycan and a number of different cell envelope components have also been described. A paracrystalline protein surface layer, commonly referred to as S-layer, is present in nearly all archaea described to date. S-layers are composed of only one or two proteins and form different lattice structures. In this review, we summarize current understanding of archaeal S-layer proteins, discussing topics such as structure, lattice type distribution among archaeal phyla and glycosylation. The hexagonal lattice type is dominant within the phylum Euryarchaeota, while in the Crenarchaeota this feature is mainly associated with specific orders. S-layers exclusive to the Crenarchaeota have also been described, which are composed of two proteins. Information regarding S-layers in the remaining archaeal phyla is limited, mainly due to organism description through only culture-independent methods. Despite the numerous applied studies using bacterial S-layers, few reports have employed archaea as a study model. As such, archaeal S-layers represent an area for exploration in both basic and applied research.

A novel bifunctional transcriptional regulator of riboflavin metabolism in Archaea

Riboflavin (vitamin B2) is the precursor of flavin mononucleotide (FMN) and flavin adenine dinucleotide, which are essential coenzymes in all free-living organisms. Riboflavin biosynthesis in many Bacteria but not in Archaea is controlled by FMN-responsive riboswitches. We identified a novel bifunctional riboflavin kinase/regulator (RbkR), which controls riboflavin biosynthesis and transport genes in major lineages of Crenarchaeota, Euryarchaeota and Thaumarchaeota. RbkR proteins are composed of the riboflavin kinase domain and a DNA-binding winged helix-turn-helix-like domain. Using comparative genomics, we predicted RbkR operator sites and reconstructed RbkR regulons in 94 archaeal genomes. While the identified RbkR operators showed significant variability between archaeal lineages, the conserved core of RbkR regulons includes riboflavin biosynthesis genes, known/predicted vitamin uptake transporters and the rbkR gene. The DNA motifs and CTP-dependent riboflavin kinase activity of two RbkR proteins were experimentally validated in vitro. The DNA binding activity of RbkR was stimulated by CTP and suppressed by FMN, a product of riboflavin kinase. The crystallographic structure of RbkR from Thermoplasma acidophilum was determined in complex with CTP and its DNA operator revealing key residues for operator and ligand recognition. Overall, this study contributes to our understanding of metabolic and regulatory networks for vitamin homeostasis in Archaea.

Integration of Metagenomic and Stable Carbon Isotope Evidence Reveals the Extent and Mechanisms of Carbon Dioxide Fixation in High-Temperature Microbial Communities

Although the biological fixation of CO₂ by chemolithoautotrophs provides a diverse suite of organic compounds utilized by chemoorganoheterotrophs as a carbon and energy source, the relative amounts of autotrophic C in chemotrophic microbial communities are not well-established. The extent and mechanisms of CO₂ fixation were evaluated across a comprehensive set of high-temperature, chemotrophic microbial communities in Yellowstone National Park by combining metagenomic and stable ¹³C isotope analyses. Fifteen geothermal sites representing three distinct habitat types (iron-oxide mats, anoxic sulfur sediments, and filamentous “streamer” communities) were investigated. Genes of the 3-hydroxypropionate/4-hydroxybutyrate, dicarboxylate/4-hydroxybutyrate, and reverse tricarboxylic acid CO₂ fixation pathways were identified in assembled genome sequence corresponding to the predominant Crenarchaeota and Aquificales observed across this habitat range. Stable ¹³C analyses of dissolved inorganic and organic C (DIC, DOC), and possible landscape C sources were used to interpret the ¹³C content of microbial community samples. Isotope mixing models showed that the minimum fractions of autotrophic C in microbial biomass were >50% in the majority of communities analyzed. The significance of CO₂ as a C source in these communities provides a foundation for understanding community assembly and succession, and metabolic linkages among early-branching thermophilic autotrophs and heterotrophs.

The distribution, diversity and function of predominant Thermoproteales in high-temperature environments of Yellowstone National Park

High-temperature environments (> 70°C) contain diverse and abundant members of the crenarchaeal order Thermoproteales. However, a comprehensive study of the distribution and function of diverse members of this group across different habitat types has not been conducted. Consequently, the goals of this study were to determine the distribution of different Thermoproteales genera across geochemically distinct geothermal habitats of Yellowstone National Park, and to identify key functional attributes of major genera that correlate with environmental parameters. Curated sequence assemblies belonging to five genera were characterized in replicate samples of 11 high-temperature communities ranging in pH from 3 to 9. Thermocladium, Vulcanisaeta and Caldivirga spp. were the primary Thermoproteales populations present in low pH (pH < 5) habitats, whereas Thermoproteus populations were found in mildly-acidic (pH 5-6) sulfur sediments, and Pyrobaculum populations were confined to higher pH (pH > 6) sulfur sediments and/or filamentous 'streamer' communities. Metabolic reconstruction and comparative genomics among assemblies show that these populations are primarily chemoorganotrophs that utilize different electron acceptors depending on geochemical conditions. The presence of potential CO₂ fixation pathways in some Thermoproteales populations appears to be linked with NiFe hydrogenases, which combined with high levels of H₂ in many sulfidic systems, may provide the energy required to fix inorganic C.

Novel, Deep-Branching Heterotrophic Bacterial Populations Recovered from Thermal Spring Metagenomes

Thermal spring ecosystems are a valuable resource for the discovery of novel hyperthermophilic Bacteria and Archaea, and harbor deeply-branching lineages that provide insight regarding the nature of early microbial life. We characterized bacterial populations in two circumneutral (pH ~8) Yellowstone National Park thermal (T ~80°C) spring filamentous “streamer” communities using random metagenomic DNA sequence to investigate the metabolic potential of these novel populations. Four de novo assemblies representing three abundant, deeply-branching bacterial phylotypes were recovered. Analysis of conserved phylogenetic marker genes indicated that two of the phylotypes represent separate groups of an uncharacterized phylum (for which we propose the candidate phylum name “Pyropristinus”). The third new phylotype falls within the proposed Calescamantes phylum. Metabolic reconstructions of the “Pyropristinus” and Calescamantes populations showed that these organisms appear to be chemoorganoheterotrophs and have the genomic potential for aerobic respiration and oxidative phosphorylation via archaeal-like V-type, and bacterial F-type ATPases, respectively. A survey of similar phylotypes (>97% nt identity) within 16S rRNA gene datasets suggest that the newly described organisms are restricted to terrestrial thermal springs ranging from 70 to 90°C and pH values of ~7–9. The characterization of these lineages is important for understanding the diversity of deeply-branching bacterial phyla, and their functional role in high-temperature circumneutral “streamer” communities.

Geomicrobiology of sublacustrine thermal vents in Yellowstone Lake: geochemical controls on microbial community structure and function

Yellowstone Lake (Yellowstone National Park, WY, USA) is a large high-altitude (2200 m), fresh-water lake, which straddles an extensive caldera and is the center of significant geothermal activity. The primary goal of this interdisciplinary study was to evaluate the microbial populations inhabiting thermal vent communities in Yellowstone Lake using 16S rRNA gene and random metagenome sequencing, and to determine how geochemical attributes of vent waters influence the distribution of specific microorganisms and their metabolic potential. Thermal vent waters and associated microbial biomass were sampled during two field seasons (2007–2008) using a remotely operated vehicle (ROV). Sublacustrine thermal vent waters (circa 50–90°C) contained elevated concentrations of numerous constituents associated with geothermal activity including dissolved hydrogen, sulfide, methane and carbon dioxide. Microorganisms associated with sulfur-rich filamentous “streamer” communities of Inflated Plain and West Thumb (pH range 5–6) were dominated by bacteria from the Aquificales, but also contained thermophilic archaea from the Crenarchaeota and Euryarchaeota. Novel groups of methanogens and members of the Korarchaeota were observed in vents from West Thumb and Elliot's Crater (pH 5–6). Conversely, metagenome sequence from Mary Bay vent sediments did not yield large assemblies, and contained diverse thermophilic and nonthermophilic bacterial relatives. Analysis of functional genes associated with the major vent populations indicated a direct linkage to high concentrations of carbon dioxide, reduced sulfur (sulfide and/or elemental S), hydrogen and methane in the deep thermal ecosystems. Our observations show that sublacustrine thermal vents in Yellowstone Lake support novel thermophilic communities, which contain microorganisms with functional attributes not found to date in terrestrial geothermal systems of YNP.

The distribution, diversity, and importance of 16S rRNA gene introns in the order Thermoproteales

Background

Intron sequences are common in 16S rRNA genes of specific thermophilic lineages of Archaea, specifically the Thermoproteales (phylum Crenarchaeota). Environmental sequencing (16S rRNA gene and metagenome) from geothermal habitats in Yellowstone National Park (YNP) has expanded the available datasets for investigating 16S rRNA gene introns. The objectives of this study were to characterize and curate archaeal 16S rRNA gene introns from high-temperature habitats, evaluate the conservation and distribution of archaeal 16S rRNA introns in geothermal systems, and determine which “universal” archaeal 16S rRNA gene primers are impacted by the presence of intron sequences.

Results

Several new introns were identified and their insertion loci were constrained to thirteen locations across the 16S rRNA gene. Many of these introns encode homing endonucleases, although some introns were short or partial sequences. Pyrobaculum, Thermoproteus, and Caldivirga 16S rRNA genes contained the most abundant and diverse intron sequences. Phylogenetic analysis of introns revealed that sequences within the same locus are distributed biogeographically. The most diverse set of introns were observed in a high-temperature, circumneutral (pH 6) sulfur sediment environment, which also contained the greatest diversity of different Thermoproteales phylotypes.

Conclusions

The widespread presence of introns in the Thermoproteales indicates a high probability of misalignments using different “universal” 16S rRNA primers employed in environmental microbial community analysis.

Reviewers

This article was reviewed by Dr. Eugene Koonin and Dr. W. Ford Doolittle.

Electronic supplementary material

The online version of this article (doi:10.1186/s13062-015-0065-6) contains supplementary material, which is available to authorized users.