# Microbial food web components, bulk metabolism, and single-cell physiology of piconeuston in surface microlayers of high-altitude lakes

Hugo Sarmento, Emilio Casamayor, J.-C. Auguet, Maria Vila-Costa, Marisol Felip, Lluís Camarero, Josep Gasol

► To cite this version:

Hugo Sarmento, Emilio Casamayor, J.-C. Auguet, Maria Vila-Costa, Marisol Felip, et al.. Microbial food web components, bulk metabolism, and single-cell physiology of piconeuston in surface microlayers of high-altitude lakes. Frontiers in Microbiology, 2015, 6, pp.361. 10.3389/fmicb.2015.00361. hal-01841336

HAL Id: hal-01841336

https://hal.umontpellier.fr/hal-01841336

Submitted on 8 Jun 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Microbial food web components, bulk metabolism, and single-cell physiology of piconeuston in surface microlayers of high-altitude lakes

OPEN ACCESS

**Edited by:**

Jérôme Comte,
Laval University, Canada

**Reviewed by:**

Michael Cunliffe,
Marine Biological Association of the
UK, UK
Craig E. Nelson,
University of Hawaii at Manoa, USA

**\*Correspondence:**

Hugo Sarmento,
Laboratory of Microbial Processes and Biodiversity, Department of
Hydrobiology, Federal University of
São Carlos, Rodovia Washington Luiz,
Sao Carlos, SP, Brazil
hsarmento@ufscar.br

†Present Address:

Jean-Christophe Auguet,
Equipe Environnement et
Microbiologie, Université de Pau et
des Pays de l'Adour, Pau, France;
Maria Vila-Costa,
Department of Environmental
Chemistry, Instituto de Diagnóstico
Ambiental y Estudios del Agua,
Consejo Superior de Investigaciones
Científicas, University of Barcelona,
Barcelona, Spain

**Specialty section:**

This article was submitted to
Aquatic Microbiology,
a section of the journal
Frontiers in Microbiology

**Received:** 28 January 2015

**Accepted:** 09 April 2015

**Published:** 05 May 2015

**Citation:**

Sarmento H, Casamayor EO, Auguet
J-C, Vila-Costa M, Felip M, Camarero
L and Gasol JM (2015) Microbial food
web components, bulk metabolism,
and single-cell physiology of
piconeuston in surface microlayers of
high-altitude lakes.
Front. Microbiol. 6:361.
doi: 10.3389/fmicb.2015.00361

Hugo Sarmento<fnref n="1" />,²*, Emilio O. Casamayor<fnref n="3" />, Jean-Christophe Auguet<fnref n="3" />†,
Maria Vila-Costa<fnref n="3" />†, Marisol Felip⁴,⁵, Lluís Camarero<fnref n="3" />,⁵ and Josep M. Gasol<fnref n="1" />

Sharp boundaries in the physical environment are usually associated with abrupt shifts in organism abundance, activity, and diversity. Aquatic surface microlayers (SML) form a steep gradient between two contrasted environments, the atmosphere and surface waters, where they regulate the gas exchange between both environments. They usually harbor an abundant and active microbial life: the neuston. Few ecosystems are subjected to such a high UVR regime as high altitude lakes during summer. Here, we measured bulk estimates of heterotrophic activity, community structure and single-cell physiological properties by flow cytometry in 19 high-altitude remote Pyrenean lakes and compared the biological processes in the SML with those in the underlying surface waters. Phototrophic picoplankton (PPP) populations, were generally present in high abundances and in those lakes containing PPP populations with phycoerythrin (PE), total PPP abundance was higher at the SML. Heterotrophic nanoflagellates (HNF) were also more abundant in the SML. Bacteria in the SML had lower leucine incorporation rates, lower percentages of "live" cells, and higher numbers of highly-respiring cells, likely resulting in a lower growth efficiency. No simple and direct linear relationships could be found between microbial abundances or activities and environmental variables, but factor analysis revealed that, despite their physical proximity, microbial life in SML and underlying waters was governed by different and independent processes. Overall, we demonstrate that piconeuston in high altitude lakes has specific features different from those of the picoplankton, and that they are highly affected by potential stressful environmental factors, such as high UVR radiation.

**Keywords:** picoplankton, leucine incorporation, bacterioplankton, flagellates, flow cytometry, NADS, CTC, ultraviolet radiation

## Introduction

The upper micrometers of the water column of aquatic ecosystems form a substantially different environment than that of the underlying surface water, despite its close proximity as reviewed by

<footnote n="1">¹Institut de Ciències del Mar, ICM-Consejo Superior de Investigaciones Científicas, Barcelona, Spain, ²Laboratory of Microbial Processes and Biodiversity, Department of Hydrobiology, Federal University of São Carlos, São Carlos, Brazil,</footnote>

<footnote n="3">³Integrative Freshwater Ecology Group, Limnological Observatory of the Pyrenees, Centre d'Estudis Avançats de Blanes, CEAB-Consejo Superior de Investigaciones Científicas, Blanes, Spain, ⁴Department of Environmental Chemistry, IDAEA-CSIC, Barcelona, Spain, ⁵Centre de Recerca d'Alta Muntanya, Universitat de Barcelona, Lleida, Spain</footnote>
Cunliffe et al. (2011). The group of organisms living at the water-air interface, the surface microlayer (SML), was named neuston (Naumann, 1917), and the smallest organisms (<2 µm) living there can be defined as piconeuston. The SML is the boundary between two totally different environments and the site where fundamental biological and chemical processes, such as gas exchanges between atmosphere and water occur (e.g., Conrad and Seiler, 1988). The SML is known to concentrate substances, which may enhance or inhibit microbial activities. Previous studies showed that the marine SML was enriched with amino acids and bacteria and, in most instances, virus-like particles (Kuznetsova et al., 2004), and that heterotrophic piconeuston had enhanced extracellular enzymatic peptide hydrolysis rates (Kuznetsova and Lee, 2001), compared to subsurface waters.

However, and at least in marine/estuarine environments, early reports indicated that heterotrophic piconeuston has lower metabolic activity than surface picoplankton, (Dietz et al., 1976). This difference in metabolism was explained as a community response to different environmental stress factors occurring preferentially in the SML, such as high incident solar UVR, and more variable water temperature and salinity, as well as potentially higher concentrations of toxic substances (Maki, 1993). In aquatic environments naturally subjected to high UVR doses, such as high altitude lakes in summer, a difference in activity should therefore be expected.

The differences in metabolism, however, might also be attributed to differences in community compositions between the SML and the surface. Auguet and Casamayor (2008), using catalyzed reporter deposition fluorescence *in situ* hybridization (CARD-FISH) and 16S rRNA gene sequence analysis, reported a higher abundance of archaeal communities in the SML of Pyrenean oligotrophic high-mountain lakes as compared to surface communities. These Archaea populations were composed mainly of Crenarchaeota, whereas surface populations were mainly comprised of Euryarchaeota. Similarly, Vila-Costa et al. (2014) found distinctive populations of both archaea and bacteria inhabiting SML and surface waters of the same Pyrenean lakes using 454 pyrosequencing, and the differences were exacerbated under atmospheric loadings that stimulated microbial activities. A less clear pattern was observed in a set of six Alpine lakes located across an altitude gradient (Hörtnagl et al., 2010), where Betaproteobacteria (enumerated by CARD-FISH) dominated in both SML and underlying water, and the differences observed among lakes were attributed to lake-specific intrinsic factors.

Living in the SML is rather challenging, mainly due to the harsh prevailing conditions resulting from summer extreme UVR levels (Sommaruga, 2001). Previous reports indicating that UVR negatively affects bacterial activity (i.e., Ruiz-González et al., 2013), HNF growth, and bacterial consumption rates (Sommaruga et al., 1999) suggest that microorganisms living in the neuston should experience heavy environmental stress. Independently of the peculiarity of SML's prokaryotic taxonomic composition described in the studies cited above, there is little information available on the *in situ* microbial food web structure (i.e., both abundance and composition of heterotrophic prokaryotes, phototrophic picoplankton [PPP], and heterotrophic nanoflagellates [HNF]) and of bacterial

single-cell activity and physiology, which could illustrate the ecological processes shaping life in the SML.

The aim of this study was to study microbial community structure, metabolism, and physiology of piconeuston of SML compared to underlying water in high mountain lakes. Our working hypothesis is that microbial communities living in the SML of high altitude lakes are subjected to environmental harshness that affects their composition, community structure, activity, and physiology in a different way than that of surface waters communities. In order to achieve this goal, we carried out a comprehensive flow-cytometry measurement of (i) microbial community structure, (ii) prokaryotic bulk and single-cell activity, and (iii) physiological status in 19 remote high altitude lakes sampled under summer high solar radiation conditions, in order to determine the variability of these parameters in the SML as compared to surface waters. To the best of our knowledge, most of the variables studied, such as detailed microbial community structure by flow cytometry and bacterial single-cell activity, had never before been measured in the SML.

# Materials and Methods

## Sampling Sites and Limnological Parameters

A set of 19 high mountain lakes from the Central Pyrenees were sampled from 17th to 24th, June 2008 at 3 depths: in the first ~400 µm of the water column, here defined as the SML; at 0.5 m depth—which we label as “surface”; and at the depth equivalent to 1.5-fold Secchi disk value, usually corresponding to the depth of the summer deep chlorophyll maximum (DCM) (Catalan et al., 2006), which ranged from 2 to 30 m depth, depending on the lake. In this report the DCM values of Chlorophyll-a (Chla) were only used to characterize the lakes according to their nutrient and trophic status. In these clear water mountain lakes Chla at the surface does not reflect the trophic status of the lake because most primary production is located at the DCM (Catalan et al., 2006). The lakes were selected in order to maximize variability in chemical and morphological characteristics and were accessed by foot as they are located in uninhabited remote locations.

SML samples were collected from the upper ~400 µm water with a nylon screen sampler (Agogué et al., 2004, 2005) near the deepest point of each lake. Surface (0.5 m depth) and deeper samples were taken using a 3-litre sampler (either Ruttner or Patalas bottles). Samples were pre-screened through a 40 µm pore-size net to remove large plankton components.

Water transparency was measured with a Secchi disk. Temperature profiles were measured *in situ* using a PT10 type thermistor. In the lab, pH was measured with an Orion instrument equipped with a probe for low ionic strength samples (Crison). Conductivity was measured with conductivimeter of the brand WTW. Total phosphorus (TP) was extracted with an oxidative persulfate digestion, and analyzed by spectrophotometry using the malachite green method (Camarero, 1994). Total nitrogen (TN) was extracted with an oxidative persulfate digestion in the autoclave, and measured by UV spectrometry. Dissolved organic carbon (DOC) was assessed on Whatman GF/F filtered lake-water with a Shimadzu TOC5000 analyzer. Chla concentrations were measured in 90% acetone extracts (1 L sample filtered on Whatman GF/F filter and
extracted by 3 min sonication in 4 mL 90% acetone) according to
Jeffrey and Humphrey (1975).

# Abundances of Microbial Food Web Components by Flow Cytometry

## Sample Collection

Microbial food web components (cell abundance of heterotrophic prokaryotes, PPP, and HNF) were quantified by flow cytometry. Three subsamples were taken for separate counts of heterotrophic prokaryotes, PPP, and heterotrophic flagellates. The last two samples were kept in water baths with *in situ* water and analyzed as soon as possible (3 h on average after collection). For heterotrophic prokaryotes counts, 4 ml of lake water were collected from each depth, fixed immediately with cold glutaraldehyde 10% (final concentration 1%), left in the dark for 10 min at *in situ* temperature, placed in liquid nitrogen and then stored at -80°C. Samples were analyzed 2 months after sampling. An aliquot of these frozen samples were used to quantify PPP and results were compared to the fresh analysis done *in situ* (see below).

## Heterotrophic Prokaryote Abundances

Stored samples were thawed at room temperature and 400 µl were mixed with DMSO-diluted SYTO-13 (Molecular Probes Inc., Eugene, OR, USA) at 2.5 µmol l⁻¹ final concentration. The mixture was left for about 10 min in the dark for complete staining and was run in a FACSCalibur (BectonDickinson) flow cytometer equipped with a 15 mW Argon-ion laser (488 nm emission). At least 30,000 events were acquired for each subsample. Heterotrophic prokaryotes were detected by their signature in a plot of side scatter (SSC) vs. FL1 (green fluorescence) following del Giorgio et al. (1996) as discussed in Gasol and del Giorgio (2000).

## Phototrophic Plankton Abundances

The fixed samples were thawed and run without addition of any stain in a FACSCalibur flow cytometer equipped with a blue laser (488 nm, 15 mW) and a red laser diode (~635 nm). Small algae were identified in plots of SSC vs. FL3 (red fluorescence of Chla), FL2 (orange fluorescence of phycoerithrine, PE) vs FL3 (Olson et al., 1993), and FL3 vs. FL4 (far red fluorescence after red-light excitation, indicative of phycocyanine, PC). Samples were run twice with different settings (different voltages for scatter and fluorescence), in order to maximize the “sampling window” and cover cells sizes below 1 µm to ca. 20 µm, thus including pico- and nanophototrophs. Data analysis was performed with the CellQuest software (BectonDickinson), and we combined the information collected with the two settings into one set of populations. Examples of cytograms with distinct populations are presented in Schiaffino et al. (2013) and Sarmento et al. (2008). While the data presented in this study corresponds to samples analyzed a few months after collection and from fixed replicates, we also analyzed fresh samples the same day of collection, ~3 h after sampling. Since we did not have the red laser for these live samples, the numbers presented are those of the fixed sampled. The fresh samples were used to confirm that no populations were negatively affected by fixation and storage.

## Small Heterotrophic Eukaryotes Abundances

In this study, “small heterotrophic eukaryotes” included pico and nano eukaryotes, and refers mainly to small HNF. Small eukaryote abundances were estimated following the protocol by Rose et al. (2004). From a stock solution of 1 mM Lysotracker Green (Molecular Probes), 1 µl was added to 99 µl of <0.2 µm MilliQ, and 3.8 µl of this diluted Lysotracker were added to 0.5 ml of the sample, ending at 75 nM Lysotracker final concentration. We analyzed the samples as in Rose et al. (2004), using a combination of light scatter (SSC) and green (FL1) and red (FL3) fluorescence. Samples were run alive, less than 3 h after sample collection, and at high (ca. 100 µl min⁻¹) speed. Concentrations were obtained from weight measurement of the volume analyzed.

## Bacterial Bulk and Single-Cell Activity

Bacterial bulk heterotrophic activity was estimated using the ³H-leucine incorporation method (Kirchman et al., 1985). Quadruplicate aliquots of 1.2 ml and 2 trichloroacetic acid (TCA) killed controls were taken immediately after sample collection. The samples were incubated with 40 nM ³H-leucine final concentration added immediately and for about 2 h in the dark in a water bath with lake water and at *in situ* temperature. The incorporation was stopped with the addition of 120 µl of cold TCA 50% to each replicate and the samples were kept frozen at -20°C until processing, which was carried out by the centrifugation method described by Smith and Azam (1992).

We used two physiological probes to test the metabolic and physiological single-cell status of prokaryotic microbes, respectively (i.e., Del Giorgio and Gasol, 2008). The abundance of respiring bacteria was determined using the probe 5-cyano-2,3-ditolyl tetrazolium chloride (CTC, Polysciences), an indicator of electron transport system respiratory activity. Highly respiring prokaryotes were considered those able to reduce CTC. Reduced CTC turns into a red fluorescent formazan that is detectable by epifluorescence and flow cytometry (Sherr et al., 1999; Sieracki et al., 1999). A stock solution of 50 mmol L⁻¹ CTC (Polysciences) was prepared daily, filtered through 0.1 µm, and kept in the dark at 4°C until use. The CTC stock solution was then added to 0.45 mL of sample (5 mmol L⁻¹ final CTC concentration) and incubated for 3 h at room temperature in the dark. The red fluorescence of CTC (FL3) and SSC were used to discriminate the CTC positive cells from other particles and an FL2-FL3 plot to exclude picoautotrophs. The percentage of CTC positive cells was calculated relative to the total bacterial counts obtained as mentioned above.

Cells with intact membranes were enumerated using the Nucleic-Acid-Double-Staining (NADS) viability protocol, which uses a combination of the cell-permanent nucleic acid strain SybrGreen I (SG1, Molecular Probes, Eugene, OR) and the cell-impermeant propidium iodine (PI, Sigma Chemical Co.) fluorescent probe. We used a 1:10 SG1 and 10 µg ml⁻¹ PI concentrations that were added to fresh samples, less than 2 h after sampling. After simultaneous addition of each stain, the samples were incubated for 20 min in the dark at room temperature and then analyzed by flow cytometry. SG1 and PI fluorescence were detected in the green (FL1) and orange-red (FL3) cytometric channels, respectively. A dot plot of red vs.
green fluorescence allowed distinction of the “live” cells (i.e., cells with intact membranes and DNA present) from the “dead” cells (i.e., with compromised membranes, (Grégori et al., 2001; Falcioni et al., 2008).

## Statistical Analyses

We computed a neuston enrichment factor as the ratio between abundance or activity in the SML divided by values measured at the surface to quantify the magnitude of the differences between layers. We tested for differences between layers using paired Student's t-tests. In order to explore possible relationships between environmental variables and microbial abundance or activity obtained for the different lakes and lake layers, we run a Factor Analysis (FA) using SPSS 20 Statistics and the stats package in R software (R Core Team, 2014). We separated physical and chemical variables and performed two different FA for each set of environmental variables. We included lake and catchment areas, altitude, maximum depth, surface temperature, and Secchi depth as physical variables. As for chemical variables, we included alkalinity, conductivity, DOC, total dissolved phosphorus (TDP), $NH_4^+$, $NO_2^-$, $NO_3^-$, $SO_4^{2-}$, $Ca^{2+}$-dissolved reactive silica (DRSi), Chla. We had measured other parameters (pH, soluble reactive phosphorus, total dissolved nitrogen, DIC, $Cl^-$, $Na^+$, $K^+$, $Mg^{2+}$) that were excluded from the FA because they were highly correlated or were a linear combination of some of the parameters included and violated the requirement of uniqueness. All variables were previously log transformed to normalize their distribution. The Kolmogorov-Smirnoff test was applied to check that there were no strong departures from normality. The abundance or activity variables at the surface and the SML, and the enrichment factor were tested for correlation against the scores obtained for the physical and chemical factors.

## Results

### In situ Physical and Chemical Conditions

The formation of a true SML requires stable conditions of wind and temperature for a certain amount of time. We sampled 19 lakes in a window of 9 days in which stable and sunny conditions occurred in the mountain. In order to detect the presence of a true SML, we first scanned our database, comparing temperature in the first millimeters of surface water (temperature profiles in Supplementary Material). We considered that there was no true SML in lakes in which water temperature was equal (at the 0.01°C resolution, the smallest reliable according to the resolution of the temperature probe used) in the lake skin and few millimeters below the surface, and we discarded these lakes from the analysis (indicated with “*” in Table 1), by considering that a true SML did not exist in those lakes. From here on, we refer only to the results obtained in the 16 selected lakes, which developed a true SML at the moment they were sampled.

Most of the lakes sampled were located at >2000 m above sea level and can be considered as high altitude lakes. This set of lakes includes a large range of maximal depths (from 5 to 113 m) and geochemical types. As an example, the pH among the 16 lakes ranged from 4.6 to 9.2, and conductivity ranged from 6 to 133 µS cm⁻¹ (Table 1). All the lakes were oligotrophic except

Muntanyó d'Àrreu, which was mesotrophic (Table 1). Average daily-integrated long-wave solar radiation during the sampling period (17th–24th July 2008) was 26.4 ± 8.3 MJ m⁻². The UV doses corresponding to this radiation were 1.4 ± 0.4 kJ m⁻² at 305 nm, 8.2 ± 2.5 kJ m⁻² at 320 nm, 15.3 ± 4.5 kJ m⁻² at 340 nm, and 20.7 ± 6.3 kJ m⁻² at 380 nm.

Chla (measured at the DCM) ranged from 0.2 µg l⁻¹ in Lake Certascan to 12.0 µg l⁻¹ in Lake Muntanyó d'Àrreu. Accordingly, water transparency ranged between 1.3 and 19.6 m of Secchi disk depth. Chla was correlated to total phosphorus (n = 17, *p* < 0.001, r² = 0.70) and nitrogen (n = 17, *p* < 0.01, r² = 0.58). Dissolved organic carbon (DOC) was relatively low, ranging from 0.2 to 1.4 mg l⁻¹ in Lake Certascan and Lake Plan, respectively.

### Abundances of Microbial Food Web Components
#### Phototrophic Components

Flow cytometry-determined abundances of PPP ranged two orders of magnitude from 2.5 × 10² to 1.8 × 10⁴ ml⁻¹, with slightly higher average abundances in the SML than in the surface waters, although not significant (Table 2). We observed a high diversity of populations in the cytograms obtained. The number of identified populations was on average 4.1 (range 1–7). In general we identified higher diversity of small phototrophic eukaryotes (red-fluorescing cells) than of phycoerythrin (PE)-rich (picocyanobacteria). None of the considered lakes had red-fluorescing (phycocyanin-rich) cyanobacteria. While all lakes had small phototrophic eukaryotes, some did not have cyanobacteria. In lakes with no populations of cyanobacteria, total PPP abundances were higher in the surface than in the SML. Lakes harboring PE-rich populations had higher total PPP abundance in the SML. Moreover, PE-rich (Synechococcus-like) cells were more abundant in the SML than in the surface, in those lakes where such populations were present (Figure 1). In most lakes, the number of different populations observed by flow cytometry (cytometric diversity) was higher at the surface than in the SML (Table 2). In a general comparison, the diversity of populations was higher at the surface but abundances were higher at the SML, although most comparisons were not significant or at the significance borderline (Table 2).

#### Heterotrophic Components

We observed significantly more heterotrophic prokaryotes and HNF in the SML than in surface waters (Table 2, Figure 2D). The total heterotrophic prokaryote:HNF ratio was relatively low (~480–490 bacteria per HNF, Table 2), and no significant differences were found between layers.

### Bulk and Single-Cell Prokaryote Heterotrophic Activities

We measured several indices of bacterial activity, some at the bulk level, such as leucine incorporation, and others at the single-cell level, such as abundance of bacterial cells presenting high levels of respiration (i.e., CTC positive). We used plots of surface against SML samples for the different parameters to discriminate
<table><caption>TABLE 1 | Location, general features, and some limnological parameters of the 21 samples from Pyrenean high-altitude lakes (DOC, dissolved organic carbon; TP, total phosphorous; TN, total nitrogen).</caption><thead><tr><th>Lake</th><th>Latitude</th><th>Longitude</th><th>Altitude (m asl)</th><th>Max. depth (m)</th><th>pH units</th><th>Conductivity (mS cm<sup>-1</sup>)</th><th>DOC (mg l<sup>-1</sup>)</th><th>TP (nM)</th><th>TN (µM)</th><th>Chl a (µg l<sup>-1</sup>)</th><th>Secchi depth (m)</th></tr></thead><tbody><tr><td>*Aixeus</td><td>42°36′40″</td><td>1°22′18″</td><td>2370</td><td>15.5</td><td>4.97</td><td>49.9</td><td>0.2</td><td>74</td><td>18</td><td>0.6</td><td>15.5</td></tr><tr><td>*Bassa de les Granotes</td><td>42°34′24″</td><td>0°58′16″</td><td>2330</td><td>5</td><td>6.45</td><td>9.8</td><td>2.3</td><td>119</td><td>26</td><td>4.6</td><td>2.4</td></tr><tr><td>Botornàs</td><td>42°35′35″</td><td>0°40′52″</td><td>2340</td><td>22</td><td>7.23</td><td>24.1</td><td>0.3</td><td>168</td><td>26</td><td>2.0</td><td>12.4</td></tr><tr><td>Certascan</td><td>42°42′40″</td><td>1°18′12″</td><td>2335</td><td>113</td><td>5.67</td><td>9.1</td><td>0.2</td><td>78</td><td>35</td><td>0.2</td><td>19.6</td></tr><tr><td>Filià</td><td>42°27′4″</td><td>0°57′12″</td><td>2140</td><td>5.5</td><td>7.79</td><td>133.3</td><td>not available</td><td>94</td><td>21</td><td>1.7</td><td>1.3</td></tr><tr><td>Gerber</td><td>42°37′50″</td><td>0°59′41″</td><td>2170</td><td>63</td><td>7.13</td><td>23.4</td><td>0.6</td><td>88</td><td>9</td><td>1.3</td><td>11.8</td></tr><tr><td>Ilbonet Perramó</td><td>42°38′34″</td><td>0°29′49″</td><td>2293</td><td>5</td><td>7.49</td><td>33.3</td><td>0.4</td><td>150</td><td>17</td><td>0.6</td><td>5.2</td></tr><tr><td>Illa</td><td>42°37′6″</td><td>0°59′37″</td><td>2452</td><td>18</td><td>6.68</td><td>13.3</td><td>0.4</td><td>70</td><td>7</td><td>1.3</td><td>12.0</td></tr><tr><td>Llauset</td><td>42°39′22″</td><td>0°41′11″</td><td>2190</td><td>90</td><td>7.52</td><td>58.7</td><td>not available</td><td>113</td><td>18</td><td>1.1</td><td>10.0</td></tr><tr><td>Llebreta</td><td>42°33′3″</td><td>0°53′25″</td><td>1620</td><td>11.5</td><td>7.55</td><td>30.6</td><td>1.1</td><td>203</td><td>not available</td><td>1.8</td><td>7.9</td></tr><tr><td>Long de Liat</td><td>42°48′24″</td><td>0°52′26″</td><td>2140</td><td>32</td><td>7.27</td><td>20.8</td><td>0.7</td><td>68</td><td>5</td><td>1.0</td><td>11.8</td></tr><tr><td>Muntanyó d'Àrreu</td><td>42°40′26″</td><td>1°00′29″</td><td>2210</td><td>29</td><td>9.18</td><td>74.3</td><td>0.8</td><td>674</td><td>9</td><td>12.0</td><td>10.1</td></tr><tr><td>Pica Palomèra</td><td>42°47′38″</td><td>0°52′8″</td><td>2308</td><td>10</td><td>4.61</td><td>29.6</td><td>not available</td><td>50</td><td>9</td><td>0.5</td><td>8.0</td></tr><tr><td>Plan</td><td>42°37′21″</td><td>0°55′51″</td><td>2188</td><td>11</td><td>7.04</td><td>16.7</td><td>1.4</td><td>63</td><td>7</td><td>0.6</td><td>9.0</td></tr><tr><td>Pòdo</td><td>42°36′11″</td><td>0°56′21″</td><td>2450</td><td>20</td><td>6.41</td><td>9.4</td><td>0.3</td><td>273</td><td>79</td><td>4.2</td><td>9.6</td></tr><tr><td>Pois</td><td>42°39′20″</td><td>0°42′27″</td><td>2056</td><td>19.5</td><td>8.19</td><td>67.2</td><td>0.2</td><td>82</td><td>14</td><td>2.9</td><td>8.6</td></tr><tr><td>*Redon</td><td>42°38′33″</td><td>0°46′13″</td><td>2240</td><td>73</td><td>6.67</td><td>10.5</td><td>0.4</td><td>n.a.</td><td>n.a.</td><td>1.5</td><td>n.a.</td></tr><tr><td>Roi</td><td>42°34′29″</td><td>0°48′12″</td><td>2310</td><td>10</td><td>7.06</td><td>25.8</td><td>0.5</td><td>126</td><td>18</td><td>5.2</td><td>5.9</td></tr><tr><td>Romedo de Dalt</td><td>42°42′22″</td><td>1°19′29″</td><td>2114</td><td>40</td><td>6.31</td><td>5.8</td><td>0.8</td><td>103</td><td>13</td><td>2.8</td><td>14.4</td></tr></tbody></table>

Chemical and Chla values measured at the depth of the deep chlorophyll maxima. (*Lakes discarded from the analysis, see text).

<table><caption>TABLE 2 | Average and range of values of the different variables measured in the SML and surface plankton of the studied lakes.</caption><thead><tr><th>Variable</th><th>SML</th><th>Surface</th><th>P (paired t-test)</th></tr></thead><tbody><tr><td colspan="4"><strong>PHOTOTROPHIC PICOPLANKTON</strong></td></tr><tr><td>Total PPP abundance</td><td>7.1 10<sup>3</sup> (0.3–18.9)</td><td>6.3 10<sup>3</sup> (0.5–12.6)</td><td>NS</td></tr><tr><td>Picocyanobacteria</td><td>2.4 10<sup>3</sup> (0–13.0)</td><td>1.3 10<sup>3</sup> (0–9.5)</td><td>0.02</td></tr><tr><td>Small eukaryotes</td><td>4.7 10<sup>3</sup> (0.3–13.7)</td><td>5.0 10<sup>3</sup> (0.2–10.8)</td><td>NS</td></tr><tr><td>Ratio small eukaryotes:cyanobacteria</td><td>150 (0–650)</td><td>238 (0–800)</td><td>0.03</td></tr><tr><td>Number of diff. PPP populations</td><td>4.1 (1–7)</td><td>4.5 (2–7)</td><td>NS</td></tr><tr><td>Number of diff. picocyanobacterial pop.</td><td>0.8 (0–2)</td><td>0.6 (0–2)</td><td>NS</td></tr><tr><td>Number of diff. small eukaryotes populations</td><td>3.3 (1–5)</td><td>3.9 (2–6)</td><td>0.04</td></tr><tr><td colspan="4"><strong>HETEROTROPHIC PICOPLANKTON</strong></td></tr><tr><td>Total heterotrophic prokaryotes</td><td>1.0 10<sup>6</sup> (0.3–2.0)</td><td>0.8 10<sup>6</sup> (0.3–1.8)</td><td>0.03</td></tr><tr><td>NADS-Live prokaryotes</td><td>5.0 10<sup>5</sup> (0.7–9.9)</td><td>4.9 10<sup>5</sup> (1.1–12.1)</td><td>NS</td></tr><tr><td>NADS-Dead prokaryotes</td><td>5.1 10<sup>5</sup> (1.8–11.2)</td><td>3.2 10<sup>5</sup> (0.8–5.4)</td><td>0.01</td></tr><tr><td>HNF abundance</td><td>3.0 10<sup>3</sup> (0.9–8.8)</td><td>2.0 10<sup>3</sup> (0.9–5.7)</td><td>0.02</td></tr><tr><td>Ratio Total heterotrophic prokaryotes:HNF</td><td>482 (63–1297)</td><td>491 (135–1122)</td><td>NS</td></tr><tr><td>CTC positive prokaryotes</td><td>3.7 10<sup>4</sup> (2.1–6.1)</td><td>2.8 10<sup>4</sup> (1.6–5.8)</td><td>0.01</td></tr><tr><td>Prokaryotic heterotrophic act. (leucine inc.)</td><td>9.5 (1.4–29.0)</td><td>27.9 (7.4–61.4)</td><td>&lt;0.01</td></tr><tr><td>Cell-specific prokaryotic activity</td><td>9.9 (3–31)</td><td>36.0 (16–64)</td><td>&lt;0.01</td></tr><tr><td>% Live prokaryotes</td><td>49 (22–70)</td><td>59 (44–82)</td><td>&lt;0.01</td></tr></tbody></table>

Abundances in cells ml⁻¹. Bacterial activity in pmol Leu ml⁻¹ h⁻¹ and cell-specific activity in amol cell⁻¹ h⁻¹. Averages and the statistical comparison of the picophytoplankton (PPP) abundance data do not include Bassa de les Granotes, an obvious outlier. All variables were log transformed to achieve normality prior to the paired t-test.

the variables that did not randomly lay around the 1:1 ratio (Figure 2).

Leucine incorporation varied in a large range, from 1.4 to 121 pmol Leu l⁻¹ d⁻¹ and was systematically higher in the

surface than in the SML with the exception of lakes Pois and Roi (Figure 2A). A similar plot for CTC indicated that the number of heterotrophic prokaryotic cells with high respiration rates was higher in the SML than at the surface, as most values were above
<figure><img src="image_2.png" /></figure>

<figure><img src="image_3.png" /></figure>

<figure><img src="image_4.png" /><figcaption>**FIGURE 1** | (A) Phototrophic small eukaryotes, (B) *Synechococcus*-like, and (C) total phototrophic picoplankton abundance in SML vs. surface waters of the different high altitude lakes studied, as determined by flow cytometry. In all panels, the dashed line indicates the 1:1 ratio.</figcaption></figure>

the 1:1 ratio (Figure 2B), Overall, there were more respiring prokaryotic cells at the SML, while they produced more in surface waters, both trends being highly significant (Table 2).

The NADS staining protocol (that measures the relationships between membrane-intact and membrane-damaged heterotrophic prokaryotic cells) showed healthier heterotrophic prokaryotes in surface plankton than at the neuston (Figure 2C). In fact, damaged cells accumulated at the SML while intact cells had similar abundances in both layers (Table 2).

## Enrichment Factors

Enrichment factors were calculated as the ratio between the values of abundance or activity in the SML divided by those in surface waters. The results of enrichment factors indicated different degrees of enrichment in the neuston for CTC positive, and bacterial, HNF and PPP abundances. On the contrary, leucine incorporation and the percentage of “live” bacterial cells were lower in the neuston (Figure 3).

## Factor Analysis

FA simplifies the interpretation of the many environmental variables measured in the 16 lakes. We performed separated FA for physical and chemical variables, and the three first factors explained 85.3–64.5% of the variance, respectively (Table 3). The first three physical factors reflected lake morphometry (Factor F1), altitude (Factor F2), and water temperature (Factor F3). Regarding chemical factors, the FA analysis revealed productivity (Factor C1) as the main factor, then DOC, silica, sulfate, and nitrate concentration (Factor C2), and finally conductivity and calcium concentration (Factor C3), which would capture the lithology as reflected in water chemistry. A fourth factor was obtained on which $NH_4^+$ (negatively) and $NO_2$ and alkalinity (positively) loaded. However, this factor was forced by a single outlier lake. It did not show any correlation with any biological parameter and therefore it is not discussed here. Most factors were easily interpretable with exception of factor C2 (DOC, silica, sulfate, and nitrate). In this factor, DOC, silica, and sulfate were the strongest variables (higher loadings), and nitrate, to a lesser extend. Silica and sulfate are usually related to the metamorphic nature of rocks in the Pyrenees (Catalan et al., 1993). This factor C2 could be interpreted as a “watershed DOC factor,” representing a gradient from a rockier watershed with low DOC production (high silica and sulfate from metamorphic rock weathering), to higher soil coverage with higher DOC production (less silica and sulfate). An alternative explanation is that the covariation of these three parameters could be related to the landscape position of the lakes (Sadro et al., 2011).

We tested for correlation of the scores of the physical factors against those of the chemical factors and no significant correlations were found, indicating that the factors were independent (and so were the variables that laid behind them). One exception was the correlation between factors F3 (water temperature) and C2 (watershed DOC), which were positively correlated ($p < 0.05$, Pearson coefficient = 0.518).

In a second step, the abundance or activity variables in the two layers (surface and SML), and the enrichment factors were
<figure><img src="image_5.png" /><figcaption>**FIGURE 2** | (A) Leucine incorporation in SML samples vs. values observed in the surface; (B) Number of heterotrophic prokaryote cells with high respiration rates (CTC+) in SML samples vs. those observed in the surface as determined by flow cytometry and the CTC staining protocol. (C) Percentage of "live" cells in the SML samples vs. the value observed at the surface as determined by flow cytometry and the NADS staining protocol.</figcaption></figure>

(D) Abundance of HNF in SML samples vs. values observed in the surface measured by flow cytometry with the Lysotracker staining protocol. In all panels, the dashed line indicates the 1:1 ratio. Insets in each panel indicate the distribution of values for surface and the SML. The central line indicates the median value, boxes indicate the lower and upper quartiles, and vertical lines indicate the 10th and 90th percentiles.

tested for correlation against the scores of the main chemical and physical factors (Table 4). In general, biological variables in the SML were more related to physical factors, especially lake morphometry (Factor F1) and temperature (Factor F3) and to the nature of the rock (Factor C3). Conversely, surface biological variables were frequently correlated with the Factor C2 (watershed DOC), but also with water temperature (Factor F3). Concerning enrichment factors, the same pattern was observed: stronger and more correlations were found with the Factor C2 (watershed DOC).

Discussion

Consistent differences were observed for picoplankton abundance, activity, and single cell physiological status between the SML and surface waters for the great majority of the 16 high mountain lakes that had developed a SML. In most lakes, PPP abundances were higher in the SML, particularly when PE-rich populations were present, whereas leucine incorporation rates in the SML were clearly below the values found in underlying waters, and these differences were probably related to the *in*
<figure><img src="image_6.png" /><figcaption>FIGURE 3 | Boxplot of abundance or activity enrichment factors in the piconeuston (i.e., SML:surface). Central line indicates median value, boxes the lower and upper quartiles, and vertical lines indicate 10th and 90th percentiles.</figcaption></figure>

situ stressful environmental conditions experienced by the neuston.

All the lakes studied had several PPP populations, and roughly half of them had picocyanobacteria. As far as we are aware, no other study had reported such a widespread presence of PPP in this type of environments, i.e., lakes located above 2000 m, although a large genetic diversity of planktonic eukaryotic microorganisms (3–40 microns sized) has been recently reported from the same Pyrenean lakes dataset (Triadó-Margarit and Casamayor, 2012).

Concerning heterotrophic prokaryotes, leucine incorporation was clearly lower in the SML compared to surface waters. A possible explanation for this observation is the direct effect of sunlight, particularly UVR, on the prokaryotic community. Several studies have demonstrated the inhibitory effect of solar radiation on substrate (thymidine, leucine, glucose, and acetate) incorporation rates on exposed prokaryotic communities, compared to dark or UVR-free incubations (Herndl et al., 1993; Sommaruga et al., 1997; Santos et al., 2012b; Ruiz-González et al., 2013). Actually, these processes are known to be by far more complex, as different bacterial groups have different sensitivities to UVR exposure, some of them being strongly inhibited (such as the SAR11 clade, at least in marine systems), some resistant (such as *Gammaproteobacteria* and *Bacteroidetes*), and some stimulated (such as *Roseobacter*) by natural sunlight (Alonso-Sáez et al., 2006; Santos et al., 2012a). The observed bulk community response could therefore be the sum of the positive and negative effects on the different prokaryotic populations.

The theoretical framework outlined by Carlson et al. (2007) postulates that increased environmental hostility (lack

<table><caption>TABLE 3 | Factor loadings and cumulative variance explained obtained in the factor analysis of chemical and physical variables of the 17 high altitude lakes studied.</caption><thead><tr><th rowspan="2">Physical variables</th><th>Factor F1</th><th>Factor F2</th><th>Factor F3</th></tr><tr><th>Lake morphometry</th><th>-Altitude</th><th>-Temperature</th></tr></thead><tbody><tr><td>Zmax</td><td>0.946</td><td></td><td></td></tr><tr><td>Lake area</td><td>0.874</td><td></td><td></td></tr><tr><td>Secchi depth</td><td>0.685</td><td></td><td></td></tr><tr><td>Watershed area</td><td></td><td>0.926</td><td></td></tr><tr><td>Altitude</td><td></td><td>-0.884</td><td></td></tr><tr><td>Water temperature</td><td></td><td></td><td>-0.911</td></tr><tr><td>Cumulative Variance</td><td>37.8%</td><td>67.4%</td><td>85.3%</td></tr></tbody><thead><tr><th rowspan="2">Chemical variables</th><th>Factor C1</th><th>Factor C2</th><th>Factor C3</th></tr><tr><th>Productivity</th><th>-Watershed DOC</th><th>Lithology</th></tr></thead><tbody><tr><td>TP</td><td>0.922</td><td></td><td></td></tr><tr><td>Chla</td><td>0.776</td><td></td><td></td></tr><tr><td>DOC</td><td></td><td>-0.769</td><td></td></tr><tr><td>DRSi</td><td></td><td>0.765</td><td></td></tr><tr><td>SO<sub>4</sub></td><td></td><td>0.678</td><td></td></tr><tr><td>NO<sub>3</sub></td><td></td><td>0.526</td><td></td></tr><tr><td>Conductivity</td><td></td><td></td><td>0.951</td></tr><tr><td>Ca<sub>2</sub></td><td></td><td></td><td>0.937</td></tr><tr><td>Cumulative Variance</td><td>21.8%</td><td>43.4%</td><td>64.5%</td></tr></tbody></table>

of resources - organic matter or nutrients, salinity, extreme temperature, UVR, extreme pH values, or a combination of several of these) should result in an increase in the proportion of the total flux of energy that is devoted to cell maintenance. Associated with this increase in cell maintenance, cell-specific respiration should increase to fuel the increased maintenance and the repair mechanisms. Thus, bacterial growth efficiency (BGE) would be expected to decrease with increasing environmental hostility. Our results provide experimental evidence that validate this hypothesis, at least in high-mountain lakes with high levels of UVR: not only bulk bacterial activity was lower in the SML (ca. four times lower on average), but the number of heterotrophic prokaryote cells with high respiration rates (CTC positives) was higher in the SML than at the surface (1.5x). We also observed a higher percentage of membrane-damaged heterotrophic prokaryotes in the SML than in the surface, probably as a consequence of the postulated environmental hostility. The number of CTC-positive cells has been found to correlate with total bacterial respiration in several previous studies (Cook and Garland, 1997; Smith, 1998; Berman et al., 2001).

It is unlikely that virus-induced lysogeny could be the main cause for the differences observed between layers, as it is known that solar radiation, particularly UV-B radiation (290–320 nm), is an important factor contributing to the decay of aquatic viruses (Suttle and Chen, 1992). Moreover, a previous study in a high mountain lake reported a negative correlation between viral abundances and solar radiation (Hofer and
<table><caption>TABLE 4 | Pearson correlation coefficient (r) of the statistically significant correlations (*p* < 0.05 or **p* < 0.01) between the variables measured in the 17 high altitude lakes and the chemical and physical scores (obtained in the factor analyses, see Table 3).</caption><thead><tr><th></th><th>Factor F1</th><th>Factor F2</th><th>Factor F3</th><th>Factor C1</th><th>Factor C2</th><th>Factor C3</th></tr><tr><th></th><th>Lake morphometry</th><th>-Altitude</th><th>-Temperature</th><th>Productivity</th><th>-Watershed DOC</th><th>Lithology</th></tr></thead><tbody><tr><td><strong>SML</strong></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Total heterotrophic prokaryotes</td><td>-0.654**</td><td></td><td></td><td></td><td></td><td>0.558*</td></tr><tr><td>% Live prokaryotes</td><td>-0.536*</td><td></td><td>-0.523*</td><td></td><td></td><td></td></tr><tr><td>Total PPP abundance</td><td></td><td></td><td>-0.652**</td><td></td><td></td><td></td></tr><tr><td>Prokaryotic heterotrophic activity (leucine incorporation)</td><td></td><td>0.540*</td><td></td><td></td><td></td><td></td></tr><tr><td><strong>SURFACE</strong></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>% Live prokaryotes</td><td>-0.498*</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>CTC positive prokaryotes</td><td></td><td>0.535*</td><td></td><td></td><td></td><td></td></tr><tr><td>Ratio Total heterotrophic prokaryotes:HNF</td><td></td><td>-0.531*</td><td></td><td></td><td>0.488*</td><td>-0.544*</td></tr><tr><td>Total heterotrophic prokaryotes</td><td></td><td></td><td>-0.628**</td><td></td><td>-0.533*</td><td></td></tr><tr><td>Total PPP abundance</td><td></td><td></td><td></td><td>0.530*</td><td>-0.540*</td><td></td></tr><tr><td>Prokaryotic heterotrophic activity (leucine incorporation)</td><td></td><td></td><td></td><td></td><td>-0.567*</td><td></td></tr><tr><td><strong>ENRICHMENT FACTOR</strong></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Total PPP abundance</td><td></td><td></td><td>-0.506*</td><td></td><td></td><td></td></tr><tr><td>Ratio Total heterotrophic prokaryotes:HNF</td><td></td><td></td><td></td><td>0.677**</td><td></td><td></td></tr><tr><td>Prokaryotic heterotrophic activity (leucine incorporation)</td><td></td><td></td><td></td><td></td><td>0.713**</td><td></td></tr></tbody></table>

Sommaruga, 2001). In addition to HNF pressure, higher number of membrane-damaged heterotrophic prokaryotes cells in the SML comparing to surface, can also be explained by the fact that many airborne bacteria can be trapped in the air-water interface (Hervàs and Casamayor, 2009).

Our observations contrast with a recent study in 6 alpine lakes situated in an altitude gradient, which could not find a consistent pattern in the comparison of bacterial production between layers, and attributed the variation observed to lake-specific communities and environmental conditions (Hörtnagl et al., 2010). However, in such study, 3 out of 4 lakes located above 1500 m a.s.l. had lower bacterial production values in the SML than at the surface. In our study, most lakes were above 2000 m a.s.l.. As UVR increases with altitude (Blumthaler et al., 1997), this might explain the discrepancies between our study and that of Hörtnagl et al. Additionally, in our case, we only sampled under very stable weather conditions, sunny and calm, something that might have facilitated the detection of a well-formed SML with a microbiota strongly affected by solar radiation.

To the best of our knowledge, this is the first report of HNF abundance in freshwater SML and our results show that the neuston holds generally more HNF than surface lake water. The enrichment factor of HNF in the SML was, in average, slightly higher (1.47) than that of bacteria (1.27). These observations are somehow unexpected, as previous studies had shown that UVR affects negatively HNF growth and grazing on bacteria (Sommaruga et al., 1999). The total heterotrophic prokaryote:HNF ratio might provide some hints on this issue, and although not significantly different, this ratio was slightly higher in the SML. In addition, absolute values (~480–490 bacteria per HNF) were relatively low compared to literature data (Gasol, 1994; Fermani et al., 2013). This might indicate relatively high grazing pressure of HNF on bacteria, which might

have to do with specific features of food web structure in high altitude lakes. *Daphnia*, and large cladocerans in general, are particularly sensitive to UVR reviewed by Sommaruga (2001). The dominance of large filter feeding cladocerans reduces HNF abundance, and affects the microbial food web (Gasol et al., 1995; Zöllner et al., 2003; Sarmento, 2012). Thus, HNF from high altitude lakes are likely to undergo low predation pressure, and this would result in higher grazing pressure on bacteria, providing an explanation for the relatively low heterotrophic prokaryote:HNF ratio observed in this study.

No simple and direct linear relationships could be found between microbial abundances or activities and raw environmental variables. Nevertheless, factors derived from the FA did correlate in some cases with biological variables. Both physical and chemical FA captured a fairly good percentage of the total variance (85.3–64.5% respectively), which means that the extracted factors successfully explained most of the variability observed in our dataset. It has to be noted that our data do not account of course for all aspects of environmental variability. Our discussion includes only those aspects captured by our measurements, and further research may seed light on other aspects not considered here.

The physical and chemical multivariate factors identified in the FA were not correlated between them, reinforcing the independency of the factors, with the exception of factors F3 (-temperature) and C2 (-watershed DOC), which were positively correlated (*p* < 0.05, *r* = 0.518). It is possible that higher temperatures promote higher DOC production in the watershed increasing allochthonous DOC inputs to the water. Two alternative hypotheses are that lakes with higher DOC are warmer because of solar radiation absorption, or that snowmelt increases stream DOC inputs into lakes.
Interestingly, the factors that influenced heterotrophic prokaryote abundance in the different layers were different. On one hand, prokaryote abundance in SML was negatively correlated with lake morphometry (factor F1). On the other hand, surface abundance was negatively correlated with factor F3 (–temperature). Larger lakes have greater fetch and turbulence, and probably present less well-defined SML and present them less often than smaller lakes. Larger lakes might thus lack some taxa well-adapted to SML conditions (environmental filtering), while surface prokaryote abundance depends on water temperature, regardless of the size of the lake.

Prokaryote abundance and activity in the surface co-varied with watershed DOC (factor C1) but not in the SML (Table 4). This suggests that surface and piconeuston microbial communities rely on different energy and/or nutrient sources. Pyrenean high mountain lakes receive considerable amounts of atmospheric deposition, namely from the Sahara desert (Mladenov et al., 2011; Barberán et al., 2014) that might be the main source of nutrients and organic matter for the piconeuston.

Another interesting outcome of the FA is the fact that altitude (factor F2) correlated negatively with prokaryotic activity in the SML. In other words, leucine incorporation in the SML decreased with altitude. UVR increases with altitude, but part of the radiation is reflected and absorbed in the first millimeters of the water column, and the UVR that reaches subsurface water is not as high as that reaching the SML (Sommaruga and Psenner, 1997), This probably explains why leucine incorporation in subsurface waters had no relationship with altitude. However, in the SML exposure levels are of such magnitude that anabolic processes are affected according to the UVR dose (i.e., directly correlated to altitude).

Finally, the prokaryotic activity enrichment factor (which was always lower than one, **Figure 3**) was correlated to factor C2 (–watershed DOC), meaning that in rocky watersheds with low DOC production, neuston and surface prokaryotic activity had relatively small differences, and these differences increased in DOC enriched watersheds. Basically, prokaryotic activity increased in DOC-rich surface waters, but not in neuston, indicating that despite the physical proximity between these two habitats, the processes shaping microbial metabolism were substantially different.

## References

* Agogué, H., Casamayor, E. O., Bourrain, M., Obernosterer, I., Joux, F., Herndl, G. J., et al. (2005). A survey on bacteria inhabiting the sea surface microlayer of coastal ecosystems. *FEMS Microbiol. Ecol.* 54, 269–280. doi: 10.1016/j.femsec.2005.04.002

* Agogué, H., Casamayor, E. O., Joux, F., Obernosterer, I., Dupuy, C., Lantoine, F., et al. (2004). Comparison of samplers for the biological characterization of the sea surface microlayer. *Limnol. Oceanogr. Methods* 2, 213–225. doi: 10.4319/lom.2004.2.213

* Alonso-Sáez, L., Gasol, J. M., Lefort, T., Hofer, J., and Sommaruga, R. (2006). Effect of natural sunlight on bacterial activity and differential sensitivity of natural bacterioplankton groups in northwestern Mediterranean coastal waters. *Appl. Environ. Microbiol.* 72, 5806–5813. doi: 10.1128/AEM.00597-06

In conclusion, we found evidence that heterotrophic prokaryotic communities in the piconeuston of high altitude lakes had lower percentage of healthy cells, lower leucine incorporation rates and higher number of cells with high respiration rates (CTC positive), likely resulting in lower BGE as compared to those of the underlying surface water, evidencing a high degree of environmental harshness. PPP were widespread in these high mountain lakes and PE-rich picocyanobacteria were particularly prominent in the SML. Piconeuston was also enriched in HNF and heterotrophic prokaryotes. Overall, we demonstrate that steep gradients in the upper millimeters of surface waters of high mountain lakes are common in stable conditions, and that different mechanisms drive microbial processes in the SML and the underlying surface water.

## Acknowledgments

We thank the reviewers and the Guest Associate Editor Jérôme Comte for valuable comments on the manuscript. The authors are thankful to Daniel Diaz de Quijano, Miguel Ballen and Berta Fueyo for help with sampling. Thanks to the authorities of the Parc de Sant Maurici i Aigüestortes for allowing sampling within this protected area. This work was supported by grants 2009SGR/1177 “Grup d’estructura i funcio de xarxes tròfiques microbianes planctòniques” and 2009SGR/361 “Grup d’Ecologia dels Canvis Ambientals, GECA” from Generalitat de Catalunya, AERBAC-2 178/2010 (to EOC) and EGALA 124/2010 (MF, LC) from the Ministerio de Medio Ambiente-Red de Parques Nacionales, and PIRENA CGL2009-13318 (EOC) and NITROPIR CGL2010-19373 (MF, LC) from MINECO. HS’s work was supported by CNPq and FAPESP (Process: 2014/14139-3).

## Supplementary Material

The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fmicb.2015.00361/abstract

Detailed temperature profiles in the high mountain lakes studied in the sampling day.

* Auguet, J. C., and Casamayor, E. O. (2008). A hotspot for cold crenarchaeota in the neuston of high mountain lakes. *Environ. Microbiol.* 10, 1080–1086. doi: 10.1111/j.1462-2920.2007.01498.x

* Barberán, A., Henley, J., Fierer, N., and Casamayor, E. O. (2014). Structure, inter-annual recurrence, and global-scale connectivity of airborne microbial communities. *Science Total Environ.* 487, 187–195. doi: 10.1016/j.scitotenv.2014.04.030

* Berman, T., Kaplan, B., Chava, S., Viner, Y., Sherr, B. F., and Sherr, E. B. (2001). Metabolically active bacteria in Lake Kinneret. *Aquat. Microb. Ecol.* 23, 213–224. doi: 10.3354/ame023213

* Blumthaler, M., Ambach, W., and Ellinger, R. (1997). Increase in solar UV radiation with altitude. *J. Photochem. Photobiol. B Biol.* 39, 130–134. doi: 10.1016/S1011-1344(96)00018-8
Camarero, L. (1994). Assay of soluble reactive phosphorus at nanomolar levels in nonsaline waters. *Limnol. Oceanogr.* 39, 707–711. doi: 10.4319/lo.1994.39.3.0707

Carlson, C. A., del Giorgio, P. A., and Herndl, G. J. (2007). Microbes and the dissipation of energy and respiration: from cells to ecosystems. *Oceanography* 20, 89–100. doi: 10.5670/oceanog.2007.52

Catalan, J., Ballesteros, E., Gacia, E., Palau, A., and Camarero, L. (1993). Chemical composition of disturbed and undisturbed high-mountain lakes in the Pyrenees: a reference for acidified sites. *Water Res.* 27, 133–141. doi: 10.1016/0043-1354(93)90203-T

Catalan, J., Camarero, L., Felip, M., Pla, S., Ventura, M., Buchaca, T., et al. (2006). High mountain lakes: extreme habitats and witnesses of environmental change. *Limnetica* 25, 551–584.

Conrad, R., and Seiler, W. (1988). Influence of the surface microlayer on the flux of nonconservative trace gases (CO, H₂, CH₄, N₂O) across the ocean-atmosphere interface. *J. Atmos. Chem.* 6, 83–94. doi: 10.1007/BF00048333

Cook, K. L., and Garland, J. L. (1997). The relationship between electron transport activity as measured by CTC Reduction and CO₂ production in mixed microbial communities. *Microb. Ecol.* 34, 237–247. doi: 10.1007/s002489900053

Cunliffe, M., Upstill-Goddard, R. C., and Murrell, J. C. (2011). Microbiology of aquatic surface microlayers. *FEMS Microbiol. Rev.* 35, 233–246. doi: 10.1111/j.1574-6976.2010.00246.x

Del Giorgio, P. A., and Gasol, J. M. (2008). "Physiological structure and single-cell activity in marine bacterioplankton," in *Microbial Ecology of the Oceans, 2nd Edn.*, ed D. L. Kirchman (New York, NY: John Wiley & Sons, Inc.), 243–298. doi: 10.1002/9780470281840.ch8

del Giorgio, P. A., Gasol, J. M., Vaqué, D., Mura, P., Agustí, S., and Duarte, C. M. (1996). Bacterioplankton community structure: protists control net production and the proportion of active bacteria in a coastal marine community. *Limnol. Oceanogr.* 41, 1169–1179. doi: 10.4319/lo.1996.41.6.1169

Dietz, A. S., Albright, L. J., and Tuominen, T. (1976). Heterotrophic activities of bacterioneuston and bacterioplankton. *Can. J. Microbiol.* 22, 1699–1709. doi: 10.1139/m76-251

Falcioni, T., Papa, S., and Gasol, J. M. (2008). Evaluating the flow-cytometric nucleic acid double-staining protocol in realistic situations of planktonic bacterial death. *Appl. Environ. Microbiol.* 74, 1767–1779. doi: 10.1128/AEM. 01668-07

Fermani, P., Diovisalvi, N., Torremorell, A., Lagomarsino, L., Zagarese, H., and Unrein, F. (2013). The microbial food web structure of a hypertrophic warm-temperate shallow lake, as affected by contrasting zooplankton assemblages. *Hydrobiologia* 714, 115–130. doi: 10.1007/s10750-013-1528-3

Gasol, J. M., and del Giorgio, P. A. (2000). Using flow cytometry for counting natural planktonic bacteria and understanding the structure of planktonic bacterial communities. *Sci. Mar.* 64, 197–224. doi: 10.3989/scimar.2000.64n2197

Gasol, J. M., Simons, A. M., and Kalff, J. (1995). Patterns in the top-down versus bottom-up regulation of heterotrophic nanoflagellates in temperate lakes. *J. Plankton Res.* 17, 1879–1903. doi: 10.1093/plankt/17.10.1879

Gasol, J. M. (1994). A framework for the assessment of top-down vs bottom-up control of heterotrophic nanoflagellate abundance. *Mar. Ecol. Prog. Ser.* 113, 291–300. doi: 10.3354/meps113291

Grégori, G., Citterio, S., Ghiani, A., Labra, M., Sgorbati, S., Brown, S., et al. (2001). Resolution of viable and membrane-compromised bacteria in freshwater and marine waters based on analytical flow cytometry and nucleic acid double staining. *Appl. Environ. Microbiol.* 67, 4662–4670. doi: 10.1128/AEM.67.10.4662-4670.2001

Herndl, G. J., Muller-Niklas, G., and Frick, J. (1993). Major role of ultraviolet-B in controlling bacterioplankton growth in the surface layer of the ocean. *Nature* 361, 717–719. doi: 10.1038/361717a0

Hervás, A., and Casamayor, E. O. (2009). High similarity between bacterioneuston and airborne bacterial community compositions in a high mountain lake area. *FEMS Microbiol. Ecol.* 67, 219–228. doi: 10.1111/j.1574-6941.2008.00617.x

Hofer, J. S., and Sommaruga, R. (2001). Seasonal dynamics of viruses in an alpine lake: Importance of filamentous forms. *Aquat. Microb. Ecol.* 26, 1–11. doi: 10.3354/ame026001

Hörtnagl, P., Perez, M. T., and Sommaruga, R. (2010). Living at the border: a community and single-cell assessment of lake bacterioneuston activity. *Limnol. Oceanogr.* 55, 1134–1144. doi: 10.4319/lo.2010.55.3.1134

Jeffrey, S. W., and Humphrey, G. F. (1975). New spectrophotometric equations for determining chlorophylls a1, b1, c1 and c2 in higher plants, algae and natural phytoplankton. *Biochem. Physiol. Pflanz.* 167, 194–204.

Kirchman, D., K'nees, E., and Hodson, R. (1985). Leucine incorporation and its potential as a measure of protein synthesis by bacteria in natural aquatic systems. *Appl. Environ. Microbiol.* 49, 599–607.

Kuznetsova, M., and Lee, C. (2001). Enhanced extracellular enzymatic peptide hydrolysis in the sea-surface microlayer. *Mar. Chem.* 73, 319–332. doi: 10.1016/S0304-4203(00)00116-X

Kuznetsova, M., Lee, C., Aller, J., and Frew, N. (2004). Enrichment of amino acids in the sea surface microlayer at coastal and open ocean sites in the North Atlantic Ocean. *Limnol. Oceanogr.* 49, 1605–1619. doi: 10.4319/lo.2004.49.5.1605

Maki, J. S. (1993). "The air-water interface as an extreme environment," in *Aquatic Microbiology: An Ecological Approach*, ed T. E. Ford (Boston; Oxford: Blackwell Scientific Publications), 409–439.

Mladenov, N., Sommaruga, R., Morales-Baquero, R., Laurion, I., Camarero, L., Dieguez, M. C., et al. (2011). Dust inputs and bacteria influence dissolved organic matter in clear alpine lakes. *Nat. Commun.* 2:405. doi: 10.1038/ncomms1411

Naumann, E. (1917). Beiträge zur kenntnis des teichnanoplanktonns. II. Über das neuston des süßwasser. *Biologisches Zentralblatt* 37, 98–106.

Olson, R. J., Zettler, E. R., and Du Rand, M. D. (1993). "Phytoplankton analysis using flow cytometry," in *Handbook of Methods in Aquatic Microbial Ecology*, eds P. F. Kemp, B. F. Sherr, E. B. Sherr and J. J. Cole (Boca Raton FL: Lewis Publishers), 175–186.

R Core Team (2014). *R: A Language and Environment for Statistical Computing*. Foundation for Statistical Computing, Vienna, Austria. Available online at: http://www.R-project.org/

Rose, J. M., Caron, D. A., Sieracki, M. E., and Poulton, N. (2004). Counting heterotrophic nanoplanctonic protists in cultures and aquatic communities by flow cytometry. *Aquat. Microb. Ecol.* 34, 263–277. doi: 10.3354/ame034263

Ruiz-González, C., Simó, R., Sommaruga, R., and Gasol, J. M. (2013). Away from darkness: a review on the effects of solar radiation on heterotrophic bacterioplankton activity. *Front. Microbiol.* 4:131. doi: 10.3389/fmicb.2013. 00131

Sadro, S., Nelson, C. E., and Melack, J. M. (2011). The influence of landscape position and catchment characteristics on aquatic biogeochemistry in high-elevation lake-chains. *Ecosystems* 15, 363–386. doi: 10.1007/s10021-011-9515-x

Santos, A. L., Baptista, I., Lopes, S., Henriques, I., Gomes, N. C. M., Almeida, A., et al. (2012a). The UV responses of bacterioneuston and bacterioplankton isolates depend on the physiological condition and involve a metabolic shift. *FEMS Microbiol. Ecol.* 80, 646–658. doi: 10.1111/j.1574-6941.2012.01336.x

Santos, A. L., Oliveira, V., Baptista, I., Henriques, I., Gomes, N. C. M., Almeida, A., et al. (2012b). Effects of UV-B Radiation on the structural and physiological diversity of bacterioneuston and bacterioplankton. *Appl. Environ. Microbiol.* 78, 2066–2069. doi: 10.1128/AEM.06344-11

Sarmento, H., Unrein, F., Isumbisho, M., Stenuite, S., Gasol, J. M., and Descy, J. P. (1998). Abundance and distribution of picoplankton in tropical, oligotrophic Lake Kivu, eastern Africa. *Freshw. Biol.* 53, 756–771. doi: 10.1111/j.1365- 2427.2007.01939.x

Sarmento, H. (2012). New paradigms in tropical limnology: the importance of the microbial food web. *Hydrobiologia* 686, 1–14. doi: 10.1007/s10750-012-1011-6

Schiaffino, M. R., Gasol, J. M., Izaguirre, I., and Unrein, F. (2013). Picoplankton abundance and cytometric group diversity along a trophic and latitudinal lake gradient. *Aquat. Microb. Ecol.* 68, 231–250. doi: 10.3354/ame01612

Sherr, B. F., del Giorgio, P., and Sherr, E. B. (1999). Estimating abundance and single-cell characteristics of respiring bacteria via the redox dye CTC. *Aquat. Microb. Ecol.* 18, 117–131. doi: 10.3354/ame018117

Sieracki, M. E., Cucci, T. L., and Nicinski, J. (1999). Flow cytometric analysis of 5-cyano-2,3-ditolyl tetrazolium chloride activity of marine bacterioplankton in dilution cultures. *Appl. Environ. Microbiol.* 65, 2409–2417.

Smith, D. C., and Azam, F. (1992). A simple, economical method for measuring bacterial protein synthesis rates in seawater using 3H-leucine. *Mar. Microb. Food Webs* 6, 107–114.

Smith, E. M. (1998). Coherence of microbial respiration rate and cell-specific bacterial activity in a coastal planktonic community. *Aquat. Microb. Ecol.* 16, 27–35. doi: 10.3354/ame016027
Sommaruga, R., Obernosterer, I., Herndl, G. J., and Psenner, R. (1997). Inhibitory effect of solar radiation on thymidine and leucine incorporation by freshwater and marine bacterioplankton. *Appl. Environ. Microbiol.* 63, 4178–4184.

Sommaruga, R., and Psenner, R. (1997). Ultraviolet radiation in a high mountain lake of the austrian alps: air and underwater measurements. *Photochem. Photobiol.* 65, 957–963. doi: 10.1111/j.1751-1097.1997.tb07954.x

Sommaruga, R., Sattler, B., Oberleiter, A., Wille, A., Wögrath-Sommaruga, S., Psenner, R., et al. (1999). An *in situ* enclosure experiment to test the solar UVB impact on plankton in a high-altitude mountain lake. *II. Effects on the microbial food web*. *J. Plankton Res.* 21, 859–876. doi: 10.1093/plankt/21.5.859

Sommaruga, R. (2001). The role of solar UV radiation in the ecology of alpine lakes. *J. Photochem. Photobiol. B Biol.* 62, 35–42. doi: 10.1016/S1011-1344(01)00154-3

Suttle, C. A., and Chen, F. (1992). Mechanisms and rates of decay of marine viruses in seawater. *Appl. Environ. Microbiol.* 58, 3721–3729.

Triadó-Margarit, X., and Casamayor, E. O. (2012). Genetic diversity of planktonic eukaryotes in high mountain lakes (Central Pyrenees, Spain). *Environ. Microbiol.* 14, 2445–2456. doi: 10.1111/j.1462-2920.2012.02797.x

Vila-Costa, M., Bartrons, M., Catalan, J., and Casamayor, E. O. (2014). Nitrogen-Cycling Genes in Epilithic biofilms of oligotrophic high-altitude lakes (Central Pyrenees, Spain). *Microb. Ecol.* 68, 60–69. doi: 10.1007/s00248-014-0417-2

Zöllner, E., Santer, B., Boersma, M., Hoppe, H.-G., and Jürgens, K. (2003). Cascading predation effects of Daphnia and copepods on microbial food web components. *Freshw. Biol.* 48, 2174–2193. doi: 10.1046/j.1365-2426.2003.01158.x

**Conflict of Interest Statement:** 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.

Copyright © 2015 Sarmento, Casamayor, Auguet, Vila-Costa, Felip, Camarero and Gasol. 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) or licensor 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.