Seasonal patterns in microbial communities inhabiting the hot springs of Tengchong, Yunnan Province, China

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.

Metagenome-assembled genomes reveal carbohydrate degradation and element metabolism of microorganisms inhabiting Tengchong hot springs, China

A hot spring is a distinctive aquatic environment that provides an excellent system to investigate microorganisms and their function in elemental cycling processes. Previous studies of terrestrial hot springs have been mostly focused on the microbial community, one special phylum or category, or genes involved in a particular metabolic step, while little is known about the overall functional metabolic profiles of microorganisms inhabiting the terrestrial hot springs. Here, we analyzed the microbial community structure and their functional genes based on metagenomic sequencing of six selected hot springs with different temperature and pH conditions. We sequenced a total of 11 samples from six hot springs and constructed 162 metagenome-assembled genomes (MAGs) with completeness above 70% and contamination lower than 10%. Crenarchaeota, Euryarchaeota and Aquificae were found to be the dominant phyla. Functional annotation revealed that bacteria encode versatile carbohydrate-active enzymes (CAZYmes) for the degradation of complex polysaccharides, while archaea tend to assimilate C1 compounds through carbon fixation. Under nitrogen-deficient conditions, there were correspondingly fewer genes involved in nitrogen metabolism, while abundant and diverse set of genes participating in sulfur metabolism, particularly those associated with sulfide oxidation and thiosulfate disproportionation. In summary, archaea and bacteria residing in the hot springs display distinct carbon metabolism fate, while sharing the common energy preference through sulfur metabolism. Overall, this research contributes to a better comprehension of biogeochemistry of terrestrial hot springs.

Resource partitioning and amino acid assimilation in a terrestrial geothermal spring

High-temperature geothermal springs host simplified microbial communities; however, the activities of individual microorganisms and their roles in the carbon cycle in nature are not well understood. Here, quantitative stable isotope probing (qSIP) was used to track the assimilation of 13C-acetate and 13C-aspartate into DNA in 74 °C sediments in Gongxiaoshe Hot Spring, Tengchong, China. This revealed a community-wide preference for aspartate and a tight coupling between aspartate incorporation into DNA and the proliferation of aspartate utilizers during labeling. Both 13C incorporation into DNA and changes in the abundance of taxa during incubations indicated strong resource partitioning and a significant phylogenetic signal for aspartate incorporation. Of the active amplicon sequence variants (ASVs) identified by qSIP, most could be matched with genomes from Gongxiaoshe Hot Spring or nearby springs with an average nucleotide similarity of 99.4%. Genomes corresponding to aspartate primary utilizers were smaller, near-universally encoded polar amino acid ABC transporters, and had codon preferences indicative of faster growth rates. The most active ASVs assimilating both substrates were not abundant, suggesting an important role for the rare biosphere in the community response to organic carbon addition. The broad incorporation of aspartate into DNA over acetate by the hot spring community may reflect dynamic cycling of cell lysis products in situ or substrates delivered during monsoon rains and may reflect N limitation.

Evidence for nontraditional mcr-containing archaea contributing to biological methanogenesis in geothermal springs

Recent discoveries of methyl-coenzyme M reductase-encoding genes (mcr) in uncultured archaea beyond traditional euryarchaeotal methanogens have reshaped our view of methanogenesis. However, whether any of these nontraditional archaea perform methanogenesis remains elusive. Here, we report field and microcosm experiments based on ¹³C-tracer labeling and genome-resolved metagenomics and metatranscriptomics, revealing that nontraditional archaea are predominant active methane producers in two geothermal springs. Archaeoglobales performed methanogenesis from methanol and may exhibit adaptability in using methylotrophic and hydrogenotrophic pathways based on temperature/substrate availability. A five-year field survey found Candidatus Nezhaarchaeota to be the predominant mcr-containing archaea inhabiting the springs; genomic inference and mcr expression under methanogenic conditions strongly suggested that this lineage mediated hydrogenotrophic methanogenesis in situ. Methanogenesis was temperature-sensitive , with a preference for methylotrophic over hydrogenotrophic pathways when incubation temperatures increased from 65° to 75°C. This study demonstrates an anoxic ecosystem wherein methanogenesis is primarily driven by archaea beyond known methanogens, highlighting diverse nontraditional mcr-containing archaea as previously unrecognized methane sources.

Temporal dynamics of geothermal microbial communities in Aotearoa-New Zealand

Microbial biogeography studies, in particular for geothermal-associated habitats, have focused on spatial patterns and/or individual sites, which have limited ability to describe the dynamics of ecosystem behaviour. Here, we report the first comprehensive temporal study of bacterial and archaeal communities from an extensive range of geothermal features in Aotearoa-New Zealand. One hundred and fifteen water column samples from 31 geothermal ecosystems were taken over a 34-month period to ascertain microbial community stability (control sites), community response to both natural and anthropogenic disturbances in the local environment (disturbed sites) and temporal variation in spring diversity across different pH values (pH 3, 5, 7, 9) all at a similar temperature of 60–70°C (pH sites). Identical methodologies were employed to measure microbial diversity via 16S rRNA gene amplicon sequencing, along with 44 physicochemical parameters from each feature, to ensure confidence in comparing samples across timeframes. Our results indicated temperature and associated groundwater physicochemistry were the most likely parameters to vary stochastically in these geothermal features, with community abundances rather than composition more readily affected by a changing environment. However, variation in pH (pH ±1) had a more significant effect on community structure than temperature (±20°C), with alpha diversity failing to adequately measure temporal microbial disparity in geothermal features outside of circumneutral conditions. While a substantial physicochemical disturbance was required to shift community structures at the phylum level, geothermal ecosystems were resilient at this broad taxonomic rank and returned to a pre-disturbed state if environmental conditions re-established. These findings highlight the diverse controls between different microbial communities within the same habitat-type, expanding our understanding of temporal dynamics in extreme ecosystems.

Spatial and temporal dynamics at an actively silicifying hydrothermal system

Steep Cone Geyser is a unique geothermal feature in Yellowstone National Park (YNP), Wyoming, actively gushing silicon-rich fluids along outflow channels possessing living and actively silicifying microbial biomats. To assess the geomicrobial dynamics occurring temporally and spatially at Steep Cone, samples were collected at discrete locations along one of Steep Cone's outflow channels for both microbial community composition and aqueous geochemistry analysis during field campaigns in 2010, 2018, 2019, and 2020. Geochemical analysis characterized Steep Cone as an oligotrophic, surface boiling, silicious, alkaline-chloride thermal feature with consistent dissolved inorganic carbon and total sulfur concentrations down the outflow channel ranging from 4.59 ± 0.11 to 4.26 ± 0.07 mM and 189.7 ± 7.2 to 204.7 ± 3.55 μM, respectively. Furthermore, geochemistry remained relatively stable temporally with consistently detectable analytes displaying a relative standard deviation <32%. A thermal gradient decrease of ~55°C was observed from the sampled hydrothermal source to the end of the sampled outflow transect (90.34°C ± 3.38 to 35.06°C ± 7.24). The thermal gradient led to temperature-driven divergence and stratification of the microbial community along the outflow channel. The hyperthermophile Thermocrinis dominates the hydrothermal source biofilm community, and the thermophiles Meiothermus and Leptococcus dominate along the outflow before finally giving way to more diverse and even microbial communities at the end of the transect. Beyond the hydrothermal source, phototrophic taxa such as Leptococcus, Chloroflexus, and Chloracidobacterium act as primary producers for the system, supporting heterotrophic growth of taxa such as Raineya, Tepidimonas, and Meiothermus. Community dynamics illustrate large changes yearly driven by abundance shifts of the dominant taxa in the system. Results indicate Steep Cone possesses dynamic outflow microbial communities despite stable geochemistry. These findings improve our understanding of thermal geomicrobiological dynamics and inform how we can interpret the silicified rock record.

Interspecific and intraspecific Taylor's laws for frog skin microbes

Amphibians are known to have an abundance of microorganisms colonizing their skin, and these symbionts often protect the host from disease. There are now many comprehensive studies on amphibian skin microbes, but the interspecific and intraspecific abundance distributions (or abundance heterogeneity) of amphibian skin microbes remain unclear. Furthermore, we have a very limited understanding of how the abundance and heterogeneity of microbial communities relate to the body size (or more specifically, skin surface area) of amphibian hosts. In this study, we evaluated the interspecific and intraspecific abundance distribution patterns of amphibian skin microbes and evaluated whether the symbiotic skin microbes of different anuran species share a fundamental heterogeneity scaling parameter. If scaling invariance exists, we hypothesize that a fundamental heterogeneity scaling value also exists. A total of 358 specimens of 10 amphibian host species were collected, and we used Type-I and III Taylor's power law expansions (TPLE) to assess amphibian skin microbial heterogeneity at the community and mixed-species population levels, respectively. The obtained results showed that, at the community scale, a high aggregation of the microbial abundance distribution on the skin barely changed with host size. In a mixed-species population (i.e., a community context), the abundance distribution pattern of mixed microbial species populations also does not change with host size and always remains highly aggregated. These findings suggest that while amphibian skin microbiomes located in different hosts may have different environmental conditions, they share a fundamental heterogeneity scaling parameter, and thus, scale invariance exists. Finally, we found that microhabitat area provided by the host skin is vital to the stability of the symbiotic microbial community.

Microbial community composition and hydrochemistry of underexplored geothermal waters in Croatia

In Croatia, a variety of geothermal springs with a wide temperature range and varied hydrochemical conditions exist, and they may harbor different niches for the distribution of microbial communities. In this study, 19 different sites, mainly located in central and eastern Croatia, were selected for primary characterization of spring hydrochemistry and microbial community composition. Using 16S rRNA gene amplicon sequencing, it was found that the bacterial communities that dominated most geothermal waters were related to Proteobacteria and Campylobacteria, while most archaeal sequences were related to Crenarchaeota. At the genus level, the prokaryotic community was highly site-specific and was often dominated by a single genus, including sites dominated by Hydrogenophilus, Sulfuricurvum, Sulfurovum, Thiofaba and Nitrospira, while the most abundant archaeal genera were affiliated to the ammonia-oxidizing archaea, Candidatus Nitrosotenuis and Candidatus Nitrososphaera. Whereas the microbial communities were overall highly location-specific, temperature, pH, ammonia, nitrate, total nitrogen, sulfate and hydrogen sulfide, as well as dissolved organic and inorganic carbon, were the abiotic factors that significantly affected microbial community composition. Furthermore, an aquifer-type effect was observed in the community composition, but there was no pronounced seasonal variability for geothermal spring communities (i.e. the community structure was mainly stable during the three seasons sampled). These results surprisingly pointed to stable and geographically unique microbial communities that were adapted to different geothermal water environments throughout Croatia. Knowing which microbial communities are present in these extreme habitats is essential for future research. They will allow us to explore further the microbial metabolisms prevailing at these geothermal sites that have high potential for biotechnological uses, as well as the establishment of the links between microbial community structure and the physicochemical environment of geothermal waters.

Primary microbial succession in the anchialine ecosystem

When new land is created, initial microbial colonization lays the foundation for further ecological succession of plant and animal communities. Primary microbial succession of new aquatic habitats formed during volcanic activity has received little attention. The anchialine ecosystem, which includes coastal ponds in young lava flows, offers an opportunity to examine this process. Here, we characterized microbial communities of anchialine habitats in Hawaii that were created during volcanic eruptions in 2018. Benthic samples from three habitats were collected ∼2 years after their formation and at later time points spanning ∼1 year. Sequence profiling (16S and 18S) of prokaryotic and eukaryotic communities was used to test whether communities were similar to those from older, established anchialine habitats, and if community structure changed over time. Results show that microbial communities from the new habitats were unlike any from established anchialine microbial communities, having higher proportions of Planctomycetota and Chloroflexi but lower proportions of green algae. Each new habitat also harbored its own unique community relative to other habitats. While community composition in each habitat underwent statistically significant changes over time, they remained distinctive from established anchialine habitats. New habitats also had highly elevated temperatures compared to other habitats. These results suggest idiosyncratic microbial consortia form during early succession of Hawaiian anchialine habitats. Future monitoring will reveal whether the early communities described here remain stable after temperatures decline and macro-organisms become more abundant, or if microbial communities will continue to change and eventually resemble those of established habitats. This work is a key first step in examining primary volcanic succession in aquatic habitats and suggests young anchialine habitats may warrant special conservation status.

Minerals Determined a Special Ecological Niche and Selectively Enriched Microbial Species from Bulk Water Communities in Hot Springs

Minerals provide physical niches and supply nutrients or serve as electron donors/acceptors for microorganism survival and growth, and thus minerals and microbes co-evolved. Yet, little is known about how sediment minerals impact microbial community assembly in hot springs and to what extent mineralogical composition influences microbial community composition and diversity. Here the influences of minerals on thermophiles in Tengchong hot springs were revealed by network analysis of field samples, as well as in-situ microcosm experiments with minerals. A molecular ecological network was constructed based on high throughput sequencing data of 16S rRNA gene, with a combination of water geochemistry and sedimentary mineralogical compositions. Six modules were identified and this highly modular network structure represents the microbial preference to different abiotic factors, consequently resulting in niche partitioning in sedimentary communities in hot springs. Diverse mineralogical compositions generated special niches for microbial species. Subsequently, the in-situ microcosm experiments with four minerals (aragonite, albite, K-feldspar, and quartz) and spring water were conducted in a silicate-hosted alkaline spring (i.e., Gmq) and a carbonate-hosted neutral hot spring (i.e., Gxs) for 70 days. Different microbial preferences were observed among different mineral types (carbonate versus silicate). Aragonite microcosms in Gmq spring enriched archaeal genera Sulfophobococcus and Aeropyrum within the order Desulfurococcales by comparison with both in-situ water and silicate microcosms. Sulfophobococcus was also accumulated in Gxs aragonite microcosms, but the contribution to overall dissimilarity is much lower than that in Gmq spring. Besides, Caldimicrobium was a bacterial genus enriched in Gxs aragonite microcosms, in contrast to in-situ water and silicate microcosms, whereas Candidatus Kryptobacter and Thermus were more abundant in silicate microcosms. The differences in microbial accumulations among different mineral types in the same spring implied that mineral chemistry may exert extra deterministic selective pressure in drawing certain species from the bulk water communities, in addition to stochastic absorption on mineral surface. Taken together, our results highlight the special niche partitioning determined by mineralogical compositions and further confirm that minerals could be used as “fishing bait” to enrich certain rare microbial species.

An emerging view of the diversity, ecology and function of Archaea in alkaline hydrothermal environments

The described diversity within the domain Archaea has recently expanded due to advances in sequencing technologies, but many habitats that likely harbor novel lineages of archaea remain understudied. Knowledge of archaea within natural and engineered hydrothermal systems, such as hot springs and engineered subsurface habitats, has been steadily increasing, but the majority of the work has focused on archaea living in acidic or circumneutral environments. The environmental pressures exerted by the combination of high temperatures and high pH likely select for divergent communities and distinct metabolic pathways from those observed in acidic or circumneutral systems. In this review, we examine what is currently known about the archaea found in thermoalkaline environments, focusing on the detection of novel lineages and knowledge of the ecology, metabolic pathways and functions of these populations and communities. We also discuss the potential of emerging multi-omics approaches, including proteomics and metabolomics, to enhance our understanding of archaea within extreme thermoalkaline systems.

The Effect of Spring Water Geochemistry on Copper Proteins in Tengchong Hot Springs, China

Copper (Cu) is an essential trace metal cofactor for a variety of proteins; however, excess Cu is toxic to most organisms. Cu homeostasis is maintained by a complex machinery of Cu binding proteins that control the uptake, transport, sequestration, and efflux of Cu ions. Despite the importance of Cu binding proteins in electron transfer, substrate oxidation, superoxide dismutation, and denitrification, little information exists about microbial Cu utilization in extreme environments, where the geochemical conditions may affect Cu bioavailability. Using metagenomic data from 9 hot springs in Tengchong, China, which range in temperature from 42°C to 96°C and in pH from 2.3 to 9, the effects of pH, temperature, and spring geochemistry on the distribution of Cu binding domains of proteins and oxidoreductases were studied. Dissolved Cu and Cu binding domains were detected across all temperature and pH gradients. Cu binding domains of cytochrome c oxidase subunits, heavy-metal-associated domains, and nitrous oxide reductase were detected at all sites. DoxB, a quinol oxidase, and other quinol oxidase subunits were the dominant Cu binding oxidoreductase subunits present at low-pH and high-temperature sites, whereas cbb ₃-type cytochrome c oxidase subunits were dominant at high-pH and high-temperature sites. Additionally, aa ₃-type cytochrome c oxidase was more prominent than cbb ₃-type cytochrome c oxidase under circumneutral-pH conditions. This suggests that the type of cytochrome c oxidase pathway and the Cu proteins employed by microbes to carry out important functions such as energy acquisition and efflux of excess Cu are affected by the physicochemical conditions of the springs.IMPORTANCE Copper is present in a variety of proteins and is required to carry out essential functions by all organisms. However, in hot spring environments, copper availability may be limited due to the high temperatures and the wide range in pH. The significance of our research is in relating the physicochemical environment to the distribution of copper proteins across hot spring environments, which provides increased understanding of primary functions and adaptions in these environments.

Temperature governs the distribution of hot spring microbial community in three hydrothermal fields, Eastern Tibetan Plateau Geothermal Belt, Western China

The eastern Tibetan Plateau geothermal belt in the southwest of China hosts a number of hot springs with a wide range of temperature and hydrogeochemical conditions, which may harbor different niches for the distribution of microbial communities. In this study, we investigated hydrochemical characteristics and microbial community composition in 16 hot springs with a temperature range of 34.6 to 88.2 °C within and across three typical hydrothermal fields (Kangding, Litang, and Batang). According to aquifer lithologic and tectonic differences, the hydrochemical compositions of hot springs displayed an apparent regional-specific pattern with distinct distributions of major and trace elements (e.g., Ca²⁺, Mg²⁺, F^-/B) and were primarily formed by water-rock interaction across the three hydrothermal fields. Nonetheless, microbial communities significantly assembled with the temperature rather than the geographic locations with distinct hydrogeological features. Low temperature (<45 °C), moderate temperature (55-70 °C) and high temperature (>70 °C) groups were identified based on their community compositions. Proteobacteria and Nitrospirae were the predominant phyla in low-temperature hot springs, while in moderate to high-temperature springs they were mainly composed of Aquificae, Deinococcus-Thermus, Thermodesulfobacteria, Thermotogae and Cyanobacteria. Variation partition analysis suggested a higher explanation of temperature (29.6%) than spatial variable (1.8%) and other geochemical variables (2.5%) on the microbial distribution. Microbial co-occurrence network showed >80% negative associations hinting a low co-existence pattern and highlighted the driving force of temperature as well as F^- or total organic carbon (TOC) for microbial interactions. Microbial dissimilarity displayed significant linear correlations with environmental (temperature) and geographic distance in Batang but only with temperature in Kangding area, which might be attributed to the regional-specific hydrogeochemistry. This study may help us to better understand the distribution of the microbial community in hot spring across different hydrothermal fields.

Comparative power law analysis for the spatial heterogeneity scaling of the hot-spring microbiomes

Spatial heterogeneity is a fundamental property of any natural ecosystems, including hot spring and human microbiomes. Two important scales that spatial heterogeneity exhibits are population and community scales, and Taylor's power law (PL) and its extensions (PLEs) offer ideal quantitative models to assess population- and community-level heterogeneities. Here we analyse 165 hot spring microbiome samples at the global scale that cover a wide range of temperatures (7.5-99°C) and pH levels (3.3-9). We explore a question of fundamental importance for measuring the spatial heterogeneity of the hot-spring microbiome and further discuss their ecological implications: How do critical environmental factors such as temperature and pH influence the scaling of community spatial heterogeneity? We are particularly interested in the existence of a universal scaling model that is independent of environmental gradients. By applying PL and PLEs, we were able to obtain such scaling parameters of the hot spring at both community and population levels, which are temperature- and pH-invariant. These findings suggest that while the hot-spring microbiomes located at different regions may have different environmental conditions, they share a fundamental heterogeneity scaling parameter, analogically similar to the gravitational acceleration on Earth, which may vary slightly depending on altitude and latitude, but is invariant overall. In contrast, similar to the physics of the Moon and Earth, which have different gravitational accelerations, the hot spring and human microbiomes can have different scaling parameters as demonstrated in this study.

High variations of methanogenic microorganisms drive full-scale anaerobic digestion process

Anaerobic digestion is one of the most successful waste management strategies worldwide, wherein microorganisms play an essential role in reducing organic pollutants and producing renewable energy. However, variations of microbial community in full-scale anaerobic digesters, particularly functional groups relevant to biogas production, remain elusive. Here, we examined microbial community in a year-long monthly time series of 3 full-scale anaerobic digesters. We observed substantial diversification in community composition, with only a few abundant OTUs (e.g. Clostridiales, Anaerolineaceae and Methanosaeta) persistently present across different samples. Similarly, there were high variations in relative abundance of methanogenic archaea and methanogenic genes, which were positively correlated (r² = 0.530, P < 0.001). Variations of methanogens explained 55.7% of biogas producing rates, much higher than the explanatory percentage of environmental parameters (16.4%). Hydrogenotrophic methanogens, especially abundant Methanomicrobiales taxa, were correlated with biogas production performance (r = 0.665, P < 0.001) and nearly all methanogenic genes (0.430 < r < 0.735, P < 0.012). Given that methanogenic archaea or genes are well established for methanogenesis, we conclude that high variations in methanogenic traits (e.g. taxa or genes) are responsible for biogas production variations in full-scale anaerobic digesters.

Seasonal microbial variation accounts for arsenic dynamics in shallow alluvial aquifer systems

Determining the temporal variation of microbial communities in groundwater systems is essential to improve our understanding of hydrochemical dynamics in aquifers, particularly as it relates to the fate of redox-sensitive contaminants like arsenic (As). Therefore, a high-resolution hydrobiogeochemical investigation was conducted in the As-affected alluvial aquifer systems of the Jianghan Plain. In two 25 m-deep monitoring wells, the seasonal variation in the composition of groundwater microbial communities was positively correlated with the change in groundwater level (R = 0.47 and 0.39 in NH03B and NH05B, respectively, P <  0.01), implying that the latter could be a primary driver of the seasonal microbial dynamics. In response to the fluctuating groundwater level, iron (Fe) reducers within the Desulfuromonadales were dominant (9.9 ± 4.7% among different sampling sites) in groundwater microbial communities during the monsoon season and associated with high concentrations of Fe(II) and As, while the predominance (16.7 ± 15.2% among different sampling sites) of iron-oxidizers the Gallionellaceae was accompanied by low Fe(II) and As in the non-monsoon season. These results suggest that microbially-mediated iron reduction/oxidation may have governed the seasonal mobilization/scavenging of As in groundwater. Our results provide new insights into mechanisms responsible for seasonal variations in groundwater As concentrations in similar aquifer systems.

Temperature-controlled thermophilic bacterial communities in hot springs of western Sichuan, China

Background

Ganzi Prefecture in Western China is situated geographically at the transition regions between Tibetan Plateau and Sichuan Basin in a highly tectonically active boundary area between the India and Eurasia plates. The region hosts various hot springs that span a wide range of temperature from 30 to 98 °C and are located at high altitude (up to 4200 m above sea level) in the region of large geothermal anomalies and active Xianshuihe slip-fault that has been active since Holocene. The site represents a biodiversity reservoir for thermophiles, yet their diversity and relationship to geochemical parameters are largely unknown. In the present work, bacterial diversity and community structure in 14 hot springs of Ganzi were investigated using Illumina MiSeq sequencing.

Results

Bacterial community compositions were evidently distinct among the 14 hot springs, and the bacterial communities in hot springs were majorly abundant in phyla Aquificae, Cyanobacteria and Proteobacteria. Both clustering and PCoA analysis suggested the existence of four bacterial community patterns in these hot springs. Temperature contributed to shaping bacterial community structure of hot springs as revealed by correlation analysis. Abundant unassigned-genus sequences detected in this study strongly implied the presence of novel genera or genetic resources in these hot springs.

Conclusion

The diversity of hot springs of Ganzi prefecture in Western Sichuan, China is evidently shaped by temperature. Interestingly disproportionally abundant unassigned-genus sequences detected in this study show indicate potential of novel genera or phylotypes. We hypothesize that frequent earthquakes and rapidly changing environment might have contributed to evolution of these potentially new lineages. Overall, this study provided first insight into the bacterial diversity of hot springs located in Western Sichuan, China and its comparison with other similar communities worldwide.

Electronic supplementary material

The online version of this article (10.1186/s12866-018-1271-z) contains supplementary material, which is available to authorized users.

Toward an Intensive Longitudinal Understanding of Activated Sludge Bacterial Assembly and Dynamics

Temporal microbial community studies have broadened our knowledge of the dynamics and correlations among microbes in both natural and artificial engineering systems. Using activated sludge as a model system, we utilized the intensive longitudinal sampling method to identify overlooked diversity and the hidden dynamics of microbes, detect cross-associations among microbes after detrending, and reveal the central microbial dynamics during sludge bulking and foaming. We discovered that the accumulative alpha diversity in activated sludge sampled daily over 392 days could be as high as 14 000 OTUs, and that the bacterial community dynamics followed a gradual succession, drifting away from the initial observed day and displaying a significant time-dependent trend. Cross-associations among bacteria were modulated after removing potential spurious correlations based on autocorrelation in microbial time series. Moreover, clusters of bacteria displaying rapid turnover were discovered during the beginning, ongoing, and fading of sludge bulking and foaming, and their physicochemical parameters are resolved. These identified groups of bacteria and their related environmental factors could potentially supply clues to form hypotheses for treating operational problems, such as sludge bulking and foaming.

Microbial biogeography of 925 geothermal springs in New Zealand

Geothermal springs are model ecosystems to investigate microbial biogeography as they represent discrete, relatively homogenous habitats, are distributed across multiple geographical scales, span broad geochemical gradients, and have reduced metazoan interactions. Here, we report the largest known consolidated study of geothermal ecosystems to determine factors that influence biogeographical patterns. We measured bacterial and archaeal community composition, 46 physicochemical parameters, and metadata from 925 geothermal springs across New Zealand (13.9–100.6 °C and pH < 1–9.7). We determined that diversity is primarily influenced by pH at temperatures <70 °C; with temperature only having a significant effect for values >70 °C. Further, community dissimilarity increases with geographic distance, with niche selection driving assembly at a localised scale. Surprisingly, two genera (Venenivibrio and Acidithiobacillus) dominated in both average relative abundance (11.2% and 11.1%, respectively) and prevalence (74.2% and 62.9%, respectively). These findings provide an unprecedented insight into ecological behaviour in geothermal springs, and a foundation to improve the characterisation of microbial biogeographical processes.

Power et al. catalogue the microbial biodiversity and physicochemistry of around 1000 hotsprings across New Zealand, providing insights into the ecological conditions that drive community assembly in these ecosystems.

Abundance and taxonomic affiliation of molybdenum transport and utilization genes in Tengchong hot springs, China

The nitrogen, sulfur and carbon cycles all rely on critical microbial transformations that are carried out by enzymes that require molybdenum (Mo) as a cofactor. Despite Mo importance in these biogeochemical cycles, little information exists about microbial Mo utilization in extreme environments where, due to geochemical conditions, bioavailable Mo may be limited. Using metagenomic data from nine hot springs in Tengchong, Yunnan Province, China, which range in temperature from 42°C to 96°C and pH from 2.3 to 9, the effects of pH, temperature and spring geochemistry on the abundance and taxonomic affiliation of genes related to Mo were studied. Dissolved Mo was only detected at sites with circumneutral pH. However, processes and organisms that require Mo were detected at all sites across all temperature and pH gradients. All sites contained xanthine dehydrogenase, formate dehydrogenase, carbon-monoxide dehydrogenase, nitrate reductase, sulfite oxidase and methionine-sulfoxide reductase despite different community compositions. This suggests that different microbial communities, resulting from different physicochemical conditions, may be performing similar metabolic functions. Furthermore, the abundance and taxonomic diversity of Mo-related annotations increased with higher concentrations of Mo. This study shows that despite geochemical conditions that can limit Mo bioavailability, microbes require Mo for a variety of processes.

Different Diversity and Distribution of Archaeal Community in the Aqueous and Oil Phases of Production Fluid From High-Temperature Petroleum Reservoirs

To get a better knowledge on how archaeal communities differ between the oil and aqueous phases and whether environmental factors promote substantial differences on microbial distributions among production wells, we analyzed archaeal communities in oil and aqueous phases from four high-temperature petroleum reservoirs (55–65°C) by using 16S rRNA gene based 454 pyrosequencing. Obvious dissimilarity of the archaeal composition between aqueous and oil phases in each independent production wells was observed, especially in production wells with higher water cut, and diversity in the oil phase was much higher than that in the corresponding aqueous phase. Statistical analysis further showed that archaeal communities in oil phases from different petroleum reservoirs tended to be more similar, but those in aqueous phases were the opposite. In the high-temperature ecosystems, temperature as an environmental factor could have significantly affected archaeal distribution, and archaeal diversity raised with the increase of temperature (p < 0.05). Our results suggest that to get a comprehensive understanding of petroleum reservoirs microbial information both in aqueous and oil phases should be taken into consideration. The microscopic habitats of oil phase, technically the dispersed minuscule water droplets in the oil could be a better habitat that containing the indigenous microorganisms.

Topographical diversity of common skin microflora and its association with skin environment type: An observational study in Chinese women

This study evaluated cutaneous microbial distribution, and microbial co-occurrence at different body sites and skin environments in Chinese women (39.6 ± 11.9 years, N = 100) during the winter season. Microbial distribution (Propionibacterium acnes, Staphylococcus aureus, Staphylococcus epidermidis, Lactobacillus, Pseudomonadaceae, and Malassezia furfur), association with biomarkers (antimicrobial peptides: LL-37, β-defensins [HBD-2, HBD-3]), and claudin-1) and skin biophysical parameters (transepidermal water loss, pH, skin scaliness and roughness, sebum and hydration levels) were also determined. Skin sites (glabella [GL], hand-back [HB], interdigital web-space [IS], antecubital fossa [AF], volar forearm [VF], back [BA]) were classified as normal, oily or dry based on two-step cluster analysis and exposed or unexposed (uncovered or covered by clothes, respectively) based on seasonal apparel. Pseudomonadaceae and Staphylococcus aureus had the highest and lowest detection rate respectively at all sites. Cluster analysis identified skin sites as ‘normal’ (HB, BA, AF, VF), ‘dry’ (IS) and ‘oily’ (GL). Bacterial alpha diversity was higher in exposed (HB, IS, and GL) compared with unexposed sites (BA, AF and VF). Co-occurrence of Staphylococcus aureus with any of the other five microorganisms was lower in dry and oily skin versus normal skin. Skin exposure, biophysical/barrier profile and biomarkers were found to be associated with bacterial distribution and co-occurrence.

Phylogenetic Analysis and Antimicrobial Profiles of Cultured Emerging Opportunistic Pathogens (Phyla Actinobacteria and Proteobacteria) Identified in Hot Springs

Hot spring water may harbour emerging waterborne opportunistic pathogens that can cause infections in humans. We have investigated the diversity and antimicrobial resistance of culturable emerging and opportunistic bacterial pathogens, in water and sediment of hot springs located in Limpopo, South Africa. Aerobic bacteria were cultured and identified using 16S ribosomal DNA (rDNA) gene sequencing. The presence of Legionella spp. was investigated using real-time polymerase chain reaction. Isolates were tested for resistance to ten antibiotics representing six different classes: β-lactam (carbenicillin), aminoglycosides (gentamycin, kanamycin, streptomycin), tetracycline, amphenicols (chloramphenicol, ceftriaxone), sulphonamides (co-trimoxazole) and quinolones (nalidixic acid, norfloxacin). Gram-positive Kocuria sp. and Arthrobacter sp. and gram-negative Cupriavidus sp., Ralstonia sp., Cronobacter sp., Tepidimonas sp., Hafnia sp. and Sphingomonas sp. were isolated, all recognised as emerging food-borne pathogens. Legionella spp. was not detected throughout the study. Isolates of Kocuria, Arthrobacter and Hafnia and an unknown species of the class Gammaproteobacteria were resistant to two antibiotics in different combinations of carbenicillin, ceftriaxone, nalidixic acid and chloramphenicol. Cronobacter sp. was sensitive to all ten antibiotics. This study suggests that hot springs are potential reservoirs for emerging opportunistic pathogens, including multiple antibiotic resistant strains, and highlights the presence of unknown populations of emerging and potential waterborne opportunistic pathogens in the environment.

Thioarsenate Formation Coupled with Anaerobic Arsenite Oxidation by a Sulfate-Reducing Bacterium Isolated from a Hot Spring

Thioarsenates are common arsenic species in sulfidic geothermal waters, yet little is known about their biogeochemical traits. In the present study, a novel sulfate-reducing bacterial strain Desulfotomaculum TC-1 was isolated from a sulfidic hot spring in Tengchong geothermal area, Yunnan Province, China. The arxA gene, encoding anaerobic arsenite oxidase, was successfully amplified from the genome of strain TC-1, indicating it has a potential ability to oxidize arsenite under anaerobic condition. In anaerobic arsenite oxidation experiments inoculated with strain TC-1, a small amount of arsenate was detected in the beginning but became undetectable over longer time. Thioarsenates (AsO4-xSx2- with x = 1-4) formed with mono-, di- and tri-thioarsenates being dominant forms. Tetrathioarsenate was only detectable at the end of the experiment. These results suggest that thermophilic microbes might be involved in the formation of thioarsenates and provide a possible explanation for the widespread distribution of thioarsenates in terrestrial geothermal environments.

Microbial communities and arsenic biogeochemistry at the outflow of an alkaline sulfide-rich hot spring

Alkaline sulfide-rich hot springs provide a unique environment for microbial community and arsenic (As) biogeochemistry. In this study, a representative alkaline sulfide-rich hot spring, Zimeiquan in the Tengchong geothermal area, was chosen to study arsenic geochemistry and microbial community using Illumina MiSeq sequencing. Over 0.26 million 16S rRNA sequence reads were obtained from 5-paired parallel water and sediment samples along the hot spring’s outflow channel. High ratios of As(V)/As_Sum (total combined arsenate and arsenite concentrations) (0.59–0.78), coupled with high sulfide (up to 5.87 mg/L), were present in the hot spring’s pools, which suggested As(III) oxidation occurred. Along the outflow channel, As_Sum increased from 5.45 to 13.86 μmol/L, and the combined sulfide and sulfate concentrations increased from 292.02 to 364.28 μmol/L. These increases were primarily attributed to thioarsenic transformation. Temperature, sulfide, As and dissolved oxygen significantly shaped the microbial communities between not only the pools and downstream samples, but also water and sediment samples. Results implied that the upstream Thermocrinis was responsible for the transformation of thioarsenic to As(III) and the downstream Thermus contributed to derived As(III) oxidation. This study improves our understanding of microbially-mediated As transformation in alkaline sulfide-rich hot springs.

Spatial isolation and environmental factors drive distinct bacterial and archaeal communities in different types of petroleum reservoirs in China

To investigate the spatial distribution of microbial communities and their drivers in petroleum reservoir environments, we performed pyrosequencing of microbial partial 16S rRNA, derived from 20 geographically separated water-flooding reservoirs, and two reservoirs that had not been flooded, in China. The results indicated that distinct underground microbial communities inhabited the different reservoirs. Compared with the bacteria, archaeal alpha-diversity was not strongly correlated with the environmental variables. The variation of the bacterial and archaeal community compositions was affected synthetically, by the mining patterns, spatial isolation, reservoir temperature, salinity and pH of the formation brine. The environmental factors explained 64.22% and 78.26% of the total variance for the bacterial and archaeal communities, respectively. Despite the diverse community compositions, shared populations (48 bacterial and 18 archaeal genera) were found and were dominant in most of the oilfields. Potential indigenous microorganisms, including Carboxydibrachium, Thermosinus, and Neptunomonas, were only detected in a reservoir that had not been flooded with water. This study indicates that: 1) the environmental variation drives distinct microbial communities in different reservoirs; 2) compared with the archaea, the bacterial communities were highly heterogeneous within and among the reservoirs; and 3) despite the community variation, some microorganisms are dominant in multiple petroleum reservoirs.

Microbial Community Structure and Arsenic Biogeochemistry in an Acid Vapor-Formed Spring in Tengchong Geothermal Area, China

Arsenic biogeochemistry has been studied extensively in acid sulfate-chloride hot springs, but not in acid sulfate hot springs with low chloride. In this study, Zhenzhuquan in Tengchong geothermal area, a representative acid sulfate hot spring with low chloride, was chosen to study arsenic geochemistry and microbial community structure using Illumina MiSeq sequencing. Over 0.3 million 16S rRNA sequence reads were obtained from 6-paired parallel water and sediment samples along its outflow channel. Arsenic oxidation occurred in the Zhenxhuquan pool, with distinctly high ratios of arsenate to total dissolved arsenic (0.73–0.86). Coupled with iron and sulfur oxidation along the outflow channel, arsenic accumulated in downstream sediments with concentrations up to 16.44 g/kg and appeared to significantly constrain their microbial community diversity. These oxidations might be correlated with the appearance of some putative functional microbial populations, such as Aquificae and Pseudomonas (arsenic oxidation), Sulfolobus (sulfur and iron oxidation), Metallosphaera and Acidicaldus (iron oxidation). Temperature, total organic carbon and dissolved oxygen significantly shaped the microbial community structure of upstream and downstream samples. In the upstream outflow channel region, most microbial populations were microaerophilic/anaerobic thermophiles and hyperthermophiles, such as Sulfolobus, Nocardia, Fervidicoccus, Delftia, and Ralstonia. In the downstream region, aerobic heterotrophic mesophiles and thermophiles were identified, including Ktedonobacteria, Acidicaldus, Chthonomonas and Sphingobacteria. A total of 72.41–95.91% unassigned-genus sequences were derived from the downstream high arsenic sediments 16S rRNA clone libraries. This study could enable us to achieve an integrated understanding on arsenic biogeochemistry in acid hot springs.

The impact of temperature on microbial diversity and AOA activity in the Tengchong Geothermal Field, China

Using a culture-independent method that combines CARD-FISH, qPCR and 16S rDNA, we investigated the abundance, community structure and diversity of microbes along a steep thermal gradient (50–90 °C) in the Tengchong Geothermal Field. We found that Bacteria and Archaea abundance changed markedly with temperature changes and that the number of cells was lowest at high temperatures (90.8 °C). Under low-temperature conditions (52.3–74.6 °C), the microbial communities were dominated by Bacteria, which accounted for 60–80% of the total number of cells. At 74.6 °C, Archaea were dominant, and at 90.8 °C, they accounted for more than 90% of the total number of cells. Additionally, the microbial communities at high temperatures (74.6–90.8 °C) were substantially simpler than those at the low-temperature sites. Only a few genera (e.g., bacterial Caldisericum, Thermotoga and Thermoanaerobacter, archaeal Vulcanisaeta and Hyperthermus) often dominated in high-temperature environments. Additionally, a positive correlation between Ammonia-Oxidizing Archaea (AOA) activity and temperature was detected. AOA activity increased from 17 to 52 pmol of NO₂⁻ per cell d⁻¹ with a temperature change from 50 to 70 °C.

Carbon source preference in chemosynthetic hot spring communities

Rates of dissolved inorganic carbon (DIC), formate, and acetate mineralization and/or assimilation were determined in 13 high-temperature (>73 °C) hot springs in Yellowstone National Park (YNP), Wyoming, in order to evaluate the relative importance of these substrates in supporting microbial metabolism. While 9 of the hot spring communities exhibited rates of DIC assimilation that were greater than those of formate and acetate assimilation, 2 exhibited rates of formate and/or acetate assimilation that exceeded those of DIC assimilation. Overall rates of DIC, formate, and acetate mineralization and assimilation were positively correlated with spring pH but showed little correlation with temperature. Communities sampled from hot springs with similar geochemistries generally exhibited similar rates of substrate transformation, as well as similar community compositions, as revealed by 16S rRNA gene-tagged sequencing. Amendment of microcosms with small (micromolar) amounts of formate suppressed DIC assimilation in short-term (<45-min) incubations, despite the presence of native DIC concentrations that exceeded those of added formate by 2 to 3 orders of magnitude. The concentration of added formate required to suppress DIC assimilation was similar to the affinity constant (K(m)) for formate transformation, as determined by community kinetic assays. These results suggest that dominant chemoautotrophs in high-temperature communities are facultatively autotrophic or mixotrophic, are adapted to fluctuating nutrient availabilities, and are capable of taking advantage of energy-rich organic substrates when they become available.

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

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

Greater temporal changes of sediment microbial community than its waterborne counterpart in Tengchong hot springs, Yunnan Province, China

Temporal variation in geochemistry can cause changes in microbial community structure and diversity. Here we studied temporal changes of microbial communities in Tengchong hot springs of Yunnan Province, China in response to geochemical variations by using microbial and geochemical data collected in January, June and August of 2011. Greater temporal variations were observed in individual taxa than at the whole community structure level. Water and sediment communities exhibited different temporal variation patterns. Water communities were largely stable across three sampling times and dominated by similar microbial lineages: Hydrogenobaculum in moderate-temperature acidic springs, Sulfolobus in high-temperature acidic springs, and Hydrogenobacter in high-temperature circumneutral to alkaline springs. Sediment communities were more diverse and responsive to changing physicochemical conditions. Most of the sediment communities in January and June were similar to those in waters. However, the August sediment community was more diverse and contained more anaerobic heterotrophs than the January and June: Desulfurella and Acidicaldus in moderate-temperature acidic springs, Ignisphaera and Desulfurococcus in high-temperature acidic springs, the candidate division OP1 and Fervidobacterium in alkaline springs, and Thermus and GAL35 in neutral springs. Temporal variations in physicochemical parameters including temperature, pH, and dissolved organic carbon may have triggered the observed microbial community shifts.

Distribution of ether lipids and composition of the archaeal community in terrestrial geothermal springs: impact of environmental variables

Archaea can respond to changes in the environment by altering the composition of their membrane lipids, for example, by modification of the abundance and composition of glycerol dialkyl glycerol tetraethers (GDGTs). Here, we investigated the abundance and proportions of polar GDGTs (P-GDGTs) and core GDGTs (C-GDGTs) sampled in different seasons from Tengchong hot springs (Yunnan, China), which encompassed a pH range of 2.5-10.1 and a temperature range of 43.7-93.6°C. The phylogenetic composition of the archaeal community (reanalysed from published work) divided the Archaea in spring sediment samples into three major groups that corresponded with spring pH: acidic, circumneutral and alkaline. Cluster analysis showed correlation between spring pH and the composition of P- and C-GDGTs and archaeal 16S rRNA genes, indicating an intimate link between resident Archaea and the distribution of P- and C-GDGTs in Tengchong hot springs. The distribution of GDGTs in Tengchong springs was also significantly affected by temperature; however, the relationship was weaker than with pH. Analysis of published datasets including samples from Tibet, Yellowstone and the US Great Basin hot springs revealed a similar relationship between pH and GDGT content. Specifically, low pH springs had higher concentrations of GDGTs with high numbers of cyclopentyl rings than neutral and alkaline springs, which is consistent with the predominance of high cyclopentyl ring-characterized Sulfolobales and Thermoplasmatales present in some of the low pH springs. Our study suggests that the resident Archaea in these hot springs are acclimated if not adapted to low pH by their genetic capacity to effect the packing density of their membranes by increasing cyclopentyl rings in GDGTs at the rank of community.

Microbial community structures of novel Icelandic hot spring systems revealed by PhyloChip G3 analysis

Microbial community profiles of recently formed hot spring systems ranging in temperatures from 57°C to 100°C and pH values from 2 to 4 in Hveragerði (Iceland) were analyzed with PhyloChip G3 technology. In total, 1173 bacterial operational taxonomic units (OTUs) spanning 576 subfamilies and 38 archaeal OTUs covering 32 subfamilies were observed. As expected, the hyperthermophilic (∼100°C) spring system exhibited both low microbial biomass and diversity when compared to thermophilic (∼ 60°C) springs. Ordination analysis revealed distinct bacterial and archaeal diversity in geographically distinct hot springs. Slight variations in temperature (from 57°C to 64°C) within the interconnected pools led to a marked fluctuation in microbial abundance and diversity. Correlation and PERMANOVA tests provided evidence that temperature was the key environmental factor responsible for microbial community dynamics, while pH, H2S, and SO2 influenced the abundance of specific microbial groups. When archaeal community composition was analyzed, the majority of detected OTUs correlated negatively with temperature, and few correlated positively with pH.