Microbial dynamics in shallow CO2 seeps system off Panarea Island (Italy)

Shallow-water hydrothermal vents are extreme environments characterized by high temperatures, low pH, and high CO2 concentrations; therefore, they are considered as suitable laboratories for studying the effect of global changes on marine microbes. We hypothesized a direct effect of vents on prokaryotic community structure and functioning in the Panarea Island’s hydrothermal system. Sampling was conducted along a 9-station transect characterized by three active emission points. The water column was stratified with a thermocline at 25 m depth and a deep chlorophyll maximum between 50 and 100 m. Prokaryotic abundance ranged from 0.2 to 1.5 × 109 cells L−1, prokaryotic carbon production from 2.4 to 75.4 ng C L−1 h−1, and exoenzymatic activities degrading proteins, phosphorylated compounds, and polysaccharides were on the order of 4–28, 2–31 and 0.2–4.16 nM h−1, respectively. While microbial abundance and production were shaped by the water column's physical structure, alkaline phosphatase and beta-glucosidase activities seemed to be enhanced by hydrothermal fluids. The 16S rRNA gene amplicon sequencing analysis identified a surface, a deep, and a vent-influenced microbial community. In terms of relative abundance members of the SAR11 group dominated the water column, alongside Synechococcus and Prochlorococcus in surface and bottom samples, respectively. Vent-influenced stations were characterized by the presence of Thiomicrorhabdus, a sulfur-oxidizer chemolithoautotroph. Overall, this study provides insights on the coupling between microbial community structure and the biogeochemical cycling of nutrients in low-pH conditions (CO2 and H2S-based), thus addressing some of the opened questions about the response of microbes to acidification.


Introduction
Marine hydrothermal systems are part of active submarine volcanic areas (Gugliandolo and Maugeri 2019). Fluids emitted from vents are characterized by high temperatures, and high concentrations of chemicals (CO 2 , H 2 S, hydrocarbons, and heavy metals), thus representing extreme habitats for marine microbes, impacting the local biogeochemical cycles (Jannasch et al. 1985;Karl 1995).
Since November 2002, when a gas emission spike was detected, Panarea's hydrothermal vent system became a critical area to study for a better understanding of hydrothermal vent dynamics (Italiano and Nuccio 1991;Maugeri et al. 2009). On a biological perspective, systems like the hydrothermal field of Panarea can be used as natural laboratories for studying the effect of ocean acidification on the marine biota. In fact, the increased concentration of carbon dioxide that lowers the seawater's pH can simulate acidification both due to increased atmospheric CO 2 partial pressure (Gonzáles-Delgado and Hernández 2018) and to leakages from CO 2 geological storage sites 97 Page 2 of 14 (e.g., Jones et al. 2015;Molari et al. 2018) with consequent solubilization in the ocean.
Several studies have investigated the impact of increased CO 2 concentration on planktonic prokaryotes under controlled and simulated conditions (Riebesell and Gattuso 2015;Boyd et al. 2018). On a general perspective, many studies have been focusing on the alteration of the organic carbon pool due to different responses of phytoplankton production to ocean acidification (e.g., Flynn et al. 2015;Bach and Taucher 2019;Taucher et al. 2021). Since organic matter is the source of C and energy for most organotrophic prokaryotes, a secondary effect on prokaryotic community structure and functioning has been hypothesized. For example, Tsiola et al. (2023) recently demonstrated that members of SAR11 can be favored over other taxa by acidified conditions, yet other studies have reported controversial results on community shifts (e.g., Liu et al. 2010 and references therein). Likewise, the effect of ocean acidification on the processing of organic matter by prokaryotes (degradation and utilization rates) has been shown to vary in different locations and trophic regimes (e.g., Sala et al. 2016;Celussi et al. 2017;Hu et al. 2021).
Utilizing a different approach, other authors have explored the possibility of exploiting the natural features of shallow hydrothermal vents to simulate future acidification scenarios (e.g., Lidbury et al. 2012;Gonzáles-Delgado and Hernández 2018;Tangherlini et al. 2021). However, a disadvantage of using natural systems of volcanic origin is the presence of other associated chemicals in hydrothermal fluids such as H 2 S and CH 4 (Gonzáles-Delgado and Hernández 2018) that can create a unique environment for the growth of specialized microorganisms. In fact, hydrothermal vent microbes are able to use the energy and chemicals released from vents to produce organic matter through chemosynthesis (Tarasov et al. 2005). An extensive body of literature indicates that low levels of S-compounds in the emitted fluids can promote the prevalence of sulfur cycling prokaryotes in the proximity of vents (Maugeri et al. 2009;González et al. 2020;Rizzo et al. 2022;). In addition to chemosynthetic microorganisms, photosynthetic bacteria have been found to play an important role in shallow hydrothermal vent systems, by converting CO 2 to organic matter contributing significantly to primary production (Gomez-Saez et al. 2016).
Studies at the shallow hydrothermal system of the Panarea Island, showed the presence of photosynthetic and chemolithotrophic primary producers both in seawater and sediments (e.g., Manini et al. 2008;Maugeri et al. 2009Maugeri et al. , 2010Karuza et al. 2012;Lentini et al. 2014;Bortoluzzi et al. 2017;Tangherlini et al. 2021;Rizzo et al. 2022 and references therein), highlighting fine-tuned adaptation and survival strategies in this peculiar environment.
Within this framework, we tested whether the vents (i.e., lower pH) influenced the composition of microbial communities and, for the first time in this environment, we combined the molecular-based data with the organic matter pool processing (i.e., degradation and utilization) in the water column.

Sampling design
Sampling was carried out between 28th and 31st July 2011 in the Panarea Island area ( Fig. 1; Supplementary Table S1) along a 9-station transect. One extra station located north of the vent field was sampled as a reference (i.e., no ventinfluenced). At each station, temperature, salinity, oxygen, and pH vertical profiles were registered using a CTD SBE9 Seabird. Seawater was collected with 10 L Niskin bottles for determining pH, the concentrations of chlorophyll a, macronutrients, dissolved inorganic carbon, particulate organic carbon, and nitrogen concentrations, as well as prokaryotic abundance, the rates of heterotrophic carbon production and extracellular enzymatic activities, and to depict the prokaryotic community structure by 16S rRNA gene amplicon sequencing (at seven stations). Samples were collected in the surface (~2 m depth, named surface) and bottom layers (collected a few meters above the seafloor, or at a depth ~45 m at the stations deeper than 100 m-stations P0 and P15-, named bottom). Additional intermediate sampling depths were selected according to the CTD profiles of temperature, salinity and pH (Supplementary Table S1). Within the transect, three active vents were detected (stations V23, V35, and V31).

Chemical and microbiological analyses
Seawater samples for determining chlorophyll a concentration were filtered through Whatman GF/F glass-fiber filters (0.7 µm nominal pore size, 47 mm diameter) and immediately frozen (-20 °C) until analysis. Pigments were extracted from the grinded filters overnight in the dark at 4 °C with 90% acetone and determined fluorometrically according to Lorenzen and Jeffrey (1980). The readings of the extracts were performed before and after acidification with two drops of HCl 1 N using a JASCO FP 6500 spectrofluorometer (450 nm excitation and 665 nm emission wavelengths). Calibration curves were prepared using pure chlorophyll a standard from spinach (Sigma-Aldrich).
Water samples for pH measurement (pH T ) were collected directly into quartz cuvettes with a 10 cm path length, overflowing with twice the cuvette volume and no head space. pH analyses were carried out immediately onboard against MilliQ water in a Varian Cary 50 single beam spectrophotometer with thermostatted cylindrical cuvette holders. A circulating bath was used to maintain the temperature at 25.0 °C. The temperature in the cuvette was controlled with a PT100 temperature sensor (precision: 0.01 °C, GMH 3750). The indicator m-cresol violet was added to the cuvette using a micropipette. The SOP6b ver 3.01 (Dickson et al. 2007) method was followed. The in situ pH, expressed on the total scale, was calculated on the basis of the in situ temperature and salinity measured by the CTD probe.
Samples for dissolved inorganic carbon concentration (DIC) analyses were collected in HCl-washed glass vials minimizing gas exchange with the atmosphere, treated with a mercuric chloride solution (0.04% by volume of a saturated aqueous solution) in order to prevent biological activity, sealed with a Teflon cap, and stored refrigerated until analyses. DIC was determined using a Shimadzu TOC-V CSH analyzer. Samples were injected into the instrument port and directly acidified with H 3 PO 4 (25%). Phosphoric acidification generated CO 2 that was carried to a nondispersive infrared detector (NDIR). The analysis showed a variation coefficient <2%. The reproducibility of the method was between 1.5 and 3%.
Seawater samples for dissolved inorganic nutrient analysis (NH 4 + , NO 2 − , NO 3 − and PO 4 3− ) were prefiltered on 0.7 µm pore size glass-fiber filters (Whatman GF/F) and then immediately frozen. Samples were analyzed at room temperature on a continuous flow autoanalyzer (Bran + Luebbe, Norderstedt, Germany), according to Koroleff and Grasshof (1983). The efficiency of the system was checked before and after sample analyses by doing replicates of internal standards. Dissolved inorganic nitrogen (DIN) was calculated as the sum of NH 4 + , NO 2 − and NO 3 − . For determining particulate organic carbon and nitrogen concentration, seawater samples were filtered through precombusted 25 mm Whatman GF/F filters. The filters were frozen at −20 °C. Before the analysis, the filters were treated with 1 N HCl to remove carbonates (Lorrain et al. 2003) and then oven-dried at 60 °C for 1 h. A CHNO-S elemental analyzer Costech mod. ECS 4010 was used to measure carbon and nitrogen according to Pella and Colombo (1973) and Sharp (1974). The analysis was performed by combustion of the sample in oxygen excess at 980 °C with Cr 2 O 3 catalyst and reduction in a column filled with reduced copper wires at 650 °C. CO 2 and N 2 were separated in a HayeSepQ divinylbenzene packed column at 70 °C, using helium as carrier gas, and measured with a thermal conductivity detector. Filter blanks were analyzed and subtracted. The C and N sample concentrations were expressed as µg C L −1 and µg N L −1 , respectively.
Samples for prokaryotic abundance (heterotrophs-HP, Synechococcus-SYN) were fixed with 0.22 µm filtered buffered formaldehyde solution (2% final concentration) and flash frozen in liquid nitrogen according to Patel et al.  (2007). Samples were processed according to Noble and Fuhrman (1998) within 1 month. Samples (0.5-1.5 mL) were filtered onto 0.02 µm pore-size Al 2 O 3 inorganic membrane filters (Anodisc, Whatman), then stained with SYBR Green I 50× (final concentration), mounted on a glass slide with antifade solution (50% glycerol, 49% PBS and 1% ascorbic acid) and stored at −20 °C. The enumeration was done by epifluorescence microscopy (Olympus BX 60 F5) at 1000× magnification under a blue (BP 450-480 nm, BA 515 nm) and a green (BP 480-550 nm, BA 590 nm) (Waterbury et al. 1979) filters set for HP and SYN, respectively. A minimum of 300 cells was counted for each filter in at least 20 randomly selected fields. All samples were run in three replicates and values represent their mean ± a standard deviation <10%.
Prokaryotic carbon production (PCP) was measured by the incorporation of 3 H-leucine (Leu) (Kirchman et al. 1985). Triplicate aliquots (1.7 mL) and two controls killed by the addition of 90 µL 100% trichloroacetic acid (TCA) were amended with a 20 nmol L −1 radiotracer and incubated at in situ temperature in the dark. Incubations were stopped with TCA (5% final concentration) after 1 h. The extraction with 5% TCA and 80% ethanol was carried out using the microcentrifugation method (Smith and Azam 1992). Activity in the samples was determined by a β-counter (TRI-CARB 2900 TR Liquid Scintillation Analyzer) after the addition of 1 mL scintillation cocktail (Ultima GoldTM MV; Packard). Carbon biomass production was estimated using the conversion factor 3.1 kg C mol −1 Leu assuming a twofold isotope dilution (Simon and Azam 1989). The standard deviation among replicates was <10% of the average. Extracellular enzymatic activities were measured following the fluorogenic substrate analogues method (Hoppe 1993) using 7-amino-4-methyl-coumarin (AMC) and 4-methyl-umbelliferone (MUF). Leucine aminopeptidase, β-glucosidase and alkaline phosphatase were assayed as the hydrolysis rate of leucine-AMC, MUF-β-d-glucoside, and MUF-phosphate, respectively. Enzyme activities were expressed in terms of the rate of MUF or AMC production. After evaluation of the saturating concentrations, hydrolysis rates were measured by incubating 2.5 mL subsamples with 200 µmol L −1 MUF-β-d-glucoside, leucine-AMC, and 50 µmol L −1 MUF-phosphate for 1 h in the dark at in situ temperature. The fluorescence increase due to AMC and MUF hydrolyzed from the model substrates was measured using a Shimadzu RF-1501 spectrofluorometer (AMC = 380 nm excitation and 440 nm emission: MUF = 365 nm excitation and 455 nm emission). Triplicate blanks without fluorogenic substrate were used to determine the natural fluorescence increase in the samples not attributable to the tested enzymes. Filtered seawater (0.2 µm) collected at every environmental pH was used to produce calibration curves (in three replicates) by adding 5 µmol L −1 standard solutions of MUF and AMC. All samples were run in three replicates and values represent their mean ± a standard deviation <10%.

DNA extraction, amplicon library preparation, and sequencing
For each sample, 1 L of seawater was filtered onto a 0.22 µm Supor membrane on a sterile filtration column, filters were flash frozen in liquid nitrogen and once in the laboratory, stored at −20 °C. DNA was extracted using the Dneasy Power Water kit (Qiagen) following the manufacturer's instructions. Extracted DNA was quantified by a Qubit Fluorometer (Thermo Fisher Scientific) and stored at −80 °C until further analysis. DNA was PCR amplified by using the primer pair 515F-Y (5ʹ-GTG YCA GCMGCC GCG GTAA-3ʹ) and 926R (5ʹ-CCG YCA ATTYMTTT RAG TTT-3ʹ), which encompass the V4 and V5 hypervariable loops of 16S rRNA genes (Parada et al. 2016). Libraries were prepared following the 16S Metagenomic Sequencing Library Preparation protocol (Illumina 2013) and run on an Illumina NovaSeq System for a read length of 2 × 250 bp at the genetic and epigenetic ARGO Open Lab Platform, Area Science Park, Trieste, Italy.

Bioinformatic pipeline
Bioinformatic analyses were carried out with QIIME2 2020.6 (Bolyen et al. 2019). Given the high number of reads produced by the NovaSeq (i.e., 46,985,351) we applied a 30% subsampling threshold prior running the bioinformatic pipeline in order to decrease the computational demand. The inspection of the rarefaction curves (Supplementary Figure  S1) ensured that the sequencing effort was more than enough to assess the microbial diversity. Based on visual inspection of quality profiles, forward and reverse raw sequences were then trimmed at 240 bp and 220 bp, respectively. Reads were then filtered, denoised, and used to infer amplicon sequence variants (ASVs) using the dada2 qiime2 plugin (Callahan et al. 2016). Taxonomy was assigned using the sklearn I Bayes taxonomy classifier (Bokulich et al. 2018) against the Silva 138.1 99% reference database (Quast et al. 2013). Reads belonging to Eukarya, mitochondria, chloroplasts, and ASVs with frequency <2 (singletons) were removed.

Statistical analysis
We used the R environment (v. 3.6.2, R Core Team 2019) for statistical analyses and data visualization. A principal component analysis (PCA) was performed using all the physical, chemical and microbiological data by using the FactoMineR 4.2.2 package in R. Because the variables used had different units, we scaled them to mean = 0 and variance = 1 before performing the analysis. The PCA results were extracted and visualized using the factoextra 4.2.2 R package (Le et al. 2008).
To further investigate the environmental drivers of the communities, we used a distance-based redundancy analysis (dbRDA) based on a Bray-Curtis distance matrix calculated on relative abundances using the capscale function in the R package vegan 2.6-2 (Oksanen et al. 2019). The analysis was applied only to surface and bottom data, i.e., where samples for DNA metabarcoding were collected. The environmental variables salinity, depth, pH, Chla, dissolved oxygen, DIC, POC, PN, temperature, and nutrients were first checked for collinearity inspecting pairwise scatterplots and correlation coefficients, according to the approached proposed by Zuur et al. (2010). The final set of environmental variables used in the dbRDA was then selected taking into account the multicollinearity and their ecological relevance. The significance of the selected dbRDA model was assessed by hypothesis testing (i.e., ANOVA).

Water column features
The area was characterized by a thermocline between 25 and 50 m depth separating warmer surface waters (>25.6 °C) from colder subsurface ones (14.3 to 16.1 °C below 50 m) ( Fig. 2a). High salinity (>37.9) surface and deep layers were separated by a lower salinity (37.5-37.8) layer between 20 and 70 m depth (Fig. 2b). CTD-derived chlorophyll a concentration profiles (Fig. 2c) highlighted the presence of a deep chlorophyll maximum from 40 to approximately 150 m (data not shown), displaying the maximum values around 70-90 m depth. These features affected the dissolved oxygen concentration, whose spatial distribution showed lower levels (195-200 μmol L −1 ) in surface waters in the upper 20 m, below which an oxygen-rich layer was present in correspondence of lower salinity and higher phytoplankton biomass, reaching the highest values of 250-260 μmol L −1 between 35 and 55 m.
The higher phytoplankton biomass in deep waters was also confirmed by the analyses of discrete samples for chlorophyll a concentration that showed the highest values (0.3 μg L −1 ) below 70 m (Fig. 3a). Along the transect, a drop in pH T (Fig. 3b) was detected in the vent-influenced stations V31, V35, and V23. The minimum value was 7.31, measured at 20 m depth at station V35. It is also noteworthy that the acidification effect of the vent plume reached the surface (pH T = 7.62) at this station, likely due to the shallow depth. Because of dilution and mixing of the vent fluid with seawater, the pH T at the third vent-associated station (V31) increased to ambient seawater levels (8.00) at middepth (54 m).
Dissolved inorganic carbon (DIC, Fig. 3c) concentrations were higher (up to 39.76 mg L −1 ) in bottom waters, especially in the vent-associated station V35. Station P19 had the lowest DIC concentration in the intermediate layer (30.11 mg L −1 ), along with maximum pH values. The concentration of inorganic nutrients (Fig. 3d, e) indicated a general oligotrophic condition. Dissolved inorganic nitrogen concentrations were highly heterogeneous along the transect (DIN, Fig. 3d) ranging from 0.7 to 2.4 µmol L −1 . Phosphate concentrations were in general between 0.01 and 0.07 µmol L −1 , but at station V35, in proximity to one of the vents, the highest value of 1.28 µmol L −1 was measured both at the surface and at 20 m depth. Particulate nitrogen (PN, Fig. 3f) and particulate organic carbon (POC, Fig. 3g) showed comparable distributions along the transect with higher values in the bottom layer than at the surface. The highest values of PN and POC were measured at station P15. POC concentration was higher than PN with a range of 30-80 µg C L −1 vs 4-14 µg N L −1 , respectively. All data generated at station P0, not represented in the contour plots, are reported in Supplementary Table S2.

Microbial dynamics in water column
HP (Fig. 4a) ranged from 0.2 to 1.5 × 10 9 cells L −1 , with the maximum value recorded at station P25 and the lowest abundance at station P15. Synechococcus was more abundant at the surface (Fig. 4b), ranging from 0.3 to 2.1 × 10 7 cells L −1 .
To assess the vent effect on microbial organic matter degradation, we measured the activity (Fig. 5) of leucineaminopeptidase (AMA), alkaline phosphatase (APA), and  Fig. 5a). Alkaline phosphatase (Fig. 5b) rates were generally low (<10 nM h −1 ) in the area, with a spike (20 nM h −1 ) of activity in venting station V35 where the highest phosphate and DIC concentrations were measured. A similar pattern was evident for BGLU (Fig. 5c), with faster rates (>3 nM h −1 ) in the surface layer of vent-affected stations (V31, V35, and V23).
Prokaryotic carbon production rates (PCP, Fig. 5d) mirrored the chlorophyll a concentration with higher values in the bottom samples collected in the southern sector of the transect (Max = 75.4 ng C L −1 h −1 ; Fig. 5d).

Relationship between environmental and biological variables
The principal component analysis (Fig. 6) of physical, chemical, and microbiological data showed a separation of samples according to depth (and the depth-related variables temperature, Chla, DO, POC, PN and PCP) along the first principal component, accounting for 37.1% of the total variance. APA, DIC, Phosphate, and pH were correlated with PC2 (Fig. 6) and explained 19.4% of the total variance.
Overall, the variables correlated with PC1 were affected by the physical structure of the water column, whereas the variables correlated with PC2 were affected by the vent emissions.

Prokaryotic diversity and community composition
16S rRNA gene analysis was performed at six stations along the transect and at the reference station P0 (Supplementary Table S1).
The relative abundance of the top 15 genera accounted for ~75% and ~66% of all reads in surface and bottom samples, respectively (Fig. 7). At the surface, the prokaryotic assemblage was dominated by SAR11 clade Ia, SAR116 and Synechococcus (average relative abundance = 18.7%, 12.4% and 12.2%, respectively). At the bottom, the highest relative abundances were attributed to SAR11 clade Ia, Prochlorococcus and Thiomicrorhabdus (on average 10.3%, 9.9% and 7.6% respectively), Thiomicrorhabdus, a sulfur oxidizer member of the Gammaproteobacteria was predominant especially at the vent-associated stations with  the highest relative abundance in V35 (24.53%) followed by station V23 (18.8%). Overall, the prokaryotic community structure of the surface differed from the bottom one, where the presence of sulfur oxidizers was typical of vent-affected stations.
The environmental data exploration highlighted collinearity between the variable depth and Chla, PN, POC, temperature, DO, and between DIC and pH. Therefore, the dbRDA analysis was performed using the variables depth, DIC, DIN and PO 4 (Fig. 8). The ANOVA on the dbRDA model highlighted depth (df: 1, SSE: 1.73810, F: 52.9201, p < 0.001) and DIC (df: 1, SSE: 0.12638, F: 3.8479, p < 0.05) as the main drivers shaping the community structure. In fact, also in this case the samples followed a main gradient (along RDA1) related to depth and depth-connected drivers with vent-influences samples further segregating from the other bottom samples (RDA2).

Effect of vents on environmental features and microbial dynamics
Investigating the functioning of microbial communities in shallow hydrothermal vents is important to understand how microbes respond to the rapid changes in nutrient availability, pH, electron donors and in general to the chemistry of such environments. The results are discussed in the light of shallow hydrothermal vent systems and ocean acidification literature, as the main vent feature we detected at ventaffected stations was an increase in DIC (Rizzo et al. 2022 and references therein). The resulting pH variations could impact microbial responses in terms of metabolic activities and microbial community structure and diversity.
In the study area (Panarea Island), deep-origin CO 2 leaks from the sea bottom into the water column (Esposito et al. 2006). The composition of the gas is stable, with 96.6-98.6% of CO 2 and 1.2-2.2% of H 2 S (Caliro et al. 2004). Thus, the vents in Panarea, as in other natural vent systems, have an impact on the chemistry of the water column, mostly affecting its pH (Dittmar and Koch 2006;Lang et al. 2006;Walker et al. 2008;Bennett et al. 2011;Gomez-Saez et al. 2016). Our results showed a drop in pH at the vent stations (V35, V23 and V31) due to elevated DIC concentrations. In addition, the area was characterized by low nutrient and chlorophyll concentrations, typical of oligotrophic Mediterranean waters (Bergamin et al. 2019). Only phosphate concentration showed a remarkable increase in the vent station V35, probably due to the turbulence induced by the seepage of gas bubbles that transported the phosphate-enriched interstitial water upwards (Molari et al. 2018). Chla, POC, and PON concentrations did not show evident differences among the venting and non-venting stations, as previously shown (Acosta Pomar and Giuffrè 1996).
Moreover, the thermocline formed a two-layered euphotic zone with the upper layer being characterized by higher temperature and the lower one characterized by higher Chla concentration (i.e., a deep chlorophyll maximum, DCM). In such a system, the thermocline had an important effect on structuring the environmental and biological parameters. For example, both prokaryotic abundance and PCP rates were higher in the DCM, thus highlighting the tight coupling between carbon fixation and prokaryotic carbon production processes (Longhurst and Harrison 1989). Furthermore, Chla and POC (and PN) were positively correlated (Spearman's ρ = 0.71, df = 22, p < 0.001), thus reflecting the increase of phytoplankton carbon biomass in the DCM (Macias et al. 2019). In contrast, APA and BGLU rates were more affected by the presence of the vent rather than by the vertical structure of the water column (Figs 5, 6).
We did not find any apparent effect of the vent on prokaryotic abundance and production. Many ocean acidification studies found that prokaryotic biomass was weakly influenced by acidification either increasing or decreasing in cell abundance (Grossart et al. 2006;Allgaier et al. 2008) although its photoautotrophic fraction (Cyanobacteria) may benefit from additional CO 2 (Hutchins et al. 2007) in fueling photosynthesis. Prokaryotic production has been found to be either stimulated (Grossart et al. 2006), inhibited (Coffin et al. 2004), or unaffected (Allgaier et al. 2008), during ocean acidification mesocosm experiments. In our study, we did not observe any increase or decrease in prokaryotic production associated with the vent, yet a dependance of organic matter utilization on phytoplankton biomass was evidenced (Spearman's ρ between PCP and Chla = 0.51, df = 22, p = 0.05). The relationship between environmental factors and the extracellular enzymatic activities in extreme environments, including hydrothermal vent systems, has not been adequately explored (Jacobson Meyers et al. 2014). It is known that enzyme activity depends on several conditions, such as pH, temperature, catalyst availability, concentration and structure of substrates (Fukuda et al. 2000;Pomeroy and Wiebe 2001;Arnosti 2011;Celussi et al. 2017). The vent effect was detected in terms of enhanced hydrolysis rates of alkaline phosphatase and beta-glucosidase extracellular activities (Fig. 5). The faster alkaline phosphatase rates detected in our study are most likely to be caused by P-rich compounds ejected by the vents (Molari et al. 2018). On the contrary higher beta-glucosidase activity values at the lowest pH are possibly related to high DIC concentration, as previously observed in a mesocosms experiment in the Mediterranean Sea (Celussi et al. 2017). pH plays an important role in determining the redox state of metal cofactor, the stability of enzymes and their performance, the organic matter substrate and structure, ultimately affecting the measured exoenzymatic rates (Gottschalk 1979). However, the direct effect remains difficult to assess because it is hard to separate the effect of hydrogen ion concentrations from other synchronizing factors, such as temperature, nutrient concentration, or other anthropogenic influences (Piontek et al. 2010(Piontek et al. , 2013. Our results showed that vents had no clear effect on AMA (Fig. 5), in contrast with the studies of Yamada and Suzumura (2010), and Grossart et al. (2006), where protease rates were found to be highly sensitive to acidification conditions, either suppressed or enhanced by decreasing pH. However, when comparing the results of this study with literature findings it has to be taken into account that most acidification studies are conducted using artificial CO 2 enrichments, often involving microbial communities not adapted to chronic exposure to low pH conditions.
In order to deepen our understanding of the consequences of ocean acidification by atmospheric CO 2 dissolution and by hydrothermal fluid enrichment, it will be important to examine specifically the influence of chemical processes on the activation and regulations of genes involved in extracellular enzymatic activities as well as in membrane transporters responsible for organic matter utilization.

Prokaryotic community structure in the Panarea Island hydrothermal system
Our results show that prokaryotic assemblages were mainly shaped by depth and depth-related drivers (Fig. 8), as already found in other shallow Mediterranean systems under summer stratification (e.g., Tinta et al. 2015). Interestingly, prokaryotic assemblages within the two separate layers were very homogenous (horizontally) not only along the transect but also compared to the 'reference' station P0, located 10 km north of the island. However, in our data set, the vent did affect the structure of microbial communities only at their proximity, possibly due to the rapid dispersion of the vent fluids by the hydrodynamic regime.
The water column community was dominated by Proteobacteria and Cyanobacteria (Fig. 7), which are recognized as major players in many marine systems (Cottrell and Kirchman 2000;Rappé et al. 2000;Morris et al. 2002). In particular, the most represented taxa in the upper warmer and low-chlorophyll layer were SAR11 and Synechococcus. SAR11 clades are made up of various subgroups with different environmental preferences (Carlson et al. 2009;Morris et al. 2012), living in specific niches. In our study, we found a prevalence of clade Ia over clade II, as often evidenced in coastal systems (Auladell et al. 2022). The high abundance of SAR11 could also be explained by the metabolic effect of proteorhodopsin (light dependent proton pump) that provide an energetic boost in organic-poor environments (Giovannoni et al. 2005). Synechococcus is known to thrive in warm surface oligotrophic waters where its accessory pigments (phycoerythrin; Everroad and Wood 2012) provide an adaptive advantage over Prochlorococcus whereas its high surface to volume ratio (Moutin et al. 2002) outcompetes phototrophic eukaryotes in nutrient uptake.
Besides a consistent relative contribution of SAR11 (ca. 10%), the lower, cooler and high(er)-chlorophyll layer was dominated by Prochlorococcus and genera belonging to the families Flavobacteriaceae (NS2, NS4, NS5, NS9) and Rhodobacteraceae (formerly known as the Roseobacter group). The small cell size of Prochlorococcus leads to a high surface-to-volume ratio that facilitates efficient nutrient acquisition and enhances light absorption, which, when combined with its unique pigmentation and broad temperature optima for growth, make it the most efficient light absorber of any photosynthetic cell, thus adapted to subsurface ocean layers (Biller et al. 2015). Both the Flavobacteriaceae and the Rhodobacteracea are often found in association with high microalgal biomass, being capable of effectively degrading phytoplankton-derived organic matter, spanning from sugars and proteins to low molecular weight compounds such as DMSP and taurine (Moran et al. 2007;Teeling et al. 2012).
Within this 'typical' coastal Mediterranean picture (Auladell et al. 2022), the venting stations were characterized by high relative abundance of the chemoautotrophic bacterial genus Thiomicrohabdus (formerly known as Thimicrospira), specifically where DIC concentration was high (bottom of stations V35 and V23, Fig. 3c). This genus is known to be a culturable sulfide-oxidizer, dominant at shallow hydrothermal vents and extreme environments (Tarasov et al. 2005;Scott et al. 2019). It uses reduced sulfur compounds (e.g., thiosulfate, tetrathionate, sulfur and sulfide) as electron donors and obtains carbon from CO 2 fixation. Cultivation-based studies revealed that sulfur compounds are essential substrates for this bacterium in hydrothermal vent systems (Jannasch et al. 1985). Culture-independent methods revealed that Thiomicrorhabdus are predominant members in hydrothermal vent communities, such as sulfide chimneys and the Lost City carbonate chimney in the Southwest Indian Ridge (Cao et al. 2014) and the Mid-Atlantic Ridge site (Wirsen et al. 1998) and in extreme environments (Scott et al. 2019). In our study, Thiomicrorhabdus was dominant in bottom samples of vent-affected stations, in agreement with previous surveys in the Panarea Island area (Manini et al. 2008;Maugeri et al. 2009;Karuza et al. 2012;Lentini et al. 2014) and in other shallow hydrothermal vents in the Mediterranean Sea (Sievert et al. 1999;Patwardhan et al. 2018Patwardhan et al. , 2021. Noteworthy, the relative abundance of Thiomicrorhabdus decreased with increasing distances from the emission points, suggesting a very rapid dilution effect.
In this natural laboratory, future studies should focus on the benthic-pelagic coupling at the vent sites to unravel the magnitude and the extent of the vent's influence on the C and S cycles. Furthermore, resolving the plume dynamics at a higher spatial resolution would allow to depict the fate of the vent-associated microbiome and its biogeochemical consequences in the oligotrophic Mediterranean. When considering natural laboratories, in order to dissect the effect of atmospheric-driven from hydrothermal-driven acidification on marine microbes and their role in biogeochemical cycles it will be important to compare vents characterized by CO 2 only and by CO 2 + H 2 S emissions. These findings could contribute to a better understanding of the role of microbes in the future ocean.

Acknowledgements
The project's ship-time was funded by Eurofleets, with work being performed as a sub-project of the Seventh Framework Programme projects "ECO2" and "RISCS", which provided subsidiary funding. The help of the captain and of the crew of the R/V Urania is kindly acknowledged. The authors would like to thank M. Kralj for nutrient analyses, C. Comici for CHN data and A. Karuza for microscopy counts. The comments of two anonymous reviewers helped improving the original version of the manuscript.
Author contributions MC, MG and CDV coordinated the fieldwork. MC, MG, and CDV provided resources. MC, FM, MG, AS, EB, VF, VM and CDV performed data curation, sequencing work and lab work. AS, EB, VF and VM conducted the bioinformatics analysis. AS, FM and MC wrote the manuscript with the help and inputs from all the authors.
Funding Open access funding provided by Istituto Nazionale di Oceanografia e di Geofisica Sperimentale within the CRUI-CARE Agreement.

Conflicts 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.
Ethical approval No ethical approval or permit for animal experimentation was required.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.