Changing Isotopic Food Webs of Two Economically Important Fish in Mediterranean Coastal Lakes with Di ﬀ erent Trophic Status

: Transitional waters are highly productive ecosystems, providing essential goods and services to the biosphere and human population. Human inﬂuence in coastal areas exposes these ecosystems to continuous internal and external disturbance. Nitrogen-loads can a ﬀ ect the composition of the resident community and the trophic relationships between and within species, including ﬁsh. Based on carbon ( δ 13 C) and nitrogen ( δ 15 N) stable isotope analyses of individuals, we explored the feeding behaviour of two ecologically and economically important omnivorous ﬁsh, the eel Anguilla anguilla and the seabream Diplodus annularis , in three neighbouring lakes characterised by di ﬀ erent trophic conditions. We found that A. anguilla showed greater generalism in the eutrophic lake due to the increased contribution of basal resources and invertebrates to its diet. By contrast, the diet of D. annularis, which was mainly based on invertebrate species, became more specialised, focusing especially on polychaetes. Our results suggest that changes in macroinvertebrate and ﬁsh community composition, coupled with anthropogenic pressure, a ﬀ ect the trophic strategies of high trophic level consumers such as A. anguilla and D. annularis . Detailed food web descriptions based on the feeding choices of isotopic trophospecies (here Isotopic Trophic Units, ITUs) enable identiﬁcation of the prey taxa crucial for the persistence of omnivorous ﬁsh stocks, thus providing useful information for their management and habitat conservation.


Introduction
Transitional waters are extremely complex ecosystems [1][2][3]. The Water Framework Directive of the European Communities (European Communities, 2000. Directive 2000/60/EC of the European Parliament and of the Council of 23 October 2000) defines them as "superficial bodies of water near the mouths of rivers which have a partially saline character due to their proximity to coastal waters, but which are substantially influenced by freshwater flows". Their high productivity provides habitats, refuge areas and food sources for a wide range of aquatic animals from resident brackish to freshwater and marine migratory species [4]. Transitional waters support important ecosystem services, including good water quality, fisheries, aquaculture and tourism, as well as agricultural activities in their watersheds [5]. Anthropic activities expose these ecosystems to continuous internal and external disturbance [2,[6][7][8], including nitrogen (N) pollution arising from agricultural and urban Isotopic Trophic Units are defined as groups of individuals with similar isotopic signatures occupying the same position in the δ 13 C-δ 15 N niche space [31].
We studied the diet of each population in detail without excluding a priori any food source in the area. We hypothesised that differences in invertebrate and fish assemblages across lakes with differing trophic status could determine patterns of variation in the trophic niche width of the two fish species depending on their specific feeding habits. Specifically, we sought to verify whether a lower abundance and diversity of species at higher trophic levels caused A. anguilla to become more generalist and D. annularis to become more specialized.

Study Area
The samplings were carried out in three neighbouring Mediterranean brackish costal lakes located on the Tyrrhenian coast of central Italy (42 • 28 00" North-12 • 51 00" East): Lake Caprolace, Lake Fogliano and Lake Sabaudia (  The lakes are affected by various forms of anthropogenic disturbance related to organic and inorganic nitrogen inputs from urban treated sewage, livestock farming and agricultural activities, which are widespread in the surrounding areas [2,3].
Lake Caprolace and Lake Fogliano (hereafter respectively LP and IP), characterised respectively by low and intermediate levels of eutrophication [3,47], are Sites of Community Importance (SCIs) located within the Circeo National Park (Lazio).
Lake Caprolace does not receive water inputs from the hinterland, while Lake Fogliano is affected by nutrient inputs from both the River Rio Martino and the livestock breeding activities practised in the surrounding areas. The annual concentration of Chlorophyll a was generally lower in Caprolace (2.1 ± 0.4 µg/L) than Fogliano (5.8 ± 1.2 µg/L) in 2006-2010.
Lake Sabaudia, the southernmost lake (hereafter HP), is affected by the highest anthropogenic pressure [3], mainly due to runoff from both the city of Sabaudia and cultivated fields in the surrounding areas as well as fishing and mussel farming. In this lake, freshwater inputs are present throughout the year.
Annual algal biomass and Chlorophyll a concentrations in this lake vary from 10.2 to 40.9 µg/L, with an average recorded value in 2006-2010 of 24.2 ± 6.15 µg/L. Further details regarding the study area can be found in Santoro et al. [2] and Jona-Lasinio et al. [3].

Field Collections
Samples of basal resources (primary producers and detritus), invertebrates and fish were collected in 4 sites per lake between April and May 2012, when primary productivity and invertebrate abundances were high. The sampling sites within each lake were selected from areas with heterogeneous physical and biotic characteristics and a range of anthropogenic impacts deriving from the surrounding areas [2,3]. The sampling sites were located at the northern and southern ends of each lake, and both on the landward and seaward sides (see also Santoro et al. [2]). Macrophytes, algae, and detritus samples were collected by hand and invertebrates by Van Veen grab (volume: 3.5 L) in three replicates per sampling site. The dominant macrophytes were Ruppia sp. and Cymodocea nodosa (Ucria) Ascherson, while the macroalgae were represented by taxa of the genera Chetomorpha, Chondria, Gracilaria, Rytiphloea and Ulva. The detritus was mostly composed of fragments of dead leaves delicately scraped to remove any epibionts and rinsed in distilled water. Phytoplankton samples were collected using a plankton net (20-µm mesh size) and concentrated by centrifugation (2000 rpm for 20 min).
Samples of fish were collected once a day for 3 days in each site. In order to collect pelagic, benthic, resident and migratory fish species, fish samples were collected using fixed weirs and fishing traps placed on the bottom. The fishing traps, made of very fine mesh (0.5 cm), were 1.5 m in diameter at the mouth and were composed of four consecutive chambers of decreasing diameter with a total length of 3.6 m. In addition to A. anguilla and D. annularis, the sampled fish community included the sand smelt Atherina boyeri (Risso, 1810), black goby Gobius niger (Linnaeus, 1758) and the mullets Chelon ramada (Risso, 1827), C. aurata (Risso, 1810), C. saliens (Risso, 1810) and C. labrosus (Risso, 1827), which are known to be prey species of A. anguilla and D. annularis [19,24,25,42,43]. Further fish samples included species belonging to the Sparidae, Scorpaenidae, Clupeidae, Cyprinodontidae, Blenniidae and Belonidae, Gobiidae, Labridae, Moronidae, Mugilidae, Soleidae and Syngnathidae families. Standard length measured in centimetres was recorded for each fish specimen. For each fish species, individuals of different sizes were collected in order to reduce the effects of size variability on isotopic signals. From the sampled fish specimens, including A. anguilla and D. annularis, samples of dorsal white muscle were taken. This tissue provides a long-term (several months) integrated indicator of food sources due to its slow turnover with respect to other tissues (e.g., liver and blood) [32].
After collection, all samples were transported to the laboratory, where specimens were sorted, counted, and identified to the lowest possible taxonomic level and processed for the stable isotope analysis.

Stable Isotope Analysis (SIA)
Samples were individually stored at −80 • C and freeze-dried for 24 h. Fish specimens were considered individually for isotopic analysis. Muscle samples were also taken from large invertebrates such as crustaceans, for which the tissue was taken from the claws, and bivalves and sea snails, whose tissue was taken from the feet [7]. When present, shells, valves and other exoskeletal parts of animals were removed under dissection microscopes in order to avoid tissue acidification before the stable isotope analysis. For small invertebrates (such as amphipods and polychaetes), the whole body was used. Samples were individually analysed. Plankton biomass was analysed as a whole due to the difficulty of obtaining sufficient biomass for isotopic analysis.
Before the stable isotope analysis, each sample was homogenised to a fine powder using a ball mill (Mini-Mill Fritsch Pulverisette 23: Fritsch Instruments, Idar-Oberstein, Germany). When necessary, samples were pre-acidified using 1M HCl according to the drop-by-drop method [48] in order to eliminate inorganic carbon and re-dried (60 • C) for 72 h to remove the remaining moisture. δ 15 N signatures were measured in un-acidified powders to prevent acidification from interfering with the nitrogen analysis [30,49].
Aliquots of 0.25 ± 0.10 mg for the animals and 2.00 ± 0.10 mg for basal resources were placed into tin capsules for C and N stable isotopic analysis (SIA). Each sample was analysed in two replicates. The analyses were carried out using a continuous flow mass spectrometer (IsoPrime100, Isoprime Ltd., Cheadle Hulme, United Kingdom) coupled with an elemental analyser (Elementar Vario Micro-Cube, Elementar Analysensysteme GmbH, Germany).
The isotopic signatures of each sample were expressed in δ units (δ 15 N; δ 13 C) as parts per thousand (% ) deviations from international standards (atmospheric N 2 for N; PD-belemnite (PDB) carbonate for C), in accordance with the formula: , where X is 13 C or 15 N and R is the corresponding ratio of heavy to light isotope for the element ( 13 C/ 12 C or 15 N/ 14 N). Outputs were standardised with the internal laboratory standard Caffeine IAEA-600 (C 8 H 10 N 4 O 2 ). Measurement errors were found to be typically smaller than 0.05% .

Data Analysis
Differences between lakes in terms of animal community composition (considering both fish and benthic invertebrates) were tested using contingency tables based on chi-square (χ 2 ) tests, Monte Carlo permutation tests and the associated Cramer's V index (a measure of the strength of association among communities; Past 3.0 software package). Specimens collected in each sampling site (and replicates) were grouped by type or taxon (respectively for basal resources and animals) for each lake.
The Shannon diversity index (Hs) of invertebrate fauna for each lake was calculated at family level considering a total abundance of the taxa collected in each lake. Given that assessing Shannon diversity is only possible at the level of equal identification of all taxa, the few individuals belonging to the Gastropoda, Oligochaeta, Nematoda, and Nemertea classes (together accounting for less than 1.5% of total fauna) were excluded from the Shannon diversity index computation. Hutcheson's diversity t-test and the associated bootstrap procedure (9999 replicates), both available in the Past 3.0 software package, were applied to Hs values to test for significant differences [51]. Hutcheson's diversity t-test is a modified version of the classic t-test and is based on comparison of Hs variances. The t statistics of Hutcheson are defined as: which follows Student's t distribution. In the equation, i and j referred to the invertebrate communities of the lakes in paired comparisons, Hs represents the Shannon diversity index and var(Hs) its variance.
The isotopic values of collected organisms were used to reconstruct the diets of the eel Anguilla anguilla and the annular seabream Diplodus annularis in each lake. The diets were estimated on the basis of the Isotopic Trophic Unit (ITU) method [31]. The isotopic signatures of single basal resources, invertebrates, and fish were represented in the bi-dimensional isotopic space ( Figure S1). This was subdivided into squares (ITUs) corresponding to 1 × 1% δ 15 N and δ 13 C values, starting from the lowest δ 13 C value in the dataset and a δ 15 N value of zero. The ITUs were thus identified and labelled ( Figure S1).
The diets of each ITU containing individuals of the two fish species were calculated by means of Bayesian Mixing Models (R software ver. 3.5.3, SIMMr package) [52] considering a Trophic Enrichment factor (TEF) of 3.4 ± 1.0% for δ 15 N and 1.0 ± 0.5% for δ 13 C [18,37,49,[53][54][55][56] and uninformative priors. These TEF values (expressed as mean ± standard deviation) are considered a robust and widely applicable assumption in the presence of multiple trophic pathways and different types of food sources [37,56]. For all SIMMr models, we ran three Markov Chain Monte Carlo chains of 300,000 iterations each with a burn-in of 200,000 and a thinning rate of 100 iterations. We assumed that all incoming food items had the same probability of being included in the consumer's diet. The model considers both variance in the isotopic signatures of the resources and uncertainty regarding the trophic enrichment of the consumer (TEF). The model results were expressed in the form of a probability distribution of plausible contribution values. The central tendency values of the distribution (mode, mean, median) allowed us to identify the most important food sources, while the upper and lower limits of the credibility ranges (CI: 50%, 75%, 95%) revealed the range of feasible contributions. The pool of food sources was selected based on the mixing model outputs in accordance with Rossi et al. [31]. Since the A. anguilla and D. annularis diet was obtained by starting from the foraging choice of each individual, the overall contribution of some food sources, important at individual level (>5%), could be relatively small (<5%) if considered at the population level (see also [31]). In order to obtain detailed information on the diet of the eel and the seabream these contributions were also considered.
Individuals other than A. anguilla and D. annularis, including basal resources and invertebrates, were excluded from ITU-consumers (but not from potential ITU food sources) before performing the Bayesian mixing models. This was done in order to correctly estimate the diet of A. anguilla and D. annularis. The set of potential ITU food sources was considered on the entire δ 13 C axis and within a given range on the δ 15 N axis, i.e., within ±3.4% (the TEF) of the value of the consumer [31]. The Bray-Curtis similarity index (BC), based on the contribution of each resource to the diet of the two fish, was also calculated in order to quantify the diet similarity among lakes [36,56]. BC is expressed as proportional similarity ranging from 0, when no common food sources are found for the compared groups, to 1, when the compared groups have the same food sources in the same proportions [36,56].
The symmetric overlap in resource use [57][58][59] was measured in accordance with the Pianka equation [59]: where the Pianka index (O jk ) represents a symmetric measure of overlap between species j and k, and p ij and p ik are the proportional contributions of any given resource i used by species j and species k. The Pianka index ranges from 0 (overlap absent) to 1 (complete overlap).
Chesson's selectivity index [60] was calculated for each food item to determine possible preferences for particular food sources among those offered: where α i is Chesson's selectivity index, m is the number of food source types, r i is the proportion of food type i in the diet and n i is the proportion of food type i in the environment. The value of α i ranges from 0 to 1, with 0 indicating complete avoidance, values above 1/m indicating preference and 1 indicating absolute preference [61]. Since consumer isotopic ratios provide an integrated measure of prey assimilated over time, we hypothesized that the composition of the taxon in each lake did not vary considerably over the course of a season. Therefore, the Chesson index based on the relationship between assimilated prey and its abundance in the environment could measure the selectivity of food products with a good approximation. χ 2 tests were performed to test for differences between lakes in terms of the relative abundance of fauna and differences between food sources in terms of their proportional contribution to the diet of the fish population in each lake. Although it is not possible to establish a theoretical expected value, a χ 2 test was performed to test for possible differences between IP and LP and between HP and LP, considering the least polluted lake as the reference value.
Differences between lakes in both the δ 15 N and δ 13 C isotopic signatures of basal resources and fauna were tested by one-way ANOVA for comparisons between normal distributions (Shapiro-Wilk normality test, p > 0.05) while the Mann-Whitney with Bonferroni correction in cases of multiple comparisons and Kruskal-Wallis tests were used if non-normality was observed (Shapiro-Wilk normality test, p-value < 0.05). Levene's test for variances was used to test for differences within and between lakes in the δ 13 C variance of primary producers. Kruskal-Wallis tests and associated Mann-Whitney pairwise comparisons were also used to compare the proportional contribution of food items to the diet of A. anguilla and D. annularis.
The niche metrics for both species in each lake were also calculated [62][63][64]. These metrics, originally proposed by Layman et al. [62] for application at community level, can be used at population level to obtain information about trophic diversity within a single population [35,63,64]. These included the ranges (highest to lowest) of δ 13 C (Carbon Range, CR) and δ 15 N (Nitrogen Range, NR) values. CR provides information about the variety of food sources exploited by the population (i.e., its trophic generalism), while NR indicates the number of trophic levels (i.e., degree of omnivory) of the population. The isotopic niche widths of both A. Anguilla and D. annularis were calculated as SEAc (Standard Ellipse Area corrected by degree of freedom) using R software ver. 3.5.3, SIBER analysis package [64,65]. The SEAc encompasses the core (about 40%) of the population's isotopic observations. This is a solid metric for comparing the isotopic niche of populations regardless of sample size and any isotopic outliers in the data [62,64]. Linkage density (L/S) was measured as the average number of feeding links (L) per ITU (S). Finally, based on the proportional contribution of each food source, the trophic niche width (TNW) of each population was measured as the diversity of resources consumed (Hs) by each population and compared among lakes. If not specified otherwise, the results are reported as mean ± standard error (s.e.).

Community Composition and Isotopic Signatures
A total of 8752 samples comprising basal resources, invertebrates and fish were collected from the three lakes, 8645 (148 taxa) of which were invertebrates and fish ( Table 1, Table S1).
Invertebrate abundance was lower in IP than the other two lakes (Table 1, paired-χ 2 test, χ 2 at least 20.46, p-value always <0.0001, Table S2). The composition of both the invertebrate and fish community also varied (contingency table, χ 2 at least 170.2, p always <0.001, Cramer's V at least 0.46, Tables S2 and S3). The abundance of some taxa, such as Decapoda and Anthozoa, decreased, while that of others (such as Amphipoda) increased with the pollution level of the lake (paired-χ 2 test, χ 2 at least 48.36, p-value always <0.0001). The number of fish taxa varied, i.e., 23 in LP, 17 in IP and 7 in HP. The relative abundance of fish differed between lakes (paired-χ 2 test, χ 2 at least 34.21 p-value always <0.0001) and was lowest in HP. The standard length of Anguilla anguilla was lower in HP (35.10 ± 3.63 cm) than the other two lakes (48.81 ± 6.07 cm in LP, 47.66 ± 2.12 cm in IP) (Mann-Whitney test with Bonferroni correction in cases of multiple comparisons, U = 13.0, p-value always <0.05). Similarly, Diplodus annularis had an average standard length of 6.25 ± 0.14 cm in HP, which was lower than LP (8.76 ± 0.31 cm) and IP (9.43 ± 0.41 cm) (Mann-Whitney test with Bonferroni correction in cases of multiple comparisons, U = 2.0, p-value <0.05). Among the basal resources, detritus showed depleted δ 13 C values, while primary producers were δ 13 C-enriched ( Figure 2). δ 13 C-enrichment was also observed in pelagic fish with specialist diets such as Atherina boyeri (δ 13 C = −15.34 ± 0.08% in LP, −17.45 ± 1.35% in IP and −16.14 ± 0.90% in HP). Since neither the mean nor the variance (σ 2 ) of δ 13 C in the primary producers differed significantly either within each lake or between lakes (one-way ANOVA and associated Levene's test for homogeneity of variances, F at least 0.1412, p-value always >0.05), we concluded that the presence of a salinity gradient within a lake could not have an effect on the isotopic variability of the baseline. δ 15 N values of primary producers increased with eutrophication (one-way ANOVA, F: 5.80, p <0.01).
The isotopic differences observed in basal resources reflected those observed in the whole community (Table 1, Figure 3; Kruskal-Wallis, Hc at least 127.1, p-value < 0.001; for δ 13 C Mann-Whitney with Bonferroni correction for multiple comparisons, U at least 45,590, p-value always <0.001 and for δ 15 N Mann-Whitney, U at least 421.57 with Bonferroni correction for multiple comparisons, p-value always <0.001, Figure 3) and in A. anguilla and D. annularis.
The δ 13 C and δ 15 N isotopic signatures of Anguilla anguilla and Diplodus annularis differed between lakes (Kruskal-Wallis, Hc at least 35.5, p-value <0.001). δ 13 C values were higher in the least polluted lake ( Table 2; Mann-Whitney with Bonferroni correction for multiple comparisons, U at least 3392.5, p-value always <0.001), while δ 15 N values increased with eutrophication ( Table 2; Mann-Whitney with Bonferroni correction for multiple comparisons, U at least 555, p-value always <0.001).
Specifically, in A. anguilla, the more generalist of the two species, δ 13 C values reflected the shift of inputs from marine to terrestrial origin passing from the least to the most eutrophic lake (Figure 2 and Tables 1 and 2).

Niche Metrics and Diet of Anguilla anguilla
The isotopic signatures and niche metrics of Anguilla anguilla varied among lakes (  Table 2; Mann-Whitney with Bonferroni correction in cases of multiple comparisons, U at least 0.1, p-value always <0.001). The Carbon Range increased with eutrophication ( Figure 4, Table 2) and the largest Nitrogen Range was observed in the eutrophic lake ( Table 2). Overall, no correlation between the body length and δ 13 C (% ) of Anguilla anguilla was observed in any lake (Pearson correlation, p >0.05).
Anguilla anguilla had 5 ITUs in LP and HP and 12 ITUs in IP, where the eel-resource ITU linkage density was highest ( Table 2). ITU-based mixing models showed no differences between lakes in terms of the overall contribution of invertebrates to the eels' diet ( Figure 5a).
By contrast, the consumption of basal resources increased and piscivory decreased with increasing levels of pollution (i.e., from LP to HP; Table 3 Figure 5a).
A. anguilla showed a generalist diet including 20 different categories of food source (Table 3, Figures 6 and 7). Some of these were common to the three lake populations (e.g., Actinopterygii, Bivalvia, Gastropoda, Decapoda and Polychaeta) but their consumption varied. The Bray-Curtis index (BC) applied to diet showed a lower similarity between the HP population and the others (76% similarity between LP and IP vs. 41% between LP and HP, and 54% between IP and HP). Specifically, in LP the diet of A. anguilla was mostly based on Actinopterygii (34.76% ± 1.90), Decapoda (27.84% ± 4.60) and Gastropoda (13.06% ± 0.40), in IP on Actinopterygii (30.65% ± 0.50) and Decapoda (26.13% ± 1.90) and in HP on Polychaeta (30.41% ± 0.60), Actinopterygii (12.18% ± 1.50), Bivalvia (10.35% ± 0.90), detritus (9.03% ± 0.50) and Decapoda (9.11% ± 1.90) (Table 3, Figure 6).  letters (a,b,c) within panels indicate differences between lakes in the contribution of food sources to the diet (χ 2 -test, p-value <0.001). Table 3. Proportional contribution (in %) of food sources to the diet of A. anguilla in each lake, obtained from ITU-based mixing models. LP: low anthropogenic pressure, IP: intermediate anthropogenic pressure, HP: high anthropogenic pressure. The contribution of each food source is reported as the mean (±s.e.). "Taxa" indicates the number of taxa belonging to the respective group in the diet of A. anguilla. The overall contribution of basal resources, invertebrates and fish is reported as the mean (%) ± s.e. Different superscript letters (a,b,c) indicate differences between lakes in the contribution of categories of food sources to the diet (χ 2 -test, p-value <0.05). For details please refer to the methods section.    The difference in resource use was associated with a difference in trophic niche width (Figures 2-6), which increased with eutrophication (TNW: 1.81, 2.06 and 2.32 in LP, IP and HP respectively), with significant differences between HP and LP (bootstrap comparison among populations, p <0.0001).

Niche Metrics and Diet of Diplodus annularis
The δ 13 C and δ 15 N of Diplodus annularis varied across lakes (Figures 2-4, Kruskal-Wallis for both δ 13 C and δ 15 N, Hc at least 35.5, p-value <0.001). In addition, no significant correlation between the body length and δ 13 C (% ) of Diplodus annularis was observed in any lake (Pearson correlation, p >0.05).
δ 15 N increased with the level of pollution (Table 2, Figure 4). The Carbon Range was highest in LP while no differences in Nitrogen Range were observed between lakes ( Table 2). The seabreams had four ITUs in LP and two ITUs in both IP and HP (Table 2). No differences in the linkage density (L/S) of the ITUs or in SEAc were observed between lakes (Table 2).
Mixing models based on the single ITU values showed similar average contributions of basal resources, invertebrates and fish in the lakes (Figure 5b; paired-χ 2 test, p-value always >0.05).
However, when taxa in each category were distinguished, diet similarity between LP and HP was 55%, while between IP and both LP and HP it was 61% (Bray-Curtis index, BC).

LP
IP HP The trophic niche was significantly narrower in HP than LP and IP ( Figure 6; TNW: 1.98 vs. 2.15 and 2.28 respectively; bootstrap comparison among populations, p <0.05). The niche overlap between D. annularis and A. anguilla decreased from the eutrophic to the unpolluted lake (O = 0.95, 0.87 and 0.84 in HP, IP and LP respectively).

Discussion
Our results indicate that anthropogenic inputs affected the composition and abundance of the lake animal community. Specifically, the diversity of fish decreased, and trophic choices of the eel Anguilla anguilla and the seabream Diplodus annularis changed, with increasing eutrophication. A. anguilla, which was more piscivorous than the seabream at low and intermediate eutrophication, increased its consumption of invertebrates and basal resources. On the other hand, the seabream, which fed more on invertebrates, increased its preference for polychaetes. The high selectivity for polychaetes in the highly eutrophicated lake could be due to the facilitated capture of these preys and their good contribution to the energy supply of their predators compared to other aquatic invertebrates [66][67][68]. Changes in the feeding choices of the two fish species resulted in increased interspecific niche overlap, suggesting that eutrophication may have strong bottom-up effects on interspecific interactions [69,70].
Previous research demonstrates that the δ 15 N values of the aquatic biota reflect anthropogenic nitrogen inputs from the surrounding terrestrial areas [7,10,71,72]. In our study, increased eutrophication was associated with higher δ 15 N in the fish community, as previously observed by Santoro et al. [2] and Jona-Lasinio et al. [3] for invertebrates and primary producers respectively. Our results indicate that in transitional waters, individual fish δ 15 N could be a useful indicator of anthropogenic N transfer along food webs [10,73], while the range of δ 13 C in the fish population could reflect the diversity of C inputs, emphasising the need for ecological monitoring in these productive ecosystems.
The distinct carbon isotopic signatures (δ 13 C) of primary producers (e.g., terrestrial vs. aquatic vegetation) allowed us to discern the origin of the organic matter contributing to the nutrient pool of water bodies [7,33,34,74]. Specifically, depleted δ 13 C organic matter values indicated the contribution to the organic matter pool of allochthonous (terrestrial) carbon, while enriched δ 13 C indicated autochthonous primary production, as also observed in other aquatic ecosystems [7,34,35,74]. Although it was not in the remit of this study to investigate the cause, carbon enrichment was also visible in pelagic and strictly specialist species such as the sand smelt, A. boyeri. In addition, neither the mean nor the variance of δ 13 C of primary producers differed significantly between sites within the same lake or between different lakes. This allowed us to exclude possible interferences in the isotopic baseline arising from environmental parameters such as the salinity and oxygenation levels of the respective area. The larger contribution of terrestrial organic matter in HP may be due to the large input of fresh water that this lake receives from the hinterland, as indicated by the low salinity generally observed in this lake compared to the other two. These results are consistent with what has been observed in similar environments [7,74,75].
Large Carbon Ranges suggest multiple carbon sources at the base of the food web [7,35,[62][63][64]. The supply of organic matter from multiple sources in the three lakes was evident in the δ 13 C values of A. anguilla, which shifted from marine to terrestrial input with the increasing pollution. Its larger Carbon Range in the eutrophic lake indicates that this species integrated both autochthonous and allochthonous carbon pathways, while in the unpolluted lake it relied mainly on the autochthonous one.
It is acknowledged that increased N-loads promote significant changes in aquatic productivity [76] that could potentially affect the composition of the prey community [77,78]. This in turn might be reflected in the feeding behaviour of consumers at all trophic levels [2,7,79].
In our study, increased N pollution resulted in altered community composition, with decreasing diversity, which seemed to affect the feeding preference and niche width of Anguilla anguilla and Diplodus annularis. The feeding regime of the two species is known to be characterised by marked generalism and trophic plasticity [19,22,42,43,80]. However, while the eel enlarged its trophic niche, feeding off multiple resources at various trophic levels in eutrophic conditions, the seabream concentrated on a small number of invertebrates. Although dietary changes are known to depend primarily on size and growth stage [81][82][83], in our study neither A. anguilla nor D. annularis showed a relationship between body size and δ 13 C signatures. The greater trophic generalism and omnivory of eels with eutrophication may be due to the different density, accessibility and availability of the prey at higher trophic levels [19,22,24,42,43] as well as to an altered presence of potential competitors for the same food sources [78]. This was also confirmed by the results of the Bayesian mixing models, which showed a shift with eutrophication in the diet of Anguilla anguilla mostly from pelagic (fish) to benthic (invertebrates) prey, regardless of their abundance. Bouchereau et al. [24,25] reported Teleosts, amphipods and decapods as the predominant prey in the diet of A. anguilla in two North-Mediterranean lagoons, and that prey selection could be linked to the activity and accessibility of the prey itself. Rosecchi [42], Pita et al. [23] and Chaouch et al. [19] indicated molluscs (bivalves and gastropods), crustaceans, polychaetes and Teleosts as the main items in the diet of D. annularis in lagoons and coastal waters. Lammens et al. [84] and Dörner et al. [80] identified A. anguilla as belonging to the piscivorous community in many European lakes. In our study, the adoption of piscivory in the least eutrophic lake allowed the eel to reduce niche overlap [80,84,85] and therefore potential competition with other fish.
However, the trophic behaviour of the two species can be expressed differently by individuals within populations [18,20,21], enabling A. anguilla and D. annularis to include several food sources in their diet even in a single area [19,25]. In this context, the individual isotopic characterisations of carbon (δ 13 C) and nitrogen (δ 15 N), coupled with diet reconstruction at the isotopic trophic unit level (ITU, as recently proposed by Rossi et al. [31]) were crucial to the detailed estimation of the diets of the two fish species with high trophic generalism and omnivory. This allowed us to consider the possibility that each individual consumer could draw on the whole spectrum of potential food sources available within each lake [31]. In this way, we are able to describe the trophic plasticity and generalism of two ecologically and economic important fish species, and hence the real variation of the diet within the same populations under a range of eutrophication conditions.

Concluding Remarks
Understanding the trophic response of the community to eutrophication and depicting the structure of food webs in coastal lakes is still problematic [14], mainly due to the extraordinary biological diversity and complexity of the potential trophic links between species [2,86,87]. In our study, the stable isotope analysis of carbon (δ 13 C) and nitrogen (δ 15 N) provided an effective approach with which to (1) track the propagation of anthropogenic nutrient inputs along food chains, (2) evaluate the relative contributions of food sources to fish diets, and (3) quantify the trophic relationships between organisms [31,37,62,[87][88][89]. Here, the diet resulting from the application of the method recently proposed by Rossi et al. [31] substantially improved our ability to understand the response of communities to increasing eutrophication, as well as its effect on the feeding behaviour and food choices of important fish species in Mediterranean coastal lakes. Together, our results confirm food web theory as a powerful approach for obtaining valuable information for the management and conservation of these complex and productive ecosystems.

Conflicts of Interest:
The authors declare no competing interests.