Metabolic characteristics of dominant microbes and key rare species from an acidic hot spring in Taiwan revealed by metagenomics

Influence of Geothermal Fumaroles in Driving the Microbial Community Dynamics and Functions of Adjacent Ecosystems

Growing evidence suggests that the hydrochemical properties of geothermal fumaroles may play a crucial role in shaping the diversity and functions of microbial communities in various environments. In the present study, the impact of geothermal furaneols on the microbial communities and their metabolic functions across the rock-soil-plant continuum was explored considering varying distances from the fumarole source. The results revealed that bacterial phylum Proteobacteria was predominant in all sample types, except in the 10 m rock sample, irrespective of the sampling distance. Archaeal phyla, such as Euryarchaeota and Crenarchaeota, were more prevalent in rock and soil samples, whereas bacterial phyla were more prevalent in plant samples. Thermoacidophilic archaeons, including Picrophilus, Ferroplasma, and Thermogymnomonas were dominant in rocks and soil samples of 1 and 5 m distances; acidophilic mesophiles, including Ferrimicrobium and Granulicella were abundant in the rhizoplane samples, whereas rhizosphere-associated microbes including Pseudomonas, Pedobacter, Rhizobium, and Novosphingobium were found dominant in the rhizosphere samples. The functional analysis highlighted the higher expression of sulfur oxidative pathways in the rock and soil samples; dark iron oxidation and nitrate/nitrogen respiratory functions in the rhizosphere samples. The findings underscore microbial adaptations across the rock-soil-plant continuum, emphasizing the intricate relationship between geothermal fumaroles and microbial communities in adjacent ecosystems. These insights offer a crucial understanding of the evolution of microbial life and highlight their pivotal roles in shaping ecosystem dynamics and functions.

pH and thiosulfate dependent microbial sulfur oxidation strategies across diverse environments

Sulfur oxidizing bacteria (SOB) play a key role in sulfur cycling in mine tailings impoundment (TI) waters, where sulfur concentrations are typically high. However, our understanding of SOB sulfur cycling via potential S oxidation pathways (sox, rdsr, and S4I) in these globally ubiquitous contexts, remains limited. Here, we identified TI water column SOB community composition, metagenomics derived metabolic repertoires, physicochemistry, and aqueous sulfur concentration and speciation in four Canadian base metal mine, circumneutral-alkaline TIs over four years (2016 - 2019). Identification and examination of genomes from nine SOB genera occurring in these TI waters revealed two pH partitioned, metabolically distinct groups, which differentially influenced acid generation and sulfur speciation. Complete sox (csox) dominant SOB (e.g., Halothiobacillus spp., Thiomonas spp.) drove acidity generation and S2O3 2- consumption via the csox pathway at lower pH (pH ~5 to ~6.5). At circumneutral pH conditions (pH ~6.5 to ~8.5), the presence of non-csox dominant SOB (hosting the incomplete sox, rdsr, and/or other S oxidation reactions; e.g. Thiobacillus spp., Sulfuriferula spp.) were associated with higher [S2O3 2-] and limited acidity generation. The S4I pathway part 1 (tsdA; S2O3 2- to S4O6 2-), was not constrained by pH, while S4I pathway part 2 (S4O6 2- disproportionation via tetH) was limited to Thiobacillus spp. and thus circumneutral pH values. Comparative analysis of low, natural (e.g., hydrothermal vents and sulfur hot springs) and high (e.g., Zn, Cu, Pb/Zn, and Ni tailings) sulfur systems literature data with these TI results, reveals a distinct TI SOB mining microbiome, characterized by elevated abundances of csox dominant SOB, likely sustained by continuous replenishment of sulfur species through tailings or mining impacted water additions. Our results indicate that under the primarily oxic conditions in these systems, S2O3 2- availability plays a key role in determining the dominant sulfur oxidation pathways and associated geochemical and physicochemical outcomes, highlighting the potential for biological management of mining impacted waters via pH and [S2O3 2-] manipulation.

Analysis of nearly 3000 archaeal genomes from terrestrial geothermal springs sheds light on interconnected biogeochemical processes

Terrestrial geothermal springs are physicochemically diverse and host abundant populations of Archaea. However, the diversity, functionality, and geological influences of these Archaea are not well understood. Here we explore the genomic diversity of Archaea in 152 metagenomes from 48 geothermal springs in Tengchong, China, collected from 2016 to 2021. Our dataset is comprised of 2949 archaeal metagenome-assembled genomes spanning 12 phyla and 392 newly identified species, which increases the known species diversity of Archaea by ~48.6%. The structures and potential functions of the archaeal communities are strongly influenced by temperature and pH, with high-temperature acidic and alkaline springs favoring archaeal abundance over Bacteria. Genome-resolved metagenomics and metatranscriptomics provide insights into the potential ecological niches of these Archaea and their potential roles in carbon, sulfur, nitrogen, and hydrogen metabolism. Furthermore, our findings illustrate the interplay of competition and cooperation among Archaea in biogeochemical cycles, possibly arising from overlapping functional niches and metabolic handoffs. Taken together, our study expands the genomic diversity of Archaea inhabiting geothermal springs and provides a foundation for more incisive study of biogeochemical processes mediated by Archaea in geothermal ecosystems.

Acidithiobacillia class members originating at sites within the Pacific Ring of Fire and other tectonically active locations and description of the novel genus 'Igneacidithiobacillus'

Recent studies have expanded the genomic contours of the Acidithiobacillia, highlighting important lacunae in our comprehension of the phylogenetic space occupied by certain lineages of the class. One such lineage is 'Igneacidithiobacillus', a novel genus-level taxon, represented by 'Igneacidithiobacillus copahuensis' VAN18-1T as its type species, along with two other uncultivated metagenome-assembled genomes (MAGs) originating from geothermally active sites across the Pacific Ring of Fire. In this study, we investigate the genetic and genomic diversity, and the distribution patterns of several uncharacterized Acidithiobacillia class strains and sequence clones, which are ascribed to the same 16S rRNA gene sequence clade. By digging deeper into this data and contributing to novel MAGs emerging from environmental studies in tectonically active locations, the description of this novel genus has been consolidated. Using state-of-the-art genomic taxonomy methods, we added to already recognized taxa, an additional four novel Candidate (Ca.) species, including 'Ca. Igneacidithiobacillus chanchocoensis' (mCHCt20-1TS), 'Igneacidithiobacillus siniensis' (S30A2T), 'Ca. Igneacidithiobacillus taupoensis' (TVZ-G3 TS), and 'Ca. Igneacidithiobacillus waiarikiensis' (TVZ-G4 TS). Analysis of published data on the isolation, enrichment, cultivation, and preliminary microbiological characterization of several of these unassigned or misassigned strains, along with the type species of the genus, plus the recoverable environmental data from metagenomic studies, allowed us to identify habitat preferences of these taxa. Commonalities and lineage-specific adaptations of the seven species of the genus were derived from pangenome analysis and comparative genomic metabolic reconstruction. The findings emerging from this study lay the groundwork for further research on the ecology, evolution, and biotechnological potential of the novel genus 'Igneacidithiobacillus'.

Deciphering microbial communities and their unique metabolic repertoire across rock-soil-plant continuum in the Dayoukeng fumarolic geothermal field of the Tatun Volcano Group

High temperature and sulfur concentrations in geothermal sulfur fumaroles host unique microbial ecosystems with niche-specific metabolic diversity and physiological functions. In this study, the microbial communities and their functionalities associated with the Dayoukeng geothermal field and the rock-soil-plant continuum were investigated to underpin the microbial modulation at different distances from the fumaroles source. At the phylum level, Bacteroidota, Planctomycetota, Armatimonadota, and Patescibacteria were abundant in plant samples; Elusimicrobiota and Desulfobacterota were in the rock samples while Nitrospirota, Micrarchaeota, and Deinococcota were dominant in the soil samples. Acidophilic thermophiles were enriched in samples within close proximity to the fumaroles, primarily at a distance of 1 m. The sulfur and iron-oxidizing acidophilic bacterial genera such as Acidothiobacillus and Sulfobacillus were abundant in the rock samples. The thermoacidophilic archaeon Acidianus and acidophilic bacteria Acidiphilium were abundant in the soil samples. Additionally, Thermosporothrix and Acidothermus were found abundant in the plant samples. The results of the functional annotation indicated that dark sulfur oxidation, iron oxidation, and hydrogen oxidation pathways were abundant in the soil samples up to 1 m from the fumaroles, while methanogenic and fermentation pathways were more prevalent in the soil samples located 10 m from the fumaroles. Interestingly, the results of this study indicated a higher microbial richness and abundance of acidophilic communities in the soils and plants compared to the rocks of the DYK fumarolic geothermal field.

Analysis and interpretation of hot springs water, biofilms, and sediment bacterial community profiling and their metabolic potential in the area of Taiwan geothermal ecosystem

Microorganisms developed a mechanism that copes with heat, acidity, and high dissolved metal concentrations that likely first evolved. The geothermal fluids emerging in the geothermal springs of Taiwan, located at a subduction zone, are still under signs of progress in the characterization of the various microbial taxonomic changes over time. However, no systematic studies have been performed to compare water, biofilms, and sediment bacterial communities and the primary driving force of dissolved and mineral substrates capable of supporting microbial metabolism. In this study, 16S rRNA gene sequencing was employed for bacterial community exploration, and their potential metabolic pathways involved from water, biofilms, and sediment samples, collected from the geothermal valley (Ti-re-ku). Metagenomic data revealed that the water samples had higher bacterial diversity and richness than biofilms and sediment samples. At the genus level, Alicyclobacillus, Thiomonas, Acidocella, Metallibacterium, Picrophilus, and Legionella were significantly abundant in the water samples. The biofilms were rich in Aciditerrimonas, Bacillus, Acidithiobacillus, and Lysinibacillus, whereas the sediment samples were abundant in Sulfobacillus. The PICRUSt2-predicted functional results revealed that heavy metal-related functions such as heavy-metal exporter system, cobalt‑zinc‑cadmium resistance, arsenical pump, high-affinity nickel-transport, and copper resistance metabolisms were significant in the water samples. Moreover, sulfur-related pathways such as thiosulfate oxidation, dissimilatory sulfate reduction, and assimilatory sulfate reduction were important in water samples, followed by biofilms and sediment. Therefore, our findings highlighted the comparative taxonomic diversity and functional composition contributions to geothermal fluid, with implications for understanding the evolution and ecological niche dimension of microbes which are the key to geothermal ecosystem function.

nifH gene expression and diversity in geothermal springs of Tengchong, China

Terrestrial hot springs have been suggested to harbor diverse diazotrophic lineages by using DNA-based nifH gene phylogenetic analysis. However, only a small amount of diazotrophs were ever confirmed to perform nitrogen fixation. In order to explore the compositions of active diazotrophic populations in hot springs, the in situ expression and diversity of nifH and 16S rRNA genes were investigated in the sediments of hot springs (pH 4.3-9.1; temperature 34-84°C) in Tengchong, China, by using high-throughput sequencing. The results showed that active diazotrophs were diverse in the studied Tengchong hot springs. The main active diazotrophs in high-temperature hot springs were affiliated with Aquificae, while those in low-temperature hot springs belonged to Cyanobacteria and Nitrospirae. Such dominance of Aquificae and Nitrospirae of diazotrophs has not been reported in other ecosystems. This suggests that hot springs may harbor unique active diazotrophs in comparison with other type of ecosystems. Furthermore, there were significant differences in the phylogenetic lineages of diazotrophs between hot springs of Tengchong and other regions, indicating that diazotrophs have geographical distribution patterns. Statistical analysis suggests that the expression and distribution of nifH gene were influenced by temperature and concentrations of ammonia and sulfur seem in Tengchong hot springs. These findings avail us to understand element cycling mediated by diazotrophs in hot spring ecosystems.

Geochemical and microbiological profiles in hydrothermal extreme acidic environments (Pisciarelli Spring, Campi Flegrei, Italy)

Although terrestrial hydrothermal systems are considered among the most fascinating environments, how their unique and extreme conditions can affect microorganisms selection and the role in biogeochemical cycles has not yet been well elucidated. A combined geochemical and microbiological exploration in waters and sediments from ten sampling points along a sharp temperature gradient (15-90 °C) within an extremely acidic hydrothermal system (Pisciarelli Spring, Campi Flegrei area, southern Italy) displayed how hydrothermal fluids influence the microbial dynamics. This area was characterized by high levels of reduced gaseous species (e.g. H2S, H2, CH4, CO), and very low pH values (<2.3). Thermodynamic calculations revealed a high microbial catabolic potential in oxidation/reduction reactions of N-, S-, and Fe-bearing species. Overall, an increase of the archaeal/bacterial abundance ratio was observed by decreasing temperature and pH values. In particular, Archaea and Bacteria were present in almost equal cell abundance (up to 1.1 × 109 and 9.3 × 108 cell/g, respectively) in the <70 °C sampling points (average pH = 2.09); on the contrary, highest temperature waters (85-90 °C; average pH = 2.26) were characterized by low abundance of archaeal cells. The high-throughput sequencing of 16S rRNA gene indicated strong differences in archaeal and bacterial communities' composition along temperature gradient. However, the microbiome in this extreme environment was mainly constituted by chemoautotrophic microorganisms that were likely involved in N-, S-, and Fe-bearing species transformations (e.g. Acidianus infernus, Ferroplasma acidarmanus, Acidithiobacillus, Sulfobacillus, Thaumarchaeota), in agreement with thermodynamic calculations.

The Evaluation of Bacterial Abundance and Functional Potentials in the Three Major Watersheds, Located in the Hot Spring Zone of the Tatun Volcano Group Basin, Taiwan

The Tatun Volcanic Group (TVG), located in northern Taiwan, is characterized by acidic hot springs where the outflow of the hot springs may affect the properties of the associated lotic water bodies. We investigated the bacterial diversity and functional profiles of the Peihuang (PHC), HuangGang (HGC), and Nanhuang Creeks (NHC) located in the TVG basin using 16S rRNA gene sequencing coupled with statistical analyses. Water samples were collected from various streams of the creeks for two months of the year. The NHC showed the highest diversity, richness, and a unique number of phyla, which was followed by the HGC. A reduced number of phyla and a lower diversity was noticed in the PHC. The NHC was found to be abundant in the genera Armatimonas, Prosthecobacter, Pirellula, and Bdellovibrio, whereas the HGC was rich in Thiomonas, Acidiphilium, Prevotella, Acidocella, Acidithiobacillus, and Metallibacterium. The PHC was abundant in Thiomonsa, Legionella, Acidocella, and Sulfuriferula. The samples did not show any strong seasonal variations with the bacterial diversity and abundance; however, the relative abundance of each sampling site varied within the sampling months. The iron transport protein- and the sulfur metabolism-related pathways were predicted to be the key functions in all the creeks, whereas the heavy metal-related functions, such as the cobalt/nickel transport protein and the cobalt-zinc-cadmium efflux system were found to be abundant in the HGC and PHC, respectively. The abundance of Bdellovibrio in the NHC, Diplorickettsia in the HGC, and Legionella in the PHC samples indicated a higher anthropogenic impact over the creek water quality. This study provides the data to understand the distinct bacterial community structure, as well as the functional potentials of the three major watersheds, and helps the knowledge of the impact of the physicochemical properties of the TVG hot springs upon the watersheds.

Comparative Metagenomic Analysis of Two Hot Springs From Ourense (Northwestern Spain) and Others Worldwide

With their circumneutral pH and their moderate temperature (66 and 68°C, respectively), As Burgas and Muiño da Veiga are two important human-use hot springs, previously studied with traditional culture methods, but never explored with a metagenomic approach. In the present study, we have performed metagenomic sequence-based analyses to compare the taxonomic composition and functional potential of these hot springs. Proteobacteria, Deinococcus-Thermus, Firmicutes, Nitrospirae, and Aquificae are the dominant phyla in both geothermal springs, but there is a significant difference in the abundance of these phyla between As Burgas and Muiño da Veiga. Phylum Proteobacteria dominates As Burgas ecosystem while Aquificae is the most abundant phylum in Muiño da Veiga. Taxonomic and functional analyses reveal that the variability in water geochemistry might be shaping the differences in the microbial communities inhabiting these geothermal springs. The content in organic compounds of As Burgas water promotes the presence of heterotrophic populations of the genera Acidovorax and Thermus, whereas the sulfate-rich water of Muiño da Veiga favors the co-dominance of genera Sulfurihydrogenibium and Thermodesulfovibrio. Differences in ammonia concentration exert a selective pressure toward the growth of nitrogen-fixing bacteria such as Thermodesulfovibrio in Muiño da Veiga. Temperature and pH are two important factors shaping hot springs microbial communities as was determined by comparative analysis with other thermal springs.

Genomic evolution of the class Acidithiobacillia: deep-branching Proteobacteria living in extreme acidic conditions

Members of the genus Acidithiobacillus, now ranked within the class Acidithiobacillia, are model bacteria for the study of chemolithotrophic energy conversion under extreme conditions. Knowledge of the genomic and taxonomic diversity of Acidithiobacillia is still limited. Here, we present a systematic analysis of nearly 100 genomes from the class sampled from a wide range of habitats. Some of these genomes are new and others have been reclassified on the basis of advanced genomic analysis, thus defining 19 Acidithiobacillia lineages ranking at different taxonomic levels. This work provides the most comprehensive classification and pangenomic analysis of this deep-branching class of Proteobacteria to date. The phylogenomic framework obtained illuminates not only the evolutionary past of this lineage, but also the molecular evolution of relevant aerobic respiratory proteins, namely the cytochrome bo₃ ubiquinol oxidases.

Metatranscriptomic analysis reveals active microbes and genes responded to short-term Cr(VI) stress

Heavy metals have been severely polluting the environment. However, the response mechanism of microbial communities to short-term heavy metals stress remains unclear. In this study, metagenomics (MG) and metatranscriptomics (MT) was performed to observe the microbial response to short-term Cr(VI) stress. MG data showed that 99.1% of species were similar in the control and Cr(VI) treated groups. However, MT data demonstrated that 83% of the microbes were active in which 58.7% increased, while the relative abundance of 41.3% decreased after short-term Cr(VI) incubation. The MT results also revealed 9% of microbes were dormant in samples. Genes associated with oxidative stress, Cr(VI) transport, resistance, and reduction, as well as genes with unknown functions were 2-10 times upregulated after Cr(VI) treatment. To further confirm the function of unknown genes, two genes (314 and 494) were selected to detect the Cr(VI) resistance and reduction ability. The results showed that these genes significantly increased the Cr(VI) remediation ability of Escherichia coli. MT results also revealed an increase in the expression of some rare genera (at least two times) after Cr(VI) treatment, indicating these rare species played a crucial role in microbial response to short-term Cr(VI) stress. In summary, MT is an efficient way to understand the role of active and dormant microbes in specific environmental conditions.

Genomic adaptations enabling Acidithiobacillus distribution across wide-ranging hot spring temperatures and pHs

BACKGROUND

Terrestrial hot spring settings span a broad spectrum of physicochemistries. Physicochemical parameters, such as pH and temperature, are key factors influencing differences in microbial composition across diverse geothermal areas. Nonetheless, analysis of hot spring pools from the Taupo Volcanic Zone (TVZ), New Zealand, revealed that some members of the bacterial genus, Acidithiobacillus, are prevalent across wide ranges of hot spring pHs and temperatures. To determine the genomic attributes of Acidithiobacillus that inhabit such diverse conditions, we assembled the genomes of 19 uncultivated hot spring Acidithiobacillus strains from six geothermal areas and compared these to 37 publicly available Acidithiobacillus genomes from various habitats.

RESULTS

Analysis of 16S rRNA gene amplicons from 138 samples revealed that Acidithiobacillus comprised on average 11.4 ± 16.8% of hot spring prokaryotic communities, with three Acidithiobacillus amplicon sequence variants (ASVs) (TVZ_G1, TVZ_G2, TVZ_G3) accounting for > 90% of Acidithiobacillus in terms of relative abundance, and occurring in 126 out of 138 samples across wide ranges of temperature (17.5-92.9 °C) and pH (1.0-7.5). We recovered 19 environmental genomes belonging to each of these three ASVs, as well as a fourth related group (TVZ_G4). Based on genome average nucleotide identities, the four groups (TVZ_G1-TVZ_G4) constitute distinct species (ANI < 96.5%) of which three are novel Acidithiobacillus species (TVZ_G2-TVZ_G4) and one belongs to Acidithiobacillus caldus (TVZ_G1). All four TVZ Acidithiobacillus groups were found in hot springs with temperatures above the previously known limit for the genus (up to 40 °C higher), likely due to significantly higher proline and GC contents than other Acidithiobacillus species, which are known to increase thermostability. Results also indicate hot spring-associated Acidithiobacillus have undergone genome streamlining, likely due to thermal adaptation. Moreover, our data suggest that Acidithiobacillus prevalence across varied hot spring pHs is supported by distinct strategies, whereby TVZ_G2-TVZ_G4 regulate pH homeostasis mostly through Na⁺/H⁺ antiporters and proton-efflux ATPases, whereas TVZ_G1 mainly relies on amino acid decarboxylases.

CONCLUSIONS

This study provides insights into the distribution of Acidithiobacillus species across diverse hot spring physichochemistries and determines genomic features and adaptations that potentially enable Acidithiobacillus species to colonize a broad range of temperatures and pHs in geothermal environments. Video Abstract.

Metagenomics revealing molecular profiling of community structure and metabolic pathways in natural hot springs of the Sikkim Himalaya

BACKGROUND

Himalaya is an ecologically pristine environment. The geo-tectonic activities have shaped various environmental niches with diverse microbial populations throughout the Himalayan biosphere region. Albeit, limited information is available in terms of molecular insights into the microbiome, including the uncultured microbes, of the Himalayan habitat. Hence, a vast majority of genomic resources are still under-explored from this region. Metagenome analysis has simplified the extensive in-depth exploration of diverse habitats. In the present study, the culture-independent whole metagenome sequencing methodology was employed for microbial diversity exploration and identification of genes involved in various metabolic pathways in two geothermal springs located at different altitudes in the Sikkim Himalaya.

RESULTS

The two hot springs, Polok and Reshi, have distinct abiotic conditions. The average temperature of Polok and Reshi was recorded to be 62 °C and 43 °C, respectively. Both the aquatic habitats have alkaline geochemistry with pH in the range of 7-8. Community profile analysis revealed genomic evidence of plentiful bacteria, with a minute fraction of the archaeal population in hot water reservoirs of Polok and Reshi hot spring. Mesophilic microbes belonging to Proteobacteria and Firmicutes phyla were predominant at both the sites. Polok exhibited an extravagant representation of Chloroflexi, Deinococcus-Thermus, Aquificae, and Thermotogae. Metabolic potential analysis depicted orthologous genes associated with sulfur, nitrogen, and methane metabolism, contributed by the microflora in the hydrothermal system. The genomic information of many novel carbohydrate-transforming enzymes was deciphered in the metagenomic description. Further, the genomic capacity of antimicrobial biomolecules and antibiotic resistance were discerned.

CONCLUSION

The study provided comprehensive molecular information about the microbial treasury as well as the metabolic features of the two geothermal sites. The thermal aquatic niches were found a potential bioresource of biocatalyst systems for biomass-processing. Overall, this study provides the whole metagenome based insights into the taxonomic and functional profiles of Polok and Reshi hot springs of the Sikkim Himalaya. The study generated a wealth of genomic data that can be explored for the discovery and characterization of novel genes encoding proteins of industrial importance.

Control of ecological outcomes through deliberate parameter changes in a model of the gut microbiome

The generalized Lotka-Volterra (gLV) equations are a mathematical proxy for ecological dynamics. We focus on a gLV model of the gut microbiome, in which the evolution of the gut microbial state is determined in part by pairwise interspecies interaction parameters that encode environmentally mediated resource competition between microbes. We develop an in silico method that controls the steady-state outcome of the system by adjusting these interaction parameters. This approach is confined to a bistable region of the gLV model. In this method, a dimensionality reduction technique called steady-state reduction (SSR) is first used to generate a two-dimensional (2D) gLV model that approximates the high-dimensional dynamics on the 2D subspace spanned by the two steady states. Then a bifurcation analysis of the 2D model analytically determines parameter modifications that drive an initial condition to a target steady state. This parameter modification of the reduced 2D model guides parameter modifications of the original high-dimensional model, resulting in a change of steady-state outcome in the high-dimensional model. This control method, called SSR-guided parameter change (SPARC), bypasses the computational challenge of directly determining parameter modifications in the original high-dimensional system. SPARC could guide the development of indirect bacteriotherapies, which seek to change microbial compositions by deliberately modifying gut environmental variables such as gut acidity or macronutrient availability.

Species Sorting and Neutral Theory Analyses Reveal Archaeal and Bacterial Communities Are Assembled Differently in Hot Springs

Although the recognition of archaea as one of the three kingdoms in the tree of life has been nearly a half-century long, the comparative investigations on their ecological adaptations with bacteria have been limited. The mechanisms of their community assembly and diversity maintenance in hot springs have not been addressed. The mechanistic study is critical not only for understanding the hot-spring microbiome structure and dynamics, but also for shedding light on their evolutionary adaptations. We applied the neutral theory model and species sorting paradigm of metacommunity theory to investigate how hot-spring microbial communities were assembled, how their diversities were maintained, and how the temperature and pH influence these mechanisms. Through rigorous statistical tests based on the neutral theory and species sorting paradigm, we found (i) According to the neutral theory, archaeal and bacterial communities are assembled differently, with stochastic neutral force playing a more significant role in archaeal communities than in bacterial communities (neutrality-rate = 52.9 vs. 15.8%, p-value < 0.05). (ii) Temperature and pH account for rather limited (<10%) variations in hot-spring microbiomes based on the species sorting paradigm. The pH has more significant influences than temperature on archaeal communities, and both pH and temperature have similarly low influences on bacterial community structure. (iii) We postulate that the differences between archaea and bacteria are likely due to the longer evolutionary history and better adaptation of archaea to host spring environments.

Diurnal Changes in Active Carbon and Nitrogen Pathways Along the Temperature Gradient in Porcelana Hot Spring Microbial Mat

Composition, carbon and nitrogen uptake, and gene transcription of microbial mat communities in Porcelana neutral hot spring (Northern Chilean Patagonia) were analyzed using metagenomics, metatranscriptomics and isotopically labeled carbon (H13CO3) and nitrogen (15NH4Cl and K15NO3) assimilation rates. The microbial mat community included 31 phyla, of which only Cyanobacteria and Chloroflexi were dominant. At 58°C both phyla co-occurred, with similar contributions in relative abundances in metagenomes and total transcriptional activity. At 66°C, filamentous anoxygenic phototrophic Chloroflexi were >90% responsible for the total transcriptional activity recovered, while Cyanobacteria contributed most metagenomics and metatranscriptomics reads at 48°C. According to such reads, phototrophy was carried out both through oxygenic photosynthesis by Cyanobacteria (mostly Mastigocladus) and anoxygenic phototrophy due mainly to Chloroflexi. Inorganic carbon assimilation through the Calvin-Benson cycle was almost exclusively due to Mastigocladus, which was the main primary producer at lower temperatures. Two other CO2 fixation pathways were active at certain times and temperatures as indicated by transcripts: 3-hydroxypropionate (3-HP) bi-cycle due to Chloroflexi and 3-hydroxypropionate-4-hydroxybutyrate (HH) cycle carried out by Thaumarchaeota. The active transcription of the genes involved in these C-fixation pathways correlated with high in situ determined carbon fixation rates. In situ measurements of ammonia assimilation and nitrogen fixation (exclusively attributed to Cyanobacteria and mostly to Mastigocladus sp.) showed these were the most important nitrogen acquisition pathways at 58 and 48°C. At 66°C ammonia oxidation genes were actively transcribed (mostly due to Thaumarchaeota). Reads indicated that denitrification was present as a nitrogen sink at all temperatures and that dissimilatory nitrate reduction to ammonia (DNRA) contributed very little. The combination of metagenomic and metatranscriptomic analysis with in situ assimilation rates, allowed the reconstruction of day and night carbon and nitrogen assimilation pathways together with the contribution of keystone microorganisms in this natural hot spring microbial mat.

Active Crossfire Between Cyanobacteria and Cyanophages in Phototrophic Mat Communities Within Hot Springs

Cyanophages are viruses with a wide distribution in aquatic ecosystems, that specifically infect Cyanobacteria. These viruses can be readily isolated from marine and fresh waters environments; however, their presence in cosmopolitan thermophilic phototrophic mats remains largely unknown. This study investigates the morphological diversity (TEM), taxonomic composition (metagenomics), and active infectivity (metatranscriptomics) of viral communities over a thermal gradient in hot spring phototrophic mats from Northern Patagonia (Chile). The mats were dominated (up to 53%) by cosmopolitan thermophilic filamentous true-branching cyanobacteria from the genus Mastigocladus, the associated viral community was predominantly composed of Caudovirales (70%), with most of the active infections driven by cyanophages (up to 90% of Caudovirales transcripts). Metagenomic assembly lead to the first full genome description of a T7-like Thermophilic Cyanophage recovered from a hot spring (Porcelana Hot Spring, Chile), with a temperature of 58°C (TC-CHP58). This could potentially represent a world-wide thermophilic lineage of podoviruses that infect cyanobacteria. In the hot spring, TC-CHP58 was active over a temperature gradient from 48 to 66°C, showing a high population variability represented by 1979 single nucleotide variants (SNVs). TC-CHP58 was associated to the Mastigocladus spp. by CRISPR spacers. Marked differences in metagenomic CRISPR loci number and spacers diversity, as well as SNVs, in the TC-CHP58 proto-spacers at different temperatures, reinforce the theory of co-evolution between natural virus populations and cyanobacterial hosts. Considering the importance of cyanobacteria in hot spring biogeochemical cycles, the description of this new cyanopodovirus lineage may have global implications for the functioning of these extreme ecosystems.

Towards long-read metagenomics: complete assembly of three novel genomes from bacteria dependent on a diazotrophic cyanobacterium in a freshwater lake co-culture

Here we report three complete bacterial genome assemblies from a PacBio shotgun metagenome of a co-culture from Upper Klamath Lake, OR. Genome annotations and culture conditions indicate these bacteria are dependent on carbon and nitrogen fixation from the cyanobacterium Aphanizomenon flos-aquae, whose genome was assembled to draft-quality. Due to their taxonomic novelty relative to previously sequenced bacteria, we have temporarily designated these bacteria as incertae sedis Hyphomonadaceae strain UKL13-1 (3,501,508 bp and 56.12% GC), incertae sedis Betaproteobacterium strain UKL13-2 (3,387,087 bp and 54.98% GC), and incertae sedis Bacteroidetes strain UKL13-3 (3,236,529 bp and 37.33% GC). Each genome consists of a single circular chromosome with no identified plasmids. When compared with binned Illumina assemblies of the same three genomes, there was ~7% discrepancy in total genome length. Gaps where Illumina assemblies broke were often due to repetitive elements. Within these missing sequences were essential genes and genes associated with a variety of functional categories. Annotated gene content reveals that both Proteobacteria are aerobic anoxygenic phototrophs, with Betaproteobacterium UKL13-2 potentially capable of phototrophic oxidation of sulfur compounds. Both proteobacterial genomes contain transporters suggesting they are scavenging fixed nitrogen from A. flos-aquae in the form of ammonium. Bacteroidetes UKL13-3 has few completely annotated biosynthetic pathways, and has a comparatively higher proportion of unannotated genes. The genomes were detected in only a few other freshwater metagenomes, suggesting that these bacteria are not ubiquitous in freshwater systems. Our results indicate that long-read sequencing is a viable method for sequencing dominant members from low-diversity microbial communities, and should be considered for environmental metagenomics when conditions meet these requirements.

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.