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

Tectonic settings influence the geochemical and microbial diversity of Peru hot springs

Tectonic processes control hot spring temperature and geochemistry, yet how this in turn shapes microbial community composition is poorly understood. Here, we present geochemical and 16 S rRNA gene sequencing data from 14 hot springs from contrasting styles of subduction along a convergent margin in the Peruvian Andes. We find that tectonic influence on hot spring temperature and geochemistry shapes microbial community composition. Hot springs in the flat-slab and back-arc regions of the subduction system had similar pH but differed in geochemistry and microbiology, with significant relationships between microbial community composition, geochemistry, and geologic setting. Flat-slab hot springs were chemically heterogeneous, had modest surface temperatures (up to 45 °C), and were dominated by members of the metabolically diverse phylum Proteobacteria. Whereas, back-arc hot springs were geochemically more homogenous, exhibited high concentrations of dissolved metals and gases, had higher surface temperatures (up to 81 °C), and host thermophilic archaeal and bacterial lineages.

Fungal Community Shift Along Steep Environmental Gradients from Geothermal Soils in Yellowstone National Park

Geothermal soils offer unique insight into the way extreme environmental factors shape communities of organisms. However, little is known about the fungi growing in these environments and in particular how localized steep abiotic gradients affect fungal diversity. We used metabarcoding to characterize soil fungi surrounding a hot spring-fed thermal creek with water up to 84 °C and pH 10 in Yellowstone National Park. We found a significant association between fungal communities and soil variable principal components, and we identify the key trends in co-varying soil variables that explain the variation in fungal community. Saprotrophic and ectomycorrhizal fungi community profiles followed, and were significantly associated with, different soil variable principal components, highlighting potential differences in the factors that structure these different fungal trophic guilds. In addition, in vitro growth experiments in four target fungal species revealed a wide range of tolerances to pH levels but not to heat. Overall, our results documenting turnover in fungal species within a few hundred meters suggest many co-varying environmental factors structure the diverse fungal communities found in the soils of Yellowstone National Park.

Microbial diversity in extreme environments

A wide array of microorganisms, including many novel, phylogenetically deeply rooted taxa, survive and thrive in extreme environments. These unique and reduced-complexity ecosystems offer a tremendous opportunity for studying the structure, function and evolution of natural microbial communities. Marker gene surveys have resolved patterns and ecological drivers of these extremophile assemblages, revealing a vast uncultured microbial diversity and the often predominance of archaea in the most extreme conditions. New omics studies have uncovered linkages between community function and environmental variables, and have enabled discovery and genomic characterization of major new lineages that substantially expand microbial diversity and change the structure of the tree of life. These efforts have significantly advanced our understanding of the diversity, ecology and evolution of microorganisms populating Earth's extreme environments, and have facilitated the exploration of microbiota and processes in more complex ecosystems.

Exploring the taxonomical and functional profile of As Burgas hot spring focusing on thermostable β-galactosidases

In the present study we investigate the microbial community inhabiting As Burgas geothermal spring, located in Ourense (Galicia, Spain). The approximately 23 Gbp of Illumina sequences generated for each replicate revealed a complex microbial community dominated by Bacteria in which Proteobacteria and Aquificae were the two prevalent phyla. An association between the two most prevalent genera, Thermus and Hydrogenobacter, was suggested by the relationship of their metabolism. The high relative abundance of sequences involved in the Calvin-Benson cycle and the reductive TCA cycle unveils the dominance of an autotrophic population. Important pathways from the nitrogen and sulfur cycle are potentially taking place in As Burgas hot spring. In the assembled reads, two complete ORFs matching GH2 beta-galactosidases were found. To assess their functional characterization, the two ORFs were cloned and overexpressed in E. coli. The pTsbg enzyme had activity towards o-Nitrophenyl-β-D-galactopyranoside (ONPG) and p-Nitrophenyl-β-D-fucopyranoside, with high thermal stability and showing maximal activity at 85 °C and pH 6, nevertheless the enzyme failed to hydrolyze lactose. The other enzyme, Tsbg, was unable to hydrolyze even ONPG or lactose. This finding highlights the challenge of finding novel active enzymes based only on their sequence.

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.

Relationship between Microorganisms Inhabiting Alkaline Siliceous Hot Spring Mat Communities and Overflowing Water

The compositions of Octopus Spring and Mushroom Spring (Yellowstone National Park, Wyoming, USA) microbial mats have been thoroughly studied, but the compositions of the effluent waters that flow above the mats have not. In this study, cells in the mats and overflowing waters of both springs were investigated at multiple sites where Synechococcus spp. are the dominant cyanobacteria (ca. 72°C to ca. 50°C), and on several dates. In addition to microscopic analyses of stained and autofluorescent cells, 16S rRNA gene sequencing was used to characterize the major taxa present and a protein-encoding gene (psaA) was sequenced and analyzed by ecotype simulation to predict species of Synechococcus The mats of both springs were similar in terms of the downstream distribution of predominant taxa detected previously. However, waters above these mats were predominated by taxa that reside in upstream mats or communities above the upper-temperature limit of the mat. A disturbance/recolonization study was performed at a site normally predominated by Synechococcus species adapted to low temperatures. After removing indigenous Synechococcus cells, Synechococcus species adapted to higher temperatures, which were predominant in the water overflowing this site, colonized the newly forming mat. Differences in recolonization under reduced and UV-screened irradiance suggested that, in addition to physical transport, environmental conditions likely select for species that are better adapted to these different conditions and can influence mat recovery. A transport model was developed and used to predict that, in Mushroom Spring, erosion predominates in the narrower and deeper upstream effluents and deposition predominates over erosion in wider and shallower downstream effluents.IMPORTANCE In flowing aquatic systems, cell erosion and deposition are important to the dispersal of cells from one location to another. Very little is known about microbial dispersal and the physical processes that underlie it. This study demonstrates its importance to colonization of downstream surfaces and especially to the recolonization and functioning of disturbed sites. Ecological systems in flowing environments are often, roughly speaking, pseudosteady, in that nutrients enter the system and by-products leave at relatively steady rates. Over time, material inputs and outputs must balance. Measurements of input fluxes (e.g., growth rates and proxies, such as photosynthesis rates) are frequent. However, erosion and deposition of cells are seldom measured and ecological significance is sometimes neglected. The importance of these parameters is immediately evident in any attempt to construct a model of long-time community behavior, as spatial ecological structure is significantly impacted and can be dominated by migration of organisms, even in small numbers.

Hot in Cold: Microbial Life in the Hottest Springs in Permafrost

Chukotka is an arctic region located in the continuous permafrost zone, but thermal springs are abundant there. In this study, for the first time, the microbial communities of the Chukotka hot springs (CHS) biofilms and sediments with temperatures 54–94 °C were investigated and analyzed by NGS sequencing of 16S rRNA gene amplicons. In microbial mats (54–75 °C), phototrophic bacteria of genus Chloroflexus dominated (up to 89% of all prokaryotes), while Aquificae were the most numerous at higher temperatures in Fe-rich sediments and filamentous “streamers” (up to 92%). The electron donors typical for Aquificae, such as H₂S and H₂, are absent or present only in trace amounts, and the prevalence of Aquificae might be connected with their ability to oxidize the ferrous iron present in CHS sediments. Armatimonadetes, Proteobacteria, Deinococcus-Thermus, Dictyoglomi, and Thermotogae, as well as uncultured bacteria (candidate divisions Oct-Spa1-106, GAL15, and OPB56), were numerous, and Cyanobacteria were present in low numbers. Archaea (less than 8% of the total community of each tested spring) belonged to Bathyarchaeota, Aigarchaeota, and Thaumarchaeota. The geographical location and the predominantly autotrophic microbial community, built on mechanisms other than the sulfur cycle-based ones, make CHS a special and unique terrestrial geothermal ecosystem.

Carbon Oxidation State in Microbial Polar Lipids Suggests Adaptation to Hot Spring Temperature and Redox Gradients

The influence of oxidation-reduction (redox) potential on the expression of biomolecules is a topic of ongoing exploration in geobiology. In this study, we investigate the novel possibility that structures and compositions of lipids produced by microbial communities are sensitive to environmental redox conditions. We extracted lipids from microbial biomass collected along the thermal and redox gradients of four alkaline hot springs in Yellowstone National Park (YNP) and investigated patterns in the average oxidation state of carbon (Z_C), a metric calculated from the chemical formulae of lipid structures. Carbon in intact polar lipids (IPLs) and their alkyl chains becomes more oxidized (higher Z_C) with increasing distance from each of the four hot spring sources. This coincides with decreased water temperature and increased concentrations of oxidized inorganic solutes, such as dissolved oxygen, sulfate, and nitrate. Carbon in IPLs is most reduced (lowest Z_C) in the hot, reduced conditions upstream, with abundance-weighted Z_C values between −1.68 and −1.56. These values increase gradually downstream to around −1.36 to −1.33 in microbial communities living between 29.0 and 38.1°C. This near-linear increase in Z_C can be attributed to a shift from ether-linked to ester-linked alkyl chains, a decrease in average aliphatic carbons per chain (nC), an increase in average degree of unsaturation per chain (nUnsat), and increased cyclization in tetraether lipids. The Z_C of lipid headgroups and backbones did not change significantly downstream. Expression of lipids with relatively reduced carbon under reduced conditions and oxidized lipids under oxidized conditions may indicate microbial adaptation across environmental gradients in temperature and electron donor/acceptor supply.

The Evolution of Reverse Gyrase Suggests a Nonhyperthermophilic Last Universal Common Ancestor

Reverse gyrase (RG) is the only protein found ubiquitously in hyperthermophilic organisms, but absent from mesophiles. As such, its simple presence or absence allows us to deduce information about the optimal growth temperature of long-extinct organisms, even as far as the last universal common ancestor of extant life (LUCA). The growth environment and gene content of the LUCA has long been a source of debate in which RG often features. In an attempt to settle this debate, we carried out an exhaustive search for RG proteins, generating the largest RG data set to date. Comprising 376 sequences, our data set allows for phylogenetic reconstructions of RG with unprecedented size and detail. These RG phylogenies are strikingly different from those of universal proteins inferred to be present in the LUCA, even when using the same set of species. Unlike such proteins, RG does not form monophyletic archaeal and bacterial clades, suggesting RG emergence after the formation of these domains, and/or significant horizontal gene transfer. Additionally, the branch lengths separating archaeal and bacterial groups are very short, inconsistent with the tempo of evolution from the time of the LUCA. Despite this, phylogenies limited to archaeal RG resolve most archaeal phyla, suggesting predominantly vertical evolution since the time of the last archaeal ancestor. In contrast, bacterial RG indicates emergence after the last bacterial ancestor followed by significant horizontal transfer. Taken together, these results suggest a nonhyperthermophilic LUCA and bacterial ancestor, with hyperthermophily emerging early in the evolution of the archaeal and bacterial domains.

Metal and organic pollutants bioremediation by extremophile microorganisms

Extremophiles comprise microorganisms that are able to grow and thrive in extreme environments, including in an acidic or alkaline pH, high or low temperatures, high concentrations of pollutants, and salts, among others. These organisms are promising for environmental biotechnology due to their unique physiological and enzymatic characteristics, which allow them to survive in harsh environments. Due to the stability and persistence of these microorganisms under adverse environmental conditions, they can be used for the bioremediation of environments contaminated with extremely recalcitrant pollutants. Here, we provide an overview of extremophiles and the role of "omics" in the field of bioremediation of environmental pollutants, including hydrocarbons, textile dyes and metals.

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.

Metagenomic Analysis of Hot Springs in Central India Reveals Hydrocarbon Degrading Thermophiles and Pathways Essential for Survival in Extreme Environments

Extreme ecosystems such as hot springs are of great interest as a source of novel extremophilic species, enzymes, metabolic functions for survival and biotechnological products. India harbors hundreds of hot springs, the majority of which are not yet explored and require comprehensive studies to unravel their unknown and untapped phylogenetic and functional diversity. The aim of this study was to perform a large-scale metagenomic analysis of three major hot springs located in central India namely, Badi Anhoni, Chhoti Anhoni, and Tattapani at two geographically distinct regions (Anhoni and Tattapani), to uncover the resident microbial community and their metabolic traits. Samples were collected from seven distinct sites of the three hot spring locations with temperature ranging from 43.5 to 98°C. The 16S rRNA gene amplicon sequencing of V3 hypervariable region and shotgun metagenome sequencing uncovered a unique taxonomic and metabolic diversity of the resident thermophilic microbial community in these hot springs. Genes associated with hydrocarbon degradation pathways, such as benzoate, xylene, toluene, and benzene were observed to be abundant in the Anhoni hot springs (43.5–55°C), dominated by Pseudomonas stutzeri and Acidovorax sp., suggesting the presence of chemoorganotrophic thermophilic community with the ability to utilize complex hydrocarbons as a source of energy. A high abundance of genes belonging to methane metabolism pathway was observed at Chhoti Anhoni hot spring, where methane is reported to constitute >80% of all the emitted gases, which was marked by the high abundance of Methylococcus capsulatus. The Tattapani hot spring, with a high-temperature range (61.5–98°C), displayed a lower microbial diversity and was primarily dominated by a nitrate-reducing archaeal species Pyrobaculum aerophilum. A higher abundance of cell metabolism pathways essential for the microbial survival in extreme conditions was observed at Tattapani. Taken together, the results of this study reveal a novel consortium of microbes, genes, and pathways associated with the hot spring environment.

Unraveling the microbial and functional diversity of Coamo thermal spring in Puerto Rico using metagenomic library generation and shotgun sequencing

In Puerto Rico, the microbial diversity of the thermal spring (ThS) in Coamo has never been studied using metagenomics. The focus of our research was to generate a metagenomic library from the ThS of Coamo, Puerto Rico and explore the microbial and functional diversity. The metagenomic library from the ThS waters was generated using direct DNA isolation. High molecular weight (40 kbp) DNA was end-repaired, electro eluted and ligated into a fosmid vector (pCCFOS1); then transduced into Escherichia coli EPI300-T1_R using T1 bacteriophages. The library consisted of approximately 6000 clones, 90% containing metagenomic DNA. Next-Generation-Sequencing technology (Illumina MiSeq) was used to process the ThS metagenome. After removing the cloning vector, 122,026 sequences with 33.10 Mbps size and 64% of G + C content were annotated and analyzed using the MG-RAST online server. Bacteria showed to be the most abundant domain (95.84%) followed by unidentified sequences (2.28%), viruses (1.67%), eukaryotes (0.15%), and archaea (0.01%). The most abundant phyla were Proteobacteria (95.03%), followed by unidentified (2.28%), unclassified from viruses (1.74%), Firmicutes (0.20%) and Actinobacteria (0.18%). The most abundant species were Escherichia coli, Polaromonas naphthalenivorans, Albidiferax ferrireducens and Acidovorax sp. Subsystem functional analysis showed that 20% of genes belong to transposable elements, 10% to clustering-based subsystems, and 8% to the production of cofactors. Functional analysis using NOG annotation showed that 82.79% of proteins are poorly characterized indicating the possibility of novel microbial functions and with potential biomedical and biotechnological applications. Metagenomic data was deposited into the NCBI database under the accession number SAMN06131862.

Metagenomics of Thermophiles with a Focus on Discovery of Novel Thermozymes

Microbial populations living in environments with temperatures above 50°C (thermophiles) have been widely studied, increasing our knowledge in the composition and function of these ecological communities. Since these populations express a broad number of heat-resistant enzymes (thermozymes), they also represent an important source for novel biocatalysts that can be potentially used in industrial processes. The integrated study of the whole-community DNA from an environment, known as metagenomics, coupled with the development of next generation sequencing (NGS) technologies, has allowed the generation of large amounts of data from thermophiles. In this review, we summarize the main approaches commonly utilized for assessing the taxonomic and functional diversity of thermophiles through metagenomics, including several bioinformatics tools and some metagenome-derived methods to isolate their thermozymes.