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₄²⁻ 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.

Introduction

Continental hydrothermal systems (CHSs) are usually formed within diverse igneous lithologies (i.e., mafic, andesitic, and felsic) and the hydrothermal fluids derived from deep subsurface sources may be mixed with marine, brine, or meteoritic water resulting in geochemically diverse environments. These diverse environments exhibit substantial differences in the availability and abundance of electron acceptors and donors (Fournier, 1989; Shock et al., 2010; Lowenstern et al., 2015; Lindsay et al., 2018; Amenabar and Boyd, 2019), which microorganisms exploit as sources of energy (Amenabar et al., 2017; Lindsay et al., 2019; Puopolo et al., 2020; Aulitto et al., 2021).

CHSs have been the subject of many studies as they are a surface manifestation of hydrothermal activities, which are easily accessible for diverse microbiological studies (Inskeep et al., 2013; Menzel et al., 2015; Colman et al., 2016; Mardanov et al., 2018; Oliverio et al., 2018; Power et al., 2018; Boyer et al., 2020; Podar et al., 2020; Crognale et al., 2022 as a few examples). These studies have revealed that hot springs support microbial communities that are exceptionally diverse and vary in their response to physical and geochemical parameters (Inskeep et al., 2013; Sharp et al., 2014; Lindsay et al., 2018; Power et al., 2018; Colman et al., 2019a,b; Podar et al., 2020). Some studies have found pH to be the primary physical parameter that influences the microbial community composition of hot springs in Yellowstone National Park (Boyd et al., 2010; Inskeep et al., 2013; Colman et al., 2016, 2019a,b), Tengchong, China (Hou et al., 2013; Xie et al., 2015), Iceland (Moreras-Marti et al., 2021), and volcanic provinces in New Zealand (Power et al., 2018). In contrast, other studies have identified temperature as the most important driver of microbial community composition in hot springs (Miller et al., 2009; Cole et al., 2013; Sharp et al., 2014; Podar et al., 2020). Further, at Yellowstone National Park, subsurface processes including phase separation and mixing with meteoritic fluids were shown to shape the ecology of hot spring communities through their influence on the availability of nutrients that support microbial metabolism (Lindsay et al., 2018; Colman et al., 2019a). In addition to temperature, pH, phase separation and mixing of fluids, Fullerton et al. (2021) found that microbial diversity analyzed from fluids and associated sediments reflects the subsurface geological structures that fluids traverse, which in turn influences carbon cycling within the subduction zone in Costa Rica.

While we have learned a great deal about hot springs, their geochemical complexity and the observed ecological differences between water and hot spring deposits (Cole et al., 2013; Colman et al., 2016) highlight the need for studies analyzing different substrates (e.g., fumarolic deposits) to understand the full range of microbial ecology in CHSs. Some studies have analyzed microbial communities from different available substrates (Glamoclija et al., 2004; Ellis et al., 2008; Benson et al., 2011; Sharp et al., 2014; Wall et al., 2015; Medrano-Santillana et al., 2017; Crognale et al., 2018; Marlow et al., 2020; Arif et al., 2021); however, to the best of our knowledge, no CHSs study has used as many different substrates from the same hydrothermal system and correlated their environmental settings and metagenomes to characterize microbial ecology. The lack of comparison among different substrates (fumarolic deposits, water and mud from mud pools) may derive from the fact that not all substrates are present at all CHSs. In this study, we use metagenomes to investigate the microbial ecology in different substrates (i.e., water, mud and fumarolic deposits) from two CHSs in the context of their environmental geochemistry. The objectives of this study were to: (1) assess the microbial diversity of different CHSs and a variety of available substrates (2) identify the metabolic potential of these microbial communities in connection to environmental geochemistry.

The study areas are Solfatara and Pisciarelli CHSs (Figure 1) located within Campi Flegrei Caldera (CFC); an 8 km nested caldera in Naples (Italy) formed by the Campania Ignimbrite (39 ka) and Neapolitan Yellow Tuff (15 ka) eruptions (De Natale et al., 2016; Rolandi et al., 2020a,b). Present-day activities at CFC are characterized by large-scale hydrothermal circulation, gaseous emissions, and intense ground deformation (Troise et al., 2019; Moretti et al., 2020; Chiodini et al., 2021). According to Caliro et al. (2007), fumaroles and mud pools at Solfatara and Pisciarelli are fed by fluids of mixed magmatic-meteoritic origin. The Solfatara crater is 0.6 km in diameter and maintains hydrothermal activity at fumaroles where gas emissions reach temperatures from 145°C up to 165°C, while mud pools have an average temperature of approximately 45°C (Glamoclija et al., 2004; Chiodini et al., 2011). Furthermore, the Solfatara environment is extremely acidic (mud pool: 1.9–2.1 pH and fumarolic deposits: 1.3–2.2 pH; Crognale et al., 2018). The second location, Pisciarelli, sits on the outside northeastern wall of the Solfatara crater. Pisciarelli holds an unstable fumarole field characterized by consistent deposit degassing, fluid emission from ephemeral vents, and boiling mud pools. This degassing activity is episodically accompanied by seismic swarms and macroscopic morphology changes such as the appearance of vigorously degassing vents, collapsing landslides, and bubbling mud (Fedele et al., 2021). Fumarolic deposits and mud pools at Pisciarelli are also very acidic (pH: 0.5–3) and temperatures of fumaroles ranges from 95°C to 110°C, while mud pools have temperatures between 84°C and 95°C (Ciniglia, 2005; Troiano et al., 2014; Poichi et al., 2019). In Solfatara-Pisciarelli CHSs, the major gases in fumarole vents are H₂O and CO₂ while the minor gases include H₂S, N₂, H₂, CH₄, He, Ar, and CO (Chiodini et al., 2010; Aiuppa et al., 2013). In both systems, fumaroles have a similar content of H₂O, CO₂, Ar, He, and N₂, however, they differ in their concentration of H₂S, H₂, and CO (Chiodini et al., 2001). Solfatara and Pisciarelli also differ in their fluid geochemistry; fluid in mud pools of Pisciarelli have been reported to have higher NH₄⁺ concentrations (508–1,026 mg L⁻¹) compared to Solfatara (<1 mg L⁻¹; Martini et al., 1991; Valentino and Stanzione, 2003; Glamoclija et al., 2004; Poichi et al., 2019) resulting in the precipitation of the minerals mascagnite [(NH₄)₂SO₄] and tschermigite [(NH₄)Al(SO₄)₂·12(H₂O)] (Poichi et al., 2019). The geochemical diversity of substrates in Solfatara-Pisciarelli CHSs provides an excellent location for evaluating the composition, structure, and functional potential of thermophilic microorganisms in different hydrothermal substrates and locations.

FIGURE 1

Figure 1. Sampling points at Solfatara-Pisciarelli hydrothermal systems. (A) Location map of Solfatara-Pisciarelli hydrothermal systems on contoured map of Italy (red dot). (B) Oblique view of Solfatara and Pisciarelli systems; red dashed lines represent faults; yellow dots are sampling points (E and D stand for Solfatara mud pool and fumaroles, respectively). (C) Sampling point of Solfatara fumarolic deposits (SF) next to the main fumarole Bocca Grande. (D) Close-up view of deposits to show heterogenicity of the material in regard to grain size and mineralogy (note different colors of deposits). (E) Solfatara mud pool (SMP) with bubbling water and warm muddy substrate near the pool. (F) Overview of the Pisciarelli location with marked sampling points. (G) Pisciarelli large pool (PLP) with bubbling hot mud (temp. 84.1°C). (H) Pisciarelli fumarolic deposits (PF) with sulfur crystals precipitate. (I) Pisciarelli small pool (PSP) sampling point. (J) Pisciarelli epilithic microbial layer (PLP-E) observed and sampled on the wall near the outflow channel of the large pool.

Materials and methods

To investigate the near-surface microbial community composition, structure and function at two CHSs, water, mud, and fumarolic deposits samples were collected from Solfatara and Pisciarelli (Figures 1A,B), in October 2012 (Supplementary Method). All samples were collected in triplicates, aseptically, using sterile Falcon tubes, Nalgene bottles, scoops, and gloves. After the collection, the samples were stored at −20°C until further processing; samples for long term storage were stored at −80°C. Temperature, pH, and redox potential (Eh) were measured in situ using a portable probe (Table 1). Gas readings were taken from the continuous Istituto Nazionale di Geofisica e Volcanologia (INGV) gas monitoring station at Solfatara and Pisciarelli (Table 1).

TABLE 1

Table 1. Field measurements for water (W), mud (M), mud outlet (MO), epilith from dry mud wall (E), and fumarolic deposits (D) samples collected from Solfatara-Pisciarelli hydrothermal systems.

Geochemical characterization

About 1 g of solid sample (fumarolic deposits and mud) was added to 20 ml of MilliQ water and extracted as described in Lezcano et al. (2019). The concentration of water-soluble anions (F⁻, Cl⁻, NO₂⁻, Br⁻, NO₃⁻, PO₄³⁻, and SO₄²⁻) in extracts and water samples was measured in triplicate using an 881 compact IC pro ion chromatography system (Metrohm, Switzerland) with a Metrosep A Supp 5-250/4.0 column. The concentration of NH₄⁺ in extracts and water samples was measured using the alkaline hypochlorite/phenol nitroprusside method, after adding sodium citrate to prevent the precipitation of calcium and magnesium salts (Solorzano, 1969). Before measuring concentrations of NH₄⁺ in water samples, samples were diluted 10-fold. Ammonium sulfate (NH₄)₂SO₄ solutions (0, 660, 1,320, 3,300, 6,600, and 13,200 ppm) were prepared and used as standard. The absorbance of each sample was measured in triplicate using an Evolution 60S UV-Vis Spectrophotometer at 640 nm wavelength.

X-ray fluorescence (XRF) was used to examine the chemical composition of mud and fumarolic deposits samples. Samples were analyzed in triplicate using a Horiba XGT-1000WR X-ray Fluorescence with an Rh tube X-ray source and elemental wt.% was determined using the XGT-1000WR software’s quantification (Supplementary Method).

Stable isotopic analysis of carbon, nitrogen, and sulfur

Carbon, nitrogen, and sulfur stable isotope and elemental concentration analyses were performed at the Earth and Planets Laboratory, Carnegie Institution for Science. For organic carbon (TOC), and organic δ¹³C measurements, samples were weighed into silver boats and fumed with 12 N HCl for 12–14 h. δ¹³C and δ¹⁵N isotopes were measured using a Thermo Scientific Delta VPlus isotope ratio gas-source mass spectrometer connected to a Carlo Erba (NA 2500) elemental analyzer (EA/IRMS) via a Conflo III interface. δ³⁴S was analyzed by the same gas-source mass spectrometer but connected to an Elementar Americas vario Micro CUBE elemental analyzer (EA-IRMS) via a Conflo III interface. Stable isotope values are reported in standard delta notation as ‰ variations relative to: Pee Dee Belemnite (PDB) for δ¹³C, atmospheric N₂ gas (AIR) for δ¹⁵N, Vienna Canyon Diablo Troctolite (V-CDT) for δ³⁴S, and with an analytical error of ±0.1‰. In-house δ¹³C and δ¹⁵N standards were also used and calibrated against international and certified standards as well. Additional δ³⁴S standards include the International Atomic Energy Agency reference materials IAEA-S-1 (δ³⁴S = −0.3‰), and IAEA S3 (δ³⁴S = −32.3‰), as well as NBS-123 (δ³⁴S = −17.09‰), NBS-127 (δ³⁴S = −21.17‰), and USGS-42 (δ³⁴S = −7.84‰). A subset of the samples was analyzed in 2021 at the EDGE Stable Isotope Laboratory at the University of California Riverside to confirm the very negative nitrogen isotope values in some of these samples. USGS25 (δ¹⁵N = −30.41‰) and USGS40 (δ¹⁵N = −4.52‰) were used as two of our calibrating standards and confirmed the values measured earlier at the Earth and Planets Laboratory.

DNA extraction, metagenomic sequencing and sequence processing

About 0.25 g of mud and fumarolic deposits was used for DNA extraction. Before DNA extraction, fumarolic deposits were powdered using a sterilized agate mortar and pestle. Water samples were filtered in the laboratory using a 0.2 μm VWR black polycarbonate filter. Filters were cut into small pieces with a sterile scalpel and used for DNA extraction. DNA was extracted using DNeasy PowerSoil kit (Qiagen Inc., Valencia, CA, United States) with modifications to the manufacturer’s instructions (Supplementary Method) and were stored at −80°C until further processing. Sequencing libraries were prepared using the Accel-NGS 2S Plus DNA Library Kit (Swift Biosciences Inc., Ann Arbor, MI, United States) according to the manufacturer’s protocols. The prepared library was sequenced on the Illumina Novaseq platform at the Oklahoma Medical Research Foundation. Metagenome sequence quality was assessed using FastQC v0.11.5¹ and reads with a mean quality score less than 25 were discarded. Sequences were quality trimmed and Illumina sequencing adapters removed using Trimmomatic v0.36 (Bolger et al., 2014) with the following parameters: SLIDINGWINDOW:4:15, LEADING:3, TRAILING:3, and MINLEN:36. After quality control, taxonomic annotation of metagenomes was performed using Kaiju v1.8.2 (Menzel et al., 2016) against the NCBI nr + euk database (accessed on February 16, 2022) with the following parameters (run mode: greedy, minimum match length: 11, minimum match score: 75, allowed mismatches: 5). To compare the effects of taxonomic profiling using metagenome vs. 16S rRNA gene sequences, the trimmed metagenome sequences were uploaded to the MG-RAST online server (Meyer et al., 2008) and then passed through the MG-RAST QC pipeline. Following the quality control on MG-RAST, 16S rRNA gene sequences present in the metagenome (metagenome-derived 16S rRNA gene sequences) were identified and taxonomy assigned. Taxonomic classification of the metagenome-derived 16S rRNA gene sequences was done using the Greengenes rRNA database (DeSantis et al., 2006) hosted in MG-RAST with a minimum cut-off identity of 60% and e-value of 5. Since contamination leads to an overestimation of diversity, especially in low biomass samples like that of Solfatara and Pisciarelli (Karstens et al., 2019), after assigning taxonomy, we removed molecular biology kit and laboratory contaminants (which hereinafter are referred to as “kitome”), low abundance taxa (<0.01%) and sequences that were unclassified at the domain level from the taxonomy table before performing downstream analyses (Supplementary Method).

Co-assembly, metagenomic binning, and MAG quality assessment

Bins were generated from 10 Solfatara and 17 Pisciarelli metagenomes. To maximize genome recovery, we co-assembled several samples based on their taxonomic composition. All 10 Solfatara metagenomes were co-assembled and all 17 Pisciarelli metagenomes were co-assembled using MEGAHIT v1.2.9 (Li et al., 2015). Contigs longer than 1,000 bp were then binned using the metaWRAP binning module (Uritskiy et al., 2018) which incorporates three binning methods: CONCOCT v1.1.0 (Alneberg et al., 2014), MaxBin2 v2.2.6 (Wu et al., 2016), and metaBAT2 v2.12.1 (Kang et al., 2019). The metaWRAP refinement module (Uritskiy et al., 2018) was used to merge results from the three binning methods using the -c 50 and -x 10 options to obtain bins with over 50% completeness and less than 10% contamination according to the CheckM tool v1.0.12 (Parks et al., 2015). Bins with >50% completeness and <10% contamination were then reassembled with SPAdes v3.13.0 (Bankevich et al., 2012) to improve the assembly quality. The contamination and completeness of resulting Metagenome Assembled Genomes (MAGs) were reassessed with CheckM (Parks et al., 2015).

Taxonomic and functional annotation of MAGs

Taxonomy was assigned to MAGs using GTDB-TK v1.7.0 (Chaumeil et al., 2020). To assess whether MAGs belong to the same species (species ANI ≥ 95%), average nucleotide identity (ANI) was calculated for each possible pair of MAGs using FastANI v1.33 (Jain et al., 2018). MAGs with taxonomic assignment similar to previously identified kitome (e.g., Corynebacterium and Staphylococcus; Sheik et al., 2018; Weyrich et al., 2019) were not included in downstream analyses. Functional annotation of MAGs, including gene prediction was done using Prokka v1.14.5 (Seemann, 2014). Predicted genes were compared against the Kyoto Encyclopedia of Genes and Genomes (KEGG) database using BlastKOALA server (Kanehisa et al., 2016) to obtain KEGG Orthology (KO) annotations.

Abundance of MAGs in metagenomes

To assess the abundance of recovered MAGs, quality trimmed sequences from each metagenome were mapped against each MAG using BBMap v38.96. Sequence counts were normalized as the number of sequences recruited per kilobase of MAG and gigabase of metagenome (RPKG). The normalized sequence counts allowed for direct comparison of genome abundance between metagenome of different depths.

Statistical analyses

All statistical analyses were performed in R v4.0.2 (R Core Team, 2020). Principal Component Analysis (PCA) was performed with environmental parameters using the “prcomp” function in the stat package with scaling enabled. PCA results were graphed using the “ggbiplot” function from the ggbiplot v0.55 package. Before using the geochemical data for statistical analyses, variables below detection limit in all samples were excluded and the dataset transformed using z-score. After removing kitome and all potential contaminants, we calculated alpha diversity using the species abundance table obtained from metagenome-derived 16S rRNA gene (Greengenes) and metagenome (NCBI nr + euk) taxonomic assignment. Alpha diversity (Shannon, observed richness and Pielou’s evenness) was calculated using the “estimate_richness” function from the phyloseq v1.32.0 package in R. Nonpareil v3.0 (Rodriguez-R et al., 2018) was then used to estimate coverage and also calculate diversity with kmer kernel and default parameters. Kruskal-Wallis test was used to investigate whether alpha diversity varied significantly between locations and substrates. The “lm()” function in R was used to perform linear regression to evaluate potential effects of environmental variables on alpha diversity. Bray–Curtis dissimilarity based on genus and species abundance was used to determine differences in microbial community composition. The Bray-Curtis distance calculated was visualized using a Non-metric Multidimensional Scaling (NMDS) plot. Analysis of Similarity (ANOSIM) was performed on the Bray–Curtis dissimilarity matrix using the vegan v2.5-6 package to evaluate the significance of microbial compositional differences between the two locations and different types of substrates. Similarity Percentage (SIMPER) analysis was performed using the vegan v2.5-6 package to identify genera that contributed to the dissimilarity between two locations and different types of substrates. Mantel test using Spearman’s correlation coefficient with 999 permutations was performed to evaluate the significance of correlation between community composition (Bray–Curtis dissimilarity) and environmental parameters (Euclidean distance).

Results and discussion

Geochemical context of Solfatara-Pisciarelli hydrothermal systems

At Solfatara-Pisciarelli CHSs (Figure 1), the crater structure that forms the rock substrate (alkaline potassic tephra and lava ranging from trachybasalt to phonotrachyte) for both locations was produced and shaped by the same volcanic events (Piochi et al., 2015). Solid substrates from both locations have similar bulk elemental composition [Supplementary Figure 1; Supplementary Table 1 and Piochi et al. (2014), Piochi et al. (2015)]. In addition, the same magma chamber provides a heat source for the fumaroles and mud pools (Chiodini et al., 2001; Valentino and Stanzione, 2003) of both systems. However, tectonic features diverge the fluids that feed Pisciarelli mud pools to pass through old marine sediments with strong organic imprints where the fluids become enriched in NH₄⁺ (Poichi et al., 2019).

PCA based on temperature, pH, Eh, and water-soluble nutrients showed further distinction between Solfatara and Pisciarelli (Figure 2; Supplementary Table 2). Pisciarelli an extremely acidic (pH: 1.5–2.5) and reducing (Eh: −628 to −485) environment has higher temperatures (74°C to 95°C) and gas concentrations including H₂S and CH₄ (Table 1) compared to Solfatara (pH: ~1; temperature: 42°C to 88.7°C; Eh: 198 and 330 mV), making Pisciarelli a more active and extreme hydrothermal environment than Solfatara. These higher temperatures, which inhibit the solubility of oxygen (Boyer et al., 2020), together with higher concentrations of reducing gases may contribute to the more reducing conditions we observed at Pisciarelli. The H₂S gas released in CHSs may be oxidized abiotically by oxygenated meteoritic waters (Nordstrom et al., 2005, 2009) or biotically by chemosynthetic microorganisms leading to the higher concentration of SO₄²⁻ we observed in mud pools (Pisciarelli: 3,326–6,208 ppm; Solfatara: 206–2,888 ppm) compared to fumarolic deposits (Pisciarelli: 118–439 ppm; Solfatara: 567–610 ppm). The concentrations of SO₄²⁻ in Pisciarelli mud pools were two times that of Solfatara mud pool possibly due to the higher concentration of H₂S gas measured at Pisciarelli. The concentration of NH₄⁺ was also higher in mud pools, particularly in Pisciarelli mud pools (Pisciarelli: 1,130–1,998 ppm; Solfatara: 1.1–36 ppm) than in fumarolic deposits (Pisciarelli: 0.8 to 38 ppm; Solfatara: 0.3 ppm) which corresponds with other studies that report higher concentrations of NH₄⁺ in Pisciarelli mud pools (Martini et al., 1991; Valentino and Stanzione, 2003; Glamoclija et al., 2004; Poichi et al., 2019). Temperature and geochemical measurements from Pisciarelli large mud pool revealed an environmental gradient; values were lower at the discharge channel (temperature: 79°C; SO₄²⁻: 134.4–193.6 ppm; NH₄⁺: 46.5–67.8 ppm) compared to the main bubbling pool (temperature: 84°C; SO₄²⁻: 3326–6,208 ppm; NH₄⁺: 1,130–1,998 ppm). The low concentrations of SO₄²⁻ and NH₄⁺ along the discharge channel may result from the removal of H₂S and NH₄⁺ possibly due to oxidation and volatilization as fluid flows into the channel.

FIGURE 2

Figure 2. Biplot based on Principal Component Analysis (PCA) using temperature, pH, Eh, water-soluble anions and cations. Variables that had values below detection limit in most samples (F⁻, Cl⁻, and PO₄³⁻) were excluded from PCA.

Discrepancies between metagenome-derived 16S rRNA gene and metagenomic profiles

In this study, a total of 7 Solfatara and 10 Pisciarelli samples were analyzed. Solfatara samples had between 54,932,932 to 210,855,916 (average = 126,261,781) quality filtered metagenome sequences whereas Pisciarelli had between 119,199,732 to 230,316,988 (average = 179,177,105) quality filtered metagenome sequences. We compared taxonomic profiles generated from metagenome against metagenome-derived 16S rRNA gene sequences and observed that metagenome sequences identified more taxa than metagenome-derived 16S rRNA gene sequences, which corresponds to other studies that detected an increased number of taxa with whole genome shotgun sequencing compared to the 16S amplicon method (Ranjan et al., 2016; Brumfield et al., 2020). For example, Ranjan et al. (2016) reported that with the same number of sequences, whole genome shotgun sequencing identified twice as many species as the 16S method. We observed that metagenome sequences identified 56 phyla while metagenome-derived 16S rRNA gene sequences identified 26 phyla. Only 3 out of the 26 phyla detected by metagenome-derived 16S rRNA gene profiling were not detected by metagenomic profiling. Metagenomic profiling, on the other hand, identified 35 phyla not identified by metagenome-derived 16S rRNA gene profiling with 12 of them being viral and eukaryotic phyla. Overall, the dominant bacterial and archaeal phyla detected were similar irrespective of profiling method. Predominant bacterial and archaeal phyla detected by both profiling methods include Crenarchaeota, Proteobacteria, Firmicutes, and Actinobacteria (Figures 3A,B). At the genus level, metagenomic profiling also identified more genera than metagenome-derived 16S rRNA gene profiling (Figures 4A,B). Of the 245 genera identified across all samples by metagenomic profiling methods only 48 genera (20%) were identified with both profiling methods. Metagenomic profiling identified 197 genera not identified by metagenome-derived 16S rRNA gene profiling including Acidilobus, Acidibacillus, Thermogymnomonas, Ampullavirus, and Bicaudavirus while genera including Acetobacterium and Caldococcus were unique to metagenome-derived 16S rRNA gene profiling. We observed that both methods identified dominant genera including Acidianus, Pyrobaculum, and Sulfobacillus.

FIGURE 3

Figure 3. Relative abundance of phyla based on (A) metagenome-derived 16S rRNA gene and (B) metagenomic profiling. Relative abundance was calculated after removing sequences that were unclassified at the domain level, kitome, singleton, species present in only one sample and low abundance species. Phyla with average abundance of <1% were grouped into Others <1%. W represents water, M for mud, MO for mud outlet, E for epilithic microbial layer from dry mud, and D for fumarolic deposits samples. Samples labeled a and b are replicates sampled from the same spot. Sample SMP-M_a1 and SMP-M_a1.2 are sequencing replicates.

FIGURE 4

Figure 4. Relative abundance of genera based on (A) metagenome-derived 16S rRNA gene and (B) metagenomic profiling. Relative abundance was calculated after removing sequences that were unclassified at the domain level, kitome, singleton, species present in only one sample and low abundance species. Genera with average abundance of <1% were grouped into Others <1%. W represents water, M for mud, MO for mud outlet, E for epilithic microbial layer from dry mud, and D for fumarolic deposits samples. Samples labeled a and b are replicates sampled from the same spot. Sample SMP-M_a1 and SMP-M_a1.2 are sequencing replicates.

Next, we analyzed the effect of both taxonomic profiling methods on diversity (Shannon diversity, observed richness, Pielou’s evenness). Before calculating the diversity, we identified and removed kitome (Supplementary Figure 3; Supplementary Tables 4, 5), low abundance taxa and sequences that were unclassified at the domain level from the taxonomy table to avoid an overestimation of diversity and a misrepresentation of community composition (Karstens et al., 2019). Our results showed that Shannon diversity and observed richness were consistently higher when the taxonomy table from metagenomic profiling was used as input compared to when we used the taxonomy table from metagenome-derived 16S rRNA gene profiling (Figures 5A,B). Overall, our results reflect that the 16S amplicon approach, which has been the most employed approach for studying CHSs’ microbiome, identifies a significantly lower number of bacterial species and completely excludes viruses and fungi which is similar to reports from other studies (Ranjan et al., 2016; Brumfield et al., 2020). Although metagenome sequencing permits the identification of more taxa, the choice between 16S rRNA gene and metagenome sequencing ultimately depends on the ecological questions and objectives of any given study. Since metagenome sequences encompass all members of the microbiome and provide high-resolution diversity analysis, which agrees with our objective to understand the full range of microbial ecology in available substrates from Solfatara-Pisciarelli CHSs, the further analyses and results hereinafter are based on the taxonomy table generated from metagenomic profiling.

FIGURE 5

Figure 5. Alpha diversity across locations and substrate types calculated using species abundance from (A) metagenome-derived 16S rRNA gene and (B) metagenomic profiling.

Microbial diversity of Solfatara-Pisciarelli hydrothermal systems

Shannon diversity in Solfatara-Pisciarelli CHSs ranged from 1.20 to 4.26 (Figure 5B; Supplementary Table 8). We observed significant variations in Shannon diversity (Kruskal–Wallis: p < 0.05) between Solfatara and Pisciarelli which were supported by Nonpareil diversity values (Supplementary Figure 4; Supplementary Table 8); on average, Shannon diversity was lower at Pisciarelli (Shannon: 1.20–3.52) than at Solfatara (Shannon: 3.55–4.26). Additionally, microbial communities were more uneven (Kruskal–Wallis: p < 0.05) at Pisciarelli (Pielou’s evenness: 0.26–0.74) than at Solfatara (Pielou’s evenness: 0.67–0.84). However, there was no significant difference in observed richness between these locations (Solfatara: 114–201 species; Pisciarelli: 102–166 species; Kruskal–Wallis: p = 0.24). The decreased Shannon diversity at Pisciarelli may have resulted from the combined effects of reducing conditions, extremely low pH, and high temperature. The addition of highly reducing conditions at Pisciarelli to an already extreme environment (high temperature and low pH) may cause further selection against certain microbial groups by hindering their growth and metabolism, thus leading to the lower Shannon diversity observed (Kimbrough et al., 2006; Seo and DeLaune, 2010; Husson, 2013; Zhang et al., 2015). Alpha diversity in substrates showed some variation (e.g., on average, mud and fumarolic deposits had higher Shannon diversity than water, and richness showed less variability in water and mud); however, these differences were not statistically significant (Kruskal–Wallis, Shannon: p = 0.97; Pielou’s evenness: p = 0.83; observed genus richness: p = 0.09) even when substrates were analyzed separately for each location which implies that species alpha diversity does not differ across substrates in the Solfatara-Pisciarelli CHSs. This result contrast with other studies comparing microbial diversity in water and solid substrates from hot springs, which showed differences in richness (Colman et al., 2016) and evenness (Cole et al., 2013) across substrates.

A comparison of diversity and physical parameters showed pH (Shannon: R² = 0.73; Nonpareil: R² = 0.43, p < 0.05) to be the primary driver of microbial diversity which is consistent with several studies that report pH as the major predictor of microbial diversity in some hydrothermal systems (Oliverio et al., 2018; Power et al., 2018). Results also showed that Eh had a significant correlation (Shannon: R² = 0.62; Nonpareil: R² = 0.43, p < 0.05) with diversity. Temperature on the other hand had minimal effect on diversity (Shannon: R² = 0.12, p = 0.09; Nonpareil: R² = 0.02, p = 0.25). This is in contrast to previous studies that reported a strong relationship between temperature and microbial diversity in hot springs (Miller et al., 2009; Sharp et al., 2014). Comparably, minimal effect of temperature on alpha diversity has been reported for hot spring communities in Yellowstone, United States, Fludir, Iceland and Tibetan Plateau, China (Wang et al., 2013; Podar et al., 2020).

Microbial community composition of Solfatara-Pisciarelli hydrothermal systems

After removing kitome and all potential contaminants, a total of 490 species (245 genera) were identified across all samples corresponding to 11 archaeal, 31 bacterial, 10 eukaryotic, and four viral phyla (Supplementary Table 6). Our results showed that overall, 33% (133) of the species detected were shared across the different types of substrates but each substrate contained unique species among which water had the least number of unique taxa and mud contained the highest (Figure 6B). NMDS plot (Figure 7) generated based on Bray–Curtis dissimilarity matrix revealed that the microbial community structure of Pisciarelli was distinct from Solfatara (ANOSIM: R = 0.63 for species; R = 0.67 for genus, p < 0.05). On average, the most abundant phylum at Pisciarelli was Crenarchaeota (46%) followed by Proteobacteria (5%) and Actinobacteria (2%), whereas Solfatara was dominated by Firmicutes (23%), Proteobacteria (21%), and Actinobacteria (12%; Figure 3B). In addition, Pisciarelli substrates had a higher number of viral sequences (4–49%) than Solfatara (1–4%; Supplementary Figure 2). Interestingly, we detected Cyanobacteria (0.5 to 4%) especially at Solfatara even in samples (SF-D_and SF-D_b) with a temperature (89°C) that exceeds the known temperature limit for photosynthesis in acidic environments (Brock, 1985; Cox et al., 2011; Hamilton et al., 2012). The Cyanobacteria we detected may have been introduced into the hydrothermal system from surrounding cooler environments (by wind or fluid circulations). However, we cannot decipher with certainty if they were environmentally introduced or if they are true members of the Solfatara-Pisciarelli hydrothermal community. We observed that at the species level, both locations shared 171 species (35%), however, there were over 100 species unique to each location (Figure 6A). Species unique to Pisciarelli include thermophiles such as Pyrolobus fumarii, Pyrobaculum aerophilum and Metallosphaera sedula. Unique species identified at Solfatara include Thermaerobacter sp. FW80, Methylacidiphilum kamchatkense and Ilumatobacter coccineus.

FIGURE 6

Figure 6. Venn diagram showing the number of unique and shared species between (A) two locations and (B) different types of substrates.

FIGURE 7

Figure 7. Nonmetric Dimensional Scaling (NMDS) ordination based on Bray–Curtis dissimilarity matrix. Bray–Curtis dissimilarity was calculated using species abundance from metagenomic profiling. The shapes define substrate types (water, mud, and deposits) and the colors represent locations (Pisciarelli and Solfatara).

Our results showed that the microbial community structure in the high temperature environments (74–95°C) of Pisciarelli was different. The source of Pisciarelli large mud pool (84°C) was dominated by Acidianus (21–68%), Pyrobaculum (11–21%) and Bicaudavirus (11–43%) but as temperature decreased away from the source (79°C), we detected a decrease in the abundance of Acidianus (20–23%) and Bicaudavirus (15–17%) but an increase in the abundance of Pyrobaculum (30–43%) in the discharge channel. Water from Pisciarelli small mud pool (89°C) was also dominated by Acidianus (48–59%), but Pyrobaculum accounted for <1% of total sequences in this mud pool. Interestingly, the viral genus, Ampullavirus accounted for over 80% of sequences in Pisciarelli epilithic microbial layer (74°C) and for 26% of sequences in both Pisciarelli fumarolic deposit (93°C) and mud from Pisciarelli small mud pool (88°C). Most of the viruses detected in Pisciarelli fumarolic deposits, epilithic microbial layer and mud from the small mud pool were known archaeal viruses (e.g., Ampullavirus, Bicaudavirus, and Betalipothrixvirus) even though archaea accounted for <3% of the total sequences in these samples. This difference in relative abundance between archaeal viruses and their known host suggests that each archaeal species may host more than one virus type or the archaeal viruses present in Solfatara-Pisciarelli CHSs may have a broader host range that includes bacteria or/and eukaryotes (Munson-McGee et al., 2018). When substrates were compared at each location, we observed that Pisciarelli substrates harbored significantly different (ANOSIM species: R = 0.64; genus: R = 0.62, p < 0.05) microbial communities. SIMPER analysis showed that the sulfur-oxidizing archaeal genus Acidianus was significantly (p < 0.05) more abundant in Pisciarelli water communities, the archaeal genus Pyrobaculum showed a higher relative abundance in mud communities, and Ampullavirus showed a higher relative abundance in fumarole deposits communities (Supplementary Table 9; Figure 4B). Overall, Pisciarelli water and mud communities were similar but distinct from communities in Pisciarelli fumarolic deposits. The distinction between communities in the mud pools and fumarolic deposits at Pisciarelli may be linked to fluids and/or may reflect the availability of oxygen. Although the entire Pisciarelli environment has high temperatures and is highly reducing, the vigorous mixing of hydrothermal gas with air may create pockets of atmospheric O₂ in pore spaces of fumarolic deposits which support aerobic communities that may be absent or in low abundance in mud pools. In contrast, high temperatures of mud pools inhibit the dissolution of oxygen thereby creating microaerobic environments that may favor thermophilic facultative and obligate anaerobes such as Acidianus and Pyrobaculum. Furthermore, the presence of Pyrobaculum specifically Pyrobaculum arsenaticum which is a strict anaerobe (Huber et al., 2000) and members of the genus Sulfuracidifex who are all obligate aerobes in Pisciarelli large pool may be an indication that conditions in this pool fluctuate between microaerobic and completely anaerobic, possibly in congruence with the intensity of hydrothermal activity.

At Solfatara, we observed variability in genus abundance across samples even between replicates (Figure 4B). For example, Acidibacillus accounted for 11% of sequences in SMP-M_a1 but was <1% in replicate mud sample, SMP-M_b. The difference between replicates may derive from the heterogeneous nature of Solfatara-Pisciarelli CHSs; however, seeing that genus abundance also varied in our sequencing replicates (SMP-M_a1 and SMP-M_a1.2), the variability in genus abundance may likely be bias introduced from sequencing. Results also showed that a large percentage of Solfatara sequences (20–36%) were unclassified at the genus level, suggesting a considerable amount of potentially novel genera at this environment. Interestingly, at Solfatara, we detected several genera that are not known thermophiles and not commonly found at pH < 2 (e.g., Rhizopus; 1.7–14%, Actinomadura; 1–9%, Terasakiella; 0.7–5.6%, and Clostridioides; 0.3–2.5%; Wang et al., 2014; Strazzulli et al., 2020). These genera were present in Solfatara water (43°C), mud (69°C), and fumarolic deposits (89°C). Furthermore, we detected these genera in most Pisciarelli samples but to a lesser degree. These microorganisms may not be native to the hydrothermal systems but likely introduced into the habitat from extensive human activity around the crater. In addition, the fluids in mud pools are a mixture of hydrothermal fluids and meteoric water which could deliver genetic material from non-native microorganisms to the hydrothermal systems. This possible introduction of non-native microorganisms into CHSs may be a potential problem for microbial ecology studies since the very nature of the CHSs makes it difficult to completely exclude all environmental contaminants. In contrast to Solfatara water and fumarolic deposits communities, Solfatara mud communities was dominated by Sulfobacillus (10–34%) and Acidibacillus (4–20%). However, unlike Pisciarelli, the differences across Solfatara substrates were not significant (ANOSIM species: R = 0.03, p = 0.50; genus: R = 0.03, p = 0.43) suggesting that variations in microbial community structure across substrates is not a general phenomenon but is specific to hydrothermal systems and possibly linked to the physical and geochemical conditions of individual systems (Wang et al., 2014). The microbial community structure observed in our study is in contrast with studies that investigated the microbial ecology of Solfatara-Pisciarelli CHSs and identified Acidithiobacillus and Metallosphaera as the dominant genus in Solfatara and Pisciarelli, respectively (Iacono et al., 2020; Strazzulli et al., 2020). Furthermore, the dominance of archaea which we observed in Pisciarelli mud pools is in contrast with a study that reported a low archaea/bacteria ratio in high temperature waters (>85°C) of Pisciarelli (Crognale et al., 2022). This contrast between our study and other Solfatara-Pisciarelli studies is likely the consequence of a microbial diversity that reflects differences in sampling year/periods from active CHSs with variable geochemistry. Additionally, the dissimilarity may be related to differences in sampling points (10 cm depth vs. 1 cm depth), sample handling, DNA extraction protocol, bias introduced by 16S rRNA gene primers or sequencing technology. An assessment of the effects of environmental variables on microbial community composition showed that pH and Eh (Mantel: ρ = 0.53 and ρ = 0.47, respectively, p < 0.05) had the strongest correlation to beta diversity which is consistent with reports from other studies that pH and redox potential shape microbial community composition (Alsop et al., 2014; Power et al., 2018). The concentrations of NH₄⁺ and SO₄²⁻ (Mantel: ρ = 0.38 and ρ = 0.23, respectively, p < 0.05) also had significant correlations to beta diversity. However, temperature had no significant (Mantel: ρ = 0.1, p = 0.14) effect on the microbial community composition of Solfatara-Pisciarelli CHSs, which coincides with temperature trends observed for alpha diversity.

Diversity, abundance and metabolic potential of Solfatara-Pisciarelli MAGs

To identify the metabolic potential of members of the microbial community and their connection to geochemistry, we generated MAGs. After co-assembly and binning, a total of 22 MAGs (10 from Solfatara and 12 from Pisciarelli) with >50% completeness and <10% contamination were recovered. Out of the 22 MAGs recovered, 16 had completeness >90% and contamination <5%. Taxonomic classification revealed that MAGs belonged to three bacterial (Firmicutes, Proteobacteria, and Aquificota) and 2 archaeal (Thermoproteota also known as Crenarchaeota and Thermoplasmatota) phyla (Supplementary Tables 10, 11), with over 40% (9 out of 22) of unclassified genomes at the genus level suggesting that our MAGs represent novel genera. MAGs belonging to the phylum Thermoproteota were most represented in Pisciarelli water and mud while Firmicutes were more abundant in all Solfatara substrates and Pisciarelli fumarolic deposits (Figure 8A), which is consistent with our results from metagenome fragments. On average, the most abundant MAGs in Pisciarelli water and mud belong to species of the Acidianus and Pyrobaculum genera and an unknown species of the family Sulfolobaceae whereas the most abundant MAGs in fumarolic deposits belong to 2 unknown species: one species of the Alicyclobacillaceae family and the other species of the Alicyclobacillales order and a third species belonging to the Acidianus genus (Supplementary Table 11). In contrast, the most abundant MAG in all Solfatara samples belong to species of the Acidianus genus. When MAGs were dereplicated based on FastANI of ≥95% we found that only 2 species (Acidianus infernus; and an unknown species of the Alicyclobacillaceae family) were shared between locations further confirming that although an overlap exists between both locations; they have different microbial community structures. Since acidic hydrothermal environments are dominated by chemolithoautotrophic communities that derive energy from inorganic compounds, we analyzed MAGs for marker genes involved in carbon fixation, sulfur, nitrogen, and methane metabolism.

FIGURE 8

Figure 8. Taxonomy and functional potential of Metagenome Assembled Genomes (MAGs) in Solfatara-Pisciarelli hydrothermal systems. (A) Normalized mean (relative) abundance of sequences assigned to each MAG at the phylum/class level in each substrate. (B) Tile plot of selected marker genes (y-axis) present in each MAG (x-axis). Green color represents presence of gene in MAG. MAGs on the x-axis are colored based on phylum/class. Each MAG is labeled either by genus level assignment (g) or by the lowest taxonomic level assigned (f; family and o; order) to the MAG. (C) Sulfur cycle and (D) Nitrogen cycle. Phyla/Class from the different locations participating in each process are indicated. Gene participating in each process are indicated in red text. Black arrows denote processes that were not found in MAGs.

Carbon Metabolism

Organic and inorganic carbon isotope (wt% and δ¹³C) values (Figures 9A,C) were very similar indicating that carbon measured at both locations is mostly organic and falls within the range typical for microorganisms with varied carbon fixation pathways (Sharp, 2007; Havig et al., 2011), which is also indicated by the functional potential of our MAGs. Genes encoding the key enzymes for Calvin-Benson-Bassham (CBB) cycle, ribulose-biphosphate carboxylase (rbc) and phosphoribulokinase (prk) were found in Firmicutes MAGs (Sulfobacillus and R501) as well as Gammaproteobacteria MAG (Acidithiobacillus; Supplementary Table 12; Figure 8B) similar to what has been reported in other studies (Caldwell et al., 2007; Cerqueira et al., 2018). Owing to the abundance of Firmicutes in Pisciarelli fumarolic deposits and Solfatara microbial communities, CBB may be the primary mode of carbon fixation in these environments. The rbc gene was also detected in an archaeal MAG (Desulfurococcus); however, in archaea, this gene encodes enzymes reported to be involved in the reductive hexulose-phosphate pathway (Kono et al., 2017). The marker gene encoding 4-hydroxybutyryl-CoA dehydratase (abfD), an enzyme in 3-hydroxypropionate/4-hydroxybutyrate (3-HP/4-HB) cycle, was present in all MAGs of the Sulfolobaceae family (Acidianus, Metallosphaera, Sulfuracidifex) which is consistent with the presence of the 3-HP/4-HB cycle in the crenarchaeal order Sulfolobales (Berg et al., 2010). Crenarchaeota especially the order Sulfolobales were abundant in Pisciarelli water and mud; hence, 3-HP/4-HB cycle may be the primary mode of carbon fixation in Pisciarelli mud pools. The abfD gene was also detected in one Firmicutes MAG (unknown species of the Sulfobacillaceae family) from Solfatara. The 3-HP/4-HB cycle to the best of our knowledge has not been reported in Sulfobacillaceae and the absence of other genes encoding enzymes involved in the 3-HP/4-HB cycle implies that this cycle may be absent in this Sulfobacillaceae MAG. One possible explanation is that the enzyme encoded by the abfD gene in this Sulfobacillaceae MAG may be involved in an alternative metabolic pathway like has been reported for Clostridium (Scherf and Buckel, 1993; Scherf et al., 1994). Another possible explanation is that the Sulfobacillaceae MAG may have acquired the abfD gene via horizontal gene transfer.

FIGURE 9

Figure 9. Carbon, nitrogen, and sulfur isotope values measured in Solfatara and Pisciarelli hydrothermal systems. (A) Total Carbon (TC) wt.%, Total Organic Carbon (TOC) wt.% and δ¹³C values. The wt.% values of TC and TOC are very similar indicating that carbon measured in the samples is mostly organic. The δ¹³C values exhibit a wide range of values. (B) Total Nitrogen wt.% and δ¹⁵N values. Pisciarelli mud and water have the highest concentrations of N. The δ¹⁵N values exhibit a large variation with values ranging from 0.41 to −34.04‰. (C) C:N ratio and δ¹³C values both point out that most of the carbon measured shows values typical for microorganisms with various carbon fixation pathways. (D) Total Sulfur wt.% and δ³⁴S exhibit values characteristic for oxidized fumarolic H₂S and sulfur with sedimentary origin.

Genes encoding the key enzymes in Wood-Ljungdahl (WL) pathway and 3-hydroxyprioponate (3-HP) bicycle were absent in all MAGs and in the entire metagenome dataset, suggesting that both pathways may not be important for carbon fixation at both locations. The marker gene encoding ATP citrate lyase (aclAB), an enzyme involved in reductive tricarboxylic acid (rTCA) cycle, was also absent in all MAGs. ATP citrate lyase catalyzes citrate cleavage in rTCA cycle but in Aquificaceae, this reaction is catalyzed by the combined action of citryl-CoA synthetase (ccsAB) and citryl-CoA lyase (ccl) (Aoshima et al., 2004a,b) which were found in our Aquificota MAG (UBA11096) from Pisciarelli metagenome indicating potential for rTCA.

Potential for methanogenesis was absent in MAGs from both locations. Although the gene encoding hydrogen disulfide reductase (hdrABC) involved in methanogenesis was found in Firmicutes (Sulfobacillus, R501, and unknown species of the Sulfobacillaceae family), Gammaproteobacteria (Acidithiobacillus), Aquificota (UBA11096), and Thermoproteota (Metallosphaera) MAGs, the gene encoding key enzyme, methyl-CoM reductase (mcr) in methanogenesis was absent. Given the absence of mcr gene, the absence of genes encoding other enzymes involved in methanogenesis especially in the Metallosphaera MAG that was 100% complete, and the absence of δ¹³C values lower than −30‰, it is likely that the hdrABC genes in these MAGS are not used in methanogenesis but are most likely used in sulfur metabolism (Mander et al., 2004; Wang et al., 2019). Potential for methane oxidation (methanol dehydrogenase; mdh1 and particulate methane monooxygenase; pmo) was also absent in all MAGs.

Sulfur metabolism

Sulfur isotope values (δ³⁴S) measured in deposits from both locations (Figure 9D) are characteristic of sulfur with sedimentary origin, but δ³⁴S values (−0.83 to 2.62‰) measured in water and mud where oxidation of H₂S is extensive, fall within the range for SO₄²⁻ formed from sulfide oxidation (−2.5 to 2.4‰; Allard et al., 1991). In acidic environments like Solfatara and Pisciarelli, abiotic oxidation of H₂S-rich gas produced from the disproportionation of volcanic SO₂ is inhibited (D’Imperio et al., 2008), thus providing a great habitat for sulfur-oxidizing microorganisms that utilize H₂S for energy generation. The abundance in H₂S is consistent with the presence of genes encoding enzymes involved in the oxidation of a variety of reduced sulfur compounds (H₂S, S⁰ or S₂O₃²⁻) at both locations. Genes encoding the enzyme, sulfide:quinone oxidoreductases (sqr) for H₂S oxidation to elemental sulfur (S⁰) were detected in all MAGs except one Thermoproteota MAG (Desulfurococcus; Figures 8B,C). The presence of sqr gene in most MAGs is indicative of the importance of H₂S oxidation to microbial communities at both locations. S⁰ formed can be reduced to H₂S by the enzyme molybdopterin sulfur reductase (sreABC) or disproportionated to H₂S, thiosulfate (S₂O₃²⁻) and sulfite (SO₃²⁻) by sulfur oxygenase reductase (SOR). The sreABC gene was found in Thermoproteota MAG (Acidianus) while SOR gene was found in Firmicutes (Sulfobacillus, R501, and Acidibacillus) Gammaproteobacteria (Acidithiobacillus) and Thermoproteota (Acidianus and Sulfuracidifex) MAGs.

Thiosulfate (S₂O₃²⁻) may not be available for microbial metabolism since it disproportionates to S⁰ and SO₃²⁻ at pH < 4 (Nordstrom et al., 2004); however, microbial communities exhibited potential to utilize S₂O₃²⁻. Thiosulfate can be oxidized to SO₄²⁻ via the Sox system (soxXYZABCD) which we detected in Gammaproteobacteria (Acidithiobacillus) and Aquificota (UBA11096) MAGs. In both MAGs, soxCD was absent suggesting that they oxidize S₂O₃²⁻ to S⁰ instead of SO₄²⁻ which is consistent with studies that identified the sox cluster without soxCD in Acidithiobacillus (Wang et al., 2019). Thiosulfate can also be oxidized to tetrathionate (S₄O₆²⁻) by thiosulfate:quinol oxidoreductase (doxDA) or reduced to SO₃²⁻ by thiosulfate:cyanide sulfur transferase (tst). The doxDA gene was found in all Sulfolobaceae MAGs (Acidianus, Metallosphaera, Sulfuracidifex) and in the Firmicutes (R501) MAG, while tst gene was found in 2 Thermoproteota (Pyrobaculum, Metallosphaera) MAGs and most Firmicutes MAGs. H₂S can be regenerated by reducing SO₄²⁻ via dissimilatory sulfate reduction (DSR). Genes encoding enzymes (dissimilatory sulfite reductase; dsrAB, adenylylsulfate reductase; aprAB, and sulfate adenylyltransferase; sat) involved in DSR were detected in one Thermoproteota (Pyrobaculum) MAG. The presence of both genes encoding enzymes involved in sulfur oxidation and reduction in some MAGs implies that microbial communities in Solfatara-Pisciarelli systems may take advantage of the high energy yield that results from coupling sulfur oxidation to sulfur reduction, which has been previously observed in other microbial community studies (Glamoclija et al., 2004; Crognale et al., 2022). Our results are consistent with studies that found the abundance of organisms capable of utilizing sulfur (Inskeep et al., 2013) and an enrichment of genes (sqr, doxDA, and SOR) involved in metabolizing reduced sulfur compounds in acidic hot springs (Colman et al., 2019a).

Nitrogen metabolism

Our results showed that nitrogen metabolism may not be an important energy-yielding pathway for microbial communities even in Pisciarelli mud pools where NH₄⁺ was abundant (Figures 8B,D). Exceptionally low δ¹⁵N values were observed at both locations (Figure 9B) and the lowest δ¹⁵N values ever measured were found in the Pisciarelli epilithic microbial layer. These extremely low δ¹⁵N values are difficult to explain in the context of metabolic potential present in our MAGs or metagenome fragments. Low and negative δ¹⁵N values are sometimes linked to microbial ammonia oxidation and/or nitrogen fixation (Tozer et al., 2005; Havig et al., 2011). Although Crognale et al. (2022) reported ammonia oxidation to be the most likely metabolic pathway in Pisciarelli mud pools due to its high energy yield, the absence of genes encoding enzymes involved in ammonia oxidation (ammonia monooxygenase; amoCAB and hydroxylamine dehydrogenase; hao) or nitrogen fixation (nifD, nifK, and nifH) in our MAGs and metagenome dataset indicates that microbial nitrogen metabolism cannot explain the geochemistry we observed. Previously observed negative δ¹⁵N values (about −10‰) have been associated with the dissolution of magmatic nitrogen in water at elevated temperatures and pressures (Labidi et al., 2020). While it is very likely, the mechanism by which magmatic nitrogen, high temperature and pressure would contribute to the extremely low δ¹⁵N values measured at Solfatara and Pisciarelli is unclear.

Evidence of denitrification which has not been linked to low and negative δ¹⁵N values was found in five MAGs (Sulfobacillus, R501, unknown species of the Alicyclobacillaceae family, Pyrobaculum, and UBA11096). However, none of the MAGs carried all the genes encoding enzymes for the complete denitrification pathway (nitrate reductase, narGH or napA; nitrite reductase, nirK; nitric oxide reductase, norB; nitrous oxide reductase, nosZ). The complete set of genes for denitrification were also absent in metagenome fragments. One Pisciarelli MAG (UBA11096) belonging to the family Aquificaceae, carried most of the denitrification genes except nosZ which is consistent with studies that found nitrous oxide (N₂O) to be the end-product of denitrification in members of Aquificaceae (Nakagawa et al., 2004). Dissimilatory nitrate reduction to ammonia via nitrite reductases (nirBD) was detected only in the Aquificaceae MAG from Pisciarelli.

Conclusion

Results from this study showed that Solfatara and Pisciarelli are lithologically similar but geochemically distinct. At Pisciarelli, we observed varying geochemistry among mud pools and fumarolic deposits, which results from subsurface fluid–rock interactions, whereas Solfatara substrates were geochemically similar. Shannon diversity was significantly different between locations but showed no significant difference across substrates. We found pH to be the most important driver of alpha diversity and microbial community composition. NMDS plot and ANOSIM showed that Solfatara’s microbial community structure exhibited no significant difference across substrates which also coincides with observed trends in geochemistry. In contrast, at Pisciarelli, microbial community structure followed trends in fluid availability rather than geochemical trends. Geochemically, water and mud from Pisciarelli large mud pool and water from Pisciarelli small mud pool clustered distinctly from fumarolic deposits, discharge channel mud and mud from Pisciarelli small pool. However, based on microbial communities, all water and mud samples clustered distinctly from fumarolic deposits that were the driest substrates analyzed; results showed that Pisciarelli fumarolic deposits were more similar to Solfatara samples than samples from Pisciarelli mud pools. Overall, the genus Acidianus dominated Pisciarelli water and Pyrobaculum was more abundant in Pisciarelli mud. Interestingly, Pisciarelli fumarolic deposits, epilithic microbial layer, and mud from small mud pool had a high number of viral sequences. Solfatara and Pisciarelli had distinct microbial communities; on average, Acidianus and Pyrobaculum dominated Pisciarelli while Sulfobacillus and Acidibacillus were dominant at Solfatara. Although microbial communities were distinct between the two locations and between Pisciarelli substrates, MAGs indicated functional redundancy. For example, the potential to reduce H₂S to S⁰ was observed in almost all MAGs. Further, MAGs showed that microbial communities inhabiting both locations took advantage of the major volcanic gases (CO₂ and H₂S) via carbon fixation and sulfur oxidation. The 3HP/4-HB pathway found in most archaeal MAGs is most likely the primary mode of carbon fixation in Pisciarelli mud pools dominated by (hyper)thermophilic archaea whereas CBB cycle is most likely the primary carbon fixation pathway in all other environments where bacteria was abundant. At both locations, sulfur cycling is represented by oxidation of H₂S-rich volcanic gas and other reduced sulfur compounds. Comparative geochemical and metagenomic analyses 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 CHSs to understand the full range of microbial ecology to avoid biased ecological assessment.

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/, Bioproject accession PRJNA889931.

Author contributions

MG, AS, MF, GN, and CT designed and organized the project. MG, RS, MP, AM collected the samples. IU, RB, MF, and MG collected the data. IU, MG, and MF contributed to the data processing and analyses. IU and MG wrote the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This research was enabled through the Alfred P. Sloan Foundation’s support of the Deep Carbon Observatory Deep Earth Carbon Degassing program (DECADE) to MG and Rutgers Faculty Program Start Up. IU was supported by Rutgers University Transform Graduate Fellowship.

Acknowledgments

We would like to thank Ying Lin at UCR for re-running samples to confirm δ¹⁵N values and Ashley E. Murphy for running samples on XRF. We thank the Office of Advanced Research Computing (OARC) at Rutgers, The State University of New Jersey for providing access to the Amarel cluster and the associated resources that have contributed to the results reported in this manuscript. We thank two reviewers for their thoughtful and constructive comments which have significantly improved this manuscript.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2023.1066406/full#supplementary-material

SUPPLEMENTARY FIGURE 1 | Ballot based on Principal Component Analysis (PCA) using bulk elemental composition of mud and fumarolic deposits.

SUPPLEMENTARY FIGURE 2 | Taxonomic profile of Solfatara and Pisciarelli microorganisms based on metagenomic profiling. Sequences not classified at the domain level and species with relative abundance less than 0.01% were grouped into unclassified and low abundance species. All kitome were grouped into contaminant. Samples collected from Pisciarelli large pool are denoted by PLP, Pisciarelli small pool by PSP, Pisciarelli fumarole by PF, Solfatara mud pool by SMP, and Solfatara fumarole by SF. W represents water, M for mud, MO for mud outlet, E for epilithic microbial layer from dry mud wall, and D for deposits samples. Samples labeled a and b are replicates sampled from the same spot. Sample SMP-M_a1 and SMP-M_a1.2 are sequencing replicates.

SUPPLEMENTARY FIGURE 3 | Relative abundance of kitome families. Kitome families with average abundance of <1% were grouped into Others <1%.

SUPPLEMENTARY FIGURE 4 | Estimation of the coverage of microbial communities in (A) Solfatara and (B) Pisciarelli. The empty circles on the curves represent community coverage estimate at the sequencing effort applied. The top and bottom horizontal red dashed indicate 100 and 95% average community coverage, respectively. The arrows are the Nonpareil diversity estimates.

Footnotes

1. ^https://www.bioinformatics.babraham.ac.uk/projects/fastqc/

References

Aiuppa, A., Tamburello, G., Di Napoli, R., Cardellini, C., Chiodini, G., Giudice, G., et al. (2013). First observations of the fumarolic gas output from a restless caldera: implications for the current period of unrest (2005–2013) at Campi Flegrei. Geochem. Geophys. 14, 4153–4169. doi: 10.1002/ggge.20261

CrossRef Full Text | Google Scholar

Allard, P., Maiorani, A., Tedesco, D., Cortecci, G., and Turi, B. (1991). Isotopic study of the origin of sulfur and carbon in Solfatara fumaroles, Campi Flegrei caldera. J. Volcanol. Geotherm. Res. 48, 139–159. doi: 10.1016/0377-0273(91)90039-3

CrossRef Full Text | Google Scholar

Alneberg, J., Bjarnason, B. S., De Bruijn, I., Schirmer, M., Quick, J., Ijaz, U. Z., et al. (2014). Binning metagenomic contigs by coverage and composition. Nat. Methods 11, 1144–1146. doi: 10.1038/nmeth.3103

PubMed Abstract | CrossRef Full Text | Google Scholar

Alsop, E. B., Boyd, E. S., and Raymond, J. (2014). Merging metagenomics and geochemistry reveals environmental controls on biological diversity and evolution. BMC Ecol. 14:16. doi: 10.1186/1472-6785-14-16

PubMed Abstract | CrossRef Full Text | Google Scholar

Amenabar, M. J., and Boyd, E. S. (2019). A review of the mechanisms of mineral-based metabolism in early earth analog rock-hosted hydrothermal ecosystems. World J. Microbiol. Biotechnol. 35:29. doi: 10.1007/s11274-019-2604-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Amenabar, M. J., Shock, E. L., Roden, E. E., Peters, J. W., and Boyd, E. S. (2017). Microbial substrate preference dictated by energy demand rather than supply. Nat. Geosci. 10, 577–581. doi: 10.1038/ngeo2978

PubMed Abstract | CrossRef Full Text | Google Scholar

Aoshima, M., Ishii, M., and Igarashi, Y. (2004a). A novel enzyme, citryl-CoA lyase, catalysing the second step of the citrate cleavage reaction in Hydrogenobacter thermophilus TK-6. Mol. Microbiol. 52, 763–770. doi: 10.1111/j.1365-2958.2004.04010.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Aoshima, M., Ishii, M., and Igarashi, Y. (2004b). A novel enzyme, citryl-CoA synthetase, catalysing the first step of the citrate cleavage reaction in Hydrogenobacter thermophilus TK-6. Mol. Microbiol. 52, 751–761. doi: 10.1111/j.1365-2958.2004.04009.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Arif, S., Willenberg, C., Dreyer, A., Nacke, H., and Hoppert, M. (2021). Sasso Pisano geothermal field environment Harbours diverse Ktedonobacteria representatives and illustrates habitat-specific adaptations. Microorganisms 9:1402. doi: 10.3390/microorganisms9071402

PubMed Abstract | CrossRef Full Text | Google Scholar

Aulitto, M., Gallo, G., Puopolo, R., Mormone, A., Limauro, D., Contursi, P., et al. (2021). Genomic insight of Alicyclobacillus Mali FL18 isolated from an arsenic-rich hot spring. Front. Microbiol. 12:669. doi: 10.3389/fmicb.2021.639697

CrossRef Full Text | Google Scholar

Bankevich, A., Nurk, S., Antipov, D., Gurevich, A. A., Dvorkin, M., Kulikov, A. S., et al. (2012). SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 19, 455–477. doi: 10.1089/cmb.2012.0021

PubMed Abstract | CrossRef Full Text | Google Scholar

Benson, C. A., Bizzoco, R. W., Lipson, D. A., and Kelley, S. T. (2011). Microbial diversity in nonsulfur, sulfur and iron geothermal steam vents. FEMS Microbiol. Ecol. 76, 74–88. doi: 10.1111/j.1574-6941.2011.01047.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Berg, I. A., Kockelkorn, D., Ramos-Vera, W. H., Say, R. F., Zarzycki, J., Hügler, M., et al. (2010). Autotrophic carbon fixation in archaea. Nat. Rev. Microbiol. 8, 447–460. doi: 10.1038/nrmicro2365

CrossRef Full Text | Google Scholar

Bolger, A. M., Lohse, M., and Usadel, B. (2014). Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120. doi: 10.1093/bioinformatics/btu170

PubMed Abstract | CrossRef Full Text | Google Scholar

Boyd, E. S., Hamilton, T. L., Spear, J. R., Lavin, M., and Peters, J. W. (2010). [FeFe]-hydrogenase in Yellowstone National Park: evidence for dispersal limitation and phylogenetic niche conservatism. ISME J. 4, 1485–1495. doi: 10.1038/ismej.2010.76

PubMed Abstract | CrossRef Full Text | Google Scholar

Boyer, G. M., Schubotz, F., Summons, R. E., Woods, J., and Shock, E. L. (2020). Carbon oxidation state in microbial polar lipids suggests adaptation to hot spring temperature and redox gradients. Front. Microbiol. 11:229. doi: 10.3389/fmicb.2020.00229

PubMed Abstract | CrossRef Full Text | Google Scholar

Brock, T. D. (1985). Life at high temperatures. Science 230, 132–138. doi: 10.1126/science.230.4722.132

CrossRef Full Text | Google Scholar

Brumfield, K. D., Huq, A., Colwell, R. R., Olds, J. L., and Leddy, M. B. (2020). Microbial resolution of whole genome shotgun and 16S amplicon metagenomic sequencing using publicly available NEON data. PLoS One 15:e0228899. doi: 10.1371/journal.pone.0228899

PubMed Abstract | CrossRef Full Text | Google Scholar

Caldwell, P. E., MacLean, M. R., and Norris, P. R. (2007). Ribulose bisphosphate carboxylase activity and a Calvin cycle gene cluster in Sulfobacillus species. Microbiology 153, 2231–2240. doi: 10.1099/mic.0.2007/006262-0

CrossRef Full Text | Google Scholar

Caliro, S., Chiodini, G., Moretti, R., Avino, R., Granieri, D., Russo, M., et al. (2007). The origin of the fumaroles of La Solfatara (Campi Flegrei, South Italy). Geochim. Cosmochim. Acta 71, 3040–3055. doi: 10.1016/j.gca.2007.04.007

CrossRef Full Text | Google Scholar

Cerqueira, T., Barroso, C., Froufe, H., Egas, C., and Bettencourt, R. (2018). Metagenomic signatures of microbial communities in deep-sea hydrothermal sediments of Azores vent fields. Microb. Ecol. 76, 387–403. doi: 10.1007/s00248-018-1144-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Chaumeil, P.-A., Mussig, A. J., Hugenholtz, P., and Parks, D. H. (2020). GTDB-Tk: a toolkit to classify genomes with the genome taxonomy database. Bioinformatics 36, 1925–1927. doi: 10.1093/bioinformatics/btz848

CrossRef Full Text | Google Scholar

Chiodini, G., Avino, R., Caliro, S., and Minopoli, C. (2011). Temperature and pressure gas geoindicators at the Solfatara fumaroles (Campi Flegrei). Ann. Geophys. 54:107245. doi: 10.4401/ag-5002

CrossRef Full Text | Google Scholar

Chiodini, G., Caliro, S., Avino, R., Bini, G., Giudicepietro, F., De Cesare, W., et al. (2021). Hydrothermal pressure-temperature control on CO₂ emissions and seismicity at Campi Flegrei (Italy). J. Volcanol. Geotherm. Res. 414:107245. doi: 10.1016/j.jvolgeores.2021.107245

CrossRef Full Text | Google Scholar

Chiodini, G., Caliro, S., Cardellini, C., Granieri, D., Avino, R., Baldini, A., et al. (2010). Long-term variations of the Campi Flegrei, Italy, volcanic system as revealed by the monitoring of hydrothermal activity. J. Geophys. Res. Solid Earth 115:B03205. doi: 10.1029/2008JB006258

CrossRef Full Text | Google Scholar

Chiodini, G., Frondini, F., Cardellini, C., Granieri, D., Marini, L., and Ventura, G. (2001). CO₂ degassing and energy release at Solfatara volcano, Campi Flegrei, Italy. J. Geophys. Res. Solid Earth 106, 16213–16221. doi: 10.1029/2001JB000246

CrossRef Full Text | Google Scholar

Ciniglia, C. V. (2005). Influences of geochemical and mineralogical constraints on algal distribution in acidic hydrothermal environments: Pisciarelli (Naples, Italy) as a model site. Arch. Hydrobiol. 162, 121–142. doi: 10.1127/0003-9136/2005/0162-0121

CrossRef Full Text | Google Scholar

Cole, J. K., Peacock, J. P., Dodsworth, J. A., Williams, A. J., Thompson, D. B., Dong, H., et al. (2013). Sediment microbial communities in great boiling spring are controlled by temperature and distinct from water communities. ISME J. 7, 718–729. doi: 10.1038/ismej.2012.157

PubMed Abstract | CrossRef Full Text | Google Scholar

Colman, D. R., Feyhl-Buska, J., Robinson, K. J., Fecteau, K. M., Xu, H., Shock, E. L., et al. (2016). Ecological differentiation in planktonic and sediment-associated chemotrophic microbial populations in Yellowstone hot springs. FEMS Microbiol. Ecol. 92:137. doi: 10.1093/femsec/fiw137

PubMed Abstract | CrossRef Full Text | Google Scholar

Colman, D. R., Lindsay, M. R., Amenabar, M. J., and Boyd, E. S. (2019a). The intersection of geology, geochemistry, and microbiology in continental hydrothermal systems. Astrobiology 19, 1505–1522. doi: 10.1089/ast.2018.2016

PubMed Abstract | CrossRef Full Text | Google Scholar

Colman, D. R., Lindsay, M. R., and Boyd, E. S. (2019b). Mixing of meteoric and geothermal fluids supports hyperdiverse chemosynthetic hydrothermal communities. Nat. Commun. 10:681. doi: 10.1038/s41467-019-08499-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Cox, A., Shock, E. L., and Havig, J. R. (2011). The transition to microbial photosynthesis in hot spring ecosystems. Chem. Geol. 280, 344–351. doi: 10.1016/j.chemgeo.2010.11.022

CrossRef Full Text | Google Scholar

Crognale, S., Venturi, S., Tassi, F., Rossetti, S., Cabassi, J., Capecchiacci, F., et al. (2022). Geochemical and microbiological profiles in hydrothermal extreme acidic environments (Pisciarelli Spring, Campi Flegrei, Italy). FEMS Microbiol. Ecol. 98:fiac088. doi: 10.1093/femsec/fiac088

PubMed Abstract | CrossRef Full Text | Google Scholar

Crognale, S., Venturi, S., Tassi, F., Rossetti, S., Rashed, H., Cabassi, J., et al. (2018). Microbiome profiling in extremely acidic soils affected by hydrothermal fluids: the case of the Solfatara crater (Campi Flegrei, southern Italy). FEMS Microbiol. Ecol. 94:fiy190. doi: 10.1093/femsec/fiy190

PubMed Abstract | CrossRef Full Text | Google Scholar

De Natale, G., Troise, C., Mark, D., Mormone, A., Piochi, M., Di Vito, M. A., et al. (2016). The Campi Flegrei deep drilling project (CFDDP): new insight on caldera structure, evolution and hazard implications for the Naples area (Southern Italy). Geochem. Geophys. 17, 4836–4847. doi: 10.1002/2015GC006183

CrossRef Full Text | Google Scholar

DeSantis, T. Z., Hugenholtz, P., Larsen, N., Rojas, M., Brodie, E. L., Keller, K., et al. (2006). Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl. Environ. Microbiol. 72, 5069–5072. doi: 10.1128/AEM.03006-05

PubMed Abstract | CrossRef Full Text | Google Scholar

D’Imperio, S., Lehr, C. R., Oduro, H., Druschel, G., Kühl, M., and McDermott, T. R. (2008). Relative importance of H2 and H2S as energy sources for primary production in geothermal springs. Appl. Environ. Microbiol. 74, 5802–5808. doi: 10.1128/AEM.00852-08

PubMed Abstract | CrossRef Full Text | Google Scholar

Ellis, D. G., Bizzoco, R. W., and Kelley, S. T. (2008). Halophilic archaea determined from geothermal steam vent aerosols. Environ. Microbiol. 10, 1582–1590. doi: 10.1111/j.1462-2920.2008.01574.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Fedele, A., Somma, R., Troise, C., Holmberg, K., De Natale, G., and Matano, F. (2021). Time-lapse landform monitoring in the Pisciarelli (Campi Flegrei-Italy) fumarole field using UAV photogrammetry. Remote Sens. 13:118. doi: 10.3390/rs13010118

CrossRef Full Text | Google Scholar

Fournier, R. O. (1989). Geochemistry and dynamics of the Yellowstone National Park hydrothermal system. Annu. Rev. Earth Planet. Sci. 17, 13–53. doi: 10.1146/annurev.ea.17.050189.000305

CrossRef Full Text | Google Scholar

Fullerton, K. M., Schrenk, M. O., Yücel, M., Manini, E., Basili, M., Rogers, T. J., et al. (2021). Effect of tectonic processes on biosphere–geosphere feedbacks across a convergent margin. Nat. Geosci. 14, 301–306. doi: 10.1038/s41561-021-00725-0

CrossRef Full Text | Google Scholar

Glamoclija, M., Garrel, L., Berthon, J., and López-García, P. (2004). Biosignatures and bacterial diversity in hydrothermal deposits of Solfatara Crater, Italy. Geomicrobiol J. 21, 529–541. doi: 10.1080/01490450490888235

CrossRef Full Text | Google Scholar

Hamilton, T., Vogl, K., Bryant, D., Boyd, E., and Peters, J. (2012). Environmental constraints defining the distribution, composition, and evolution of chlorophototrophs in thermal features of Yellowstone National Park. Geobiology 10, 236–249. doi: 10.1111/j.1472-4669.2011.00296.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Havig, J. R., Raymond, J., Meyer-Dombard, D. A. R., Zolotova, N., and Shock, E. L. (2011). Merging isotopes and community genomics in a siliceous sinter-depositing hot spring. J. Geophys. Res. 116:G01005. doi: 10.1029/2010JG001415

CrossRef Full Text | Google Scholar

Hou, W., Wang, S., Dong, H., Jiang, H., Briggs, B. R., Peacock, J. P., et al. (2013). A comprehensive census of microbial diversity in hot springs of Tengchong, Yunnan Province China using 16S rRNA gene pyrosequencing. PLoS One 8:e53350. doi: 10.1371/journal.pone.0053350

PubMed Abstract | CrossRef Full Text | Google Scholar

Huber, R., Sacher, M., Vollmann, A., Huber, H., and Rose, D. (2000). Respiration of arsenate and selenate by hyperthermophilic archaea. Syst. Appl. Microbiol. 23, 305–314. doi: 10.1016/S0723-2020(00)80058-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Husson, O. (2013). Redox potential (eh) and pH as drivers of soil/plant/microorganism systems: a transdisciplinary overview pointing to integrative opportunities for agronomy. Plant Soil 362, 389–417. doi: 10.1007/s11104-012-1429-7

CrossRef Full Text | Google Scholar

Iacono, R., Cobucci-Ponzano, B., De Lise, F., Curci, N., Maurelli, L., Moracci, M., et al. (2020). Spatial metagenomics of three geothermal sites in Pisciarelli hot spring focusing on the biochemical resources of the microbial consortia. Molecules 25:4023. doi: 10.3390/molecules25174023

PubMed Abstract | CrossRef Full Text | Google Scholar

Inskeep, W. P., Jay, Z. J., Tringe, S. G., Herrgard, M., and Rusch, D. B. (2013). The YNP metagenome project: environmental parameters responsible for microbial distribution in the Yellowstone geothermal ecosystem. Front. Microbiol. 4:67. doi: 10.3389/fmicb.2013.00067

PubMed Abstract | CrossRef Full Text | Google Scholar

Jain, C., Rodriguez-R, L. M., Phillippy, A. M., Konstantinidis, K. T., and Aluru, S. (2018). High throughput ANI analysis of 90K prokaryotic genomes reveals clear species boundaries. Nat. Commun. 9:5114. doi: 10.1038/s41467-018-07641-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Kanehisa, M., Sato, Y., and Morishima, K. (2016). BlastKOALA and GhostKOALA: KEGG tools for functional characterization of genome and metagenome sequences. J. Mol. Biol. 428, 726–731. doi: 10.1016/j.jmb.2015.11.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Kang, D. D., Li, F., Kirton, E., Thomas, A., Egan, R., An, H., et al. (2019). MetaBAT 2: an adaptive binning algorithm for robust and efficient genome reconstruction from metagenome assemblies. PeerJ 7:e7359. doi: 10.7717/peerj.7359

PubMed Abstract | CrossRef Full Text | Google Scholar

Karstens, L., Asquith, M., Davin, S., Fair, D., Gregory, W. T., Wolfe, A. J., et al. (2019). Controlling for contaminants in low-biomass 16S rRNA gene sequencing experiments. mSystems 4:e00290-19. doi: 10.1128/mSystems.00290-19

PubMed Abstract | CrossRef Full Text | Google Scholar

Kimbrough, D. E., Kouame, Y., Moheban, P., and Springthorpe, S. (2006). The effect of electrolysis and oxidation–reduction potential on microbial survival, growth, and disinfection. Int. J. Environ. Pollut. 27, 211–221. doi: 10.1504/IJEP.2006.010464

CrossRef Full Text | Google Scholar

Kono, T., Mehrotra, S., Endo, C., Kizu, N., Matusda, M., Kimura, H., et al. (2017). A RuBisCO-mediated carbon metabolic pathway in methanogenic archaea. Nat. Commun. 8:14007. doi: 10.1038/ncomms14007

PubMed Abstract | CrossRef Full Text | Google Scholar

Labidi, J., Barry, P. H., Bekaert, D. V., Broadley, M. W., Marty, B., Giunta, T., et al. (2020). Hydrothermal 15 N 15 N abundances constrain the origins of mantle nitrogen. Nature 580, 367–371. doi: 10.1038/s41586-020-2173-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Lezcano, M. Á., Moreno-Paz, M., Carrizo, D., Prieto-Ballesteros, O., Fernández-Martínez, M. Á., Sánchez-García, L., et al. (2019). Biomarker profiling of microbial mats in the geothermal band of Cerro Caliente, Deception Island (Antarctica): life at the edge of heat and cold. Astrobiology 19, 1490–1504. doi: 10.1089/ast.2018.2004

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, D., Liu, C.-M., Luo, R., Sadakane, K., and Lam, T.-W. (2015). MEGAHIT: an ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph. Bioinformatics 31, 1674–1676. doi: 10.1093/bioinformatics/btv033

PubMed Abstract | CrossRef Full Text | Google Scholar

Lindsay, M. R., Amenabar, M. J., Fecteau, K. M., Debes, R. V., Fernandes Martins, M. C., Fristad, K. E., et al. (2018). Subsurface processes influence oxidant availability and chemoautotrophic hydrogen metabolism in Yellowstone hot springs. Geobiology 16, 674–692. doi: 10.1111/gbi.12308

PubMed Abstract | CrossRef Full Text | Google Scholar

Lindsay, M. R., Colman, D. R., Amenabar, M. J., Fristad, K. E., Fecteau, K. M., Debes, R. V., et al. (2019). Probing the geological source and biological fate of hydrogen in Yellowstone hot springs. Environ. Microbiol. 21, 3816–3830. doi: 10.1111/1462-2920.14730

PubMed Abstract | CrossRef Full Text | Google Scholar

Lowenstern, J. B., Bergfeld, D., Evans, W. C., and Hunt, A. G. (2015). Origins of geothermal gases at Yellowstone. J. Volcanol. Geotherm. Res. 302, 87–101. doi: 10.1016/j.jvolgeores.2015.06.010

CrossRef Full Text | Google Scholar

Mander, G. J., Pierik, A. J., Huber, H., and Hedderich, R. (2004). Two distinct heterodisulfide reductase-like enzymes in the sulfate-reducing archaeon Archaeoglobus profundus. Eur. J. Biochem. 271, 1106–1116. doi: 10.1111/j.1432-1033.2004.04013.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Mardanov, A. V., Gumerov, V. M., Beletsky, A. V., and Ravin, N. V. (2018). Microbial diversity in acidic thermal pools in the Uzon caldera, Kamchatka. Antonie Van Leeuwenhoek 111, 35–43. doi: 10.1007/s10482-017-0924-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Marlow, J. J., Colocci, I., Jungbluth, S. P., Weber, N. M., Gartman, A., and Kallmeyer, J. (2020). Mapping metabolic activity at single cell resolution in intact volcanic fumarole sediment. FEMS Microbiol. Lett. 367:fnaa031. doi: 10.1093/femsle/fnaa031

PubMed Abstract | CrossRef Full Text | Google Scholar

Martini, M., Giannini, L., Buccianti, A., Prati, F., Legittimo, P. C., Iozzelli, P., et al. (1991). 1980–1990: ten years of geochemical investigation at Phlegrean fields (Italy). J. Volcanol. Geotherm. Res. 48, 161–171. doi: 10.1016/0377-0273(91)90040-7

CrossRef Full Text | Google Scholar

Medrano-Santillana, M., Souza-Brito, E. M., Duran, R., Gutierrez-Corona, F., and Reyna-López, G. E. (2017). Bacterial diversity in fumarole environments of the Paricutín volcano, Michoacán (Mexico). Extremophiles 21, 499–511. doi: 10.1007/s00792-017-0920-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Menzel, P., Gudbergsdóttir, S. R., Rike, A. G., Lin, L., Zhang, Q., Contursi, P., et al. (2015). Comparative metagenomics of eight geographically remote terrestrial hot springs. Microb. Ecol. 70, 411–424. doi: 10.1007/s00248-015-0576-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Menzel, P., Ng, K. L., and Krogh, A. (2016). Fast and sensitive taxonomic classification for metagenomics with kaiju. Nat. Commun. 7:11257. doi: 10.1038/ncomms11257

PubMed Abstract | CrossRef Full Text | Google Scholar

Meyer, F., Paarmann, D., D’Souza, M., Olson, R., Glass, E. M., Kubal, M., et al. (2008). The metagenomics RAST server–a public resource for the automatic phylogenetic and functional analysis of metagenomes. BMC Bioinformatics 9:386. doi: 10.1186/1471-2105-9-386

PubMed Abstract | CrossRef Full Text | Google Scholar

Miller, S. R., Strong, A. L., Jones, K. L., and Ungerer, M. C. (2009). Bar-coded pyrosequencing reveals shared bacterial community properties along the temperature gradients of two alkaline hot springs in Yellowstone National Park. Appl. Environ. Microbiol. 75, 4565–4572. doi: 10.1128/AEM.02792-08

PubMed Abstract | CrossRef Full Text | Google Scholar

Moreras-Marti, A., Fox-Powell, M., Zerkle, A. L., Stueeken, E., Gazquez, F., Brand, H. E., et al. (2021). Volcanic controls on the microbial habitability of Mars-analogue hydrothermal environments. Geobiology 19, 489–509. doi: 10.1111/gbi.12459

PubMed Abstract | CrossRef Full Text | Google Scholar

Moretti, R., De Natale, G., and Troise, C. (2020). Hydrothermal Versus Magmatic: Geochemical Views and Clues Into the Unrest Dilemma at Campi FlegreiVesuvius, Campi Flegrei, and Campanian Volcanism. Amsterdam, Netherlands: Elsevier, 371–406.

Google Scholar

Munson-McGee, J. H., Peng, S., Dewerff, S., Stepanauskas, R., Whitaker, R. J., Weitz, J. S., et al. (2018). A virus or more in (nearly) every cell: ubiquitous networks of virus–host interactions in extreme environments. ISME J. 12, 1706–1714. doi: 10.1038/s41396-018-0071-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Nakagawa, S., Nakamura, S., Inagaki, F., Takai, K., Shirai, N., and Sako, Y. (2004). Hydrogenivirga caldilitoris gen. Nov., sp. nov., a novel extremely thermophilic, hydrogen-and sulfur-oxidizing bacterium from a coastal hydrothermal field. Int. J. Syst. Evol. Microbiol. 54, 2079–2084. doi: 10.1099/ijs.0.03031-0

CrossRef Full Text | Google Scholar

Nordstrom, D., Ball, J., and McCleskey, R. (2004). “Oxidation reactions for reduced Fe, As, and S in thermal outflows of” in Water-Rock Interaction. eds. R. Wanty and R. Seal (London, UK: Taylor and Francis Group), 59–62.

Google Scholar

Nordstrom, D. K., Ball, J. W., and McCleskey, R. B. (2005). “Ground water to surface water: geochemistry of thermal outflows in Yellowstone National Park” in Geothermal Biology and Geochemistry in Yellowstone National Park. eds. W. Inskeep and T. R. McDermott (Bozeman, MT: Montana State University), 73–94.

Google Scholar

Nordstrom, D. K., McCleskey, R. B., and Ball, J. W. (2009). Sulfur geochemistry of hydrothermal waters in Yellowstone National Park: IV acid–sulfate waters. Appl. Geochem. 24, 191–207. doi: 10.1016/j.apgeochem.2008.11.019

CrossRef Full Text | Google Scholar

Oliverio, A. M., Power, J. F., Washburne, A., Cary, S. C., Stott, M. B., and Fierer, N. (2018). The ecology and diversity of microbial eukaryotes in geothermal springs. ISME J. 12, 1918–1928. doi: 10.1038/s41396-018-0104-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Parks, D. H., Imelfort, M., Skennerton, C. T., Hugenholtz, P., and Tyson, G. W. (2015). CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res. 25, 1043–1055. doi: 10.1101/gr.186072.114

PubMed Abstract | CrossRef Full Text | Google Scholar

Piochi, M., Kilburn, C., Di Vito, M., Mormone, A., Tramelli, A., Troise, C., et al. (2014). The volcanic and geothermally active Campi Flegrei caldera: an integrated multidisciplinary image of its buried structure. Int. J. Earth Sci. 103, 401–421. doi: 10.1007/s00531-013-0972-7

CrossRef Full Text | Google Scholar

Piochi, M., Mormone, A., Balassone, G., Strauss, H., Troise, C., and De Natale, G. (2015). Native sulfur, sulfates and sulfides from the active Campi Flegrei volcano (southern Italy): genetic environments and degassing dynamics revealed by mineralogy and isotope geochemistry. J. Volcanol. Geotherm. Res. 304, 180–193. doi: 10.1016/j.jvolgeores.2015.08.017

CrossRef Full Text | Google Scholar

Podar, P. T., Yang, Z., Björnsdóttir, S. H., and Podar, M. (2020). Comparative analysis of microbial diversity across temperature gradients in hot springs from yellowstone and Iceland. Front. Microbiol. 11:1625. doi: 10.3389/fmicb.2020.01625

PubMed Abstract | CrossRef Full Text | Google Scholar

Poichi, M. M. A., Strauss, H., and Balassone, G. (2019). The acid sulfate zone and the mineral alteration styles of the Roman Puteoli (Neapolitan area, Italy): clues on fluid fracturing progression at the Campi Flegrei volcano. Solid Earth 10, 1809–1831. doi: 10.5194/se-10-1809-2019

CrossRef Full Text | Google Scholar

Power, J. F., Carere, C. R., Lee, C. K., Wakerley, G. L., Evans, D. W., Button, M., et al. (2018). Microbial biogeography of 925 geothermal springs in New Zealand. Nat. Commun. 9:2876. doi: 10.1038/s41467-018-05020-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Puopolo, R., Gallo, G., Mormone, A., Limauro, D., Contursi, P., Piochi, M., et al. (2020). Identification of a new heavy-metal-resistant strain of Geobacillus stearothermophilus isolated from a hydrothermally active volcanic area in southern Italy. Int. J. Environ. Res. Public Health 17:2678. doi: 10.3390/ijerph17082678

PubMed Abstract | CrossRef Full Text | Google Scholar

R Core Team (2020). R: A language and environment for statistical computing. R Foundation for Statistical Computing. Vienna, Austria. Available at: https://www.R-project.org/

Google Scholar

Ranjan, R., Rani, A., Metwally, A., McGee, H. S., and Perkins, D. L. (2016). Analysis of the microbiome: advantages of whole genome shotgun versus 16S amplicon sequencing. Biochem. Biophys. Res. Commun. 469, 967–977. doi: 10.1016/j.bbrc.2015.12.083

PubMed Abstract | CrossRef Full Text | Google Scholar

Rodriguez-R, L. M., Gunturu, S., Tiedje, J. M., Cole, J. R., and Konstantinidis, K. T. (2018). Nonpareil 3: fast estimation of metagenomic coverage and sequence diversity. mSystems 3:e00039-18. doi: 10.1128/mSystems.00039-18

PubMed Abstract | CrossRef Full Text | Google Scholar

Rolandi, G., De Natale, G., Kilburn, C. R., Troise, C., Somma, R., Di Lascio, M., et al. (2020a). The 39 ka Campanian Ignimbrite Eruption: New Data on Source Area in the Campanian Plain Vesuvius, Campi Flegrei, and Campanian Volcanism. Amsterdam, Netherlands: Elsevier, 175–205.

Google Scholar

Rolandi, G., Di Lascio, M., and Rolandi, R. (2020b). The Neapolitan Yellow Tuff Eruption as the Source of the Campi Flegrei Caldera Vesuvius, Campi Flegrei, and Campanian Volcanism. Amsterdam, Netherlands: Elsevier, 273–296.

Google Scholar

Scherf, U., and Buckel, W. (1993). Purification and properties of an iron-sulfur and FAD-containing 4-hydroxybutyryl-CoA dehydratase/vinylacetyl-CoA delta 3-delta 2-isomerase from clostridium aminobutyricum. Eur. J. Biochem. 215, 421–429. doi: 10.1111/j.1432-1033.1993.tb18049.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Scherf, U., Söhling, B., Gottschalk, G., Linder, D., and Buckel, W. (1994). Succinate-ethanol fermentation in clostridium kluyveri: purification and characterisation of 4-hydroxybutyryl-CoA dehydratase/vinylacetyl-CoA Δ3-Δ2-isomerase. Arch. Microbiol. 161, 239–245. doi: 10.1007/BF00248699

PubMed Abstract | CrossRef Full Text | Google Scholar

Seemann, T. (2014). Prokka: rapid prokaryotic genome annotation. Bioinformatics 30, 2068–2069. doi: 10.1093/bioinformatics/btu153

PubMed Abstract | CrossRef Full Text | Google Scholar

Seo, D. C., and DeLaune, R. D. (2010). Effect of redox conditions on bacterial and fungal biomass and carbon dioxide production in Louisiana coastal swamp forest sediment. Sci. Total Environ. 408, 3623–3631. doi: 10.1016/j.scitotenv.2010.04.043

PubMed Abstract | CrossRef Full Text | Google Scholar

Sharp, Z. (2007). Principles of Stable Isotope Geochemistry. Pearson Education Inc., Upper Saddle River, NJ

Google Scholar

Sharp, C. E., Brady, A. L., Sharp, G. H., Grasby, S. E., Stott, M. B., and Dunfield, P. F. (2014). Humboldt’s spa: microbial diversity is controlled by temperature in geothermal environments. ISME J. 8, 1166–1174. doi: 10.1038/ismej.2013.237

PubMed Abstract | CrossRef Full Text | Google Scholar

Sheik, C. S., Reese, B. K., Twing, K. I., Sylvan, J. B., Grim, S. L., Schrenk, M. O., et al. (2018). Identification and removal of contaminant sequences from ribosomal gene databases: lessons from the census of deep life. Front. Microbiol. 9:840. doi: 10.3389/fmicb.2018.00840

PubMed Abstract | CrossRef Full Text | Google Scholar

Shock, E. L., Holland, M., Amend, J. P., Osburn, G., and Fischer, T. P. (2010). Quantifying inorganic sources of geochemical energy in hydrothermal ecosystems, Yellowstone National Park, USA. Geochim. Cosmochim. Acta 74, 4005–4043. doi: 10.1016/j.gca.2009.08.036

CrossRef Full Text | Google Scholar

Solorzano, L. (1969). Determination of ammonia in natural waters by the phenol hypochlorite method. Limnol. Oceanogr. 14, 799–801.

Google Scholar

Strazzulli, A., Cobucci-Ponzano, B., Iacono, R., Giglio, R., Maurelli, L., Curci, N., et al. (2020). Discovery of hyperstable carbohydrate-active enzymes through metagenomics of extreme environments. FEBS J. 287, 1116–1137. doi: 10.1111/febs.15080

PubMed Abstract | CrossRef Full Text | Google Scholar

Tozer, W. C., Hackell, D., Miers, D., and Silvester, W. (2005). Extreme isotopic depletion of nitrogen in New Zealand lithophytes and epiphytes; the result of diffusive uptake of atmospheric ammonia? Oecologia 144, 628–635. doi: 10.1007/s00442-005-0098-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Troiano, A., Di Giuseppe, M. G., Patella, D., Troise, C., and De Natale, G. (2014). Electromagnetic outline of the Solfatara–Pisciarelli hydrothermal system, Campi Flegrei (southern Italy). J. Volcanol. Geotherm. Res. 277, 9–21. doi: 10.1016/j.jvolgeores.2014.03.005

CrossRef Full Text | Google Scholar

Troise, C., De Natale, G., Schiavone, R., Somma, R., and Moretti, R. (2019). The Campi Flegrei caldera unrest: discriminating magma intrusions from hydrothermal effects and implications for possible evolution. Earth-Sci. Rev. 188, 108–122. doi: 10.1016/j.earscirev.2018.11.007

CrossRef Full Text | Google Scholar

Uritskiy, G. V., DiRuggiero, J., and Taylor, J. (2018). MetaWRAP—a flexible pipeline for genome-resolved metagenomic data analysis. Microbiome 6:158. doi: 10.1186/s40168-018-0541-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Valentino, G., and Stanzione, D. (2003). Source processes of the thermal waters from the Phlegraean fields (Naples, Italy) by means of the study of selected minor and trace elements distribution. Chem. Geol. 194, 245–274. doi: 10.1016/S0009-2541(02)00196-1

CrossRef Full Text | Google Scholar

Wall, K., Cornell, J., Bizzoco, R. W., and Kelley, S. T. (2015). Biodiversity hot spot on a hot spot: novel extremophile diversity in Hawaiian fumaroles. MicrobiologyOpen 4, 267–281. doi: 10.1002/mbo3.236

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, S., Dong, H., Hou, W., Jiang, H., Huang, Q., Briggs, B. R., et al. (2014). Greater temporal changes of sediment microbial community than its waterborne counterpart in Tengchong hot springs, Yunnan Province, China. Sci. Rep. 4:7479. doi: 10.1038/srep07479

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, S., Hou, W., Dong, H., Jiang, H., Huang, L., Wu, G., et al. (2013). Control of temperature on microbial community structure in hot springs of the Tibetan plateau. PLoS One 8:e62901. doi: 10.1371/journal.pone.0062901

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, R., Lin, J.-Q., Liu, X.-M., Pang, X., Zhang, C.-J., Yang, C.-L., et al. (2019). Sulfur oxidation in the acidophilic autotrophic Acidithiobacillus spp. Front. Microbiol. 9:3290. doi: 10.3389/fmicb.2018.03290

PubMed Abstract | CrossRef Full Text | Google Scholar

Weyrich, L. S., Farrer, A. G., Eisenhofer, R., Arriola, L. A., Young, J., Selway, C. A., et al. (2019). Laboratory contamination over time during low-biomass sample analysis. Mol. Ecol. Resour. 19, 982–996. doi: 10.1111/1755-0998.13011

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, Y.-W., Simmons, B. A., and Singer, S. W. (2016). MaxBin 2.0: an automated binning algorithm to recover genomes from multiple metagenomic datasets. Bioinformatics 32, 605–607. doi: 10.1093/bioinformatics/btv638

PubMed Abstract | CrossRef Full Text | Google Scholar

Xie, W., Zhang, C. L., Wang, J., Chen, Y., Zhu, Y., de la Torre, J. R., et al. (2015). Distribution of ether lipids and composition of the archaeal community in terrestrial geothermal springs: impact of environmental variables. Environ. Microbiol. 17, 1600–1614. doi: 10.1111/1462-2920.12595

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Z., Wang, H., Zhou, J., Li, H., He, Z., Van Nostrand, J. D., et al. (2015). Redox potential and microbial functional gene diversity in wetland sediments under simulated warming conditions: implications for phosphorus mobilization. Hydrobiologia 743, 221–235. doi: 10.1007/s10750-014-2039-6

CrossRef Full Text | Google Scholar