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

Tectonic and geological setting influence hot spring microbiology

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

Distribution and abundance of tetraether lipid cyclization genes in terrestrial hot springs reflect pH

Many Archaea produce membrane-spanning lipids that enable life in extreme environments. These isoprenoid glycerol dibiphytanyl glycerol tetraethers (GDGTs) may contain up to eight cyclopentyl and one cyclohexyl ring, where higher degrees of cyclization are associated with more acidic, hotter or energy-limited conditions. Recently, the genes encoding GDGT ring synthases, grsAB, were identified in two Sulfolobaceae; however, the distribution and abundance of grs homologs across environments inhabited by these and related organisms remain a mystery. To address this, we examined the distribution of grs homologs in relation to environmental temperature and pH, from thermal springs across Earth, where sequences derive from metagenomes, metatranscriptomes, single-cell and cultivar genomes. The abundance of grs homologs shows a strong negative correlation to pH, but a weak positive correlation to temperature. Archaeal genomes and metagenome-assembled genomes (MAGs) that carry two or more grs copies are more abundant in low pH springs. We also find grs in 12 archaeal classes, with the most representatives in Thermoproteia, followed by MAGs of the uncultured Korarchaeia, Bathyarchaeia and Hadarchaeia, while several Nitrososphaeria encodes >3 copies. Our findings highlight the key role of grs-catalysed lipid cyclization in archaeal diversification across hot and acidic environments.

Comparative metagenomics at Solfatara and Pisciarelli hydrothermal systems in Italy reveal that ecological differences across substrates are not ubiquitous

Introduction

Continental hydrothermal systems (CHSs) are geochemically complex, and they support microbial communities that vary across substrates. However, our understanding of these variations across the complete range of substrates in CHS is limited because many previous studies have focused predominantly on aqueous settings.

Methods

Here we used metagenomes in the context of their environmental geochemistry to investigate the ecology of different substrates (i.e., water, mud and fumarolic deposits) from Solfatara and Pisciarelli.

Results and Discussion

Results indicate that both locations are lithologically similar with distinct fluid geochemistry. In particular, all substrates from Solfatara have similar chemistry whereas Pisciarelli substrates have varying chemistry; with water and mud from bubbling pools exhibiting high SO₄ ^2- and NH₄ ⁺ concentrations. Species alpha diversity was found to be different between locations but not across substrates, and pH was shown to be the most important driver of both diversity and microbial community composition. Based on cluster analysis, microbial community structure differed significantly between Pisciarelli substrates but not between Solfatara substrates. Pisciarelli mud pools, were dominated by (hyper)thermophilic archaea, and on average, bacteria dominated Pisciarelli fumarolic deposits and all investigated Solfatara environments. Carbon fixation and sulfur oxidation were the most important metabolic pathways fueled by volcanic outgassing at both locations. Together, results demonstrate that ecological differences across substrates are not a widespread phenomenon but specific to the system. Therefore, this study demonstrates the importance of analyzing different substrates of a CHS to understand the full range of microbial ecology to avoid biased ecological assessments.

Archaeal lipids trace ecology and evolution of marine ammonia-oxidizing archaea

Archaeal membrane lipids are widely used for paleotemperature reconstructions, yet these molecular fossils also bear rich information about ecology and evolution of marine ammonia-oxidizing archaea (AOA). Here we identified thermal and nonthermal behaviors of archaeal glycerol dialkyl glycerol tetraethers (GDGTs) by comparing the GDGT-based temperature index (TEX₈₆) to the ratio of GDGTs with two and three cyclopentane rings (GDGT-2/GDGT-3). Thermal-dependent biosynthesis should increase TEX₈₆ and decrease GDGT-2/GDGT-3 when the ambient temperature increases. This presumed temperature-dependent (PTD) trend is observed in GDGTs derived from cultures of thermophilic and mesophilic AOA. The distribution of GDGTs in suspended particulate matter (SPM) and sediments collected from above the pycnocline-shallow water samples-also follows the PTD trend. These similar GDGT distributions between AOA cultures and shallow water environmental samples reflect shallow ecotypes of marine AOA. While there are currently no cultures of deep AOA clades, GDGTs derived from deep water SPM and marine sediment samples exhibit nonthermal behavior deviating from the PTD trend. The presence of deep AOA increases the GDGT-2/GDGT-3 ratio and distorts the temperature-controlled correlation between GDGT-2/GDGT-3 and TEX₈₆. We then used Gaussian mixture models to statistically characterize these diagnostic patterns of modern AOA ecology from paleo-GDGT records to infer the evolution of marine AOA from the Mid-Mesozoic to the present. Long-term GDGT-2/GDGT-3 trends suggest a suppression of today's deep water marine AOA during the Mesozoic-early Cenozoic greenhouse climates. Our analysis provides invaluable insights into the evolutionary timeline and the expansion of AOA niches associated with major oceanographic and climate changes.

Comparative Genomics of Thaumarchaeota From Deep-Sea Sponges Reveal Their Niche Adaptation

Thaumarchaeota account for a large portion of microbial symbionts in deep-sea sponges and are even dominant in some cases. In this study, we investigated three new sponge-associated Thaumarchaeota from the deep West Pacific Ocean. Thaumarchaeota were found to be the most dominant phylum in this sponge by both prokaryotic 16S rRNA amplicons and metagenomic sequencing. Fifty-seven published Thaumarchaeota genomes from sponges and other habitats were included for genomic comparison. Similar to shallow sponge-associated Thaumarchaeota, those Thaumarchaeota in deep-sea sponges have extended genome sizes and lower coding density compared with their free-living lineages. Thaumarchaeota in deep-sea sponges were specifically enriched in genes related to stress adapting, symbiotic adhesion and stability, host–microbe interaction and protein transportation. The genes involved in defense mechanisms, such as the restriction-modification system, clustered regularly interspaced short palindromic repeats (CRISPR)/Cas system, and toxin-antitoxin system were commonly enriched in both shallow and deep sponge-associated Thaumarchaeota. Our study demonstrates the significant effects of both depth and symbiosis on forming genomic characteristics of Thaumarchaeota, and provides novel insights into their niche adaptation in deep-sea sponges.

Multiple environmental parameters impact lipid cyclization in Sulfolobus acidocaldarius

Adaptation of lipid membrane composition is an important component of archaeal homeostatic response. Historically, the number of cyclopentyl and cyclohexyl rings in the glycerol dibiphytanyl glycerol tetraether (GDGT) Archaeal lipids has been linked to variation in environmental temperature. However, recent work with GDGT-making archaea highlight the roles of other factors, such as pH or energy availability, in influencing the degree of GDGT cyclization. To better understand the role of multiple variables in a consistent experimental framework and organism, we cultivated the model Crenarchaeon Sulfolobus acidocaldarius DSM639 at different combinations of temperature, pH, oxygen flux, or agitation speed. We quantified responses in growth rate, biomass yield, and core lipid compositions, specifically the degree of core GDGT cyclization. The degree of GDGT cyclization correlated with growth rate under most conditions. The results suggest the degree of cyclization in archaeal lipids records a universal response to energy availability at the cellular level, both in thermoacidophiles, and in other recent findings in the mesoneutrophilic Thaumarchaea. Although we isolated the effects of key individual parameters, there remains a need for multi-factor experiments (e.g., pH + temperature + redox) in order to more robustly establish a framework to better understand homeostatic membrane responses.

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

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

The Intersection of Geology, Geochemistry, and Microbiology in Continental Hydrothermal Systems

Decompressional boiling of ascending hydrothermal waters and separation into a vapor (gas) and a liquid phase drive extensive variation in the geochemical composition of hot spring waters. Yet little is known of how the process of phase separation influences the distribution of microbial metabolisms in springs. Here, we determined the variation in protein coding genes in 51 metagenomes from chemosynthetic hot spring communities that span geochemical gradients in Yellowstone National Park. The 51 metagenomes could be divided into 5 distinct groups that correspond to low and high temperatures and acidic and circumneutral/alkaline springs. A fifth group primarily comprised metagenomes from springs with moderate acidity and that are influenced by elevated volcanic gas input. Protein homologs putatively involved in the oxidation of sulfur compounds, a process that leads to acidification of spring waters, in addition to those involved in the reduction of sulfur compounds were enriched in metagenomes from acidic springs sourced by vapor phase gases. Metagenomes from springs with evidence for elevated volcanic gas input were enriched in protein homologs putatively involved in oxidation of those gases, including hydrogen and methane. Finally, metagenomes from circumneutral/alkaline springs sourced by liquid phase waters were enriched in protein homologs putatively involved in heterotrophy and respiration of oxidized nitrogen compounds and oxygen. These results indicate that the geological process of phase separation shapes the ecology of thermophilic communities through its influence on the availability of nutrients in the form of gases, solutes, and minerals. Microbial acidification of hot spring waters further influences the kinetic and thermodynamic stabilities of nutrients and their bioavailability. These data therefore provide an important framework to understand how geological processes have shaped the evolutionary history of chemosynthetic thermophiles and how these organisms, in turn, have shaped their geochemical environments.

Global Microbiome Diversity Scaling in Hot Springs With DAR (Diversity-Area Relationship) Profiles

The spatial distribution of biodiversity (i.e., the biogeography) of the hot-spring microbiome is critical for understanding the microbial ecosystems in hot springs. We investigated the microbiome diversity scaling (changes) over space by analyzing the diversity-area relationship (DAR), which is an extension to classic SAR (species-area relationship) law in biogeography. We built DAR models for archaea and bacteria with 16S-rRNA sequencing datasets from 165 hot springs globally. From the DAR models, we sketch out the biogeographic maps of hot-spring microbiomes by constructing: (i) DAR profile-measuring the archaea or bacteria diversity scaling over space (areas); (ii) PDO (pair-wise diversity overlap or similarity) profile-estimating the PDO between two hot springs; (iii) MAD (maximal accrual diversity) profile-predicting the global MAD; (iv) LRD/LGD (ratio of local diversity to regional or global diversity) profile. We further investigated the differences between archaea and bacteria in their biogeographic maps. For example, the comparison of DAR-profile maps revealed that the archaea diversity is more heterogeneous (i.e., more diverse) or scaling faster than the bacterial diversity does in terms of species numbers (species richness), but is less heterogeneous (i.e., less diverse) or scaling slower than bacteria when the diversity (Hill numbers) were weighted in favor of more abundant dominant species. When the diversity is weighted equally in terms of species abundances, archaea, and bacteria are equally heterogeneous over space or scaling at the same rate. Finally, unified DAR models (maps) were built with the combined datasets of archaea and bacteria.

Mixing of meteoric and geothermal fluids supports hyperdiverse chemosynthetic hydrothermal communities

Little is known of how mixing of meteoric and geothermal fluids supports biodiversity in non-photosynthetic ecosystems. Here, we use metagenomic sequencing to investigate a chemosynthetic microbial community in a hot spring (SJ3) of Yellowstone National Park that exhibits geochemistry consistent with mixing of a reduced volcanic gas-influenced end member with an oxidized near-surface meteoric end member. SJ3 hosts an exceptionally diverse community with representatives from ~50% of known higher-order archaeal and bacterial lineages, including several divergent deep-branching lineages. A comparison of functional potential with other available chemosynthetic community metagenomes reveals similarly high diversity and functional potentials (i.e., incorporation of electron donors supplied by volcanic gases) in springs sourced by mixed fluids. Further, numerous closely related SJ3 populations harbor differentiated metabolisms that may function to minimize niche overlap, further increasing endemic diversity. We suggest that dynamic mixing of waters generated by subsurface and near-surface geological processes may play a key role in the generation and maintenance of chemosynthetic biodiversity in hydrothermal and other similar environments.

Chemosynthetic microbial communities in hydrothermal environments receiving meteoric and geothermal fluids are understudied. Here, Colman et al. use metagenomics to study one such community from a hot spring at Yellowstone National Park, revealing exceptional biodiversity and unique functional potential.

Localized high abundance of Marine Group II archaea in the subtropical Pearl River Estuary: implications for their niche adaptation

Marine Group II archaea are widely distributed in global oceans and dominate the total archaeal community within the upper euphotic zone of temperate waters. However, factors controlling the distribution of MGII are poorly delineated and the physiology and ecological functions of these still-uncultured organisms remain elusive. In this study, we investigated the planktonic MGII associated with particles and in free-living forms in the Pearl River Estuary (PRE) over a 10-month period. We detected high abundance of particle-associated MGII in PRE (up to ∼10⁸ 16S rRNA gene copies/l), which was around 10-fold higher than the free-living MGII in the same region, and an order of magnitude higher than previously reported in other marine environments. 10‰ salinity appeared to be a threshold value for these MGII because MGII abundance decreased sharply below it. Above 10‰ salinity, the abundance of MGII on the particles was positively correlated with phototrophs and MGII in the surface water was negatively correlated with irradiance. However, the abundances of those free-living MGII showed positive correlations with salinity and temperature, suggesting the different physiological characteristics between particle-attached and free-living MGIIs. A nearly completely assembled metagenome, MGIIa_P, was recovered using metagenome binning methods. Compared with the other two MGII genomes from surface ocean, MGIIa_P contained higher proportions of glycoside hydrolases, indicating the ability of MGIIa_P to hydrolyse glycosidic bonds in complex sugars in PRE. MGIIa_P is the first assembled MGII metagenome containing a catalase gene, which might be involved in scavenging reactive oxygen species generated by the abundant phototrophs in the eutrophic PRE. Our study presented the widespread and high abundance of MGII in the water columns of PRE, and characterized the determinant abiotic factors affecting their distribution. Their association with heterotrophs, preference for particles and resourceful metabolic traits indicate MGII might play a significant role in metabolising organic matters in the PRE and other temperate estuarine systems.

Environmental factors shaping the archaeal community structure and ether lipid distribution in a subtropic river and estuary, China

Archaea are widespread and abundant in aquatic and terrestrial habitats and play fundamental roles in global biogeochemical cycles. Archaeal lipids, such as isoprenoid glycerol diakyl glycerol tetraethers (iGDGTs), are important biomarkers tracing changes in archaeal community structure and biogeochemical processes in nature. However, the linkage between the archaeal populations and the GDGT distribution in the natural environment is poorly examined, which hindered the application and interpretation of GDGT-based climate or environmental proxies. We addressed this question by investigating changes in archaeal lipid composition and community structure in the context of environmental variables along the subtropical Jiulong River Watershed (JRW) and Jiulong River Estuary (JRE) in southern China. The results showed that both the archaeal cells and the polar GDGTs (P-GDGTs) in the JRW and JRE were mostly autochthonous rather than exogenous input from surrounding soils. We further found that only five (Methanobacteriales, Ca. Bathyarchaeota, Marine Benthic Groups A (MBGA), Marine Benthic Groups B (MBGB), and Marine Benthic Groups D (MBGD)) out of sixteen lineages showed significant impacts on the composition of P-GDGTs, suggesting the significant contribution of those archaea to the changes of P-GDGT compositions. Salinity and total phosphorus (TP) showed significant impact on the distribution of both genetic and P-GDGTs compositions of archaea; whereas, sand and silt contents only had significant impact on the P-GDGTs. MBGD archaea, which occur widely in marine sediments, showed positive correlations with P-TEX₈₆ in the JRW and JRE, suggesting that uncultivated MBGD might also contribute to the variations in TEX₈₆ signals in marine sediments. This study provided insight into the sources of P-GDGTs and the factors controlling their distributions in river-dominated continental margins, which has relevance to applications of GDGT-based proxies in paleoclimate studies.

Geobiological feedbacks and the evolution of thermoacidophiles

Oxygen-dependent microbial oxidation of sulfur compounds leads to the acidification of natural waters. How acidophiles and their acidic habitats evolved, however, is largely unknown. Using 16S rRNA gene abundance and composition data from 72 hot springs in Yellowstone National Park, Wyoming, we show that hyperacidic (pH<3.0) hydrothermal ecosystems are dominated by a limited number of archaeal lineages with an inferred ability to respire O₂. Phylogenomic analyses of 584 existing archaeal genomes revealed that hyperacidophiles evolved independently multiple times within the Archaea, each coincident with the emergence of the ability to respire O₂, and that these events likely occurred in the recent evolutionary past. Comparative genomic analyses indicated that archaeal thermoacidophiles from independent lineages are enriched in similar protein-coding genes, consistent with convergent evolution aided by horizontal gene transfer. Because the generation of acidic environments and their successful habitation characteristically require O₂, these results suggest that thermoacidophilic Archaea and the acidity of their habitats co-evolved after the evolution of oxygenic photosynthesis. Moreover, it is likely that dissolved O₂ concentrations in thermal waters likely did not reach levels capable of sustaining aerobic thermoacidophiles and their acidifying activity until ~0.8 Ga, when present day atmospheric levels were reached, a time period that is supported by our estimation of divergence times for archaeal thermoacidophilic clades.

Microbial community analysis of pH 4 thermal springs in Yellowstone National Park

The pH of the majority of thermal springs in Yellowstone National Park (YNP) is from 1 to 3 and 6 to 10; relatively few springs (~5%) have a pH range of 4-5. We used 16S rRNA gene pyrosequencing to investigate microbial communities sampled from four pH 4 thermal springs collected from four regions of YNP that differed in their fluid temperature and geochemistry. Our results revealed that the composition of bacterial communities varied among the sites, despite sharing similar pH values. The taxonomic composition and metabolic functional potential of the site with the lowest temperature (55 °C), a thermal spring from the Seven Mile Hole (SMH) area, were further investigated using shotgun metagenome sequencing. The taxonomic classification, based on 372 Mbp of unassembled metagenomic reads, indicated that this community included a high proportion of Chloroflexi, Bacteroidetes, Proteobacteria, and Firmicutes. Functional comparison with other YNP thermal spring metagenomes indicated that the SMH metagenome was enriched in genes related to energy production and conversion, transcription, and carbohydrate transport. Analysis of genes involved in nitrogen metabolism revealed assimilatory and dissimilatory nitrate reduction pathways, whereas the majority of genes involved in sulfur metabolism were related to the reduction of sulfate to adenylylsulfate, sulfite, and H₂S. Given that pH 4 thermal springs are relatively less common in YNP and thermal areas worldwide, they may harbor novel microbiota and the communities that inhabit them deserve further investigation.

Ecological differentiation in planktonic and sediment-associated chemotrophic microbial populations in Yellowstone hot springs

Chemosynthetic sediment and planktonic community composition and sizes, aqueous geochemistry and sediment mineralogy were determined in 15 non-photosynthetic hot springs in Yellowstone National Park (YNP). These data were used to evaluate the hypothesis that differences in the availability of dissolved or mineral substrates in the bulk fluids or sediments within springs coincides with ecologically differentiated microbial communities and their populations. Planktonic and sediment-associated communities exhibited differing ecological characteristics including community sizes, evenness and richness. pH and temperature influenced microbial community composition among springs, but within-spring partitioning of taxa into sediment or planktonic communities was widespread, statistically supported (P < 0.05) and could be best explained by the inferred metabolic strategies of the partitioned taxa. Microaerophilic genera of the Aquificales predominated in many of the planktonic communities. In contrast, taxa capable of mineral-based metabolism such as S(o) oxidation/reduction or Fe-oxide reduction predominated in sediment communities. These results indicate that ecological differentiation within thermal spring habitats is common across a range of spring geochemistry and is influenced by the availability of dissolved nutrients and minerals that can be used in metabolism.