Skip to main content

ORIGINAL RESEARCH article

Front. Microbiol., 11 August 2023
Sec. Extreme Microbiology
This article is part of the Research Topic Rising Stars in Hydrothermal Vents and Cold Seeps: 2021 View all 10 articles

Mapping the microbial diversity associated with different geochemical regimes in the shallow-water hydrothermal vents of the Aeolian archipelago, Italy

  • 1Department of Biology, University of Naples “Federico II”, Naples, Italy
  • 2Blue Marine Foundation, London, United Kingdom
  • 3School of Marine and Atmospheric Sciences, Stony Brook, NY, United States
  • 4Department of Biochemistry and Microbiology, Rutgers University, New Brunswick, NJ, United States
  • 5Department of Marine and Coastal Science, Rutgers University, New Brunswick, NJ, United States
  • 6Istituto per le Risorse Biologiche e Biotecnologiche Marine, Consiglio Nazionale Delle Ricerche, CNR-IRBIM, Ancona, Italy
  • 7Earth-Life Science Institute, Tokyo Institute of Technology, Ookayama, Tokyo, Japan
  • 8Marine Chemistry and Geochemistry Department–Woods Hole Oceanographic Institution, Woods Hole, MA, United States

Shallow-water hydrothermal vents are unique marine environments ubiquitous along the coast of volcanically active regions of the planet. In contrast to their deep-sea counterparts, primary production at shallow-water vents relies on both photoautotrophy and chemoautotrophy. Such processes are supported by a range of geochemical regimes driven by different geological settings. The Aeolian archipelago, located in the southern Tyrrhenian sea, is characterized by intense hydrothermal activity and harbors some of the best sampled shallow-water vents of the Mediterranean Sea. Despite this, the correlation between microbial diversity, geochemical regimes and geological settings of the different volcanic islands of the archipelago is largely unknown. Here, we report the microbial diversity associated with six distinct shallow-water hydrothermal vents of the Aeolian Islands using a combination of 16S rRNA amplicon sequencing along with physicochemical and geochemical measurements. Samples were collected from biofilms, fluids and sediments from shallow vents on the islands of Lipari, Panarea, Salina, and Vulcano. Two new shallow vent locations are described here for the first time. Our results show the presence of diverse microbial communities consistent in their composition with the local geochemical regimes. The shallow water vents of the Aeolian Islands harbor highly diverse microbial community and should be included in future conservation efforts.

1. Introduction

The vast and diverse metabolic repertoire of microorganisms enables their ubiquitous distribution on our planet, making them key mediators of biogeochemical cycles on a planetary scale (Falkowski et al., 2008). Notwithstanding their ubiquity, their distribution and the environmental factors shaping it, are not fully resolved. However, significant correlations between microbial community composition and environmental geochemistry have been established in diverse environments (Jorgensen et al., 2012; Alsop et al., 2014; Liu et al., 2014; Yamamoto et al., 2019). Geothermal environments are ideal systems to study the complex relationships between community composition and geochemistry due to the wide geochemical gradients present across proximal geographic distances, and relatively constrained microbial communities (Power et al., 2018). Additionally, geothermal ecosystems harbor diverse communities of extremophiles that are an untapped source of potential biotechnological compounds (Sayed et al., 2020). Hence, there is an increasing necessity to characterize in detail the microbial diversity in these ecosystems and to prioritize their conservation.

Marine hydrothermal vents are dynamic systems where water is discharged through the Earth’s crust in areas where there is enough heat to drive fluid circulation, e.g., plate boundaries and mid-ocean ridges (Price and Giovannelli, 2017; Brovarone et al., 2020). The interaction between the seawater with the hot, deep-subsurface lithologies, promotes the leaching of inorganic reduced species, such as volatiles and metals (McDermott et al., 2018). Upon arrival at the surface, the mixing of the reduced hydrothermal fluids (enriched in electron donors) with the oxidized sea water (enriched in electron acceptors), creates a thermodynamic disequilibrium that fuels diverse microbial consortia and sustains complex food webs (Dahle et al., 2015; Price et al., 2015; Reed et al., 2015). In contrast with their deep-sea counterparts, shallow-water hydrothermal vents (SWHV) are closer to the surface, thus they are strongly influenced by solar energy. In these shallow systems, the primary production is reliant upon a mixture of phototrophy and chemolithotrophy, making them high-energy environments (Tarasov et al., 2005; Gomez-Saez et al., 2017). Additionally, the shallower depths and the resultant decrease in pressure allow for the presence of a generally abundant free gas phase (Bastianoni et al., 2022). The proximity to landmasses promotes the input of terrigenous organic carbon to the system, while tidal waves and tides make SWHV transient ecosystems subject to strong dynamic forcing (Giovannelli et al., 2013; Yücel et al., 2013; Price and Giovannelli, 2017).

The shallow-water hydrothermal vents of the Aeolian archipelago are some of the most studied marine shallow vents in the world. The archipelago, located north of Sicily, is composed of seven volcanic islands (Stromboli, Vulcano, Panarea, Salina, Lipari, Alicudi, and Filicudi) and several sea-mounts. The archipelago lies between the Tyrrhenian sea back-arc basin and the Calabrian fore-arc. While most of its volcanic activity developed during the Quaternary period, it continues to the present day on some of the islands (Vulcano and Stromboli) in the form of volcanic eruptions and secondary geothermal manifestations (Gamberi et al., 1997). The geomorphism of the region is very dynamic, marked by the presence of rock types belonging mainly to the calc-alkaline, shoshonitic, and potassium associations (De Astis et al., 2000; Peccerillo and Turco, 2004). The islands of Lipari, Salina, and Vulcano, are aligned on an NNW–SSE striking fault, while a NNE–SSW to NE–SW striking fault affects the islands of Stromboli and Panarea (De Astis et al., 2003; Favallim et al., 2005). Most of the described marine geothermal activity is present along the coast of the islands of Vulcano and Panarea, where the majority of microbiological studies on the shallow-water vents have been focused so far. Despite this, additional sites of geothermal activities have been described on the other islands (Cardellini et al., 2003; Inguaggiato et al., 2012, 2013).

Previous studies have mapped the thermodynamic chemical reaction landscape of SWHVs of Vulcano and Panarea islands, showing that a wide range of exergonic chemical reactions are capable of sustaining diverse metabolic strategies (Amend et al., 2003; Rogers and Amend, 2005, 2006; Rogers et al., 2007; Price et al., 2015). Accordingly, through metagenomics analysis, amplicon sequencing, and DGGE analysis, previous studies reported complex and diverse microbial communities inhabiting these ecosystems (Gugliandolo and Italiano, 1999; Simmons and Norris, 2002; Rogers and Amend, 2005, 2006; Rusch et al., 2005; Lentini et al., 2007, 2014; Manini et al., 2008; White et al., 2008; Maugeri et al., 2009, 2010, 2013; Gugliandolo et al., 2012, 2015; Placido et al., 2015; Antranikian et al., 2017; Fagorzi et al., 2019; Gugliandolo and Maugeri, 2019; Sciutteri et al., 2022). These communities were found to be enriched in sulfur-oxidizing chemolithoautotrophic bacteria at higher temperatures and at lower temperatures by heterotrophic and photoautotrophic-based life-styles. Furthermore, thermophilic microorganisms have been isolated from several vents (Hafenbradl et al., 1996; Deckert et al., 1998; Lentini et al., 2007; Gugliandolo et al., 2012). However, only a few studies tried to couple in more detail both the geochemistry of the shallow vents and the microbial communities inhabiting them. In addition, these studies were conducted exclusively on the SWHV of Vulcano and Panarea, with the exception of one study, based on culture-dependent approaches that included communities from Lipari and Stromboli (Gugliandolo and Italiano, 1999). Currently, there is a lack of a more comprehensive survey of the microbial communities associated with the diversity of shallow-water hydrothermal vent ecosystems of the Aeolian archipelago.

To address this, and to understand how geochemically different shallow-water hydrothermal vents shape the microbial community composition, we used a combination of 16S rRNA gene amplicon sequencing and geochemical approaches to survey six shallow-water hydrothermal vents located in four different islands of the Archipelago, including two new locations that have never been characterized in the literature before. Our results show that shallow-water hydrothermal vents of the Aeolian islands are ecosystems that support diverse and complex microbial communities, providing important ecological roles to the surrounding marine environment. For this reason, they should be included in current and future conservation efforts.

2. Materials and methods

2.1. Sampling and site description

The samples were retrieved on the AEO19 expedition to the islands of the Aeolian archipelago: Vulcano, Panarea, Salina, and Lipari in the summer of 2019. Six different locations were sampled: Vulcano - Levante Bay (LB), Vulcano - Levante Port (PL); Panarea - Bottaro (BO), Panarea - Baia Calcara (BC), Salina - Pollara (PO), and Lipari - Pietra del bagno (PB) (Figure 1), representing different shallow-water geothermal environments. All the samples were retrieved by SCUBA divers, using sterile falcon tubes and syringes. To accomplish a broader characterization of the different sites, different types of samples were collected: sediments, hydrothermal fluids, and biofilms. In the sites Levante Bay, Pollara, Baia Calcara, and Pietra del Bagno, abundant white microbial mats were observed near the venting sites. The shallow vents across the sampled locations were sulfide-rich, with the exception of the Pietra del Bagno site, where the presence of conspicuous orange precipitates suggests the presence of iron-rich hydrothermal fluids. Intense degassing was observed in the Bottaro venting site. All samples were stored either at +4°C or − 20°C in accordance to the specific protocols for later processing in the lab. Subsamples used for molecular analysis and organic matter determination were stored frozen at −20°C, while samples for the determination of the major ions were filtered through a 0.22 μm filter and stored at +4°C. Temperature was measured in situ using a waterproof thermocouple, while pH, oxidation–reduction potential, conductivity, and dissolved oxygen were measured using a multi-parameter sensor probe (Hanna- multiparameter, model HI98194) on the boat immediately after the dive, with a fluid sample collected with syringes at the seafloor.

FIGURE 1
www.frontiersin.org

Figure 1. Map of the volcanic Aeolian island archipelago, located north of Sicily. Sample sites are shown in red. Major fault systems are illustrated here, as well as the gradient of the different slab depths (100–250 Km).

2.2. Sedimentary organic matter

The determination of the concentration of pigments (chlorophyll-a and phaeopigments) was carried out according to Lorenzen and Downs (1986). Briefly, the samples were supplemented with 90% of acetone and further incubated in the dark, at 4°C, for 12 h. After incubation, in order to remove the sediments, the samples were centrifuged and the concentration of chlorophyll-a and phaeopigments was estimated using a spectrofluorometer (at wavelengths of 430 nm and 630 nm), before and after being supplemented with 200 μL of 0.1 N HCl, respectively. Pigment concentrations were normalized to sediment dry weight (Manini et al., 2003).

Total protein concentrations were determined according to Hartree (1972), and further modified for analysis of sediment samples from Fabiano et al. (1995). The total carbohydrate concentrations were determined according to Gerchakov and Hatcher (1972), and expressed as glucose equivalents. Total lipids were extracted from the sediments by a direct elution with chloroform and methyl alcohol (1:1 v/v), and further determined according to Bligh and Dyer (1959). The readings were performed spectrophotometrically. The concentrations of carbohydrates, proteins and lipids were converted to carbon equivalents using the conversion factors, 0.40, 0.49 and 0.75 μg C × μg−1, respectively, and further normalized to sediment dry weight (Fabiano et al., 1995). The biopolymeric carbon (BPC), calculated through the sum of the concentration of carbohydrates, proteins and lipids, represents the fraction of the total organic carbon potentially available to heterotrophs and benthic consumers in the environment. All measurements are reported as an average of three replicates.

2.3. Geochemical analysis

The concentrations of major cations (sodium-Na+, potassium-K+, magnesium-Mg2+ and calcium-Ca+) and anions (chlorine-Cl, sulfate-SO42− and bromide-Br) of the hydrothermal fluids were measured using ion chromatography (ECO, IC Metrohm) equipped with conductivity detectors in two independent measurements. Calibration curves for the ion species of interest were run in the range of 0.1 and 10 ppm with R ≥ 0.999. In order to reduce the amount of Cl and Na+ entering the system, before the analysis, all samples were filtered (0.22 μm) and diluted to guarantee a conductivity of 600 μS/cm. All dilutions were made using 18.2 MΩ/cm type I water, which was also used as a blank for blank subtractions. Anions were run using a 3.2 mM Na2CO3 + 1 mM NaHCO3 mobile phase on a Metrosep A Supp 5 column equipped with a 0.15 M ortho-Phosphoric acid suppressor. The flow of the anionic eluent was 0.7 mL min-1 for 30 min, effective for good resolution of all anions in the standard solution. Cations were run using a 2.5 mM HNO3 + 0.5 mM (COOH)2 × 2H2O mobile phase on a Metrosep C4 column. The flow of the cationic eluent was 0.9 mL min-1 with a total separation of 35 min, effective for good resolution of all cations in the standard solution. Data acquisition and analysis were carried out through MagIC Net 3.3 software. Additionally, calibration curves were carried out using certified CPA chem external standards for each of the anions and cations analyzed. Detection limits for each of the ions analyzed correspond to 95% confidence levels, and were as follows: Mg2+ (0.08 mM), Na+ (0.3 mM), Ca2+ (0.3 mM), K+ (0.7 mM), SO42− (0.03 mM), Br (0.02 mM), Cl (0.1 mM).

2.4. DNA extraction and amplicon sequencing

DNA extraction was performed using the DNeasy PowerSoil Kit (QIAGEN), following the manufacturer’s instructions, with small modifications. The modifications of the protocol included a further elution step to better recover the DNA from the column. For samples where the PowerSoil kit protocol resulted in low quantity of DNA, a modified phenol-chloroform extraction method was employed (Cordone et al., 2005, 2011; Zanfardino et al., 2016), specifically adapted for the conditions found in shallow-water hydrothermal vents, such as the higher presence of clays and anions (Giovannelli et al., 2013, 2016). Briefly, 0.8 g of sediment was suspended in 850 μL of extraction buffer (100 mM Na2HPO4 pH 8, 100 mM tris–HCL pH 8, 100 mM EDTA pH 8, 1.5 M NaCl, 1% CTAB) and further supplemented with 100 μL of Lysozyme (100 mg/mL), followed by two incubation periods at 37°C for 30 min, the latter after the addition of 5 μL of Proteinase K (1 mg/mL). Subsequently, 50 μL of 20% SDS was added to the samples and incubated 1 h at 65°C, with regular mixing. The sediments were removed through consecutive centrifugations (2 × 14,000 g) after the addition of a solution of phenol:chloroform:isoamyl alcohol (25:24:1/24:1). The supernatant was collected and supplemented with Na-acetate (0.1 vol) and isopropanol 100%, and incubated at room temperature for 12 h. The precipitated DNA was then washed with 70% ice-cold ethanol, and re-suspended in 50 μL of Tris–HCL. DNA was visualized with agarose gel electrophoresis and quantified spectrophotometrically and spectrofluorimetrically (NanoDrop and Qubit). Obtained DNA was sequenced at the Integrated Microbiome Resource (IMR)1 using primers targeting the V4-V5 of the 16S rRNA (515FB = GTGYCAGCMGCCGCGGTAA 926R = CCGYCAATTYMTTTRAGTTT515FB), using Illumina MiSeq technology.

2.5. Bioinformatic and statistical analysis

The raw sequence data was processed using the DADA2 package (Callahan et al., 2017). The analysis of the quality profile was carried out after removing primers and adapters. Sequences with an average call quality for each base between 20 and 40, were retained for downstream analysis. DADA2, through estimation of the error profile for each possible transition, identifies amplicon sequencing variants (ASVs). Taxonomy was assigned to the ASVs by matching with the SILVA Database, release 138.12 (Quast et al., 2013). Taxonomic assignments, together with the variant abundance table, were used to calculate diversity indices and investigate microbial diversity with the Phyloseq package (McMurdie and Holmes, 2013), as previously reported (Cordone et al., 2022). Briefly, the ASVs count table and the assigned taxonomy were combined into a phyloseq object. Mitochondria, Chloroplast, and Eukaryotic sequences were removed. Then, groups related to human pathogens and common DNA extraction contaminants (Sheik et al., 2018) were removed. The remaining reads represented 96.97% of the original reads, with 149,351 reads assigned to 2,145 ASVs preserved for downstream analysis. Diversity analyzes were carried out using the Phyloseq package (McMurdie and Holmes, 2013) with the relative abundance values set to 100%. Top abundance ASVs and genera were further defined as having accumulative abundances above 0.1%. Additionally, the alpha diversity was investigated using Shannon diversity index. The beta diversity was investigated using the UniFrac diversity index (Unweighted and Weighted) as implemented in the vegan package (Oksanen et al., 2018). Subsequently, the resultant similarity matrix was plotted using non-metric multidimensional scaling techniques. The resulting ordination was used to investigate the correlation between the environmental and geochemical variables using environmental fitting (envfit and ordisurf functions in vegan). All statistical investigations, data processing and plotting were carried out in the R statistical software version 4.1.2 (R Core Team, 2021), using the vegan (Oksanen et al., 2018), ggplot2 (Wickham, 2011), and the packages supra-referred. All the sequences analyzed in this study are available through the European Nucleotide Archive (ENA) under project accession PRJEB54611 under the ENA Umbrella project CoEvolve PRJEB55081. A complete R script containing all the steps to reproduce our analysis is available at https://github.com/giovannellilab/Barosa_et_al_Aeolian_shallow_vent_diversity with DOI: https://doi.org/10.5281/zenodo.7480614 together with all the environmental and geochemical data.

3. Results

3.1. Physico-chemical parameters and sedimentary organic matter content

The physico-chemical parameters of the hydrothermal sites sampled in the present study are presented in Table 1. Overall, the physico-chemical parameters measured varied between the sites. The shallow-water vent sediment temperatures at the sampled venting orifice ranged from 28°C in Pollara to 90°C in Baia Calcara. The salinity ranged from 33.0 psu in Levante Bay to 41.5 psu in Pietra del Bagno. Dissolved oxygen (DO) values ranged from 11% in Levante Bay to 57.3% in Pietra del Bagno. Additionally, the oxidation–reduction potential of the different sites was negative for Levante Bay, Bottaro, and Levante Port, while Baia Calcara, Bottaro, and Pietra del Bagno had positive reduction potentials. The pH was slightly acidic for every sample, with the exception of Bottaro, which had circumneutral pH (8).

TABLE 1
www.frontiersin.org

Table 1. Location (GPS coordinates), depth, and the physico-chemical parameters measured on the venting fluids at the different sampled locations.

The sampled sediments showed a wide range of sedimentary organic matter content, with concentrations of biopolymeric organic carbon (BPC) ranging from 99.44 μg C g−1 (± 2.22) in Pollara to 3499.0 μg C g−1 (± 2740.5) in Bottaro. The higher values of BPC in Bottaro, as well as the high variability between the measured replicates, was hypothesized to be related to the presence of a biofilm ingrained in the coarse sediment. Total protein and carbohydrates ranged from 492.88 μg C g−1 (± 39.5) and 39.23 μg C g−1 (± 4.64) in Levante Bay to 6547.72 μg C g−1 (± 599.9) and 1756.92 μg C g−1 (± 279.88) in Bottaro, respectively. Total lipids ranged from 10.17 μgC g−1 (± 3.02) in Pollara to 222.48 μg C g-1 (± 65.91) in Levante Port. Protein to carbohydrate ratios ranged from 1.1 in Pietra del Bagno to 15.2 in Baia Calcara. Chlorophyll-a concentration values ranged from 0.01 μg g−1 (± 0.003) in Levante Bay to 6.9 μg g−1 (± 0.126) in Levante Port. The highest chlorophyll-a concentrations were measured at the moderate temperature sites, Bottaro and Levante Port, while the lowest values were measured at the high temperature sites, Levante Bay and Baia Calcara (Figure 2).

FIGURE 2
www.frontiersin.org

Figure 2. Content and quality of the organic matter measured in the sampled sediments. (A) Biopolymeric organic carbon (BCP); (B) Chlorophyll-a pigments content; (C) Proteins; (D) Carbohydrates; (E) Lipids. All sites are organized by increasing temperatures and are expressed in μg g−1 of dry sediment. Error bars represent standard deviation (n = 3).

3.2. Geochemical analysis

The concentration of the major elements in the hydrothermal fluids were measured in the forms of major anions and cations (Figure 3; Table 2). Chloride concentrations ranged from 535.9 (% RDS 0.38) to 701.9 mM (% RDS 0.38), compared to 618.0 mM in seawater. Bottaro, Pollara and Pietra del Bagno have the highest Cl values. Sulfate concentrations ranged from 22.71 mM (% RDS 0.69) in Baia Calcara to 27.89 mM (% RDS 0.69) in Pietra del Bagno. Sulfate was depleted in the hydrothermal fluids of all the samples compared to seawater concentrations (32.9 mM). Additionally, Br concentrations ranged from 0.74 mM (% RDS 0.59) in Levante Bay to 0.91 mM (% RDS 0.59) in Bottaro, compared to sea water concentrations (0.94 mM). As for the major cations, Mg2+ concentrations ranged from 49.04 mM (% RDS 3.8) in Levante Bay to 62.36 mM (% RDS 3.8) in Bottaro, compared to seawater concentrations of 58.3 mM. Bottaro and Pollara were the only sites where fluids were enriched in Mg2+, relative to seawater. The concentrations of K+, Ca2+ and Na+ ranged from 11.18 mM (% RDS 11.75) in Pollara to 13.23 mM (% RDS 11.75) in Levante Bay, from 8.55 mM (% RDS 10.37) in Pietra del Bagno to 17.87 mM (% RDS 10.37) in Bottaro, and from 445.9 mM (% RDS 7.4) in Levante Port to 551.6 mM (% RDS 7.4) in Bottaro, compared to the seawater concentrations of 11.0 mM, 11.1 mM and 547.0 mM, respectively. Additionally, we combined our data from previous publish data on the geochemistry of the fluids sampled from the shallow water hydrothermal vents of Vulcano, Panarea, and Stromboli (Italiano and Nuccio, 1991; Sedwick and Stüben, 1996; Gugliandolo et al., 2006; Rogers et al., 2007; Sieland, 2009; Tassi et al., 2009; Inguaggiato et al., 2013; Price et al., 2015; Kürzinger, 2019). Fluid data from the present study falls within a mixing trend resembling Type 3 fluids previously described by Price et al., 2015.

FIGURE 3
www.frontiersin.org

Figure 3. Concentrations of the major anions and cations measured in the present study (bold colors), as well as from previous reports from the Literature (opaque colors): (A) SO42− vs. Mg2+, with the three visible Mg2+ end-member extrapolations; (B) K+ vs. Mg2+ with the three fluid mixing trends observed; (C) Na+ vs. Mg2+; (D) Cl vs. Mg2+; (E) SO42− vs. Ca2+ (vector illustrating end-member extrapolation); (F) Br/Cl ratio vs. Cl, where it is possible to infer the presence of halite in solution; (G) Br vs. Cl and extrapolations to zero; (H) Cl vs. Na+; (I) ternary plot with the influence of albitization on the hydrothermal fluids. The literature data was retrieved from Italiano and Nuccio (1991), Sedwick and Stüben (1996), Gugliandolo et al. (2006), Rogers et al. (2007), Sieland (2009), Tassi et al. (2009), Inguaggiato et al. (2013), Price et al. (2015), and Kürzinger (2019).

TABLE 2
www.frontiersin.org

Table 2. Concentrations of the major anions and cations (mM) of the hydrothermal fluids from the Aeolian Shallow-water vents.

3.3. Microbial diversity

The bacterial diversity of the shallow-water hydrothermal vents of the Aeolian archipelago was investigated using the V4-V5 variable region of the 16S rRNA gene. After guaranteeing high-quality data (quality check and data filtering steps), a total of 149,351 reads were obtained, and 2,145 ASVs identified. Regarding the alpha diversity of the different sample types, the biofilm samples had lower alpha diversity measure (3) in comparison with the fluid and sediment samples (4.3 and 4.2, respectively). In terms of variability, the fluid samples showed less sample variability (the difference between the lowest diversity value points to the highest) in relation to the biofilm and sediment samples, which could be related to the differences in the matrix of the different sample types (Figure 4A). For instance, deep geothermal fluids impose a stronger bottleneck in terms of possible redox couples and ecological niches to the microbial communities, compared to the sediment samples. Additionally, the temperature was found to be a major driver of diversity between samples (Pearson correlation test, R2 = 0.41), where higher temperature sites showed a decrease in alpha diversity (Figure 4B).

FIGURE 4
www.frontiersin.org

Figure 4. Shannon Alpha diversity indexes. (A) Alpha diversity measures in relation to the different sample types (Biofilm, Fluid, Sediment). (B) Alpha diversity measures in relation to the temperature, with the different Islands sampled, as well the different sample types.

The high-temperature sites, Baia Calcara and Levante Bay, characterized by temperatures above 80°C, were dominated by members belonging to Campylobacterota (formerly Epsilonproteobacteria), Proteobacteria, Aquificota, Bacteroidota, Chloroflexi and Acidobacteriota phylum (Figure 5). Among Campylobacterota, the order Campylobacterales was present in high abundances in all the high-temperature sites. Within Campylobacteriales, the genus Nitratifactor was only present at higher abundances (44%) in the biofilm sample of Baia Calcara, followed by Sulfurimonas (average abundance of 6.08% between the biofilm and sediment samples of Baia Calcara), Campylobacter (average abundance of 9% between the fluid and sediment sample of Levante Bay), and Sulfurovum genera (6.55% in the fluid sample of Levante Bay). The phylum Proteobacteria was composed entirely by class Gammaproteobacteria in the Baia Calcara site (average abundance of 12% between sediments and biofilm samples), and by Gammaproteobacteria and Alphaproteobacteria in the Levante Bay site (average abundance of 15% in both sample types and 8.5% in the fluid samples, respectively). Within Gammaproteobacteria, the orders Chromatiales and Thiomicrospirales were present in the highest abundances (average abundance of 24.6% at Levante Bay between the fluid and sediment samples, and 13.3% between the sediments and biofilm samples of Baia Calcara and the fluid samples of Levante Bay, respectively), with ASVs classified only to the genus Candidatus thiobios and Thiomicrorhabdus. The order Nitrococcales was present in lower abundances (average abundance of 6.1% between both sampling sites and all sample types, with the exception of the sediments of Levante Bay), and composed by the genus Sulfurivirga. Among Alphaproteobacteria, all ASVs were classified into the order Rhodobacterales in both the fluid and sediment samples of Levante Bay (average abundance of 8.5%), however, it was not possible to resolve them down to the genus level. The phylum Aquificota was composed entirely of members belonging to the order Hydrogenothermales, and the genus Hydrogenothermus, and only observed at the sediment samples of Baia Calcara (25.09%). Within Bacteroidota, only the order Flavobacteriales was identified in both the sediment and fluid samples at the Levante Bay site (average abundance of 9.3%), however, no ASVs were classified to the genus level. Acidobacteriota and Chloroflexi, albeit present in lower abundances (5.14% in sediment samples of Baia Calcara and 5.51% in fluid samples of Levante Bay, respectively), had as most abundant, members belonging to the order Thermoanaerobaculares, genus subgroup 23, and to the class Anaerolineae, respectively (Table 3).

FIGURE 5
www.frontiersin.org

Figure 5. Phylum level abundance and distribution of the 16S rRNA diversity in all of the sample sites and sample types at the different islands of the Aeolian Archipelago. For simplicity, only abundances above 5% at each site are illustrated. (A) Biofilm samples; (B) Fluid samples; (C) Sediment samples.

TABLE 3
www.frontiersin.org

Table 3. Taxonomic affiliation of the most abundant ASVs through SILVA database (to the lowest possible taxonomic resolution) and Ezbiocloud closest relative.

Moderate-temperature sites, such as Bottaro, Pietra del Bagno, and Levante Port, characterized by temperatures ranging from 40°C to 53°C, presented a higher diversity of bacterial phyla. These sites were dominated by members belonging to Campylobacterota, Proteobacteria, Bacteroidota, Calditrichota, Chloroflexi, Cyanobacteria, as well as the Archaeal phylum Crenarchaeota and Altiarchaeota at low abundances. Among Campylobacteriota, the order Campylobacterales was only identified at the sediments and biofilm samples of the Levante Port site, with the highest abundant ASVs belonging to the family Arcobacteraceae (37.92% at the sediment samples), without classification to the genus level. Additionally, it was also possible to classify Campylobacterales into the genera Thiovulum and Sulfurimonas in the biofilm samples (23.40 and 19.40%, respectively), and Sulfurovum in the sediment samples (6.51%). Within Proteobacteria, the order Rhodobacterales (21.96% in biofilm sample at Bottaro and 5.70% in biofilm sample of Pietra del Bagno) was the most abundant group, however, it was not possible to classify them down to the genus level. Additionally, the family Mariprofundaceae and the assigned genus Mariprofundus, was the second most abundant genus within Proteobacteria (15.03%) and only present at the sediment samples of Pietra del Bagno. While in lower abundances, it was possible to identify the family Thioprofundae, composed by the genus Thioprofundum (7.47% in the sediment samples of Levante Port), and the family Colwelliaceae, where all the ASVs were assigned to the genus Colwellia (11.54% in biofilm sample of Pietra del Bagno). Moreover, it was also possible to identify the family Thioglobaceae, with the ASVs assigned to the clade SUP05 cluster (8.03% in the biofilm samples of Pietra del Bagno). Within the Bacteroidota, the order Flavobactereales constituted most of the ASVs (26.17% in the biofilm sample of Bottaro); however, classification could not be resolved to the genus level. Lower abundance ASVs belonging to Bacteroidota were also assigned to the order Ignobacteriales and the clade LheB3-7 (11.52% in sediment samples of Levante Port). The phyla Aquificota, similarly to the high-temperature sites, was composed entirely of ASVs belonging to the genus Hydrogenothermus, and only present at the fluid sample of Bottaro (16.57%). Among Calditrichota, the genus Caldithrix represented all the ASVs (17.32% in the sediment sample of Pietra del Bagno and 8.06% in the fluid sample of Bottaro). Additionally, the phototrophic bacterial phyla Cyanobacteria, was composed by the genus Synechococcus of the order Synechococcales in only the fluid samples (6.03 and 5.55% in Levante Port and Bottaro, respectively). The Archaeal phyla were detected in low abundances. Archaeal ASVs were classified into the genera Candidatus_Nitrosopumilus and Candidatus_Altiarchaeum (6.06% in the biofilm sample of Pietra del Bagno and 5.76% in the fluid sample of Levante Port, respectively). Interestingly, some lower abundance ASVs could not be classified up to the phylum level (8.59% in the sediment sample of Pietra del Bagno and 5.67% in the biofilm sample of Bottaro, respectively).

The lowest temperature site, Pollara, at the island of Salina, marked by temperatures of 28°C, was dominated by Cyanobacteria, Proteobacteria, Bacteroidota, and Firmicutes. The majority of the ASVs at lower temperatures were assigned to the phylum Cyanobacteria and the species Synechococcus sp. CC9902 (20.15% in the fluid sample). The order Alteromonadales with the genus Pseudoalteromonas (7.53% in the fluid sample), and the order Steroidobacterales, with the genus Woesia (9.15% in the sediment sample) were the most abundant within Proteobacteria. The phyla Bacteroidetes was composed of the orders Flavobacteriales and Rhodothermales (5.83% in the fluid sample and 5.06% in the sediment sample, respectively).

To understand the role of the environmental parameters and the geochemistry in the structuring of the microbial communities inhabiting the SWHV of the Aeolian islands, we investigated the beta-diversity using both quantitative and qualitative approaches, more specifically, using weighted and unweighted UNIFRAC dissimilarity indexes, respectively, followed by environmental fitting. The unweighted UniFrac nMDS analysis shows a statistically significant separation between the communities of the islands of Lipari and Salina, and Vulcano and Panarea (ADONIS, p = 0.0016, n = 6). Additionally, the same separation was found between the composition of the communities of the biofilm samples, compared to the other sample types analyzed in the present study (ADONIS, p = 0.0075, n = 6). Through linear vector fitting versus the nMDS ordination, we observe which are the main parameters that better explain the microbial diversity of our samples (Figure 6). The group of Vulcano and Panarea can be explained by the temperature, and the group of Salina and Lipari, by the major cations and anions, Mg2+, Na+, and Br. Using the weighted UniFrac, which weights the abundance of the taxa present, we found the same statistically significant separations between the same island groups (ADONIS, p = 0.0075, n = 6), and between biofilms and the other sample types analyzed in the present study (ADONIS, p = 0.0475, n = 6). However, we found differences regarding the main parameters explaining the abundance of microbial communities with vector fitting. For instance, here the biofilm samples of Levante Port and Baia Calcara were mainly explained by temperature and longitude. Additionally, the microbial diversity of Bottaro and Levante Port, in contrast with the unweighted version, are also explained by the major anions and cations, with the addition of Cl.

FIGURE 6
www.frontiersin.org

Figure 6. Non-metric multidimensional scaling (nMDS, stress 0.075) plot of the 16S gene amplicon microbial diversity based both the weighted and unweighted UniFrac measures, overlaid with environmental vector fitting, illustrating vectors with significant p values (p < 0.01).

4. Discussion

Shallow-water hydrothermal vents are marine environments associated with volcanic activity and tectonically active margins. Due to their unique characteristics, such as the presence of light, proximity to landmasses, and dynamic geochemistry, they harbor a highly diverse microbial consortia (Tarasov et al., 2005; Price and Giovannelli, 2017). The Aeolian archipelago is one of the most studied shallow-water hydrothermal systems in the world (Gugliandolo and Italiano, 1999; Simmons and Norris, 2002; Rogers and Amend, 2005; Rusch et al., 2005; Lentini et al., 2007, 2014; Manini et al., 2008; Maugeri et al., 2009, 2010, 2013; Gugliandolo et al., 2015; Gugliandolo and Maugeri, 2019). However, the role of the diverse geochemical environments found in the different islands in controlling the microbial community structure, is still poorly understood, and usually investigated in a handful of sites. To address this, we surveyed 6 different shallow water vents located in four different islands, combining 16S rRNA amplicon sequencing with geochemical and physico-chemical data. To the best of our knowledge, two of the investigated shallow vents (Pietra del Bagno and Pollara) were never reported before in the literature.

4.1. Geochemistry of the Aeolian shallow-water hydrothermal vents

Hydrothermal systems are characterized by the convective circulation of water that percolates into the crust, gets heated, and then is discharged at the sea floor. During this journey, the fluids are chemically and physically transformed by subsurface processes, such as temperature-dependent water-rock reactions, leaching of transition metals and volatiles from the host lithologies, and phase-separation events (Fournier, 2007; German and Seyfried, 2014). These processes, although occurring in a dynamic continuum, are generally separated into three distinct areas: the recharge zone, the reaction zone, and the upflow zone. Due to the inherent characteristics of the locations within the subsurface where they take place, these areas are marked by the occurrence of different processes. The recharge zone is characterized by the fixation of alkali metals, boron, and the fixation of Mg2+ associated with the mineralization of smectite and chlorite. Deep in the system (reaction zone), anorthite is transformed to albite (albitization), resulting in the exchange of Na+ and Si4+ for Ca+ into the fluids, as well as metal mobilization. Additionally, the high temperatures and pressures present at depth allow phase separation events, separating the low weight ions (migrate to the vapor phase) from the heavier weight ions (remain in the fluids). Cl, by remaining in the fluid phase, modifies the salinity content of the hydrothermal fluids. In the upflow zone, further low-temperature water-rock reactions can occur, such as the precipitation of anhydrite and other minerals (Berndt et al., 1988; Tivey, 2007; Alt, 2013; German and Seyfried, 2014). During ascent, the fluids can either be channeled, feeding directly from deep subsurface reservoirs, or diffuse, with significant mixing of end member fluids with seawater prior to discharge. All together these processes entail different ratios of mixing with seawater, contributing to additional alterations to the fluids before they are discharged at the seafloor where they are collected (Price and Giovannelli, 2017). By analyzing the enrichment and depletion of major ions and cations in vent fluids, we can have a deeper understanding of the main processes occurring in the subsurface. These geochemical processes define the geochemical space available for microbial communities to thrive in these ecosystems, and are fundamental in investigating the microbial diversity of shallow-water hydrothermal vents (Price et al., 2015; Price and Giovannelli, 2017).

In the present study, we leverage the information available in the already published literature on the geochemistry of hydrothermal fluids of the shallow vents of Vulcano, Panarea, and Stromboli (Figure 3) combined with new data. Furthermore, we add for the first time geochemical data for the shallow vents of Lipari and Salina (Pietra del Bagno and Pollara, respectively). It is important to note that subsurface processes at shallow-water hydrothermal vents hosted in volcanic arc settings have not been characterized to the same level of detail as Mid-Ocean Ridge hosted deep-sea systems. For this reason, our study relies on previous characterization of Panarea fluids by Price et al. (2015) and conceptually extended to the other shallow vents investigated. By combining published and new geochemical data we identified three mixing trends (Type 1, Type 2, and Type 3) already described for Panarea by Price et al. (2015).

The observed different trends reveal that these hydrothermal fluids evolved differently, and therefore were subjected to different subsurface processes within the Panarea hydrothermal system. Data from Vulcano, Lipari, and Salina shallow vents shows that fluids from diverse islands fall within the Type 3 trend (Figure 3). These fluids do not show enrichments of major anions and cations, but instead, show depletions of the major species. The evolution of Type 3 fluids results in acidic fluids, with concentrations of major ions similar to sea-water concentrations, whereas type 1 and type 2 fluids are the result of deep recirculating brine and focused upflow of shallower brines, respectively. This suggests that the difference between these fluid types is directly related to subsurface processes, such as phase-separation deep in the system, the scrubbing of exsoluted gasses that percolate up the crust, as well as the mixing of the diffuse flow fluids with sea-water (Price et al., 2015). Based on our data, it appears that the evolution of Type 3 fluids might be a common process in the shallow water vents of the Aeolian Archipelago. Additionally, and as evident in Price et al. (2015), most of the fluids, including the ones in the present study, do not extrapolate to zero in both SO42− and Mg2+. This can be related to seawater mixing as the fluids approach the surface, due to their diffuse nature. However, in the Bottaro site, the fluids are enriched in Mg2+. It was hypothesized that Mg2+ enrichments in Panarea could be related to the leaching of host rocks due to acidification, by means of argillic alteration (Reeves et al., 2011). Since volcanic arc hydrothermal systems are not exposed to the same pressures compared to their deep-sea counterparts (Price and Giovannelli, 2017), the presence of free gas phases may also be contributing to the depletion of Na+ and Cl observed relative to sea water (Figure 3), altering the salinities of the fluids as they percolate through the hydrothermal system. All the fluid samples collected in our study had Br/Cl ratios lower than seawater, suggesting the presence of dissolved halite in solution. Albitization, which occurs deep in the system and results in the release of Ca+ in exchange for Na+ (Alt, 2013), does not seem to be a major process affecting the evolution of our fluid samples, compared to other fluids recovered from Panarea and Vulcano (Figure 3).

While the geochemical composition of the fluids is an important parameter in controlling the diversity of microbial communities, the availability of organic carbon for heterotrophic consumption is fundamental in driving community composition. In hydrothermal systems organic carbon can either be delivered through the subsurface fluids as dissolved organic carbon produced or recycled at depth (Fullerton et al., 2021), or from the water column, either as sinking organic matter particles, lateral advections of surface primary productivity, and terrigenous run offs (Gomez-Saez et al., 2015; Price and Giovannelli, 2017; LaRowe et al., 2020). Additionally, it could be directly produced in situ by a diverse community of phototrophic and chemolithotrophic microorganisms (Price and Giovannelli, 2017). The organic carbon delivered to the system is accumulated in the sediments, and analyzing the sedimentary organic matter content can provide an estimate of the organic load on the vent system. In our data the sedimentary organic matter, measured as biopolymeric organic carbon (BPC) was highly variable between the different sites, and did not show linear relationships with temperature, in contrast to what has been observed in a previous study on a shallow-water hydrothermal vent system (Giovannelli et al., 2013; Figure 2A). The lowest amounts of BPC was measured in the site with the lowest temperatures (Pollara, 28°C and Pietra del Bagno, 45°C), and could be related to a to hydrodynamic processes, perturbing the system in these specific locations (Bao et al., 2018), as well as to the fact that these sites are further away from the coast, and thus representing lower amounts of terrigenous inputs to the system. The high concentrations found in Bottaro are instead related to the thick microbial biofilms attached to the sand and gravel. The concentration of chlorophyll-a (chl-a) pigments was lower at the higher temperatures sites (>80°C; Figure 2B) in line with previous observations (Giovannelli et al., 2013). The presence of Chl-a in the sediments can be either related to input of fresh organic carbon from the water column or to the presence of phototrophs in the sediments and biofilms of the vent. The presence of known phototrophs in the sediment and biofilm samples (see below) suggests that the high chlorophyll concentrations found in the same sites are the results of in situ production, possibly taking advantage of the nutrients and CO2 released by the shallow vents. Additionally, the chl-a values were higher at the Levante Port and Bottaro sites, which had similar depths (average 6.4 meters), suggesting that beside temperature, the depth of the shallow vent could play a role in controlling the presence and abundance of photosynthetic organisms.

4.2. Microbial diversity

Previous studies conducted at shallow-water hydrothermal vents have demonstrated that a highly diverse microbial consortia inhabit these ecosystems (Hoaki et al., 1995; Sievert et al., 2000; Manini et al., 2008; Maugeri et al., 2009, 2010; Zhang et al., 2012; Price et al., 2013; Kerfahi et al., 2014; Lentini et al., 2014; Dávila-Ramos et al., 2015; Gugliandolo et al., 2015; Gomez-Saez et al., 2017). In the present study, we found that the microbial communities of the Aeolian Islands shallow-water vents are dominated by members belonging to diverse phyla, such as Proteobacteria, Campylobacteriota, Bacteroidota, Cyanobacteria, Chloroflexi, among others present at lower abundances. The major lineages found here have been reported in other shallow-water vents (Maugeri et al., 2010; Zhang et al., 2012; Kerfahi et al., 2014; Lentini et al., 2014; Dávila-Ramos et al., 2015; Gugliandolo et al., 2015; Gomez-Saez et al., 2017; Price and Giovannelli, 2017), illustrating their important ecological role in the functioning of the shallow-water vent ecosystem.

Our data shows that temperature is a major driver explaining the microbial diversity (alpha-diversity, Pearson correlation test, R2 = 0.41, p < 0.05). High-temperature sites (>80°C) were characterized by lower diversity in comparison to lower and moderate temperatures sites (28°C to 50°C; Figure 4B). Temperature has been shown to be one of the major factors controlling microbial diversity in a number of diverse geothermal ecosystems, including shallow-water vent environments (Giovannelli et al., 2013; Gugliandolo et al., 2015), deep-sea hydrothermal vents (Dick et al., 2013; Ding et al., 2017), and hot springs (Chan et al., 2017; Fullerton et al., 2021).

High-temperature locations, such as Baia Calcara and Levante Bay sites in the islands of Vulcano and Panarea, respectively, were characterized by temperatures above 80°C. The microbial communities inhabiting these high-temperature shallow vents were dominated by groups related to known chemolithoautotrophic, thermophilic bacteria, harnessing energy through the oxidation of reduced sulfur species and hydrogen. These types of metabolic strategies have been widely reported in shallow-water, deep-sea hydrothermal vents, as well as hot-springs (Tarasov et al., 2005; Yamamoto and Takai, 2011; Inskeep et al., 2013; Price and Giovannelli, 2017). The Baia Calcara and Levante Bay sites had as most abundant, ASVs related to the genera Nitratifactor, Thiomicroarhdbus, Hydrogenothermus, Candidatus_thiobios, and the family Rhodobactereceae. The biofilm community at the Baia Calcara site was mostly composed of members belonging to Nitratifactor and Thiomicroarhdbus. Nitratifactor has been previously described in shallow and deep-sea hydrothermal systems (Zhang et al., 2012; Zeng et al., 2021). The closest sequence to our ASVs was recovered from a shallow-water hydrothermal system in Kueishan Island, Taiwan (unpublished data). Interestingly, this group was not present in the sediment samples of the same site. Additionally, the genus Thiomicroarhdbus was also present in high abundances in the biofilm samples, and in lower abundance also in the sediment sample. It was possible to assign this ASVs to the species Galenea microaerophila (99.73% similarity) (Giovannelli et al., 2012). Initially described from a shallow-water vent from Paleochori Bay, on the island of Milos, Galenea microaerophila grows chemolithoautotrophically, through the oxidation of thiosulfate as electron donor, and oxygen (5%) as electron acceptor (Giovannelli et al., 2012). The sediment sample at Baia Calcara was mostly composed by the genus Hydrogenothermus, of the phylum Aquificota. This ancestral and deep-branching bacterial phylum is highly prevalent in geothermal environments (Hedlund et al., 2015; Giovannelli et al., 2017), and it is capable of fixing CO2 into biomass using energy derived from H2S and H2 with either elemental sulfur, nitrate or oxygen (at low partial pressures) as electron acceptor. The dominant ASV in our dataset was related to the species Hydrogenothermus marinus (99.2% similarity), a hydrogen-oxidizing, chemolithoautotrophic, thermophilic bacterium (Stöhr et al., 2001) originally isolated from sediment samples from Vulcano island.

At the fluid and sediment samples at Levante Bay, the other high temperature site, the family Rhodobactereceae and the genus Candidatus thiobios were present in highest abundances. The family Rhodobactereceae is widespread in marine environments, and are capable of growing through the utilization of multiple organic, as well as inorganic compounds, contributing to their vast ecological diversity (Pohlner et al., 2019). In our study, it could also be identified in moderate to lower temperature sites. The ASVs classified as Rhodobactereceae were further related to the species Actibacterium pelagium (100% similarity), which have been isolated from the seawater of the Mediterranean Sea, Pacific and Atlantic Oceans, growing heterotrophically on a variety of organic carbon sources (Guo et al., 2017). Even though members belonging to Rhodobactereceae have been previously identified in other hydrothermal systems (Forget et al., 2010; Zhang et al., 2012; Pohlner et al., 2019), to our knowledge, this is the first time this species has been reported in a hydrothermal vent environment. The genus Candidatus thiobios is associated with ectosymbiotic bacteria growing on the marie ciliate Zoothamnium niveum, commonly found in sulfide-rich habits, including hydrothermal vents (Rinke et al., 2006; Bright et al., 2014). The presence of high abundance of this genus could be related to the presence of the Zoothamnium niveum inhabiting this specific site.

Both sites (Baia Calcara and Levante Bay), included lower abundances of sequences assigned to the genus Sulfurimonas, Sulfurivirga, Sulforovum, and Campylobacter, all belonging to the Campylobacterota phylum (formerly known as Epsilonproteobacteria). Members of the Sulfurimonas genus are widespread in sulfidic environments and are capable of using multiple reduced sulfur species as electron donors, such as thiosulfate, hydrogen sulfide, and elemental sulfur (Han and Perner, 2015). The genus Sulfurivirga could be assigned to the species Sulfurivirga caldiculari (99.73% similarity), a thiosulfate-oxidizing, thermophilic, chemolithoautotrophic bacterium, described in a shallow-water vent system in Japan (Takai et al., 2006). Similarly, Sulfurovum is also associated with thiosulfate oxidation and chemolithoautotrophy (Inagaki et al., 2004). Campylobacter, within Campylobacteriota, is ubiquitous in hydrothermal environments, including shallow-water hydrothermal vents (Rinke et al., 2006; Anderson et al., 2011; Giovannelli and Vetriani, 2017; Wang et al., 2017). These genera have been associated with high-temperature high sulfidic locations in other shallow-water, as well as deep-sea hydrothermal systems (Forget et al., 2010; Zhang et al., 2012; Giovannelli et al., 2013; Gugliandolo et al., 2015; Gomez-Saez et al., 2017; Price and Giovannelli, 2017). Members of the Campylobacterota are present in most sites presenting higher temperatures (Levante Bay, Levante Port, Bottaro, Baia Calcara and Pietra del Bagno) and dominate the biofilm communities of Levante Port and Baia Calcara. Hydrothermal vents are often enriched in reduced sulfur species (German and Seyfried, 2014; Price and Giovannelli, 2017). The presence of sulfur-related metabolisms, coupled to inorganic carbon fixation, highlights the importance of these types of metabolic strategies for the function of hydrothermal systems, as well as the importance of these ecosystems to mediate the biogeochemical cycles of sulfur and carbon. Additionally, these high temperature locations had the highest protein/carbohydrate ratios (15.2 in Baia Calcara and 12.6 in Levante Port), suggesting high in situ productivity. Since at these locations, recovered ASVs are related to known chemolithoautotrophic groups, indicates that they are major contributors to primary productivity in high temperature shallow water hydrothermal vents. Similarly, Gomez-Saez et al. (2017), reported that chemoautotrophy constitutes an important contribution to fresh organic matter production, accounting for up to 65% of autotrophic carbon fixation in the shallow-water hydrothermal vents of Dominica.

At moderate to low temperatures an increase in phototrophic microbial lineages was observed, also reflected in the quantities of BPC and Chlorophyll-a. The prevalent photosynthetic group observed in this study was assigned to the genus Synechococcus and to the species Synechococcus sp. CC9902 (99.47% similarity), belonging to the phylum Cyanobacteria. The presence of phototrophs in shallow-water hydrothermal systems is well documented, as well as a key feature of these systems (Tarasov et al., 1990, 1999, 2005; Sievert et al., 2000; Zhang et al., 2012; Maugeri et al., 2013; Price et al., 2013; Lentini et al., 2014; Dávila-Ramos et al., 2015; Gugliandolo et al., 2015; Gomez-Saez et al., 2017; Price and Giovannelli, 2017). Here, Synechococcus sp. CC9902 was present in lower abundances at moderate temperatures in Bottaro and Levante Port, and high abundances in Pollara, the lowest temperature site. Synechococcus sp. CC9902 is widespread in the oceans and is described as having increased plasticity to diverse saline concentrations and light intensities (Ludwig and Bryant, 2012; Kim et al., 2018). Furthermore, the presence of this taxon has already been described in shallow-water hydrothermal vent environments (Maugeri et al., 2013). However, in other shallow-water vents systems, the specific lineages within Cyanobacteria differ. This suggests that in shallow-water hydrothermal vents the photoautotrophic niche space can be occupied by diverse lineages of phototrophic taxa. Notwithstanding, it is important to note that, even though the presence of Cyanobacteria in shallow water hydrothermal vents is well documented, the possibility of mixing with background seawater cannot be discarded. For instance, it has been shown that fluid mixing with seawater could be the result of the convective cells in porous media (Yücel et al., 2013).

The shallow water vent Pietra del Bagno, Lipari, due to the presence of orange precipitates surrounding the venting site, was deemed to be iron-rich. The sediment community was mainly composed by the family Anaerolineares of the phylum Chloroflexi, and the genus Mariprofundus belonging to Proteobacteria. Members belonging to Chloroflexi are ecologically and physiologically diverse, allowing them to colonize a wide variety of environments (Ward et al., 2018; Vuillemin et al., 2020). Despite being originally associated with anoxygenic phototrophy, they are also known to employ other metabolic strategies (McIlroy et al., 2017). The family Anaerolineares is associated with hydrocarbon oxidation and is known to have important roles in the functioning of anaerobic digestion systems (McIlroy et al., 2017). Nonetheless, they have been reported in deep-sea and shallow-water hydrothermal vents (Rajasabapathy et al., 2018; McGonigle et al., 2020; Wang et al., 2020). The closest relative to our sequences was reported from Fe-rich hydrothermal sediments of two submarine volcanoes of the Tonga volcanic arc (Forget et al., 2010). Interestingly, this bacterial family has also been observed to be present in high abundances in Fe-enriched locations (Hoshino et al., 2016; Ward et al., 2017). Other members of Chloroflexi, such as Ardenticatena maritima, grow by dissimilatory iron-reduction (Kawaichi et al., 2013). Thus, it is possible that Anaerolineares can couple the oxidation of organic matter to iron-reduction, making them major players in the biomineralization of organic matter in shallow-water vent environments.

The genus Mariprofundus of Zetaproteobacteria are well-known iron(II)-oxidizers that couple the iron (II) oxidation to the reduction of oxygen or nitrate (Chiu et al., 2017; McAllister et al., 2019). Members belonging to the phyla Zetaproteobacteria are widespread in hydrothermal Fe-rich environments (Emerson and Moyer, 2002; Emerson et al., 2010; Forget et al., 2010; Gomez-Saez et al., 2017; McAllister et al., 2019). Within Mariprofundus, it was possible to assign to the species Mariprofundus aestuarium (99.73% similarity), previously isolated from redox-stratified water columns, and capable of growing through the oxidation of iron (II) (Chiu et al., 2017). Additionally, albeit in lower abundances, it was possible to identify the genus Caldithrix. More specifically, we could assign our sequences to the species Caldithrix abyssi (98.94% similarity). This species is a known chemolithoautotroph, capable of growing anaerobically through the reduction of nitrate, using hydrogen as an electron donor (Miroshnichenko et al., 2003). Our geochemical analysis did not identify nitrate in the fluids at this site (data not shown); however, it was recently found that an environmentally assembled genome (MAG) of this species codes for a homolog of an extracellular respiratory system found in Geobacter, associated with iron reduction (Fortney et al., 2018). Given that this terminal electron acceptor was not tested in the original isolate (Miroshnichenko et al., 2003), the coupling of hydrogen oxidation to iron reduction could potentially suggest a metabolic strategy employed by Caldithrix. Notably, previous studies conducted at Fe-rich shallow-water vents have reported similar communities to the ones observed at this site (Handley et al., 2010; Meyer-Dombard et al., 2013; Gomez-Saez et al., 2017). Thus, these lineages can be considered to be important in driving the iron cycling in the shallow-water vent ecosystem, making important contributions to the biogeochemical cycles of iron. Interestingly, it was not possible to identify iron oxidizers in the biofilm community in high abundance at the same site. Conversely, the genus Collwelia was present in high abundances, followed by the clade SUP05. Members belonging to Collwelia are psychrophilic, heterotrophic bacteria, ubiquitous in cold marine environments (Techtmann et al., 2016), including sea ice (Zhang et al., 2008), polar sediments (Wang et al., 2015), and deep-sea trenches (Nogi et al., 2004). Specifically, we identified the species Colwellia rossensis (100% similarity), isolated from coastal Antarctic sea ice (Bowman et al., 1998). Colwellia species have been isolated from marine environments, and are characteristically psychrophiles (Bowman et al., 1998). The presence of ASVs related to this species at higher temperatures suggests a bigger temperature gradient to which this group is adapted to. The clade SUP05 is composed of known sulfur-oxidizing bacteria, also described as having a heterotrophic-based metabolism (Spietz et al., 2019; Morris and Spietz, 2022). At this location, the protein to carbohydrate ratios is 1.13. As reported by Pusceddu et al. (2003), since proteins are readily mobilized, as compared to carbohydrates, low ratios suggest the presence of aged organic matter (not freshly produced), and probably indicates external input as a potential source of organic matter to the biofilm communities.

Additionally, in moderate temperatures sites, it was possible to identify in higher abundance the families Flavobactereceae and Arcobactereceae, in the biofilms of Bottaro and Levante Port, and the genus IheB3-7 at the sediment sample of Levante Port. The family Flavobacteriaceae of the Phylum Bacteroidetes, ubiquitous in marine environments, are known to degrade complex hydrocarbons, with large numbers of glycosyl hydrolases and peptidases encoded in the genomes of this bacterial group (Gavriilidou et al., 2020). Moreover, the family Arcobactereceae, also widely distributed in marine environments (Venâncio et al., 2022), has been reported in shallow vent hydrothermal environments (Maugeri et al., 2009; Gugliandolo et al., 2015; Gomez-Saez et al., 2017).

At the lowest temperature site (Pollara), Cyanobacteria was the most abundant group, and with lower abundances, the genera Pseudoalteromonas, Woesia, NS5 marine group, and Anoxybacillus. These microbial groups are associated with heterotrophic metabolisms, and given the protein/carbohydrate ratios at this location of 3.64, as well as Chla-a concentration of 1.0 μg g−1, suggests that these groups rely on the export of organic carbon from the water column into the vent system, coupled with the usage of aged organic matter derived from external input, and in situ production from phototrophic organisms. Pseudoalteromonas, Woesia, and NS5 marine group are widely distributed in marine environments (Holmström and Kjelleberg, 1999; Du et al., 2016; Priest et al., 2022). Additionally, Anoxybacillus was also reported in deep-sea hydrothermal vent sediments (Cheng et al., 2021). Here, it was possible to assign Anoxybacillus to the species Anoxybacillus rapiensis (100% similarity), a thermophilic, anaerobic, heterotrophic bacterium, isolated from diverse hot springs in Bulgaria (Derekova et al., 2007).

While temperature is one of the key variables controlling the microbial diversity in our dataset, other environmental parameters contribute significantly in shaping microbial communities of the different sample types collected (fluids, sediments and biofilms) across the geochemically diverse vents sampled. This has been recently demonstrated in hot springs (Fullerton et al., 2021), where a combination of the geological settings, the aqueous geochemistry and trace element availability provided strong explanatory power of the beta diversity turnover among the sampled springs. In this work unweighted and weighted UniFrac based multivariate analysis (Figure 6) revealed the presence of three statistically significant groups. The first was composed of the biofilm samples, where no significant correlations with the measured parameters were found to significantly explain the diversity. The biofilm composition and architecture are found to be influenced by both deterministic, as well as stochastic effects. For instance, previous studies have shown that biofilm thickness drives the biofilm community and diversity (Suarez et al., 2019), combined with the age of the biofilm and other environmental variables. The second major group, encompassing the samples of fluids and sediments from the islands Vulcano and Panarea, could be explained by the temperature. As discussed above, these locations were marked by moderate to high temperatures, with ASVs associated with microbial lineages related to the oxidation of hydrogen or reduced sulfur species. In hydrothermal systems fluid ascent is facilitated by the permeability of lithospheric fault systems, allowing a more focused channeling of the fluids to the surface (Pearce et al., 2019). Moreover, hydrothermal systems further away from fault systems are likely to show weak thermal anomalies (Taillefer et al., 2018). Consistent with this, the site with the lowest temperature, Pollara in Salina, is the sampling site further away from known fault systems in the area. Interestingly, we found that temperature and Mg2+ follow opposite correlation patterns. Magnesium is depleted in end-member hydrothermal fluids, resulting from water-rock reactions during the circulation of the fluids in the hydrothermal system (Alt, 2013). The presence of higher quantities of Mg2+ in low temperature sites, can entail less degree of water-reactions, which could be a function of the weak thermal anomalies found in the island of Salina and Lipari. Consistent with this, these islands present less geothermal activity compared to the other island samples in this study.

In the third group, encompassing the sites from the islands of Salina and Lipari, the microbial diversity was explained by the major anions and cations Mg2+, Br, and Na+. Despite the geochemical analysis encompassing our samples in combination with the literature available placed our samples within Type 3 fluids, small variations of the major anions and cations show statistically significant correlation with the microbial distribution in both these islands. Diverse subsurface processes could result in similar concentrations of the major anions and cations, but result in a differential enrichment of minor elements, such as trace elements (Rollinson and Pease, 2021). Recent studies have shown that trace elements play key roles in shaping microbial niche occupation, as well as have been used to improve the cultivation of extremophilic organisms from natural hot springs (Meyer-Dombard et al., 2012; Shafiee et al., 2021).

Similarly to the unweighted UniFrac, with the weighted version, which weights the abundance of the taxa present, we found the same statistically significant separation between the communities of Vulcano/Panarea and Salina/Lipari, as well as the biofilm samples and the other samples types. However, conversely to the unweighted UniFrac, we found the temperature to be correlated with the abundance of taxa of Baia Calcara and Levante Port biofilm communities. These were dominated by thermophilic communities belonging to Nitratifractor and Thiomicrorhabdus. Therefore, both these parameters can act as strong selective pressures, driving not the structuring of the biofilm communities of Baia Calcara and Levante Port, but selecting for thermophilic, chemolithoautotrophic groups. The correlation with longitude could be related with the location of the islands within the Aeolian volcanic arc. It has been previously demonstrated that the magmatic composition of the Aeolian arc exhibits within-islands variations (Peccerillo et al., 2013). These differences could entail differential inputs of magmatic volatiles into the hydrothermal fluids, which in turn could be influencing the microbial communities inhabiting the shallow-water hydrothermal vents. It was also possible to find the major elements correlating with the abundance of the communities of the islands of Lipara and Salina, and with weighted UniFrac, the communities of Vulcano and Salina.

In the Aeolian Archipelago it is possible to find a highly diverse set of marine habitats. These habitats, in turn, harbor an abundant number of species, some of which with endangered status (Álvarez et al., 2019). For this reason, over the years, there have been increasing efforts to enhance the conservation measures on the Aeolian Islands, through the practices of sustainable fishing, tourism, and community education (Álvarez et al., 2019). The presence of high hydrothermal activity should further augment the need for conservation practices in these islands, for different reasons. One is that the microbial communities inhabiting these vents play important ecological roles for the entire shallow marine ecosystem (Tarasov et al., 2005; Price and Giovannelli, 2017), potentially sustaining primary productivity in surrounding ecosystems (Price and Giovannelli, 2017). Given their key roles as primary producers, these microbial communities sustain complex food webs on the surrounding islands. Secondly, extreme environments are often looked upon for the isolation of microorganisms that produce industrially and medically relevant metabolites, such as novel enzymes and antibiotics (Littlechild, 2015; Sayed et al., 2020; Petrillo et al., 2021). For instance, novel enzyme-producing Bacillus have been isolated from the Island of Panarea, whereas novel hydrolases have been identified on environmental DNA libraries from the sediments of the shallow-vents of the island of Vulcano (Lentini et al., 2007; Gugliandolo et al., 2012; Placido et al., 2015). Thirdly, in recent years there have been efforts to apply conservation status to deep-sea hydrothermal vents (Van Dover, 2012; Van Dover et al., 2018). We argue that shallow-water hydrothermal vents play an equally important role in the functioning of marine ecosystems, and should be explicitly acknowledged in marine spatial planning and conservation efforts, not only because of their uniqueness, but because they represent an accessible window for us to appreciate the intricate relationships between the geosphere and biosphere.

5. Conclusion

In conclusion, our study shows that the Aeolian archipelago shallow water hydrothermal vents harbors a highly diverse microbial consortia. In line with previous studies conducted in these unique ecosystems, we found the presence of ASVs related to organisms employing diverse metabolic strategies, ranging from chemolithoautotrophy, heterotrophy and phototrophy. Furthermore, we found those same groups to be correlated strongly with temperature, where we could observe a transition from higher temperatures (mainly ASVs related to chemolithoautotrophic groups) to moderate and lower temperatures (mainly ASVs related to phototrophic and heterotrophic groups). This transition is commonly observed in geothermal environments, and can be considered as a definite characteristic of these systems. Furthermore, through the analysis of BPC, we found that different types of organic matter are available to the resident communities, from terrigenous input exported to the system, as well as in situ production by chemolithoautotrophic and phototrophic groups. This supports the literature on how shallow water hydrothermal vents are considered high energy environments. Through beta diversity analysis, we found the presence of three significant groups, composed of the biofilm samples, and the sediment and fluid samples of the islands of Vulcano and Panarea, and Lipari and Salina. The group of Vulcano and Panarea was correlated with temperature, as these locations were marked by moderate to high temperatures, and they were dominated by ASVs related to chemolithoautotrophic groups. As for the group of Lipari and Salina, the distribution of the microbial communities was correlated with major ions. This suggests that subsurface water rock reactions have an influence on the concentration of the major ions, which in turn, can influence the microbial communities inhabiting these shallow water hydrothermal vents. Additionally, we found that the longitude, together with the temperature, was a significant vector in the weighted UNIFRAC. This correlation could be due to the location of the islands within the Aeolian volcanic arc and the different geothermal activity within each island.

Additionally, our compilation of fluid data from previous studies showed the same trend noted by Price et al. (2015), wherein in these hydrothermal systems, one can find 3 types of fluid evolution, showing differential enrichments and depletions of the major anions and cations, and therefore hinting at different subsurface processes taking place in the different islands. Given that the Aeolian archipelago shallow water hydrothermal vents, and the microbial communities inhabiting them, have important ecological roles for the surrounding marine environment, we think it is important their inclusion into present and future conservation efforts. Finally, shallow-water vents and the volcanic arc lithologies at which they occur are still understudied, compared to the basaltic, mid-ocean ridge lithologies from their deep-sea counter-parts. Our study is the first stage of a deeper characterization of the Aeolian archipelago shallow-water systems, and future studies will shed light on the ecological role of the communities in these systems, as well as on the local biogeochemical cycles.

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.ebi.ac.uk/ena, PRJEB54611, https://zenodo.org/record/7480614.

Author contributions

BB carried out data analysis, geochemical experiments and wrote the first draft of the manuscript. BB, AF, MS, MG, AB, MCO, and LD produced the geochemical and microbiological data. MS and AF carried out Dada2 and Phyloseq and statistical analysis. AF, AC, MC, GB, RP, MI, CV, and DG planned the expedition and collected the samples. RP, CV, AC, and DG supervised the study. All authors contributed equally to the final version of the manuscript.

Funding

This work is partially funded by the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation program (grant agreement No. 948972) to DG. Partial funding was from the NASA Habitable Worlds program under grant 80NSSC20K0228 to REP.

Acknowledgments

The authors wish to thank the Blue Marine Foundation (https://www.bluemarinefoundation.com/) for support in organizing the sampling expedition and for the use of the zodiac. They also thank the Salina diving center for providing the diving logistics necessary to carry out the sampling and Luca Magni for assistance with the dives and location access.

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.

Footnotes

References

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

Alt, J. C. (2013). “Subseafloor processes in Mid-Ocean ridge Hydrothennal systems” in Geophysical monograph series. eds. S. E. Humphris, R. A. Zierenberg, L. S. Mullineaux, and R. E. Thomson (Washington, DC: American Geophysical Union), 85–114.

Google Scholar

Álvarez, H., Perry, A. L., García, S., Blanco, J., and Aguilar, R. (2019). Towards the creation of a marine protected area in the Aeolian Islands. MarXiv 2019:b9dqc. doi: 10.31230/osf.io/b9dqc

CrossRef Full Text | Google Scholar

Amend, J. P., Rogers, K. L., Shock, E. L., Gurrieri, S., and Inguaggiato, S. (2003). Energetics of chemolithoautotrophy in the hydrothermal system of Vulcano Island, southern Italy. Geobiology 1, 37–58. doi: 10.1046/j.1472-4669.2003.00006.x

CrossRef Full Text | Google Scholar

Anderson, I., Sikorski, J., Zeytun, A., Nolan, M., Lapidus, A., Lucas, S., et al. (2011). Complete genome sequence of Nitratifractor salsuginis type strain (E9I37-1T). Stand. Genomic Sci. 4, 322–330. doi: 10.4056/sigs.1844518

PubMed Abstract | CrossRef Full Text | Google Scholar

Antranikian, G., Suleiman, M., Schäfers, C., Adams, M. W. W., Bartolucci, S., Blamey, J. M., et al. (2017). Diversity of bacteria and archaea from two shallow marine hydrothermal vents from Vulcano Island. Extremophiles 21, 733–742. doi: 10.1007/s00792-017-0938-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Bao, R., van der Voort, T. S., Zhao, M., Guo, X., Montluçon, D. B., McIntyre, C., et al. (2018). Influence of hydrodynamic processes on the fate of sedimentary organic matter on continental margins. Glob. Biogeochem. Cycles 32, 1420–1432. doi: 10.1029/2018GB005921

PubMed Abstract | CrossRef Full Text | Google Scholar

Bastianoni, A., Cascone, M., De Moor, J. M., Barry, P. H., Price, R., Cordone, A., et al. (2022). The missing carbon budget puzzle piece: Shallow-water hydrothermal vents contribution to global CO2 fluxes.

Google Scholar

Berndt, M. E., Seewald, J. S., and Seyfried, W. E. (1988). Hydrothermal alteration processes at Mid Ocean ridges:Constrains from diabase alteration experiments, hot-spring fluids and composition of the oceanic crust. Can. Mineral. 26, 787–804.

Google Scholar

Bligh, E. G., and Dyer, W. J. (1959). A rapid method of Total lipid extraction and purification. Can. J. Biochem. Physiol. 37, 911–917. doi: 10.1139/o59-099

PubMed Abstract | CrossRef Full Text | Google Scholar

Bowman, J. P., Lewis, T. E., Staley, J. T., and McMeekin, T. A. (1998). Colwellia demingiae sp. nov., Colwellia hornerae sp. nov., Colwellia rossensis sp. nov. and Colwellia psychrotropica sp. nov.: psychrophilic Antarctic species with the ability to synthesize docosahexaenoic acid (22: 63). Int. J. Syst. Bacteriol. 48, 1171–1180. doi: 10.1099/00207713-48-4-1171

CrossRef Full Text | Google Scholar

Bright, M., Espada-Hinojosa, S., Lagkouvardos, I., and Volland, J.-M. (2014). The giant ciliate Zoothamnium niveum and its thiotrophic epibiont Candidatus Thiobios zoothamnicoli: a model system to study interspecies cooperation. Front. Microbiol. 5:145. doi: 10.3389/fmicb.2014.00145

PubMed Abstract | CrossRef Full Text | Google Scholar

Brovarone, A. V., Butch, C. J., Ciappa, A., Cleaves, H. J., Elmaleh, A., Faccenda, M., et al. (2020). Let there be water: how hydration/dehydration reactions accompany key earth and life processes#. Am. Mineral. 105, 1152–1160. doi: 10.2138/am-2020-7380

CrossRef Full Text | Google Scholar

Callahan, B. J., McMurdie, P. J., and Holmes, S. P. (2017). Exact sequence variants should replace operational taxonomic units in marker-gene data analysis. ISME J. 11, 2639–2643. doi: 10.1038/ismej.2017.119

PubMed Abstract | CrossRef Full Text | Google Scholar

Cardellini, C., Chiodini, G., and Frondini, F. (2003). Application of stochastic simulation to CO2 flux from soil: mapping and quantification of gas release: simulation of soil CO2 flux. J. Geophys. Res. Solid Earth 108:2165. doi: 10.1029/2002JB002165

CrossRef Full Text | Google Scholar

Chan, C. S., Chan, K.-G., Ee, R., Hong, K.-W., Urbieta, M. S., Donati, E. R., et al. (2017). Effects of physiochemical factors on prokaryotic biodiversity in Malaysian Circumneutral Hot Springs. Front. Microbiol. 8:1252. doi: 10.3389/fmicb.2017.01252

PubMed Abstract | CrossRef Full Text | Google Scholar

Cheng, J.-H., Wang, Y., Zhang, X.-Y., Sun, M.-L., Zhang, X., Song, X.-Y., et al. (2021). Characterization and diversity analysis of the extracellular proteases of thermophilic Anoxybacillus caldiproteolyticus 1A02591 from Deep-Sea hydrothermal vent sediment. Front. Microbiol. 12:643508. doi: 10.3389/fmicb.2021.643508

PubMed Abstract | CrossRef Full Text | Google Scholar

Chiu, B. K., Kato, S., McAllister, S. M., Field, E. K., and Chan, C. S. (2017). Novel pelagic Iron-oxidizing Zetaproteobacteria from the Chesapeake Bay Oxic–anoxic transition zone. Front. Microbiol. 8:1280. doi: 10.3389/fmicb.2017.01280

CrossRef Full Text | Google Scholar

Cordone, A., D’Errico, G., Magliulo, M., Bolinesi, F., Selci, M., Basili, M., et al. (2022). Bacterioplankton diversity and distribution in relation to phytoplankton community structure in the Ross Sea surface waters. Front. Microbiol. 13:722900. doi: 10.3389/fmicb.2022.722900

PubMed Abstract | CrossRef Full Text | Google Scholar

Cordone, A., Lucchini, S., Felice, M., and Ricca, E. (2011). Direct and indirect control of Lrp on LEE pathogenicity genes of Citrobacter rodentium. FEMS Microbiol. Lett. 325, 64–70. doi: 10.1111/j.1574-6968.2011.02411.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Cordone, A., Mauriello, E. M. F., Pickard, D. J., Dougan, G., De Felice, M., and Ricca, E. (2005). The lrp gene and its role in type I fimbriation in Citrobacter rodentium. J. Bacteriol. 187, 7009–7017. doi: 10.1128/JB.187.20.7009-7017.2005

PubMed Abstract | CrossRef Full Text | Google Scholar

Dahle, H., Økland, I., Thorseth, I. H., Pederesen, R. B., and Steen, I. H. (2015). Energy landscapes shape microbial communities in hydrothermal systems on the Arctic Mid-Ocean ridge. ISME J. 9, 1593–1606. doi: 10.1038/ismej.2014.247

PubMed Abstract | CrossRef Full Text | Google Scholar

Dávila-Ramos, S., Estradas-Romero, A., Prol-Ledesma, R. M., and Juárez-López, K. (2015). Bacterial populations (first record) at two shallow hydrothermal vents of the Mexican Pacific West coast. Geomicrobiol J. 32, 657–665. doi: 10.1080/01490451.2014.980526

CrossRef Full Text | Google Scholar

De Astis, G., Peccerillo, A., Kempton, P. D., La Volpe, L., and Wu, T. W. (2000). Transition from calc-alkaline to potassium-rich magmatism in subduction environments: geochemical and Sr, Nd, Pb isotopic constraints from the island of Vulcano (Aeolian arc). Contrib. Mineral. Petrol. 139, 684–703. doi: 10.1007/s004100000172

CrossRef Full Text | Google Scholar

De Astis, G., Ventura, G., and Vilardo, G. (2003). Geodynamic significance of the Aeolian volcanism (Southern Tyrrhenian Sea, Italy) in light of structural, seismological, and geochemical data: geodynamics of the Aeolian volcanism. Tectonics 22:1040. doi: 10.1029/2003TC001506

CrossRef Full Text | Google Scholar

Deckert, G., Warren, P. V., Gaasterland, T., Young, W. G., Lenox, A. L., Graham, D. E., et al. (1998). The complete genome of the hyperthermophilic bacterium Aquifex aeolicus. Nature 392, 353–358. doi: 10.1038/32831

PubMed Abstract | CrossRef Full Text | Google Scholar

Derekova, A., Sjøholm, C., Mandeva, R., and Kambourova, M. (2007). Anoxybacillus rupiensis sp. Nov., a novel thermophilic bacterium isolated from Rupi basin (Bulgaria). Extremophiles 11, 577–583. doi: 10.1007/s00792-007-0071-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Dick, G. J., Anantharaman, K., Baker, B. J., Li, M., Reed, D. C., and Sheik, C. S. (2013). The microbiology of deep-sea hydrothermal vent plumes: ecological and biogeographic linkages to seafloor and water column habitats. Front. Microbiol. 4:124. doi: 10.3389/fmicb.2013.00124

PubMed Abstract | CrossRef Full Text | Google Scholar

Ding, J., Zhang, Y., Wang, H., Jian, H., Leng, H., and Xiao, X. (2017). Microbial community structure of deep-sea hydrothermal vents on the ultraslow spreading southwest Indian ridge. Front. Microbiol. 8:1012. doi: 10.3389/fmicb.2017.01012

PubMed Abstract | CrossRef Full Text | Google Scholar

Du, Z.-J., Wang, Z.-J., Zhao, J.-X., and Chen, G.-J. (2016). Woeseia oceani gen. Nov., sp. nov., a chemoheterotrophic member of the order Chromatiales, and proposal of Woeseiaceae fam. Nov. Int. J. Syst. Evol. Microbiol. 66, 107–112. doi: 10.1099/ijsem.0.000683

PubMed Abstract | CrossRef Full Text | Google Scholar

Emerson, D., Fleming, E. J., and McBeth, J. M. (2010). Iron-oxidizing Bacteria: an environmental and genomic perspective. Annu. Rev. Microbiol. 64, 561–583. doi: 10.1146/annurev.micro.112408.134208

PubMed Abstract | CrossRef Full Text | Google Scholar

Emerson, D., and Moyer, C. L. (2002). Neutrophilic Fe-oxidizing Bacteria are abundant at the Loihi seamount hydrothermal vents and play a major role in Fe oxide deposition. Appl. Environ. Microbiol. 68, 3085–3093. doi: 10.1128/AEM.68.6.3085-3093.2002

PubMed Abstract | CrossRef Full Text | Google Scholar

Fabiano, M., Danovaro, R., and Fraschetti, S. (1995). A three-year time series of elemental and biochemical composition of organic matter in subtidal sandy sediments of the Ligurian Sea (northwestern Mediterranean). Cont. Shelf Res. 15, 1453–1469. doi: 10.1016/0278-4343(94)00088-5

CrossRef Full Text | Google Scholar

Fagorzi, C., Del Duca, S., Venturi, S., Chiellini, C., Bacci, G., Fani, R., et al. (2019). Bacterial communities from extreme environments: Vulcano Island. Diversity 11:140. doi: 10.3390/d11080140

PubMed Abstract | CrossRef Full Text | Google Scholar

Falkowski, P. G., Fenchel, T., and Delong, E. F. (2008). The microbial engines that drive Earth’s biogeochemical cycles. Science 320, 1034–1039. doi: 10.1126/science.1153213

PubMed Abstract | CrossRef Full Text | Google Scholar

Favallim, M., Karátson, D., Mazzuoli, R., Pareschi, M. T., and Ventura, G. (2005). Volcanic geomorphology and tectonics of the Aeolian archipelago (southern Italy) based on integrated DEM data. Bull. Volcanol. 68, 157–170. doi: 10.1007/s00445-005-0429-3

CrossRef Full Text | Google Scholar

Forget, N. L., Murdock, S. A., and Juniper, S. K. (2010). Bacterial diversity in Fe-rich hydrothermal sediments at two South Tonga arc submarine volcanoes: bacterial diversity in Fe oxide samples. Geobiology 8, 417–432. doi: 10.1111/j.1472-4669.2010.00247.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Fortney, N. W., He, S., Converse, B. J., Boyd, E. S., and Roden, E. E. (2018). Investigating the composition and metabolic potential of microbial communities in chocolate pots Hot Springs. Front. Microbiol. 9:2075. doi: 10.3389/fmicb.2018.02075

PubMed Abstract | CrossRef Full Text | Google Scholar

Fournier, R. O. (2007). “Hydrothermal systems and volcano geochemistry” in Volcano deformation. ed. D. Dzurisin (Berlin, Heidelberg: Springer Berlin Heidelberg), 323–341.

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

Gamberi, F., Marani, M., and Savelli, C. (1997). Tectonic, volcanic and hydrothermal features of a submarine portion of the Aeolian arc (Tyrrhenian Sea). Mar. Geol. 140, 167–181. doi: 10.1016/S0025-3227(97)00020-0

CrossRef Full Text | Google Scholar

Gavriilidou, A., Gutleben, J., Versluis, D., Forgiarini, F., van Passel, M. W. J., Ingham, C. J., et al. (2020). Comparative genomic analysis of Flavobacteriaceae: insights into carbohydrate metabolism, gliding motility and secondary metabolite biosynthesis. BMC Genomics 21:569. doi: 10.1186/s12864-020-06971-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Gerchakov, S. M., and Hatcher, P. G. (1972). Improved technique for analysis of carbohydrates in sediments. 1. Limnol. Oceanogr. 17, 938–943. doi: 10.4319/lo.1972.17.6.0938

PubMed Abstract | CrossRef Full Text | Google Scholar

German, C. R., and Seyfried, W. E. (2014). “Hydrothermal processes” in Treatise on geochemistry. ed. K. K. Turekian (Amsterdam, Netherlands: Elsevier), 191–233.

Google Scholar

Giovannelli, D., d’Errico, G., Fiorentino, F., Fattorini, D., Regoli, F., Angeletti, L., et al. (2016). Diversity and distribution of prokaryotes within a shallow-water pockmark Field. Front. Microbiol. 7:941. doi: 10.3389/fmicb.2016.00941

PubMed Abstract | CrossRef Full Text | Google Scholar

Giovannelli, D., d’Errico, G., Manini, E., Yakimov, M., and Vetriani, C. (2013). Diversity and phylogenetic analyses of bacteria from a shallow-water hydrothermal vent in Milos island (Greece). Front. Microbiol. 4:184. doi: 10.3389/fmicb.2013.00184

PubMed Abstract | CrossRef Full Text | Google Scholar

Giovannelli, D., Grosche, A., Starovoytov, V., Yakimov, M., Manini, E., and Vetriani, C. (2012). Galenea microaerophila gen. Nov., sp. nov., a mesophilic, microaerophilic, chemosynthetic, thiosulfate-oxidizing bacterium isolated from a shallow-water hydrothermal vent. Int. J. Syst. Evol. Microbiol. 62, 3060–3066. doi: 10.1099/ijs.0.040808-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Giovannelli, D., Sievert, S. M., Hügler, M., Markert, S., Becher, D., Schweder, T., et al. (2017). Insight into the evolution of microbial metabolism from the deep-branching bacterium, Thermovibrio ammonificans. Elife 6:e18990. doi: 10.7554/eLife.18990

PubMed Abstract | CrossRef Full Text | Google Scholar

Giovannelli, D., and Vetriani, C. (2017). From extreme environments to human pathogens: an evolutionary journey. Biochemist 39, 4–9. doi: 10.1042/BIO03906004

CrossRef Full Text | Google Scholar

Gomez-Saez, G. V., Pop Ristova, P., Sievert, S. M., Elvert, M., Hinrichs, K.-U., and Bühring, S. I. (2017). Relative importance of chemoautotrophy for primary production in a light exposed marine shallow hydrothermal system. Front. Microbiol. 8:702. doi: 10.3389/fmicb.2017.00702

PubMed Abstract | CrossRef Full Text | Google Scholar

Gomez-Saez, G. V., Riedel, T., Niggemann, J., Pichler, T., Dittmar, T., and Bühring, S. I. (2015). Interaction between iron and dissolved organic matter in a marine shallow hydrothermal system off Dominica Island (Lesser Antilles). Mar. Chem. 177, 677–686. doi: 10.1016/j.marchem.2015.10.003

CrossRef Full Text | Google Scholar

Gugliandolo, C., and Italiano, F. (1999). Submarine hydrothermal vents of the Aeolian Islands: relationship between microbial communities and thermal fluids. Geomicrobiol J. 16, 105–117. doi: 10.1080/014904599270794

CrossRef Full Text | Google Scholar

Gugliandolo, C., Italiano, F., and Maugeri, T. L. (2006). The submarine hydrothermal system of Panarea (southern Italy): Biogeochemical processes at the thermal fluids-sea bottom interface.

Google Scholar

Gugliandolo, C., Lentini, V., Bunk, B., Overmann, J., Italiano, F., and Maugeri, T. L. (2015). Changes in prokaryotic community composition accompanying a pronounced temperature shift of a shallow marine thermal brine pool (Panarea Island, Italy). Extremophiles 19, 547–559. doi: 10.1007/s00792-015-0737-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Gugliandolo, C., Lentini, V., Spanò, A., and Maugeri, T. L. (2012). New bacilli from shallow hydrothermal vents of Panarea Island (Italy) and their biotechnological potential: Bacilli from shallow vents of Panarea Island (Italy). J. Appl. Microbiol. 112, 1102–1112. doi: 10.1111/j.1365-2672.2012.05272.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Gugliandolo, C., and Maugeri, T. L. (2019). Phylogenetic diversity of Archaea in shallow hydrothermal vents of Eolian Islands, Italy. Diversity 11:156. doi: 10.3390/d11090156

CrossRef Full Text | Google Scholar

Guo, L.-L., Wu, Y.-H., Xu, X.-W., Huang, C.-J., Xu, Y.-Y., Cheng, H., et al. (2017). Actibacterium pelagium sp. nov., a novel alphaproteobacterium, and emended description of the genus Actibacterium. Int. J. Syst. Evol. Microbiol. 67, 5080–5086. doi: 10.1099/ijsem.0.002417

CrossRef Full Text | Google Scholar

Hafenbradl, D., Keller, M., Dirmeier, R., Rachel, R., Roßnagel, P., Burggraf, S., et al. (1996). Ferroglobus placidus gen. Nov., sp. nov., a novel hyperthermophilic archaeum that oxidizes Fe2+ at neutral pH under anoxic conditions. Arch. Microbiol. 166, 308–314. doi: 10.1007/s002030050388

PubMed Abstract | CrossRef Full Text | Google Scholar

Han, Y., and Perner, M. (2015). The globally widespread genus Sulfurimonas: versatile energy metabolisms and adaptations to redox clines. Front. Microbiol. 6:989. doi: 10.3389/fmicb.2015.00989

PubMed Abstract | CrossRef Full Text | Google Scholar

Handley, K. M., Boothman, C., Mills, R. A., Pancost, R. D., and Lloyd, J. R. (2010). Functional diversity of bacteria in a ferruginous hydrothermal sediment. ISME J. 4, 1193–1205. doi: 10.1038/ismej.2010.38

PubMed Abstract | CrossRef Full Text | Google Scholar

Hartree, E. F. (1972). Determination of protein: a modification of the Lowry method that gives a linear photometric response. Anal. Biochem. 48, 422–427. doi: 10.1016/0003-2697(72)90094-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Hedlund, B. P., Reysenbach, A.-L., Huang, L., Ong, J. C., Liu, Z., Dodsworth, J. A., et al. (2015). Isolation of diverse members of the Aquificales from geothermal springs in Tengchong, China. Front. Microbiol. 6:157. doi: 10.3389/fmicb.2015.00157

PubMed Abstract | CrossRef Full Text | Google Scholar

Hoaki, T., Nishijima, M., Miyashita, H., and Maruyama, T. (1995). Dense Community of Hyperthermophilic Sulfur-Dependent Heterotrophs in a geothermally heated shallow submarine biotope near Kodakara-Jima Island, Kagoshima, Japan. Appl. Environ. Microbiol. 61, 1931–1937. doi: 10.1128/aem.61.5.1931-1937.1995

PubMed Abstract | CrossRef Full Text | Google Scholar

Holmström, C., and Kjelleberg, S. (1999). Marine Pseudoalteromonas species are associated with higher organisms and produce biologically active extracellular agents. FEMS Microbiol. Ecol. 30, 285–293. doi: 10.1016/S0168-6496(99)00063-X

PubMed Abstract | CrossRef Full Text | Google Scholar

Hoshino, T., Kuratomi, T., Morono, Y., Hori, T., Oiwane, H., Kiyokawa, S., et al. (2016). Ecophysiology of Zetaproteobacteria associated with shallow hydrothermal Iron-Oxyhydroxide deposits in Nagahama Bay of Satsuma Iwo-Jima, Japan. Front. Microbiol. 6:1554. doi: 10.3389/fmicb.2015.01554

PubMed Abstract | CrossRef Full Text | Google Scholar

Inagaki, F., Takai, K., Nealson, K. H., and Horikoshi, K. (2004). Sulfurovum lithotrophicum gen. Nov., sp. nov., a novel sulfur-oxidizing chemolithoautotroph within the ε-Proteobacteria isolated from Okinawa trough hydrothermal sediments. Int. J. Syst. Evol. Microbiol. 54, 1477–1482. doi: 10.1099/ijs.0.03042-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Inguaggiato, S., Calderone, L., Inguaggiato, C., Mazot, A., Morici, S., and Vita, F. (2012). Long-time variation of soil CO2 fluxes at the summit crater of Vulcano (Italy). Bull. Volcanol. 74, 1859–1863. doi: 10.1007/s00445-012-0637-6

CrossRef Full Text | Google Scholar

Inguaggiato, S., Jácome Paz, M. P., Mazot, A., Delgado Granados, H., Inguaggiato, C., and Vita, F. (2013). CO2 output discharged from Stromboli Island (Italy). Chem. Geol. 339, 52–60. doi: 10.1016/j.chemgeo.2012.10.008

CrossRef Full Text | Google Scholar

Inskeep, W. J., Jay, Z. G., Tringe, S. J., Herrgård, M. B., and Rusch, D., YNP Metagenome Project Steering Committee and Working Group Members (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

Italiano, F., and Nuccio, P. M. (1991). Geochemical investigations of submarine volcanic exhalations to the east of Panarea, Aeolian Islands, Italy. Clin. Lab. Stand. Inst. 46, 125–141. doi: 10.1016/0377-0273(91)90079-F

CrossRef Full Text | Google Scholar

Jorgensen, S. L., Hannisdal, B., Lanzen, A., Baumberger, T., Flesland, K., Fonseca, R., et al. (2012). Correlating microbial community profiles with geochemical data in highly stratified sediments from the Arctic Mid-Ocean ridge. Proc. Natl. Acad. Sci. 109, E2846–E2855. doi: 10.1073/pnas.1207574109

PubMed Abstract | CrossRef Full Text | Google Scholar

Kawaichi, S., Ito, N., Kamikawa, R., Sugawara, T., Yoshida, T., and Sako, Y. (2013). Ardenticatena maritima gen. Nov., sp. nov., a ferric iron- and nitrate-reducing bacterium of the phylum ‘Chloroflexi’ isolated from an iron-rich coastal hydrothermal field, and description of Ardenticatenia classis nov. Int. J. Syst. Evol. Microbiol. 63, 2992–3002. doi: 10.1099/ijs.0.046532-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Kerfahi, D., Hall-Spencer, J. M., Tripathi, B. M., Milazzo, M., Lee, J., and Adams, J. M. (2014). Shallow water marine sediment bacterial community shifts along a natural CO2 gradient in the Mediterranean Sea off Vulcano, Italy. Microb. Ecol. 67, 819–828. doi: 10.1007/s00248-014-0368-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, Y., Jeon, J., Kwak, M. S., Kim, G. H., Koh, I., and Rho, M. (2018). Photosynthetic functions of Synechococcus in the ocean microbiomes of diverse salinity and seasons. PLoS One 13:e0190266. doi: 10.1371/journal.pone.0190266

PubMed Abstract | CrossRef Full Text | Google Scholar

Kürzinger, V.. (2019). Determination and differentiation of the hydrothermal precipitates of Panarea, Italy. Freiberg Online Geoscience, p. 54.

Google Scholar

LaRowe, D. E., Arndt, S., Bradley, J. A., Estes, E. R., Hoarfrost, A., Lang, S. Q., et al. (2020). The fate of organic carbon in marine sediments–new insights from recent data and analysis. Earth-Sci. Rev. 204:103146. doi: 10.1016/j.earscirev.2020.103146

CrossRef Full Text | Google Scholar

Lentini, V., Gugliandolo, C., Bunk, B., Overmann, J., and Maugeri, T. L. (2014). Diversity of prokaryotic Community at a Shallow Marine Hydrothermal Site Elucidated by Illumina sequencing technology. Curr. Microbiol. 69, 457–466. doi: 10.1007/s00284-014-0609-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Lentini, V., Gugliandolo, C., and Maugeri, T. L. (2007). Identification of enzyme-producing thermophilic bacilli isolated from marine vents of Aeolian Islands (Italy). Ann. Microbiol. 57, 355–361. doi: 10.1007/BF03175073

CrossRef Full Text | Google Scholar

Littlechild, J. A. (2015). Enzymes from extreme environments and their industrial applications. Front. Bioeng. Biotechnol. 3:161. doi: 10.3389/fbioe.2015.00161

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, J., Hua, Z.-S., Chen, L.-X., Kuang, J.-L., Li, S.-J., Shu, W.-S., et al. (2014). Correlating microbial diversity patterns with geochemistry in an extreme and heterogeneous environment of mine tailings. Appl. Environ. Microbiol. 80, 3677–3686. doi: 10.1128/AEM.00294-14

PubMed Abstract | CrossRef Full Text | Google Scholar

Lorenzen, C. J., and Downs, J. N. (1986). The specific absorption coefficients of chlorophyllide and pheophorbide in 90% acetone, and comments on the fluorometric determination of chlorophyll and pheopigments1. Limnol. Oceanogr. 31, 449–452. doi: 10.4319/lo.1986.31.2.0449

CrossRef Full Text | Google Scholar

Ludwig, M., and Bryant, D. A. (2012). Synechococcus sp. strain PCC 7002 transcriptome: acclimation to temperature, salinity, oxidative stress, and Mixotrophic growth conditions. Front. Microbiol. 3:354. doi: 10.3389/fmicb.2012.00354

PubMed Abstract | CrossRef Full Text | Google Scholar

Manini, E., Fiordelmondo, C., Gambi, C., Pusceddu, A., and Danovaro, R. (2003). Benthic microbial loop functioning in coastal lagoons: a comparative approach. Oceanol. Acta 26, 27–38. doi: 10.1016/S0399-1784(02)01227-6

CrossRef Full Text | Google Scholar

Manini, E., Luna, G. M., Corinaldesi, C., Zeppilli, D., Bortoluzzi, G., Caramanna, G., et al. (2008). Prokaryote diversity and virus abundance in shallow hydrothermal vents of the Mediterranean Sea (Panarea Island) and the Pacific Ocean (North Sulawesi-Indonesia). Microb. Ecol. 55, 626–639. doi: 10.1007/s00248-007-9306-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Maugeri, T. L., Lentini, V., Gugliandolo, C., Cousin, S., and Stackebrandt, E. (2010). Microbial diversity at a hot, Shallow-Sea hydrothermal vent in the southern Tyrrhenian Sea (Italy). Geomicrobiol J. 27, 380–390. doi: 10.1080/01490450903451518

CrossRef Full Text | Google Scholar

Maugeri, T. L., Lentini, V., Gugliandolo, C., Italiano, F., Cousin, S., and Stackebrandt, E. (2009). Bacterial and archaeal populations at two shallow hydrothermal vents off Panarea Island (Eolian Islands, Italy). Extremophiles 13, 199–212. doi: 10.1007/s00792-008-0210-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Maugeri, T. L., Lentini, V., Spanò, A., and Gugliandolo, C. (2013). Abundance and diversity of Picocyanobacteria in shallow hydrothermal vents of Panarea Island (Italy). Geomicrobiol J. 30, 93–99. doi: 10.1080/01490451.2011.653088

CrossRef Full Text | Google Scholar

McAllister, S. M., Moore, R. M., Gartman, A., Luther, G. W., Emerson, D., and Chan, C. S. (2019). The Fe(II)-oxidizing Zetaproteobacteria: historical, ecological and genomic perspectives. FEMS Microbiol. Ecol. 95:fiz015. doi: 10.1093/femsec/fiz015

PubMed Abstract | CrossRef Full Text | Google Scholar

McDermott, J. M., Sylva, S. P., Ono, S., German, C. R., and Seewald, J. S. (2018). Geochemistry of fluids from Earth’s deepest ridge-crest hot-springs: Piccard hydrothermal field, mid-Cayman rise. Geochim. Cosmochim. Acta 228, 95–118. doi: 10.1016/j.gca.2018.01.021

CrossRef Full Text | Google Scholar

McGonigle, J. M., Lang, S. Q., and Brazelton, W. J. (2020). Genomic evidence for Formate metabolism by Chloroflexi as the key to unlocking deep carbon in lost City microbial ecosystems. Appl. Environ. Microbiol. 86, e02583–e02519. doi: 10.1128/AEM.02583-19

CrossRef Full Text | Google Scholar

McIlroy, S. J., Kirkegaard, R. H., Dueholm, M. S., Fernando, E., Karst, S. M., Albertsen, M., et al. (2017). Culture-independent analyses reveal novel Anaerolineaceae as abundant primary fermenters in anaerobic digesters treating waste activated sludge. Front. Microbiol. 8:1134. doi: 10.3389/fmicb.2017.01134

PubMed Abstract | CrossRef Full Text | Google Scholar

McMurdie, P. J., and Holmes, S. (2013). Phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS One 8:e61217. doi: 10.1371/journal.pone.0061217

PubMed Abstract | CrossRef Full Text | Google Scholar

Meyer-Dombard, D. R., Amend, J. P., and Osburn, M. R. (2013). Microbial diversity and potential for arsenic and iron biogeochemical cycling at an arsenic rich, shallow-sea hydrothermal vent (Tutum Bay, Papua New Guinea). Chem. Geol. 348, 37–47. doi: 10.1016/j.chemgeo.2012.02.024

CrossRef Full Text | Google Scholar

Meyer-Dombard, D. R., Shock, E. L., and Amend, J. P. (2012). Effects of trace element concentrations on culturing thermophiles. Extremophiles 16, 317–331. doi: 10.1007/s00792-012-0432-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Miroshnichenko, M. L., Kostrikina, N. A., Chernyh, N. A., Pimenov, N. V., Tourova, T. P., Antipov, A. N., et al. (2003). Caldithrix abyssi gen. Nov., sp. nov., a nitrate-reducing, thermophilic, anaerobic bacterium isolated from a mid-Atlantic ridge hydrothermal vent, represents a novel bacterial lineage. Int. J. Syst. Evol. Microbiol. 53, 323–329. doi: 10.1099/ijs.0.02390-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Morris, R. M., and Spietz, R. L. (2022). The physiology and biogeochemistry of SUP05. Annu. Rev. Mar. Sci. 14, 261–275. doi: 10.1146/annurev-marine-010419-010814

PubMed Abstract | CrossRef Full Text | Google Scholar

Nogi, Y., Hosoya, S., Kato, C., and Horikoshi, K. (2004). Colwellia piezophila sp. nov., a novel piezophilic species from deep-sea sediments of the Japan trench. Int. J. Syst. Evol. Microbiol. 54, 1627–1631. doi: 10.1099/ijs.0.03049-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Oksanen, J., Blanchet, F. G., Friendly, M., Kindt, R., Legendre, P., McGlinn, D., et al. (2018). Vegan: Community ecology package. Available at: https://CRAN.R-project.org/package=vegan (Accessed January 23, 2019).

Google Scholar

Pearce, R. K., Sánchez de la Muela, A., Moorkamp, M., Hammond, J. O. S., Mitchell, T. M., Cembrano, J., et al. (2019). Interaction between hydrothermal fluids and fault systems in the in the southern Andes revealed by magnetotelluric and seismic data. Geophysics 2019:1143. doi: 10.1002/essoar.10501143.1

CrossRef Full Text | Google Scholar

Peccerillo, A., De Astis, G., Faraone, D., Forni, F., and Frezzotti, M. L. (2013). Compositional variations of magmas in the Aeolian arc: implications for petrogenesis and geodynamics. Geol. Soc. Lond. Mem. 37, 491–510. doi: 10.1144/M37.15

CrossRef Full Text | Google Scholar

Peccerillo, A., and Turco, E. (2004). Petrological and geochemical variations of Plio-quaternary volcanism in the Tyrrhenian Sea area: regional distribution of magma types, petrogenesis and geodynamic implications. Per. Mineral 73, 231–251. doi: 10.1016/S1871-644X(06)80016-1

CrossRef Full Text | Google Scholar

Petrillo, C., Castaldi, S., Lanzilli, M., Selci, M., Cordone, A., Giovannelli, D., et al. (2021). Genomic and physiological characterization of Bacilli isolated from salt-pans with plant growth promoting features. Front. Microbiol. 12:2517. doi: 10.3389/fmicb.2021.715678

PubMed Abstract | CrossRef Full Text | Google Scholar

Placido, A., Hai, T., Ferrer, M., Chernikova, T. N., Distaso, M., Armstrong, D., et al. (2015). Diversity of hydrolases from hydrothermal vent sediments of the Levante Bay, Vulcano Island (Aeolian archipelago) identified by activity-based metagenomics and biochemical characterization of new esterases and an arabinopyranosidase. Appl. Microbiol. Biotechnol. 99, 10031–10046. doi: 10.1007/s00253-015-6873-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Pohlner, M., Dlugosch, L., Wemheuer, B., Mills, H., Engelen, B., and Reese, B. K. (2019). The majority of active Rhodobacteraceae in marine sediments belong to uncultured genera: a molecular approach to link their distribution to environmental conditions. Front. Microbiol. 10:659. doi: 10.3389/fmicb.2019.00659

PubMed Abstract | CrossRef Full Text | Google Scholar

Power, J. F., Carere, C. R., Lee, C. K., Wakerley, G. L. J., 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

Price, R. E., and Giovannelli, D. (2017). “A review of the geochemistry and microbiology of marine shallow-water hydrothermal vents” in Reference module in earth systems and environmental sciences ed. S. A. Elias (Amsterdam, Netherlands: Elsevier).

Google Scholar

Price, R. E., LaRowe, D. E., Italiano, F., Savov, I., Pichler, T., and Amend, J. P. (2015). Subsurface hydrothermal processes and the bioenergetics of chemolithoautotrophy at the shallow-sea vents off Panarea Island (Italy). Chem. Geol. 407–408, 21–45. doi: 10.1016/j.chemgeo.2015.04.011

CrossRef Full Text | Google Scholar

Price, R. E., Lesniewski, R., Nitzsche, K. S., Meyerdierks, A., Saltikov, C., Pichler, T., et al. (2013). Archaeal and bacterial diversity in an arsenic-rich shallow-sea hydrothermal system undergoing phase separation. Front. Microbiol. 4:158. doi: 10.3389/fmicb.2013.00158

PubMed Abstract | CrossRef Full Text | Google Scholar

Priest, T., Heins, A., Harder, J., Amann, R., and Fuchs, B. M. (2022). Niche partitioning of the ubiquitous and ecologically relevant NS5 marine group. ISME J. 16, 1570–1582. doi: 10.1038/s41396-022-01209-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Pusceddu, A., Dell’Anno, A., Danovaro, R., Manini, E., Sara, G., and Fabiano, M. (2003). Enzymatically hydrolyzable protein and carbohydrate sedimentary pools as indicators of the trophic state of detritus sink systems: a case study in a Mediterranean coastal lagoon. Estuaries 26, 641–650. doi: 10.1007/BF02711976

CrossRef Full Text | Google Scholar

Quast, C., Pruesse, E., Yilmaz, P., Gerken, J., Schweer, T., Yarza, P., et al. (2013). The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 41, D590–D596. doi: 10.1093/nar/gks1219

PubMed Abstract | CrossRef Full Text | Google Scholar

R Core Team. (2021). R: A language and environment for statistical computing, 2015. Vienna, Austria: R Core Team.

Google Scholar

Rajasabapathy, R., Mohandass, C., Bettencourt, R., Colaço, A., Goulart, J., and Meena, R. M. (2018). Bacterial diversity at a shallow-water hydrothermal vent (Espalamaca) in Azores Island. Curr. Sci. 115:2110. doi: 10.18520/cs/v115/i11/2110-2121

CrossRef Full Text | Google Scholar

Reed, D. C., Breier, J. A., Jiang, H., Anantharaman, K., Klausmeier, C. A., Toner, B. M., et al. (2015). Predicting the response of the deep-ocean microbiome to geochemical perturbations by hydrothermal vents. ISME J. 9, 1857–1869. doi: 10.1038/ismej.2015.4

PubMed Abstract | CrossRef Full Text | Google Scholar

Reeves, E. P., Seewald, J. S., Saccocia, P., Bach, W., Craddock, P. R., Shanks, W. C., et al. (2011). Geochemistry of hydrothermal fluids from the PACMANUS, northeast Pual and Vienna woods hydrothermal fields, Manus Basin. Papua New Guinea. Geochim. Cosmochim. Acta 75, 1088–1123. doi: 10.1016/j.gca.2010.11.008

CrossRef Full Text | Google Scholar

Rinke, C., Schmitz-Esser, S., Stoecker, K., Nussbaumer, A. D., Molnár, D. A., Vanura, K., et al. (2006). “Candidatus Thiobios zoothamnicoli,” an Ectosymbiotic bacterium covering the Giant marine ciliate Zoothamnium niveum. Appl. Environ. Microbiol. 72, 2014–2021. doi: 10.1128/AEM.72.3.2014-2021.2006

PubMed Abstract | CrossRef Full Text | Google Scholar

Rogers, K. L., and Amend, J. P. (2005). Archaeal diversity and geochemical energy yields in a geothermal well on Vulcano Island, Italy. Geobiology 3, 319–332. doi: 10.1111/j.1472-4669.2006.00064.x

CrossRef Full Text | Google Scholar

Rogers, K. L., and Amend, J. P. (2006). Energetics of potential heterotrophic metabolisms in the marine hydrothermal system of Vulcano Island, Italy. Geochim. Cosmochim. Acta 70, 6180–6200. doi: 10.1016/j.gca.2006.08.046

CrossRef Full Text | Google Scholar

Rogers, K. L., Amend, J. P., and Gurrieri, S. (2007). Temporal changes in fluid chemistry and energy profiles in the Vulcano Island hydrothermal system. Astrobiology 7, 905–932. doi: 10.1089/ast.2007.0128

PubMed Abstract | CrossRef Full Text | Google Scholar

Rollinson, H., and Pease, V. (2021). Using geochemical data to understand geological processes. Cambridge, UK: Cambridge University Press.

Google Scholar

Rusch, A., Walpersdorf, E., deBeer, D., Gurrieri, S., and Amend, J. P. (2005). Microbial communities near the oxic/anoxic interface in the hydrothermal system of Vulcano Island, Italy. Chem. Geol. 224, 169–182. doi: 10.1016/j.chemgeo.2005.07.026

CrossRef Full Text | Google Scholar

Sayed, A. M., Hassan, M. H. A., Alhadrami, H. A., Hassan, H. M., Goodfellow, M., and Rateb, M. E. (2020). Extreme environments: microbiology leading to specialized metabolites. J. Appl. Microbiol. 128, 630–657. doi: 10.1111/jam.14386

PubMed Abstract | CrossRef Full Text | Google Scholar

Sciutteri, V., Smedile, F., Vizzini, S., Mazzola, A., and Vetriani, C. (2022). Microbial biofilms along a geochemical gradient at the shallow-water hydrothermal system of Vulcano Island, Mediterranean Sea. Front. Virol. 13:840205. doi: 10.3389/fmicb.2022.840205

PubMed Abstract | CrossRef Full Text | Google Scholar

Sedwick, P., and Stüben, D. (1996). Chemistry of shallow submarine warm springs in an arc-volcanic setting: Vulcano Island, Aeolian archipelago, Italy. Mar. Chem. 53, 147–161. doi: 10.1016/0304-4203(96)00020-5

CrossRef Full Text | Google Scholar

Shafiee, R. T., Diver, P. J., Snow, J. T., Zhang, Q., and Rickaby, R. E. M. (2021). Marine ammonia-oxidising archaea and bacteria occupy distinct iron and copper niches. ISME Commun. 1:1. doi: 10.1038/s43705-021-00001-7

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

Sieland, R. (2009). Chemical and isotopic investigations of submarine fluid discharges from Panarea, Aeolian Islands, Italy. Freiberg Online Geol. :21:887. doi: 10.23689/fidgeo-887

CrossRef Full Text | Google Scholar

Sievert, S. M., Kuever, J., and Muyzer, G. (2000). Identification of 16S ribosomal DNA-defined bacterial populations at a shallow submarine hydrothermal vent near Milos Island (Greece). Appl. Environ. Microbiol. 66, 3102–3109. doi: 10.1128/AEM.66.7.3102-3109.2000

PubMed Abstract | CrossRef Full Text | Google Scholar

Simmons, S., and Norris, P. (2002). Acidophiles of saline water at thermal vents of Vulcano, Italy. Extremophiles 6, 201–207. doi: 10.1007/s007920100242

PubMed Abstract | CrossRef Full Text | Google Scholar

Spietz, R. L., Lundeen, R. A., Zhao, X., Nicastro, D., Ingalls, A. E., and Morris, R. M. (2019). Heterotrophic carbon metabolism and energy acquisition in Candidatus Thioglobus singularis strain PS1, a member of the SUP05 clade of marine Gammaproteobacteria. Environ. Microbiol. 21, 2391–2401. doi: 10.1111/1462-2920.14623

PubMed Abstract | CrossRef Full Text | Google Scholar

Stöhr, R., Waberski, A., Völker, H., Tindall, B. J., and Thomm, M. (2001). Hydrogenothermus marinus gen. Nov., sp. nov., a novel thermophilic hydrogen-oxidizing bacterium, recognition of Calderobacterium hydrogenophilum as a member of the genus Hydrogenobacter and proposal of the reclassification of Hydrogenobacter acidophilus as Hydrogenobaculum acidophilum gen. Nov., comb. nov., in the phylum ‘Hydrogenobacter/Aquifex. Int. J. Syst. Evol. Microbiol. 51, 1853–1862. doi: 10.1099/00207713-51-5-1853

PubMed Abstract | CrossRef Full Text | Google Scholar

Suarez, C., Piculell, M., Modin, O., Langenheder, S., Persson, F., and Hermansson, M. (2019). Thickness determines microbial community structure and function in nitrifying biofilms via deterministic assembly. Sci. Rep. 9:5110. doi: 10.1038/s41598-019-41542-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Taillefer, A., Guillou-Frottier, L., Soliva, R., Magri, F., Lopez, S., Courrioux, G., et al. (2018). Topographic and faults control of hydrothermal circulation along dormant faults in an Orogen. Geochem. Geophys. Geosystems 19, 4972–4995. doi: 10.1029/2018GC007965

CrossRef Full Text | Google Scholar

Takai, K., Miyazaki, M., Nunoura, T., Hirayama, H., Oida, H., Furushima, Y., et al. (2006). Sulfurivirga caldicuralii gen. Nov., sp. nov., a novel microaerobic, thermophilic, thiosulfate-oxidizing chemolithoautotroph, isolated from a shallow marine hydrothermal system occurring in a coral reef. Jpn. Int. J. Syst. Evol. Microbiol. 56, 1921–1929. doi: 10.1099/ijs.0.64297-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Tarasov, V. G., Gebruk, A. V., Mironov, A. N., and Moskalev, L. I. (2005). Deep-sea and shallow-water hydrothermal vent communities: two different phenomena? Chem. Geol. 224, 5–39. doi: 10.1016/j.chemgeo.2005.07.021

CrossRef Full Text | Google Scholar

Tarasov, V. G., Gebruk, A. V., Shulkin, V. M., Kamenev, G. M., Fadeev, V. I., Kosmynin, V. N., et al. (1999). Effect of shallow-water hydrothermal venting on the biota of Matupi harbour (Rabaul caldera, New Britain Island, Papua New Guinea). Cont. Shelf Res. 19, 79–116. doi: 10.1016/S0278-4343(98)00073-9

CrossRef Full Text | Google Scholar

Tarasov, V. G., Propp, M. V., Propp, L. N., Zhirmunsky, A. V., Namsakakv, B. B., Gorlenko, V. M., et al. (1990). Shallow-water Gasohydrothermal vents of Ushishir volcano and the ecosystem of Kraternaya bight (the Kurile Islands). Mar. Ecol. 11, 1–23. doi: 10.1111/j.1439-0485.1990.tb00225.x

CrossRef Full Text | Google Scholar

Tassi, F., Capaccioni, B., Caramanna, G., Cinti, D., Montegrossi, G., Pizzino, L., et al. (2009). Low-pH waters discharging from submarine vents at Panarea Island (Aeolian Islands, southern Italy) after the 2002 gas blast: origin of hydrothermal fluids and implications for volcanic surveillance. Appl. Geochem. 24, 246–254. doi: 10.1016/j.apgeochem.2008.11.015

CrossRef Full Text | Google Scholar

Techtmann, S. M., Fitzgerald, K. S., Stelling, S. C., Joyner, D. C., Uttukar, S. M., Harris, A. P., et al. (2016). Colwellia psychrerythraea strains from distant Deep Sea basins show adaptation to local conditions. Front. Environ. Sci. 4:3. doi: 10.3389/fenvs.2016.00033

CrossRef Full Text | Google Scholar

Tivey, M. K. (2007). Generation of seafloor hydrothermal vent fluids and associated mineral deposits. Oceanography 20, 50–65. doi: 10.5670/oceanog.2007.80

CrossRef Full Text | Google Scholar

Van Dover, C. (2012). Hydrothermal vent ecosystems and conservation. Oceanography 25, 313–316. doi: 10.5670/oceanog.2012.36

PubMed Abstract | CrossRef Full Text | Google Scholar

Van Dover, C. L., Arnaud-Haond, S., Gianni, M., Helmreich, S., Huber, J. A., Jaeckel, A. L., et al. (2018). Scientific rationale and international obligations for protection of active hydrothermal vent ecosystems from deep-sea mining. Mar. Policy 90, 20–28. doi: 10.1016/j.marpol.2018.01.020

CrossRef Full Text | Google Scholar

Venâncio, I., Luís, Â., Domingues, F., Oleastro, M., Pereira, L., and Ferreira, S. (2022). The prevalence of Arcobacteraceae in aquatic environments: a systematic review and Meta-analysis. Pathogens 11:244. doi: 10.3390/pathogens11020244

PubMed Abstract | CrossRef Full Text | Google Scholar

Vuillemin, A., Kerrigan, Z., D’Hondt, S., and Orsi, W. D. (2020). Exploring the abundance, metabolic potential and gene expression of subseafloor Chloroflexi in million-year-old oxic and anoxic abyssal clay. FEMS Microbiol. Ecol. 96:fiaa223. doi: 10.1093/femsec/fiaa223

CrossRef Full Text | Google Scholar

Wang, L., Cheung, M. K., Liu, R., Wong, C. K., Kwan, H. S., and Hwang, J.-S. (2017). Diversity of Total bacterial communities and chemoautotrophic populations in sulfur-rich sediments of shallow-water hydrothermal vents off Kueishan Island, Taiwan. Microb. Ecol. 73, 571–582. doi: 10.1007/s00248-016-0898-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, W., Li, Z., Zeng, L., Dong, C., and Shao, Z. (2020). The oxidation of hydrocarbons by diverse heterotrophic and mixotrophic bacteria that inhabit deep-sea hydrothermal ecosystems. ISME J. 14, 1994–2006. doi: 10.1038/s41396-020-0662-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, F.-Q., Lin, X.-Z., Chen, G.-J., and Du, Z.-J. (2015). Colwellia arctica sp. nov., isolated from Arctic marine sediment. Antonie Van Leeuwenhoek 107, 723–729. doi: 10.1007/s10482-014-0366-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Ward, L. M., Hemp, J., Shih, P. M., McGlynn, S. E., and Fischer, W. W. (2018). Evolution of Phototrophy in the Chloroflexi phylum driven by horizontal gene transfer. Front. Microbiol. 9:260. doi: 10.3389/fmicb.2018.00260

PubMed Abstract | CrossRef Full Text | Google Scholar

Ward, L. M., Idei, A., Terajima, S., Kakegawa, T., Fischer, W. W., and McGlynn, S. E. (2017). Microbial diversity and iron oxidation at Okuoku-hachikurou Onsen, a Japanese hot spring analog of Precambrian iron formations. Geobiology 15, 817–835. doi: 10.1111/gbi.12266

PubMed Abstract | CrossRef Full Text | Google Scholar

White, J. R., Escobar-Paramo, P., Mongodin, E. F., Nelson, K. E., and DiRuggiero, J. (2008). Extensive genome rearrangements and multiple horizontal gene transfers in a population of Pyrococcus Isolates from Vulcano Island, Italy. Appl. Environ. Microbiol. 74, 6447–6451. doi: 10.1128/AEM.01024-08

PubMed Abstract | CrossRef Full Text | Google Scholar

Wickham, H. (2011). ggplot2. Wiley Interdiscip. Rev. Comput. Stat. 3, 180–185. doi: 10.1002/wics.147

CrossRef Full Text | Google Scholar

Yamamoto, K., Hackley, K. C., Kelly, W. R., Panno, S. V., Sekiguchi, Y., Sanford, R. A., et al. (2019). Diversity and geochemical community assembly processes of the living rare biosphere in a sand-and-gravel aquifer ecosystem in the Midwestern United States. Sci. Rep. 9:13484. doi: 10.1038/s41598-019-49996-z

CrossRef Full Text | Google Scholar

Yamamoto, M., and Takai, K. (2011). Sulfur metabolisms in epsilon- and gamma-Proteobacteria in Deep-Sea hydrothermal fields. Front. Microbiol. 2:192. doi: 10.3389/fmicb.2011.00192

PubMed Abstract | CrossRef Full Text | Google Scholar

Yücel, M., Sievert, S. M., Vetriani, C., Foustoukos, D. I., Giovannelli, D., and Le Bris, N. (2013). Eco-geochemical dynamics of a shallow-water hydrothermal vent system at Milos Island, Aegean Sea (eastern Mediterranean). Chem. Geol. 356, 11–20. doi: 10.1016/j.chemgeo.2013.07.020

CrossRef Full Text | Google Scholar

Zanfardino, A., Migliardi, A., D’Alonzo, D., Lombardi, A., Varcamonti, M., and Cordone, A. (2016). Inactivation of MSMEG_0412 gene drastically affects surface related properties of Mycobacterium smegmatis. BMC Microbiol. 16:267. doi: 10.1186/s12866-016-0888-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Zeng, X., Alain, K., and Shao, Z. (2021). Microorganisms from deep-sea hydrothermal vents. Mar. Life Sci. Technol. 3, 204–230. doi: 10.1007/s42995-020-00086-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, D.-C., Yu, Y., Xin, Y.-H., Liu, H.-C., Zhou, P.-J., and Zhou, Y.-G. (2008). Colwellia polaris sp. nov., a psychrotolerant bacterium isolated from Arctic Sea ice. Int. J. Syst. Evol. Microbiol. 58, 1931–1934. doi: 10.1099/ijs.0.65596-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y., Zhao, Z., Chen, C.-T. A., Tang, K., Su, J., and Jiao, N. (2012). Sulfur metabolizing microbes dominate microbial communities in andesite-hosted Shallow-Sea hydrothermal systems. PLoS One 7:e44593. doi: 10.1371/journal.pone.0044593

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: shallow-water vents, Aeolian archipelago, 16S rRNA amplicon sequencing, microbial diversity, marine protected areas, hydrothermal vents

Citation: Barosa B, Ferrillo A, Selci M, Giardina M, Bastianoni A, Correggia M, di Iorio L, Bernardi G, Cascone M, Capuozzo R, Intoccia M, Price R, Vetriani C, Cordone A and Giovannelli D (2023) Mapping the microbial diversity associated with different geochemical regimes in the shallow-water hydrothermal vents of the Aeolian archipelago, Italy. Front. Microbiol. 14:1134114. doi: 10.3389/fmicb.2023.1134114

Received: 29 December 2022; Accepted: 25 July 2023;
Published: 11 August 2023.

Edited by:

Brett J. Baker, The University of Texas at Austin, United States

Reviewed by:

Florence Schubotz, University of Bremen, Germany
Craig Everroad, National Aeronautics and Space Administration, United States

Copyright © 2023 Barosa, Ferrillo, Selci, Giardina, Bastianoni, Correggia, di Iorio, Bernardi, Cascone, Capuozzo, Intoccia, Price, Vetriani, Cordone and Giovannelli. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Donato Giovannelli, ZG9uYXRvLmdpb3Zhbm5lbGxpQHVuaW5hLml0

Present address: Alessandra Ferrillo, Norwegian University of Life Sciences, Aas, Norway

These authors share first authorship

Disclaimer: 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.