DOM Biological Lability in an Estuarine System in Two Contrasting Periods

: Estuarine processes play a key role in determining the amount and quality of land-derived dissolved organic matter (DOM) reaching the oceans. Microbial-mediated reactions can affect the concentration, quality, and bioavailability of DOM within an estuary. In this study, we investigated biological DOM removal in a small estuary and its variability in two contrasting seasons (spring and autumn) characterized by natural differences in the concentration and quality of the riverine DOM. Two incubation experiments were carried out using natural DOM and heterotrophic prokaryotes community collected at the estuary in March and September. Dissolved organic carbon (DOC) concentration, DOM ﬂuorescence, and the heterotrophic prokaryotes abundance (HPA) showed marked differences between the two seasons. These parameters were followed through time for up to two months. Despite the marked differences in the initial conditions, the DOC removal rates were surprisingly similar in the two periods (16 µ M DOC month − 1 in March and 18 µ M DOC month − 1 in September), with the biggest removal in the ﬁrst 48 h. The trend of ﬂuorescent DOM (FDOM) during the incubation showed marked differences between the two periods. In March, the net removal of all the FDOM components was observed consistently with the decrease in DOC; whereas, in September, the net production of humic-like substances was observed.


Introduction
A large amount of dissolved organic matter (DOM) is continuously exported from the land to the oceans through rivers and estuaries. It has been estimated that 250 Tg of dissolved organic carbon (DOC) are delivered yearly to the ocean from the rivers [1], even if recent studies suggest that this amount might have been underestimated [2]. Estuaries represent only the 0.3% of the Earth's ocean-covered surface and they are a direct link between fresh and saline waters. The delivery of organic matter from rivers to the coastal ocean via estuaries is recognized as an important component of the global C budget [3][4][5][6]. Within the estuaries, a variety of biogeochemical reactions determine the quantity and quality of DOM reaching the oceans. These reactions include flocculation, adsorption on suspended matter, photochemical and microbial processes [2,[7][8][9], and can affect both the concentration and quality of the DOM transported to the coastal area [10].
In estuaries, microbial communities are exposed to quick changes in the environmental conditions. Mixing of seawater and freshwater, transportation of dissolved or suspended organic and inorganic material, as well as hydrological modifications due to climatic conditions or to human actions affect the estuarine biogeochemical properties, inducing specific patterns of microbial abundance, diversity, and activity in these areas [11,12]. [11,12]. Changing salinity is a considerable physiological stress factor for freshw ganisms and affects the performance and functioning of microbial communities [ ferences in microbial abundances have been therefore observed along salinity gra various sites, such as the Rhone [14], Columbia [15], and Adour estuaries [16]. Th DOM by the estuarine microbial community can be affected by several enviro parameters other than salinity, such as temperature, discharge (changing the re time of the water), nutrient availability, and the amount and quality of the riverin transported to the estuary. Apple et al. [17] suggested that bacterial carbon met can be strongly affected by local environmental conditions, and seasonal differ estuarine DOC uptake has been observed [17][18][19].
With this study we aim at investigating the biological removal and transform riverine DOM in different seasons in an estuarine system. The Arno River estuar was chosen for this study because its DOM temporal dynamic is well known and clear seasonal cycle of DOM, both in terms of concentration and quality, and of trophic prokaryotes abundance (HPA) [20]. Two incubation experiments were car using natural DOM and heterotrophic prokaryotes community collected at the es two contrasting seasons, spring (March) and autumn (September). In spring, acco the predominance of terrestrial DOM and the lower DOC concentration observe river [20], we expect a lower percentage of labile DOM, which combined with temperature might lead to lower removal rates than in autumn. To test this hyp the concentration of DOC and the HPA were followed through time in order to inv the DOM removal rates and heterotrophic prokaryotes growth efficiencies. Cha DOM quality were evaluated by measuring its fluorescence (FDOM), which give tion about the main groups of fluorophores present in its pool.

Samples Collection and Incubation Experiment Setup
Surface water samples were collected in the Arno River estuary in March a tember 2015 at intermediate salinities ( Figure 1, Table 1) into 10 L, acid washed, bonate bottles (Nalgene). Temperature and salinity were measured in situ by usin able Hanna 9033 probe (Hanna Instruments Inc., Woonsocket, RI, USA).

March
Septembe Temperature in situ (°C) 11.5 22.2  The incubation experiments were carried out in 2 L, acid washed, polycarbonate bottles (Nalgene) in triplicates. For each replicate, 1.8 L of the estuarine water was filtered through a 0.2 µm pore size filter (Whatman Polycap, 6705-3602 capsules) and inoculated with 0.2 L of the same water, filtered through a 1.2 µm pore size filter (Sartorius, 17,593 K), in order to add the local microbial community. Both filters were extensively rinsed with Milli-Q water prior to their use to avoid DOM contamination. The 1.2 µm filter was chosen because this pore size allows the majority of the heterotrophic prokaryotes to remain in the water used for the inoculum, whereas it should retain most of the autotrophic organisms and the bacterivores [21]. In order to simulate as closely as possible the natural conditions, the incubations were carried out at the in situ temperature (Table 1), and the bottles were kept in the dark to avoid any autotrophic activity not removed by filtration. A headspace of ≈15% of the volume was left empty to keep enough oxygen in the system.
Subsamples for analyses were collected immediately after the inoculum (T 0 ), and at different times during the incubation up to 2 months. The 2 L bottles were open and mixed before each subsampling, which was carried out in sterile conditions to avoid any contamination from ambient microbes.

Analytical Procedures for DOM and Heterotrophic Prokaryotes Abundance
Heterotrophic prokaryotes (HP) were enumerated by flow cytometry. Samples fixed with either glutaraldhyde alone or a mix of paraformaldehyde (PF, 1%) and glutaraldehyde (GL, 0.05%), were stained with SYBR Green (Invitrogen Milan, Italy) and analyzed with an Accuri C6 or a FACSVerse flow cytometer (both from BD BioSciences) equipped with standard laser and filter sets. HP were discriminated from other particles and background by their high green fluorescence, small size, low light scatter, and low red fluorescence, as previously described [22,23].
The samples for DOM analyses were collected into 60 mL, acid washed, polycarbonate bottles (Nalgene) and measured immediately. Total organic carbon (TOC) concentration was measured by high temperature catalytic oxidation using a Shimadzu TOC-Vcsn analyzer following the method reported by Santinelli et al. [24]. From 3 to 5 replicate injections were performed until the analytical precision was lower than 1 µM, and the measurement reliability was assessed twice daily by comparison of data with DOC consensus reference material [25] (CRM Batch #13 Lot 05-13, nominal concentration: 41-44 µM; measured concentration: 42.3 ± 0.9 µM).
To avoid any contamination, subsamples were not filtered; DOC concentration was therefore calculated by subtracting HP carbon biomass from TOC concentration. The HP biomass was calculated assuming a conversion factor of 20 fg C cell −1 [26] and ranged between 0.11 and 1.94 µM C with values <1 µM C in 80% of the samples.
The variation in biomass and in carbon concentration were used to estimate the heterotrophic prokaryotes growth efficiency (HPGE) as follows: Fluorescence excitation-emission matrices (EEMs) were measured using the Aqualog spectrofluorometer (Horiba) following the method reported by Retelletti Brogi et al. [27]. By using the drEEM toolbox [28], the EEMs were elaborated in order to remove and interpolate the Rayleigh and Raman scatter peaks, normalize the fluorescence intensities to Raman units, and carry out parallel factor analysis (PARAFAC). The model was validated by visual inspection of the residuals, split-half analysis (Supplementary Materials Figure S1), and percentage of explained variance (99.5%).

Statistics
The Kruskal-Wallis test (R software) was used in order to test if the variations in the parameters during the incubations were statistically significant. This test was chosen since is a nonparametric test and does not need any distributional assumption [29]. Differences were considered significant for p < 0.005.

Biological Removal of DOC
As expected, the initial conditions of the two experiments were biologically, chemically, and physically different ( Table 1).
It is noteworthy that in September, DOC concentration was ≈1.7 times higher than in March, whereas HPA was ≈2 times lower. These results are in good agreement with the DOM and HPA annual cycle in the Arno River [20]. The accumulation of DOM and the low HPA can be attributed to the presence of recalcitrant DOM, which cannot be consumed by the HP on the short temporal scale. The top-down control (viral lysis or grazing) would also reduce the HPA, leading to DOM accumulation. With this experiment, we investigated if the biological lability of DOM was different in the two periods and if the DOM removed was used for biomass synthesis with the same efficiency. DOC concentration showed an exponential decrease during the two months in both incubations ( Figure 2). It decreased from 179 ± 7 to 147 ± 1 µM in March and from 316 ± 2 to 279 ± 1 µM in September, suggesting an overall DOC removal rate of 16.0 µM month −1 (0.5 µM day −1 ) in March and 18.5 µM month −1 (0.6 µM day −1 ) in September. The removal rates were very similar, despite the difference in the initial DOC concentrations, HPA, and temperature (Table 1), and similar to those reported in a previous study carried out at the Arno River estuary in October 2012 (20 µM month −1 , [30]

Statistics
The Kruskal-Wallis test (R software) was used in order to test if the variations in the parameters during the incubations were statistically significant. This test was chosen since is a nonparametric test and does not need any distributional assumption [29]. Differences were considered significant for p < 0.005.

Biological Removal of DOC
As expected, the initial conditions of the two experiments were biologically, chemically, and physically different ( Table 1).
It is noteworthy that in September, DOC concentration was ≈1.7 times higher than in March, whereas HPA was ≈2 times lower. These results are in good agreement with the DOM and HPA annual cycle in the Arno River [20]. The accumulation of DOM and the low HPA can be attributed to the presence of recalcitrant DOM, which cannot be consumed by the HP on the short temporal scale. The top-down control (viral lysis or grazing) would also reduce the HPA, leading to DOM accumulation. With this experiment, we investigated if the biological lability of DOM was different in the two periods and if the DOM removed was used for biomass synthesis with the same efficiency. DOC concentration showed an exponential decrease during the two months in both incubations ( Figure  2). It decreased from 179 ± 7 to 147 ± 1 µ M in March and from 316 ± 2 to 279 ± 1 µ M in September, suggesting an overall DOC removal rate of 16.0 µ M month −1 (0.5 µ M day −1 ) in March and 18.5 µ M month −1 (0.6 µ M day −1 ) in September. The removal rates were very similar, despite the difference in the initial DOC concentrations, HPA, and temperature (Table 1), and similar to those reported in a previous study carried out at the Arno River estuary in October 2012 (20 µ M month −1 , [30]). Even though the amount of DOC removed in 48 h is different in the two periods, it is interesting to notice that its percentage, with respect to the initial concentration, is very similar (9% in March and 8% in September). These results indicate that ≈8-9% of DOC was labile on a 48 h temporal scale during both periods. In contrast, considering the whole Even though the amount of DOC removed in 48 h is different in the two periods, it is interesting to notice that its percentage, with respect to the initial concentration, is very similar (9% in March and 8% in September). These results indicate that ≈8-9% of DOC was labile on a 48 h temporal scale during both periods. In contrast, considering the whole incubations (i.e., two months), the percentage of the removed DOC was slightly different; 18% of the initial DOC was removed in March and 12% in September.
Our results are in good agreement with incubation experiments carried out in other estuaries [18,31,32]. Raymond and Bauer [18] observed no variation in DOC concentration in the first 24 h and no seasonal variability during incubation experiments, carried out with water from the York River estuary in different seasons; they also found 10% removal of the initial DOC concentration within five days of experiment. Avery et al. [31] reported the results of incubation experiments with water from the Cape Fear River estuary in three different seasons (February, April, and July) and observed similar results in the three periods with 9 ± 4.5% DOC removal during three months. Similar results were also reported by Moran and Sheldon [32] for incubation experiments carried out in five estuaries in the southeastern USA, where 3-12% of DOC was removed in 1-2 months. The environmental conditions can be totally different between estuaries and seasons and it is not easy to explain this general agreement among different studies around the globe. We can hypothesize that this 8-9% is the measurable LDOC, and that there is a fraction which is labile on a shorter temporal scale and therefore not measurable. This highly labile fraction can have different concentration among environments and seasons. However, it is not possible to exclude that the results might be also affected by the experimental conditions, which were similar between our experiment and the ones cited.
HPA showed a similar trend in both March and September: after a lag phase of 24 h, it rapidly increased to reach a maximum after 48 h (Figure 3), in correspondence with the largest DOC decrease. After reaching the maximum of abundance, HPA decreased to reach a minimum at day 10 (March) and day 7 (September). Then, it started to increase again and reached a second maximum after one month in March and 15 days in September, and then it slowly decreased (Figure 3). The marked decrease in HPA observed between 2 and 7 days can be attributed to the inability of the initial HP community to use the DOC that was left. However, we cannot exclude other factors like viral lysis [33] or an incomplete removal of grazers. The increase in abundance observed afterward may suggest the growth of a second community not dominant in the first week, which was able to grow on the remaining DOC or which was less sensitive to grazing and viral lysis. incubations (i.e., two months), the percentage of the removed DOC was slightly different; 18% of the initial DOC was removed in March and 12% in September.
Our results are in good agreement with incubation experiments carried out in other estuaries [18,31,32]. Raymond and Bauer [18] observed no variation in DOC concentration in the first 24 h and no seasonal variability during incubation experiments, carried out with water from the York River estuary in different seasons; they also found 10% removal of the initial DOC concentration within five days of experiment. Avery et al. [31] reported the results of incubation experiments with water from the Cape Fear River estuary in three different seasons (February, April, and July) and observed similar results in the three periods with 9 ± 4.5% DOC removal during three months. Similar results were also reported by Moran and Sheldon [32] for incubation experiments carried out in five estuaries in the southeastern USA, where 3-12% of DOC was removed in 1-2 months. The environmental conditions can be totally different between estuaries and seasons and it is not easy to explain this general agreement among different studies around the globe. We can hypothesize that this 8-9% is the measurable LDOC, and that there is a fraction which is labile on a shorter temporal scale and therefore not measurable. This highly labile fraction can have different concentration among environments and seasons. However, it is not possible to exclude that the results might be also affected by the experimental conditions, which were similar between our experiment and the ones cited.
HPA showed a similar trend in both March and September: after a lag phase of 24 h, it rapidly increased to reach a maximum after 48 h (Figure 3), in correspondence with the largest DOC decrease. After reaching the maximum of abundance, HPA decreased to reach a minimum at day 10 (March) and day 7 (September). Then, it started to increase again and reached a second maximum after one month in March and 15 days in September, and then it slowly decreased (Figure 3). The marked decrease in HPA observed between 2 and 7 days can be attributed to the inability of the initial HP community to use the DOC that was left. However, we cannot exclude other factors like viral lysis [33] or an incomplete removal of grazers. The increase in abundance observed afterward may suggest the growth of a second community not dominant in the first week, which was able to grow on the remaining DOC or which was less sensitive to grazing and viral lysis. These results were surprising because a higher DOC removal was expected in September when initial DOC concentration and temperature were higher. In March, the lower temperature and DOC concentration were expected to reduce the HP growth. These results were surprising because a higher DOC removal was expected in September when initial DOC concentration and temperature were higher. In March, the lower temperature and DOC concentration were expected to reduce the HP growth.

Heterotrophic Prokaryotes Growth Efficiency and CO 2 Production Estimates
Comparing DOC removal rates and the variations in HPA, both experiments can be divided in two periods that seemed to be characterized by different HP communities: in the first period (0-48 h) HP grew fast on the labile DOC (LDOC, removal rates of 9-13 µM day −1 ) while after about 1 week they grew slower on the semilabile DOC (SLDOC, mineralization rates of 0.3-0.4 µM day −1 ).
The DOC removed by HP can follow 3 metabolic pathways: it can be (1) incorporated into biomass, (2) mineralized into CO 2 through respiration, or (3) released as DOC in a different form. HPGE gives an estimate of the percentage of the removed DOC transformed in biomass.
The HPGE estimated in the first period (0-48 h) was 10 ± 0.21% in March and 5 ± 0.10% in September, whereas in the second period (7-10 days to 1 month), it was 7 ± 0.15% in March and 8 ± 0.12% in September. The lower HPGE within the first 48 h in September suggests that the HP community was less efficient in transforming DOC into biomass, and most of the DOC removed (~95%) was respired to CO 2 even with warmer conditions (T~22 • C) than in April. Although these results might seem surprising, they are in agreement with the results of Apple et al. [17], who investigated the temperature dependency of bacterioplankton carbon metabolism and showed an inverse relationship between the temperature and HPGE.
Assuming negligible the DOC released by the HP, the amount of CO 2 produced by respiration was estimated as follows: This calculation suggests that 14.4 to 24.5 µM of CO 2 were released in 48 h by the consumption of LDOC from heterotrophic prokaryotes, leading to the production of 7-12 µM of CO 2 per day. These calculations may, however, be affected by the experimental conditions.

DOM Quality Affects Its Mineralization and Vice Versa
The PARAFAC analysis of the EEMs validated a four component model ( Figures S1 and S2). According to their excitation and emission maxima and to the comparison with similar components within the OpenFluor database [34], these components were attributed to different groups of fluorophores. C1 showed the excitation and emission maxima at 310 and 410 nm, respectively, and its signal was attributed to microbial humic-like compounds [20,[35][36][37]. C2 showed the excitation peak at 360 nm, and the emission maximum at 460 and its signal can be attributed to fulvic-like fluorescence [20,38]. The excitation and emission maxima of C3 and C4 were 290/340 and 280/330 nm, which are typical of protein-like compounds fluorescence [20,35,37]. Because both C3 and C4 represent protein-like compounds, hereinafter they are reported together as C3 + C4 in order to show the overall changes in protein-like fluorescence during the incubations.
In order to evaluate the differences in DOM quality between the two periods, the EEMs were compared for the two experiments at T 0 ( Figure 4).
About 30% of the fluorescence intensity was due to the microbial humic-like component in both experiments (29 ± 0.1% and 28 ± 0.1% in March and September, respectively). The percentage of fulvic-like component in March (36 ± 0.1%) was twice that in September (18 ± 0.2%), whereas the percentage of protein-like fluorescence was markedly higher in September (54 ± 0.1%) than in March (35 ± 0.1%). These results agree well with those reported within the river [20].
Fluorescence and DOC data clearly indicate a marked difference in the DOM in the two periods; in March, DOM has a low concentration and a similar percentage of humic, fulvic, and protein-like FDOM, although the microbial humic FDOM was slightly reduced with respect to the other components. In September, DOM concentration is twice than in March, and protein-like substances dominated the fluorescence signal. It is possible that changes in DOM pool, together with the different environmental conditions (i.e., winter vs. autumn), can favor HP communities with different growth efficiencies or that according to the type of DOM, the HP community can use the available DOM mainly for biomass synthesis (March) or respiration (September).
tributed to different groups of fluorophores. C1 showed the excitation and emission maxima at 310 and 410 nm, respectively, and its signal was attributed to microbial humic-like compounds [20,[35][36][37]. C2 showed the excitation peak at 360 nm, and the emission maximum at 460 and its signal can be attributed to fulvic-like fluorescence [20,38]. The excitation and emission maxima of C3 and C4 were 290/340 and 280/330 nm, which are typical of protein-like compounds fluorescence [20,35,37]. Because both C3 and C4 represent protein-like compounds, hereinafter they are reported together as C3 + C4 in order to show the overall changes in protein-like fluorescence during the incubations.
In order to evaluate the differences in DOM quality between the two periods, the EEMs were compared for the two experiments at T0 (Figure 4). The fluorescence intensities of the components were followed during the incubations in order to investigate changes in DOM quality due to HP activity ( Figure 5). All of the fluorescence components showed no significant variations in the first 48 h of the experiment, suggesting that the DOC removed in this time frame was not fluorescent. In contrast, large changes in fluorescence were observed afterwards.
In both experiments, between 2 and 7/10 days all of the components showed a decrease in intensity, in correspondence with the decrease in HP abundance. From day 7/10 to the end of the experiments, FDOM showed a different trend in the two incubations ( Figure 5). In March, all of the components showed a significant decrease until the end of the experiment; in September, the microbial humic-like component (C1) markedly increased, reaching a fluorescence intensity >1.5 times that measured at T 0 ; the fulvic-like component (C2) followed a trend similar to that of the HPA, increasing between 7 and 15 days and then decreasing again; the protein-like components (C3 + C4) slowly decreased towards the end of the experiment.
The overall effect of the HP community on FDOM during the experiments can be seen by the subtraction of the EEMs measured after two months from the EEMs measured at T 0 ( Figure 6). In March, the HP reduced all of the FDOM components, whereas in September they mainly released microbial humic-like compounds and removed proteinlike substances ( Figure 6).
The percentage of each component on the total fluorescence was calculated at the end of the experiment and compared with that calculated at T 0 (Figure 7). This calculation showed that, in March, the relative weight of the components on the composition of FDOM did not change over the course of the incubation, meaning that all of the components were proportionally removed. After two months of incubation, the percentage of the components were 31 ± 0.1%, 39 ± 0.2%, and 31 ± 0.1% for C1, C2 and C3 + C4, respectively, which are indeed not significantly different from those calculated at T 0 . In September, a change in the relative FDOM composition was observed between the beginning and the end of the experiment. After two months, the microbial humic-like component (C1) increased from 28 ± 0.1% to 39 ± 0.05%. The protein-like fraction (C3 + C4) decreased from 54 ± 0.1% to 45 ± 0.1%, whereas the fulvic-like component (C2) remained unchanged (16 ± 0.1% after two months). In both experiments, between 2 and 7/10 days all of the components showed a decrease in intensity, in correspondence with the decrease in HP abundance. From day 7/10 to the end of the experiments, FDOM showed a different trend in the two incubations ( Figure 5). In March, all of the components showed a significant decrease until the end of the experiment; in September, the microbial humic-like component (C1) markedly increased, reaching a fluorescence intensity >1.5 times that measured at T0; the fulvic-like component (C2) followed a trend similar to that of the HPA, increasing between 7 and 15 days and then decreasing again; the protein-like components (C3 + C4) slowly decreased towards the end of the experiment.
The overall effect of the HP community on FDOM during the experiments can be seen by the subtraction of the EEMs measured after two months from the EEMs measured at T0 (Figure 6). In March, the HP reduced all of the FDOM components, whereas in September they mainly released microbial humic-like compounds and removed protein-like substances ( Figure 6). The percentage of each component on the total fluorescence was calculated at the end of the experiment and compared with that calculated at T0 (Figure 7). This calculation showed that, in March, the relative weight of the components on the composition of FDOM did not change over the course of the incubation, meaning that all of the compo- tively, which are indeed not significantly different from those calculated at T0. In September, a change in the relative FDOM composition was observed between the beginning and the end of the experiment. After two months, the microbial humic-like component (C1) increased from 28 ± 0.1% to 39 ± 0.05%. The protein-like fraction (C3 + C4) decreased from 54 ± 0.1% to 45 ± 0.1%, whereas the fulvic-like component (C2) remained unchanged (16 ± 0.1% after two months). The increase in humic-like fluorescence due to the HP activity is well documented in the literature [39][40][41][42][43], and supports that humic-like substances do not represent only terrestrial DOM. The increase in humic-like fluorescence due to the HP activity is well documented in the literature [39][40][41][42][43], and supports that humic-like substances do not represent only terrestrial DOM.
The different FDOM trends observed in the two incubations can be attributed to (i) a different composition of the microbial community, (ii) differences in the quality of DOM, or, most probably, by a combination of these two factors. It was reported that biological degradation of DOM is carried out by phylogenetically diverse communities, whose composition has been shown to be affected by the quality and quantity of the available DOM [44,45]. Different groups of bacteria have different capability to use different size-fractions of DOM. Low molecular weight compounds (LMW < 600 Da) can be taken up easily across the cell membrane [46], whereas molecules bigger than 600 Da require extracellular enzymes that not all of the bacterial taxa can produce, as observed by Berlemont and Martiny [47]. It is known that microbial communities can also release FDOM [48]. Goto et al. [41] observed the production of refractory humic-like FDOM by three strains of marine bacteria growing on a labile substrate (glucose). According to their results, the three strains produced different DOM with different efficiency, and proposed specific mechanisms of humic-like FDOM production according to the strain. A recent study highlighted how the microbial community can release humic-like refractory DOM when growing on phytoplankton-derived DOM [40]. This observation supports the results of this study since the increase in humic-like fluorescence was observed in September, after the spring/summer phytoplankton blooms, when indeed a much higher percentage of protein-like fluorescence was present than in March, suggesting a higher abundance of autotrophic DOM. The similar DOC removal rates in September with respect to March, despite the much higher DOC concentration, can be explained by the HP producing refractory humic-like compounds as shown by both Kinsey et al. [40] and Goto et al. [41]. Some studies, combining PARAFAC results with more detailed molecular analysis (FT-ICR-MS) showed that FDOM is a good tracer for the bulk DOM pool [49]. However, FDOM still represents only a fraction of the DOM and a more detailed molecular characterization of all the compounds, including the nonfluorescent ones, could give more insight into its changes and its biologically labile fraction [50].

Summary and Conclusions
The results of this study show similar DOC removal rates in the Arno River estuary in spring and autumn, despite the marked difference in DOC initial concentration, HPA, and environmental conditions (i.e., temperature). In both seasons, 8-9% of the initial DOC was removed after 48 h, suggesting that it is labile on the short temporal scale. The estimate of the HPGE growth efficiency highlighted that the HP community can use the removed DOC differently (i.e., growth vs. respiration). A higher percentage of DOC was used for biomass synthesis in March than in September. Fluorescence data suggest a different quality of the initial DOM, with a higher percentage of protein-like compounds and a lower percentage of terrestrial fulvic-like compounds in September than in March. These differences might stimulate the growth of different HP communities with different growth efficiencies. In turn, the HP community affects the FDOM differently during the incubations. In September, the HP removed the protein-like compounds and released a large amount of humic-like compounds that may inhibit a further DOM uptake. These results show that, in absence of grazers the HP community is able to remove protein-like DOM, while the accumulation of protein-rich DOM coincident with a low HPA might be attributed to a top-down control of HP by grazers or by viruses. In March, the HP removed all of the FDOM components consistently with the DOC decrease without apparently altering the relative fluorescence intensity of the components in the DOM pool. The results of this study point to a high complexity of estuarine systems that needs to be investigated further. The similar DOC removal rates, together with the low growth efficiency of September's HP community, despite DOM being more concentrated and enriched in protein, and despite the higher temperature, was an unexpected result and raises intriguing questions about the cross influences of microbial diversity, physical constrains, and DOM properties. A detailed study on the composition of the HP community at different time of the year, together with more detailed incubation experiments covering all seasons, is needed to get more insights into these processes. Moreover, a molecular characterization of the DOM during the experiments would help in the interpretation of these results and on predicting the biological lability of DOM in different conditions. Supplementary Materials: The following are available online at https://www.mdpi.com/2077-131 2/9/2/172/s1, Figure S1: Results from the PARAFAC split-validation test made by using drEEM 4.0 Matlab toolbox, Figure

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.