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

Underground fires shape the structure of microbial communities and select for thermophilic bacteria through a temperature gradient

A detailed diversity analysis of the prokaryotic and fungal communities in soil impacted by an underground fire located in the Trans-Mexican volcanic belt, Mexico, is described. Microbial diversity data obtained from soils at different depths and temperatures (27 °C, 42 °C, 50 ºC and 54 ºC) were analyzed, and Firmicutes increased in abundance as the temperature augmented, and Proteobacteria mainly decreased in abundance at high temperatures compared to unaffected soils. The fungal phylum Ascomycota was the most abundant, with no significant changes. A clear reduction in the richness of both prokaryotic and eukaryotic operational taxonomic units (OTUs) was observed in the affected soils. At the genus level, Bacillus species were the most abundant among bacteria, while Aspergillus, Penicillium, and Mortierella were dominant fungal genera at higher temperatures. Interestingly, the physicochemical parameters of the affected soils modified organic matter, which was indirectly correlated with the presence of some microbial taxa. Likewise, we obtained 308 soil bacterial isolates from both control and affected soils. Among these, the taxa from the phyla Actinobacteria and Firmicutes demonstrated the highest thermotolerance in the affected soils. Our findings shed light on the impact of underground fires on the structure of microbial communities, favoring an abundance of thermotolerant microbes.

Recovery of nearly 3,000 archaeal genomes from 152 terrestrial geothermal spring metagenomes

Terrestrial geothermal springs, reminiscent of early Earth conditions, host diverse and abundant populations of Archaea. In this study, we reconstructed 2,949 metagenome-assembled genomes (MAGs) from 152 metagenomes collected over six years from 48 geothermal springs in Tengchong, China. Among these MAGs, 1,431 (49%) were classified as high-quality, while 1,518 (51%) were considered as medium-quality. Phylogenomic analysis revealed that these MAGs spanned 12 phyla, 27 classes, 67 orders, 147 families, 265 genera, and 475 species. Notably, 575 (19%) MAGs represented new taxa at various taxonomic levels, and 2,075 (70%) lacked nomenclature and effective descriptions. The most abundant phyla of archaeal genomes were Thermoproteota, Thermoplasmatota, and Micrarchaeota. The DRTY, ZMQ, and ZZQ geothermal springs were predominated by Archaea, particularly by Thermoproteia and Thermoplasmata. These draft genomes provide new data for studying species diversity and function within terrestrial geothermal spring archaeal communities, thus contributing to the conservation and utilization of thermophilic and hyperthermophilic microbial resources.

Subsurface microbial community structure shifts along the geological features of the Central American Volcanic Arc

Subduction of the Cocos and Nazca oceanic plates beneath the Caribbean plate drives the upward movement of deep fluids enriched in carbon, nitrogen, sulfur, and iron along the Central American Volcanic Arc (CAVA). These compounds fuel diverse subsurface microbial communities that in turn alter the distribution, redox state, and isotopic composition of these compounds. Microbial community structure and functions vary according to deep fluid delivery across the arc, but less is known about how microbial communities differ along the axis of a convergent margin as geological features (e.g., extent of volcanism and subduction geometry) shift. Here, we investigate changes in bacterial 16S rRNA gene amplicons and geochemical analysis of deeply-sourced seeps along the southern CAVA, where subduction of the Cocos Ridge alters the geological setting. We find shifts in community composition along the convergent margin, with communities in similar geological settings clustering together independently of the proximity of sample sites. Microbial community composition correlates with geological variables such as host rock type, maturity of hydrothermal fluid and slab depth along different segments of the CAVA. This reveals tight coupling between deep Earth processes and subsurface microbial activity, controlling community distribution, structure and composition along a convergent margin.

Bacterial diversity along the geothermal gradients: insights from the high-altitude Himalayan hot spring habitats of Sikkim

Geothermal habitats present a unique opportunity to study microbial adaptation to varying temperature conditions. In such environments, distinct temperature gradients foster diverse microbial communities, each adapted to its optimal niche. However, the complex dynamics of bacterial populations in across these gradients high-altitude hot springs remain largely unexplored. We hypothesize that temperature is a primary driver of microbial diversity, and bacterial richness peaks at intermediate temperatures. To investigate this, we analysed bacterial diversity using 16S rRNA amplicon sequencing across three temperature regions: hot region of 56-65 °C (hot spring), warm region of 35-37 °C (path carrying hot spring water to the river), and cold region of 4-7 °C (river basin). Our findings showed that Bacillota was the most abundant phylum (45.51 %), followed by Pseudomonadota (32.81 %) and Actinomycetota (7.2 %). Bacillota and Chloroflexota flourished in the hot and warm regions, while Pseudomonadota thrived in cooler areas. Core microbiome analysis indicated that species richness was highest in the warm region, declining in both cold and hot regions. Interestingly, an anomaly was observed with Staphylococcus, which was more abundant in cases where ponds were used for bathing and recreation. In contrast, Clostridium was mostly found in cold regions, likely due to its viability in soil and ability to remain dormant as a spore-forming bacterium. The warm region showed the highest bacterial diversity, while richness decreased in both cold and hot regions. This highlights the temperature-dependent nature of microbial communities, with optimal diversity in moderate thermal conditions. The study offers new insights into microbial dynamics in high-altitude geothermal systems.

Aquificae overcomes competition by archaeal thermophiles, and crowding by bacterial mesophiles, to dominate the boiling vent-water of a Trans-Himalayan sulfur-borax spring

Trans-Himalayan hot spring waters rich in boron, chlorine, sodium and sulfur (but poor in calcium and silicon) are known based on PCR-amplified 16S rRNA gene sequence data to harbor high diversities of infiltrating bacterial mesophiles. Yet, little is known about the community structure and functions, primary productivity, mutual interactions, and thermal adaptations of the microorganisms present in the steaming waters discharged by these geochemically peculiar spring systems. We revealed these aspects of a bacteria-dominated microbiome (microbial cell density ~8.5 × 104 mL-1; live:dead cell ratio 1.7) thriving in the boiling (85°C) fluid vented by a sulfur-borax spring called Lotus Pond, situated at 4436 m above the mean sea-level, in the Puga valley of eastern Ladakh, on the Changthang plateau. Assembly, annotation, and population-binning of >15-GB metagenomic sequence illuminated the numeral predominance of Aquificae. While members of this phylum accounted for 80% of all 16S rRNA-encoding reads within the metagenomic dataset, 14% of such reads were attributed to Proteobacteria. Post assembly, only 25% of all protein-coding genes identified were attributable to Aquificae, whereas 41% was ascribed to Proteobacteria. Annotation of metagenomic reads encoding 16S rRNAs, and/or PCR-amplified 16S rRNA genes, identified 163 bacterial genera, out of which 66 had been detected in past investigations of Lotus Pond's vent-water via 16S amplicon sequencing. Among these 66, Fervidobacterium, Halomonas, Hydrogenobacter, Paracoccus, Sulfurihydrogenibium, Tepidimonas, Thermus and Thiofaba (or their close phylogenomic relatives) were presently detected as metagenome-assembled genomes (MAGs). Remarkably, the Hydrogenobacter related MAG alone accounted for ~56% of the entire metagenome, even though only 15 out of the 66 genera consistently present in Lotus Pond's vent-water have strains growing in the laboratory at >45°C, reflecting the continued existence of the mesophiles in the ecosystem. Furthermore, the metagenome was replete with genes crucial for thermal adaptation in the context of Lotus Pond's geochemistry and topography. In terms of sequence similarity, a majority of those genes were attributable to phylogenetic relatives of mesophilic bacteria, while functionally they rendered functions such as encoding heat shock proteins, molecular chaperones, and chaperonin complexes; proteins controlling/modulating/inhibiting DNA gyrase; universal stress proteins; methionine sulfoxide reductases; fatty acid desaturases; different toxin-antitoxin systems; enzymes protecting against oxidative damage; proteins conferring flagellar structure/function, chemotaxis, cell adhesion/aggregation, biofilm formation, and quorum sensing. The Lotus Pond Aquificae not only dominated the microbiome numerically but also acted potentially as the main primary producers of the ecosystem, with chemolithotrophic sulfur oxidation (Sox) being the fundamental bioenergetic mechanism, and reductive tricarboxylic acid (rTCA) cycle the predominant carbon fixation pathway. The Lotus Pond metagenome contained several genes directly or indirectly related to virulence functions, biosynthesis of secondary metabolites including antibiotics, antibiotic resistance, and multi-drug efflux pumping. A large proportion of these genes being attributable to Aquificae, and Proteobacteria (very few were ascribed to Archaea), it could be worth exploring in the future whether antibiosis helped the Aquificae overcome niche overlap with other thermophiles (especially those belonging to Archaea), besides exacerbating the bioenergetic costs of thermal endurance for the mesophilic intruders of the ecosystem.

Metagenomics insight into Puga geothermal geyser located in Himalayan Geothermal Belt (Trans-Himalayan Plateau) Ladakh, India

Puga geothermal geyser and surrounding area, located in the Himalayan Geothermal Belt of the Trans-Himalayan Plateau in Ladakh, India, are very geographically isolated and considered pristine and free of anthropogenic activities. In this study, we have conducted the first metagenomic investigation of the microbes in and around the geyser. The whole genome sequencing analysis showed the presence of a total of 44.8%, 39.7% and 41.4% bacterial phyla in the PugW, PugS, and PugSo samples respectively, 8.6% of archaeal phyla (in all the samples), unclassified (derived from other sequences, PugW: 27.6%, PugS: 27.6%, and PugSo: 15.5%) and unclassified (derived from bacteria, PugW: 12%, PugS: 13.8%, and PugSo: 13.8%). The majority of archaeal sequences were linked to Euryarchaeota (2.84%) while the majority of the bacterial communities that predominated in most geothermal locations were linked to Pseudomonadota (67.14%) and Bacteroidota (12.52%). The abundant bacterial strains at the species level included Dechloromonas aromatica, Acinetobacter baumannii, and Arcobacter butzleri, in all the samples while the most abundant archaeal species were Methanosaeta thermophile, Methanoregula boonei, and Methanosarcina berkeri. Further, this geothermal geyser metagenome has a large number of unique sequences linked to unidentified and unclassified lineages, suggesting a potential source for novel species of microbes and their products. The present study which only examined one of the many geothermal geysers and springs in the Puga geothermal area, should be regarded as a preliminary investigation of the microbiota that live in the geothermal springs on these remote areas. These findings suggest that further investigations should be undertaken to characterize the ecosystems of the Puga geothermal area, which serve as a repository for unidentified microbial lineages.

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.

A genome and gene catalog of the aquatic microbiomes of the Tibetan Plateau

The Tibetan Plateau supplies water to nearly 2 billion people in Asia, but climate change poses threats to its aquatic microbial resources. Here, we construct the Tibetan Plateau Microbial Catalog by sequencing 498 metagenomes from six water ecosystems (saline lakes, freshwater lakes, rivers, hot springs, wetlands and glaciers). Our catalog expands knowledge of regional genomic diversity by presenting 32,355 metagenome-assembled genomes that de-replicated into 10,723 representative genome-based species, of which 88% were unannotated. The catalog contains nearly 300 million non-redundant gene clusters, of which 15% novel, and 73,864 biosynthetic gene clusters, of which 50% novel, thus expanding known functional diversity. Using these data, we investigate the Tibetan Plateau aquatic microbiome's biogeography along a distance of 2,500 km and >5 km in altitude. Microbial compositional similarity and the shared gene count with the Tibetan Plateau microbiome decline along with distance and altitude difference, suggesting a dispersal pattern. The Tibetan Plateau Microbial Catalog stands as a substantial repository for high-altitude aquatic microbiome resources, providing potential for discovering novel lineages and functions, and bridging knowledge gaps in microbiome biogeography.

Microbial Diversity and Activity of Biofilms from Geothermal Springs in Croatia

Hot spring biofilms are stable, highly complex microbial structures. They form at dynamic redox and light gradients and are composed of microorganisms adapted to the extreme temperatures and fluctuating geochemical conditions of geothermal environments. In Croatia, a large number of poorly investigated geothermal springs host biofilm communities. Here, we investigated the microbial community composition of biofilms collected over several seasons at 12 geothermal springs and wells. We found biofilm microbial communities to be temporally stable and highly dominated by Cyanobacteria in all but one high-temperature sampling site (Bizovac well). Of the physiochemical parameters recorded, temperature had the strongest influence on biofilm microbial community composition. Besides Cyanobacteria, the biofilms were mainly inhabited by Chloroflexota, Gammaproteobacteria, and Bacteroidota. In a series of incubations with Cyanobacteria-dominated biofilms from Tuhelj spring and Chloroflexota- and Pseudomonadota-dominated biofilms from Bizovac well, we stimulated either chemoorganotrophic or chemolithotrophic community members, to determine the fraction of microorganisms dependent on organic carbon (in situ predominantly produced via photosynthesis) versus energy derived from geochemical redox gradients (here simulated by addition of thiosulfate). We found surprisingly similar levels of activity in response to all substrates in these two distinct biofilm communities, and observed microbial community composition and hot spring geochemistry to be poor predictors of microbial activity in the study systems.

Effects of hydrogeochemistry on the microbial ecology of terrestrial hot springs

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

Characterization of Hsp17, a Novel Small Heat Shock Protein, in Sphingomonas melonis TY under Heat Stress

Bacteria are constantly exposed to a variety of environmental stresses. Temperature is considered one of the most important environmental factors affecting microbial growth and survival. As ubiquitous environmental microorganisms, Sphingomonas species play essential roles in the biodegradation of organic contaminants, plant protection, and environmental remediation. Understanding the mechanism by which they respond to heat shock will help further improve cell resistance by applying synthetic biological strategies. Here, we assessed the transcriptomic and proteomic responses of Sphingomonas melonis TY to heat shock and found that stressful conditions caused significant changes in functional genes related to protein synthesis at the transcriptional level. The most notable changes observed were increases in the transcription (1,857-fold) and protein expression (11-fold) of Hsp17, which belongs to the small heat shock protein family, and the function of Hsp17 in heat stress was further investigated in this study. We found that the deletion of hsp17 reduced the capacity of the cells to tolerate high temperatures, whereas the overexpression of hsp17 significantly enhanced the ability of the cells to withstand high temperatures. Moreover, the heterologous expression of hsp17 in Escherichia coli DH5α conferred to the bacterium the ability to resist heat stress. Interestingly, its cells were elongated and formed connected cells following the increase in temperature, while hsp17 overexpression restored their normal morphology under high temperature. In general, these results indicate that the novel small heat shock protein Hsp17 greatly contributes to maintaining cell viability and morphology under stress conditions. IMPORTANCE Temperature is generally considered the most important factor affecting metabolic functions and the survival of microbes. As molecular chaperones, small heat shock proteins can prevent damaged protein aggregation during abiotic stress, especially heat stress. Sphingomonas species are widely distributed in nature, and they can frequently be found in various extreme environments. However, the role of small heat shock proteins in Sphingomonas under high-temperature stress has not been elucidated. This study greatly enhances our understanding of a novel identified protein, Hsp17, in S. melonis TY in terms of its ability to resist heat stress and maintain cell morphology under high temperature, leading to a broader understanding of how microbes adapt to environmental extremes. Furthermore, our study will provide potential heat resistance elements for further enhancing cellular resistance as well as the synthetic biological applications of Sphingomonas.

Community- and genome-based evidence for a shaping influence of redox potential on bacterial protein evolution

Despite deep interest in how environments shape microbial communities, whether redox conditions influence the sequence composition of genomes is not well known. We predicted that the carbon oxidation state (ZC) of protein sequences would be positively correlated with redox potential (Eh). To test this prediction, we used taxonomic classifications for 68 publicly available 16S rRNA gene sequence data sets to estimate the abundances of archaeal and bacterial genomes in river & seawater, lake & pond, geothermal, hyperalkaline, groundwater, sediment, and soil environments. Locally, ZC of community reference proteomes (i.e., all the protein sequences in each genome, weighted by taxonomic abundances but not by protein abundances) is positively correlated with Eh corrected to pH 7 (Eh7) for the majority of data sets for bacterial communities in each type of environment, and global-scale correlations are positive for bacterial communities in all environments. In contrast, archaeal communities show approximately equal frequencies of positive and negative correlations in individual data sets, and a positive pan-environmental correlation for archaea only emerges after limiting the analysis to samples with reported oxygen concentrations. These results provide empirical evidence that geochemistry modulates genome evolution and may have distinct effects on bacteria and archaea. IMPORTANCE The identification of environmental factors that influence the elemental composition of proteins has implications for understanding microbial evolution and biogeography. Millions of years of genome evolution may provide a route for protein sequences to attain incomplete equilibrium with their chemical environment. We developed new tests of this chemical adaptation hypothesis by analyzing trends of the carbon oxidation state of community reference proteomes for microbial communities in local- and global-scale redox gradients. The results provide evidence for widespread environmental shaping of the elemental composition of protein sequences at the community level and establish a rationale for using thermodynamic models as a window into geochemical effects on microbial community assembly and evolution.

High speciation rate of niche specialists in hot springs

Ecological and evolutionary processes simultaneously regulate microbial diversity, but the evolutionary processes and their driving forces remain largely unexplored. Here we investigated the ecological and evolutionary characteristics of microbiota in hot springs spanning a broad temperature range (54.8-80 °C) by sequencing the 16S rRNA genes. Our results demonstrated that niche specialists and niche generalists are embedded in a complex interaction of ecological and evolutionary dynamics. On the thermal tolerance niche axis, thermal (T) sensitive (at a specific temperature) versus T-resistant (at least in five temperatures) species were characterized by different niche breadth, community abundance and dispersal potential, consequently differing in potential evolutionary trajectory. The niche-specialized T-sensitive species experienced strong temperature barriers, leading to completely species shift and high fitness but low abundant communities at each temperature ("home niche"), and such trade-offs thus reinforced peak performance, as evidenced by high speciation across temperatures and increasing diversification potential with temperature. In contrast, T-resistant species are advantageous of niche expansion but with poor local performance, as shown by wide niche breadth with high extinction, indicating these niche generalists are "jack-of-all-trades, master-of-none". Despite of such differences, the T-sensitive and T-resistant species are evolutionarily interacted. Specifically, the continuous transition from T-sensitive to T-resistant species insured the exclusion probability of T-resistant species at a relatively constant level across temperatures. The co-evolution and co-adaptation of T-sensitive and T-resistant species were in line with the red queen theory. Collectively, our findings demonstrate that high speciation of niche specialists could alleviate the environmental-filtering-induced negative effect on diversity.

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

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

Effects of Mineral on Taxonomic and Functional Structures of Microbial Community in Tengchong Hot Springs via in-situ cultivation

Diverse mineralogical compositions occur in hot spring sediments, but the impact of minerals on the diversity and structure of microbial communities remains poorly elucidated. In this study, different mineral particles with various chemistries (i.e., hematite, biotite, K-feldspar, quartz, muscovite, aragonite, serpentine, olivine, barite, apatite, and pyrite) were incubated for ten days in two Tengchong hot springs, one alkaline (pH ~ 8.34) with a high temperature (~ 82.8 °C) (Gumingquan, short as GMQ) and one acidic (pH ~ 3.63) with a relatively low temperature (~ 43.3 °C) (Wenguangting, short as WGT), to determine the impacts of minerals on the microbial communities taxonomic and functional diversities. Results showed that the mineral-associated bacterial taxa differed from those of the bulk sediment samples in the two hot springs. The relative abundance of Proteobacteria, Euryarchaeota, and Acidobacteria increased in all minerals, indicating that these microorganisms are apt to colonize on solid surfaces. The α-diversity indices of the microbial communities on the mineral surfaces in the WGT were higher than those from the bulk sediment samples (p < 0.05), which may be caused by the stochastically adhering process on the mineral surface during 10-day incubation, different from the microbial community in sediment which has experienced long-term environmental and ecological screening. Chemoheterotrophy increased with minerals incubation, which was high in most cultured minerals (the relative contents were 5.8 - 21.4%). Most notably, the sulfate respiration bacteria (mainly related to Desulfobulbaceae and Syntrophaceae) associated with aragonite in the acidic hot spring significantly differed from other minerals, possibly due to the pH buffering effect of aragonite providing more favorable conditions for their survival and proliferation. By comparison, aragonite cultured in the alkaline hot spring highly enriched denitrifying bacteria and may have promoted the nitrogen cycle within the system. Collectively, we speculated that diverse microbes stochastically adhered on the surface of minerals in the water flows, and the physicochemical properties of minerals drove the enrichment of certain microbial communities and functional groups during the short-term incubation. Taken together, these findings thereby provide novel insights into mechanisms of community assembly and element cycling in the terrestrial hydrothermal system associated with hot springs.

Distinct distribution patterns of the abundant and rare bacteria in high plateau hot spring sediments

The diversity and distribution patterns of the abundant and rare microbial sub-communities in hot spring ecosystems and their assembly mechanisms are poorly understood. The present study investigated the diversity and distribution patterns of the total, abundant, conditionally rare, and always rare taxa in the low- and moderate-temperature hot spring sediments on the Tibetan Plateau based on high-throughput 16S rRNA gene sequencing, and explored their major environmental drivers. The diversity of these four bacterial taxa showed no significant change between the low-temperature and moderate-temperature hot spring sediments, whereas the bacterial compositions were obviously different. Stochasticity dominated the bacterial sub-community assemblages, while heterogeneous selection also played an important role in shaping the abundant and conditionally rare taxa between the low-temperature and moderate-temperature hot spring sediments. No significant difference in the topological properties of co-occurrence networks was found between the conditionally rare and abundant taxa, and the connections between the paired operational taxonomic units (OTUs) were almost positive. The diversity of the total, abundant, and conditionally rare taxa was governed by the salinity of hot spring sediments, while that of the always rare taxa was determined by the content of S element. In contrast, temperature had significant direct effect on the composition of the total, abundant, and conditionally rare taxa, but relatively weak influence on that of the always rare taxa. Besides, salinity was another major environmental factor driving the composition of the abundant and rare sub-communities in the hot spring sediments. These results reveal the assembly processes and major environmental drivers that shaped different bacterial sub-communities in the hot spring sediments on the Tibetan Plateau, and indicate the importance of conditionally rare taxa in constructing bacterial communities. These findings enhance the current understanding of the ecological mechanisms maintaining the ecosystem stability and services in extreme environment.

High Variation in Protist Diversity and Community Composition in Surface Sediment of Hot Springs in Himalayan Geothermal Belt, China

Hot springs are some of the most special environments on Earth. Many prokaryotic and eukaryotic microbes have been found to live in this environment. The Himalayan geothermal belt (HGB) has numerous hot springs spread across the area. Comprehensive research using molecular techniques to investigate eukaryotic microorganisms is still lacking; investigating the composition and diversity of eukaryotic microorganisms such as protists in the hot spring ecosystems will not only provide critical information on the adaptations of protists to extreme conditions, but could also give valuable contributions to the global knowledge of biogeographic diversity. In this study, we used high-throughput sequencing to illuminate the diversity and composition pattern of protist communities in 41 geothermal springs across the HGB on the Tibetan Plateau. A total of 1238 amplicon sequence variants (ASVs) of protists were identified in the hot springs of the HGB. In general, Cercozoa was the phylum with the highest richness, and Bacillariophyta was the phylum with the highest relative abundance in protists. Based on the occurrence of protist ASVs, most of them are rare. A high variation in protist diversity was found in the hot springs of the HGB. The high variation in protist diversity may be due to the different in environmental conditions of these hot springs. Temperature, salinity, and pH are the most important environmental factors that affect the protist communities in the surface sediments of the hot springs in the HGB. In summary, this study provides the first comprehensive study of the composition and diversity of protists in the hot springs of the HGB and facilitates our understanding of the adaptation of protists in these extreme habitats.

Associations between temperature variations and tourist arrivals: analysis based on Baidu Index of hot-spring tourism in 44 cities in China

Hot-spring tourism refers to entertainment, health preservation, commercial conferences, and other leisure activities at hot-spring locations. This tourism often shows periodic variability, which may be due to seasonal temperature variations. However, studies on the effects of temperature variations on tourist arrivals at hot springs are limited. Therefore, this study aimed to evaluate this relationship in 31 provincial capital cities and 13 s-tier cities in China. Using the Baidu Index, we obtained data for tourist arrivals to hot springs in each city and constructed a generalised additive model to explore the associations between temperature variations and tourist arrivals. We also analysed the statistical significance of the estimated effects during different seasons to explore potential effect modification. A 1 °C increase in temperature was associated with a 1.81% (95% confidence interval (CI): 1.69-1.93) decrease in daily tourist arrivals for hot-spring tourism. Significant positive associations between the abovementioned factors were observed in summer (2.18% change, 95% CI: 1.32-3.04). The effect of temperature on the volume of tourist arrivals may last for approximately 2 months. Robustness analysis confirmed the data reliability. The results indicate that significant relationships exist between temperature variations and hot-spring tourism arrivals, which vary seasonally. This study has significant implications for travel agencies to effectively manage tourist visits to hot spring locations.

Bacterial community structure in geothermal springs on the northern edge of Qinghai-Tibet plateau

Introduction:

In order to reveal the composition of the subsurface hydrothermal bacterial community in the zones of magmatic tectonics and their response to heat storage environments.

Methods:

In this study, we performed hydrochemical analysis and regional sequencing of the 16S rRNA microbial V4-V5 region in 7 Pleistocene and Lower Neogene hot water samples from the Gonghe basin.

Results:

Two geothermal hot spring reservoirs in the study area were found to be alkaline reducing environments with a mean temperature of 24.83°C and 69.28°C, respectively, and the major type of hydrochemistry was SO4-Cl·Na. The composition and structure of microorganisms in both types of geologic thermal storage were primarily controlled by temperature, reducing environment intensity, and hydrogeochemical processes. Only 195 ASVs were shared across different temperature environments, and the dominant bacterial genera in recent samples from temperate hot springs were Thermus and Hydrogenobacter, with both genera being typical of thermophiles. The correlation analysis showed that the overall level of relative abundance of the subsurface hot spring relied on a high temperature and a slightly alkaline reducing environment. Nearly all of the top 4 species in the abundance level (53.99% of total abundance) were positively correlated with temperature and pH, whereas they were negatively correlated with ORP (oxidation–reduction potential), nitrate, and bromine ions.

Discussion:

In general, the composition of bacteria in the groundwater in the study area was sensitive to the response of the thermal storage environment and also showed a relationship with geochemical processes, such as gypsum dissolution, mineral oxidation, etc.

Relationships between fluid mixing, biodiversity, and chemosynthetic primary productivity in Yellowstone hot springs

The factors that influence biodiversity and productivity of hydrothermal ecosystems are not well understood. Here we investigate the relationship between fluid mixing, biodiversity, and chemosynthetic primary productivity in three co-localized hot springs (RSW, RSN, and RSE) in Yellowstone National Park that have different geochemistry. All three springs are sourced by reduced hydrothermal fluid, but RSE and RSN receive input of vapour phase gas and oxidized groundwaters, with input of both being substantially higher in RSN. Metagenomic sequencing revealed that communities in RSN were more biodiverse than those of RSE and RSW in all dimensions evaluated. Microcosm activity assays indicate that rates of dissolved inorganic carbon (DIC) uptake were also higher in RSN than in RSE and RSW. Together, these results suggest that increased mixing of reduced volcanic fluid with oxidized fluids generates additional niche space capable of supporting increasingly biodiverse communities that are more productive. These results provide insight into the factors that generate and maintain chemosynthetic biodiversity in hydrothermal systems and that influence the distribution, abundance, and diversity of microbial life in communities supported by chemosynthesis. These factors may also extend to other ecosystems not supported by photosynthesis, including the vast subterranean biosphere and biospheres beneath ice sheets and glaciers.

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.

Post-monsoon seasonal variation of prokaryotic diversity in solfataric soil from the North Sikkim hot spring

The solfataric soil sediments of the hot springs of Sikkim located at Yume Samdung and Lachen valley were studied for deciphering the bacterial diversity. The main aim here is to present a comparative study and generate a baseline data on the post-monsoon seasonal variation for the months of October and December, analyzed through 16S rRNA V3-V4 amplicon sequencing. The results have shown that there is not much variation at phylum level in the month of October in all the three hot springs such as New Yume Samdung (NYS), Old Yume Samdung (OYS), and Tarum (TAR) hot spring. The abundant phyla mainly present were Firmicutes, followed by Proteobacteria, Actinobacteria, and Bacteroidetes. Similarly, in the month of December, Firmicutes, Proteobacteria, Actinobacteria, and Bacteroidetes were prevalent; however, the percent relative abundance of these phyla in the month of December is relatively less. Besides this decrease in percent abundance, it was interestingly seen that relatively more phyla were found contributing towards the bacterial diversity in the month of December. Similar to phylum level, at genus level, there was not much variation seen among various prevalent genera of the three studied hot springs in both months. The major genera prevalent in both months among all the three hot springs were followed by Bacillus, Desulfotomaculum, Lactobacillus, and Paenibacillus. A similar trend was also seen at gene level that relative abundance of various genera was higher in the month of October but more genera were found to be contributing towards bacterial diversity in the month of December. Few distinct genera were found to be more abundant in the month of December such as Rhodopirellula and Blastopirellula. The results may conclude that there is not much variation in the abundance and type of bacterial communities during the post-monsoon season in the month of October and December. However, this may be assumed that there is the accumulation or increase in the bacterial communities during the winter (relatively higher temperature among hot springs) and may favor few mesophilic and more thermophilic communities as well.

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.

Unique Geothermal Chemistry Shapes Microbial Communities on Mt. Erebus, Antarctica

Mt. Erebus, Antarctica, is the world's southernmost active volcano and is unique in its isolation from other major active volcanic systems and its distinctive geothermal systems. Using 16S rRNA gene amplicon sequencing and physicochemical analyses, we compared samples collected at two contrasting high-temperature (50°C-65°C) sites on Mt. Erebus: Tramway Ridge, a weather-protected high biomass site, and Western Crater, an extremely exposed low biomass site. Samples were collected along three thermal gradients, one from Western Crater and two within Tramway Ridge, which allowed an examination of the heterogeneity present at Tramway Ridge. We found distinct soil compositions between the two sites, and to a lesser extent within Tramway Ridge, correlated with disparate microbial communities. Notably, pH, not temperature, showed the strongest correlation with these differences. The abundance profiles of several microbial groups were different between the two sites; class Nitrososphaeria amplicon sequence variants (ASVs) dominated the community profiles at Tramway Ridge, whereas Acidobacteriotal ASVs were only found at Western Crater. A co-occurrence network, paired with physicochemical analyses, allowed for finer scale analysis of parameters correlated with differential abundance profiles, with various parameters (total carbon, total nitrogen, soil moisture, soil conductivity, sulfur, phosphorous, and iron) showing significant correlations. ASVs assigned to Chloroflexi classes Ktedonobacteria and Chloroflexia were detected at both sites. Based on the known metabolic capabilities of previously studied members of these groups, we predict that chemolithotrophy is a common strategy in this system. These analyses highlight the importance of conducting broader-scale metagenomics and cultivation efforts at Mt. Erebus to better understand this unique environment.

Comparative Gut Microbiome in Trachypithecus leucocephalus and Other Primates in Guangxi, China, Based on Metagenome Sequencing

The Trachypithecus leucocephalus (white-headed langur) is a highly endangered, karst-endemic primate species, inhabiting the karst limestone forest in Guangxi, Southwest China. How white-headed langurs adapted to karst limestone and special dietary remains unclear. It is the first time to study the correlation between the gut microbiome of primates and special dietary, and environment in Guangxi. In the study, 150 fecal samples are collected from nine primates in Guangxi, China. Metagenomic sequencing is used to analyze and compare the gut microbiome composition and diversity between white-headed langurs and other primates. Our results indicate that white-headed langurs has a higher diversity of microbiome than other primates, and the key microbiome are phylum Firmicutes, class Clostridia, family Lachnospiraceae, and genera Clostridiates and Ruminococcus, which are related to the digestion and degradation of cellulose. Ten genera are significantly more abundant in white-headed langurs and François’ langur than in other primates, most of which are high-temperature microbiome. Functional analysis reveals that energy synthesis-related pathways and sugar metabolism-related pathways are less abundant in white-headed langurs and François’ langur than in other primates. This phenomenon could be an adaptation mechanism of leaf-eating primates to low-energy diet. The gut microbiome of white-headed langurs is related to diet and karst limestone environment. This study could serve as a reference to design conservation breeding, manage conservation units, and determine conservation priorities.

Identification of groundwater microbial communities and their connection to the hydrochemical environment in southern Laizhou Bay, China

The microbial community plays an important role in the biogeochemical cycle in coastal groundwater ecosystems. However, the composition and controlling factors of the microbial community in coastal closed groundwater systems (CCGSs) with high salinity have rarely been studied. Here, we investigated and analyzed the hydrochemical characteristics and microbial community composition of seven brine samples with high total dissolved solid (TDS) values ranging from 74.5 to 132.3 g/L within and across three coastal saltworks (Yangkou, Hanting, and Changyi) in southern Laizhou Bay (SLB). The bacterial diversity was independent of salinity. Compared with those of low-salinity groundwater, the diversity of the microbial community in brine was lower, but the richness was slightly higher. There was a significant correlation between the microbial community diversity and groundwater sources, which indicated that the microbial communities were affected by groundwater sources. A comparison of the microbial community compositions of the three saltworks showed that the Hanting and Changyi saltworks had similar microbial communities due to their similar sampling depths. In addition, the main force shaping the differences in the microbial communities in both coastal open groundwater systems (COGSs) and CCGSs was identified as the hydraulic connection with the seawater controlled by hydrogeological conditions formed throughout geological history. This study can help to elucidate the biogeochemical processes in coastal aquifers.

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.

Co-occurrence of planktonic bacteria and archaea affects their biogeographic patterns in China's coastal wetlands

Planktonic bacteria and archaea play a key role in maintaining ecological functions in aquatic ecosystems; however, their biogeographic patterns and underlying mechanisms have not been well known in coastal wetlands including multiple types and at a large space scale. Therefore, planktonic bacteria and archaea and related environmental factors were investigated in twenty-one wetlands along China's coast to understand the above concerns. The results indicated that planktonic bacteria had different biogeographic pattern from planktonic archaea, and both patterns were not dependent on the wetland's types. Deterministic selection shapes the former's community structure, whereas stochastic processes regulate the latter's, being consistent with the fact that planktonic archaea have a larger niche breadth than planktonic bacteria. Planktonic bacteria and archaea co-occur, and their co-occurrence rather than salinity is more important in shaping their community structure although salinity is found to be a main environmental deterministic factor in the coastal wetland waters. This study highlights the role of planktonic bacteria-archaea co-occurrence on their biogeographic patterns, and thus provides a new insight into studying underlying mechanisms of microbial biogeography in coastal wetlands.

Distribution of Hydrogen-Producing Bacteria in Tibetan Hot Springs, China

Investigating the distribution of hydrogen-producing bacteria (HPB) is of great significance to understanding the source of biological hydrogen production in geothermal environments. Here, we explored the compositions of HPB populations in the sediments of hot springs from the Daggyai, Quzhuomu, Quseyongba, and Moluojiang geothermal zones on the Tibetan Plateau, with the use of Illumina MiSeq high-throughput sequencing of 16S rRNA genes and hydA genes. In the present study, the hydA genes were successfully amplified from the hot springs with a temperature of 46–87°C. The hydA gene phylogenetic analysis showed that the top three phyla of the HPB populations were Bacteroidetes (14.48%), Spirochaetes (14.12%), and Thermotogae (10.45%), while Proteobacteria were absent in the top 10 of the HPB populations, although Proteobacteria were dominant in the 16S rRNA gene sequences. Canonical correspondence analysis results indicate that the HPB community structure in the studied Tibetan hot springs was correlated with various environmental factors, such as temperature, pH, and elevation. The HPB community structure also showed a spatial distribution pattern; samples from the same area showed similar community structures. Furthermore, one HPB isolate affiliated with Firmicutes was obtained and demonstrated the capacity of hydrogen production. These results are important for us to understand the distribution and function of HPB in hot springs.

Temperature and microbial interactions drive the deterministic assembly processes in sediments of hot springs

Terrestrial geothermal ecosystems, as a representative of extreme environments, exhibit a variety of geochemical gradients, and their microbes are thought to be under high stress through environmental selection. However, it is still unclear how stochasticity and biotic interactions contribute to the microbial community assembly in hot springs. Here, we investigated the assembly processes and co-occurrence patterns of microbiota (i.e. bacteria and archaea) in both water and sediments sampled from fifteen hot springs in the Tengchong area, Southwestern of China, using 16S rRNA gene sequencing combined with multivariate ecological and statistical methods. These hot springs harbored more specialists than non-geothermal ecosystems, which are well-adapted to the extreme conditions, as shown by extremely high nearest-taxon index (NTI) and narrower niche width. Habitat differentiation led to the differences in microbial diversity, species-interactions, and community assembly between water and sediment communities. The sediment community showed stronger phylogenetic clustering and was primarily governed by heterogeneous selection, while undominated stochastic processes and dispersal limitation were the major assembly processes in the water community. Temperature and ferrous iron were the major factors mediating the balance of stochastic and deterministic assembly processes in sediment communities, as evidenced by how divergences in temperature and ferrous iron increased the proportion of determinism. Microbial interactions in sediments contributed to deterministic community assembly, as indicated by more complex associations and greater responsiveness to environmental change than water community. These findings uncover the ecological processes underlying microbial communities in hot springs, and provide potential insight into understanding the mechanism to maintain microbial diversity in extreme biospheres.

Distinct assembly processes shape bacterial communities along unsaturated, groundwater fluctuated, and saturated zones

The subsurface soil environment through the unsaturated (vadose) zone and saturated (below groundwater table) zone is one of the most active layers in the Earth's surface with biogeochemical interactions. Geochemical variables and geographic distance are key driving forces shaping the distribution of soil microbial communities, but our understandings are mainly limited to surface soil or shallow unsaturated zone (1-3 m beneath the ground). In this study, soil and sediment samples were collected from the unsaturated zone, through groundwater fluctuated zone, to saturated zone (up to 20 m) to unravel the assembly processes mediating vertical bacterial community succession across these three zones. Our results suggested both geochemical niches and bacterial diversity had different vertical patterns in each zone. With increased depth, pH increased and nutrient levels (C, N, P, K) and bacterial diversity declined in the unsaturated zone, and nutrients and bacterial diversity remained low levels after reaching the fluctuated and saturated zones. Nutrients were the key drivers shaping bacterial variation in the unsaturated zone, but limited nutrients and only 'depth' significantly explained the variations in the fluctuated zone and saturated zone, respectively. The co-occurrence network supported a more species co-existence pattern in the unsaturated zone than that in the other two zones. Due to the geochemical variations across three zones, the assembly of phylogenetically more clustered communities was observed through deterministic processes (e.g., 55% homogenizing selection) in the unsaturated zone, but the stochastic process (e.g., 50%-70% dispersal limitation) was more important in the fluctuated and saturated zones. These findings together suggested that the vertical distribution of soil bacterial community assembly was zone-specific and shaped by the degree of deterministic vs. stochastic processes. Our results provide a novel insight into the microbial community assembly across three different ecosystems in the Earth's critical zone and shed a light on subsurface biogeochemical processes.

History of petroleum disturbance triggering the depth-resolved assembly process of microbial communities in the vadose zone

Biogeochemical gradient forms in vadose zone, yet little is known about the assembly processes of microbial communities in this zone under petroleum disturbance. This study collected vadose zone soils at three sites with 0, 5, and 30 years of petroleum contamination to unravel the vertical microbial community successions and their assembly mechanisms. The results showed that petroleum hydrocarbons exhibited higher concentrations at the long-term contaminated site, showing negative impacts on some soil properties, retarding in the surface soils and decreasing along soil depth. Cultivable fraction of heterotrophic bacteria and microbial α-diversity decreased along depth in vadose zones with short-term/no contamination history, but exhibited an opposite trend with long-term contamination history. Petroleum contamination intensified the vertical heterogeneity of microbial communities based on the contamination time. Microbial co-occurrence network revealed the lowest species co-occurrence pattern at the long-term contaminated site. The distance-decay patterns and null model analysis together suggested distinct assembly mechanisms at three sites, where dispersal limitation (42-45%) was higher and variable and homogenizing selections were lower (37-38%) in vadose zones under petroleum disturbance than those in the uncontaminated vadose zone. Our findings help to better understand the subsurface biogeochemical cycles and bioremediation of petroleum-contaminated vadose zones.