IntroductionTop

Transitional environments (TEs) are unpredictable aquatic systems characterized by large environmental fluctuations of physicochemical and geochemical parameters which may influence benthic assemblages in terms of their distribution and structure (Barnes et al. 2008Barnes N., Bamber R.N., Moncrieff C.B., et al. 2008. Meiofauna in closed coastal saline lagoons in the United Kingdom: Structure and biodiversity of the nematode assemblage. Estuar. Coast. Shelf. Sci. 79: 328-340., Frontalini and Coccioni 2011Frontalini F., Semprucci F., Coccioni R., et al. 2011. On the quantitative distribution and community structure of the meio and macrofaunal communities in the coastal area of the Central Adriatic Sea (Italy). Environ. Monit. Assess. 180: 325-344.). The Adriatic coastal area hosts a large number of TEs that differ in their environmental features. These habitats range from the largest and most studied Venice Lagoon, to wetlands, estuaries, embayments, and coastal ponds, all of which differ in extent and connection with the sea. For millennia, Adriatic TEs have been modified to meet human needs and exploited for fish farming. At present, like many other coastal ecosystems, they are under severe stress due to human activity and climate change (Airoldi and Beck 2007Airoldi L., Beck M.W. 2007. Loss, status and trends for coastal marine habitats of Europe. Oceanogr. Mar. Biol. Annu. Rev. 45: 345-405.).

Free-living nematodes are the most diverse and numerically dominant metazoans in aquatic ecosystems (Balsamo et al. 2010Balsamo M., Albertelli G., Ceccherelli V.U., et al. 2010. Meiofauna of the Adriatic Sea: current state of knowledge and future perspective. Chem. Ecol. 26: 45-63., Appeltans et al. 2012Appeltans W. et al. 2012. The Magnitude of Global Marine Species Diversity. Curr. Biol. 22: 2189-2202.). This phylum is particularly important for studying TEs due to its ability to persist during environmental perturbations which may reduce or eliminate macrofauna (Hendelberg and Jensen 1993Hendelberg M., Jensen P. 1993. Vertical distribution of the nematode fauna in coastal sediment influenced by seasonal hypoxia in the bottom water. Ophelia 37: 83-94.). Nematodes are sensitive indicators of natural changes in salinity (Forster 1998Forster S.J. 1998. Osmotic stress tolerance and osmoregulation of intertidal and subtidal nematodes. J. Exp. Mar. Biol. Ecol. 224: 109-125., Barnes et al. 2008Barnes N., Bamber R.N., Moncrieff C.B., et al. 2008. Meiofauna in closed coastal saline lagoons in the United Kingdom: Structure and biodiversity of the nematode assemblage. Estuar. Coast. Shelf. Sci. 79: 328-340.), organic matter (OM) (Schratzberger and Warwick 1998Schratzberger M., Warwick R.M. 1998. Effects of the intensity and frequency of organic enrichment on two estuarine nematode communities. Mar. Ecol. Prog. Ser. 164: 83-94., Armenteros et al. 2009Armenteros M., Ruiz-Abierno A., Fernández-Garce R., et al. 2009. Biodiversity patterns of free-living marine nematodes in a tropical bay: Cienfuegos, Caribbean Sea. Estuar. Coast. Shelf. Sci. 85: 179-189.) and grain size (Vanaverbeke et al. 2002Vanaverbeke J., Gheskiere T., Steyaert M., et al. 2002. Nematode assemblages from subtidal sandbanks in the southern bight of the north sea: effect of small sedimentological differences. J. Sea. Res. 48: 197-207., Semprucci et al. 2010aSemprucci F., Colantoni P., Baldelli G., et al. 2010a. The distribution of meiofauna on back-reef sandy platforms in the Maldives (Indian Ocean). Mar. Ecol. 31: 592-607.). Furthermore, the ecological and practical advantages associated with the use of nematodes in benthic studies are good reasons for utilizing them as a bioindicator group in monitoring programmes (Balsamo et al. 2012Balsamo M., Semprucci F., Frontalini F., et al. 2012. Meiofauna as a tool for marine ecosystem biomonitoring. In: Cruzado A. (ed.), Marine Ecosystems. In Tech Publisher: pp. 77-104., Semprucci and Balsamo 2012Semprucci F., Balsamo M. 2012. Key role of free-living nematodes in the marine ecosystem. In: Boeri F., Jordan A.C. (eds). Nematodes: Morphology, Functions and Management Strategies. NOVA Science Publishers, Inc. Hauppauge, New York, pp. 109-134.). Previous meiofaunal studies in TEs have mainly been carried out at the higher taxon level, and only a few studies focusing on enclosed water bodies have been performed at species level identification (see e.g. Barnes et al. 2008Barnes N., Bamber R.N., Moncrieff C.B., et al. 2008. Meiofauna in closed coastal saline lagoons in the United Kingdom: Structure and biodiversity of the nematode assemblage. Estuar. Coast. Shelf. Sci. 79: 328-340. and references therein).

Currently, knowledge of the meiofauna in southern Italian TEs is limited to the Lesina Lagoon (Gambi et al. 2003Gambi C., Totti C., Manini E. 2003. Impact of organic loads and environmental gradients on microphytobenthos and meiofaunal distribution in a coastal lagoon. Chem. Ecol. 19: 207-223., Fabbrocini et al. 2005Fabbrocini A., Guarino A., Scirocco T., et al. 2005. Integrated biomonitoring assessment of the Lesina Lagoon (Southern Adriatic Coast, Italy): preliminary results. Chem. Ecol. 21: 479-489., Pusceddu et al. 2007Pusceddu A., Gambi C., Manini E., et al. 2007. Trophic state, eco system efficiency and biodiversity of transitional aquatic ecosystems: analysis of environmental quality based on different benthic indicators. Chem. Ecol. 23: 505-515., Frontalini et al. 2010Frontalini F., Coccioni R., Bucci C. 2010. Benthic foraminiferal assemblages and trace element contents from the lagoons of Orbetello and Lesina. Environ. Monit. Assess. 170: 245-260.). Indeed, the only information available on Lake Varano concerns macrobenthic and benthic foraminiferal assemblages (Specchiulli et al. 2010Specchiulli A., Renzi M., Scirocco T., et al. 2010. Comparative study based on sediment characteristics and macrobenthic communities in two Italian lagoons. Environ. Monit. Assess. 160: 237-256., Frontalini et al. 2013Frontalini F., Margaritelli G., Francescangeli F., et al. 2013. Benthic foraminiferal assemblages and biotopes in a coastal lake: the case study of Lake Varano (southern Italy). Acta. Protozool. 52: 147-161., 2014Frontalini F., Semprucci F., Armynot du Châtelet E., et al. 2014. Biodiversity trends of the meiofauna and foraminifera assemblages of Lake Varano (southern Italy). Proc. Biol. Soc. Wash. 127.).

This investigation documents the taxonomic and functional diversities of the nematode assemblages in the poorly known TE of Lake Varano (southern Adriatic Sea) and evaluates the ecological quality status (EQS) of this TE, thereby opening up new potential perspectives for the conservation and management of coastal areas. In particular, the following questions are addressed: (a) Are there changes in taxonomic and functional nematode diversity in response to the environmental parameters? (b) Are there significant differences in the nematode assemblages that may be imputable to anthropogenic disturbance?

Materials and MethodsTop

Study area

Located on the north side of the Gargano National Park (Foggia, Italy), Lake Varano is a coastal lagoon comprising an area of ca. 65 km² with an average water depth of ca. 4 m (Spagnoli et al. 2002Spagnoli F., Specchiulli A., Scirocco T., et al. 2002. The Lago di Varano: Hydrologic characteristics and sediment composition. Mar. Ecol. 23: 384-394.). The lake is connected to the Adriatic Sea via two artificial channels (Capoiale and Varano) that are located at the two ends of a coastal dune (Fig. 1). From a hydrological point of view, the lake is influenced both by exchanges with the sea through the two channels and by freshwater inputs from the catchment area through two tributaries (the Antonino and San Francesco Canals). The surface sediments are very heterogeneous, with the sandy fraction prevailing near the baymouth bar (Frontalini et al. 2013Frontalini F., Margaritelli G., Francescangeli F., et al. 2013. Benthic foraminiferal assemblages and biotopes in a coastal lake: the case study of Lake Varano (southern Italy). Acta. Protozool. 52: 147-161.). The lagoon has been intensively used for the farming of mussels (Mytilus galloprovincialis) and other aquaculture products (Beneduce et al. 2010Beneduce L., Vernile A., Spano G., et al. 2010. Occurrence of Vibrio vulnificus in mussel farms from the Varano lagoon environment. Lett. Appl. Microbiol. 51: 443-449.). Although the surrounding area is characterized by only a few small towns and no major industrial activity, the relatively high concentration of nitrogen species probably reflects the impact of agricultural activities (Specchiulli et al. 2008Specchiulli A., Focardi S., Renzi M., et al. 2008. Environmental heterogeneity patterns and assessment of trophic levels in two Mediterranean lagoons: Orbetello and Varano, Italy. Sci. Tot. Environ. 402: 285-298.). Though high nutrient loads and dinoflagellate blooms might suggest eutrophic conditions, Lake Varano has been regarded as an oligo-mesotrophic ecosystem (Specchiulli et al. 2008Specchiulli A., Focardi S., Renzi M., et al. 2008. Environmental heterogeneity patterns and assessment of trophic levels in two Mediterranean lagoons: Orbetello and Varano, Italy. Sci. Tot. Environ. 402: 285-298., Belmonte et al. 2011Belmonte G., Scirocco T., Denitto F. 2011. Zooplankton composition in Lake Varano (Adriatic Sea coast, Italy). Ital. J. Zool. 78: 370-378., Frontalini et al. 2013Frontalini F., Margaritelli G., Francescangeli F., et al. 2013. Benthic foraminiferal assemblages and biotopes in a coastal lake: the case study of Lake Varano (southern Italy). Acta. Protozool. 52: 147-161.).

Sample collection

Samples from 21 stations were collected using a modified model of the Van Veen grab in late March 2012. This modified model permits the insertion of a plexiglas corer from the central part of the grab, enabling portions of almost completely undisturbed sediments to be captured. A set of sub-samples of the sediments was taken from each station for the purposes of sedimentological, geochemical and nematode analyses (two and three replicates for geochemical and nematode analyses, respectively). All of the analyses were carried out on the upper surface sediment layer (0-2 cm). During the sampling process, physico-chemical parameters of the bottom water, such as temperature, pH, salinity, oxidation reduction potential (ORP) and dissolved oxygen (DO, expressed as mg/l) were measured using a conductivity, temperature and depth (CTD) probe. In particular, these environmental parameters selected as independent variables are well known to affect the structure, abundance and diversity of the nematode assemblage (Heip et al. 1985Heip C., Vincx M., Vranken G. 1985. The ecology of marine nematodes. Oceanogr. Mar. Biol. Annu. Rev. 23: 399-489., Semprucci and Balsamo 2012Semprucci F., Balsamo M. 2012. Key role of free-living nematodes in the marine ecosystem. In: Boeri F., Jordan A.C. (eds). Nematodes: Morphology, Functions and Management Strategies. NOVA Science Publishers, Inc. Hauppauge, New York, pp. 109-134.).

Environmental parameter analyses

Grain size analyses were conducted at Lille 1 University (France) following the principle of diffraction described in Trentesaux et al. (2001)Trentesaux A., Recourt P., Bout-Roumazeilles V., et al. 2001. Carbonate grain-size distribution in hemipelagic sediments from a laser particle sizer. J. Sediment. Res. 71: 858-862.. Three main fractions were considered: clay (<2 µm), silt (2 to 63 µm) and sand (63 to 2000 µm).

The quantity and quality of the OM were estimated by analysing the total carbon (TC), the total organic carbon (TOC), the total sulfur (TS), the total hydrogen (TH) and the total nitrogen (TN) content in the sediment (for details, see Armynot du Châtelet et al. 2013Armynot du Châtelet E., Bout-Roumazeilles V., Coccioni R., et al. 2013. Environmental control on shell structure and composition of agglutinated foraminifera along a proximal-distal transect in the Marmara Sea. Mar. Geol. 335: 114-128.).

Concentrations of 37 trace elements were obtained at Activation Laboratories Ltd. (Ontario, Canada, http://www.actlabs.com), where a fraction of ~0.5g was digested with aqua regia and then analysed using Inductively Coupled Plasma Optical Emission Spectrometry. Only concentrations of As, Cd, Co, Cr, Cu, Ni, Pb, Zn and Mn were considered. The enrichment of some chemical elements whose higher concentrations may be toxic was calculated using the Pollution Load Index (PLI) following Martins et al. (2013)Martins V.A, Frontalini F., Figueira R.C.L., et al. 2013. Assessment of the health quality of Ria de Aveiro (Portugal): heavy metals and benthic foraminifera. Mar. Pollut. Bull. 70: 18-33..

Nematode assemblage analyses

Samples utilized for the nematode analysis were treated with a 7% MgCl₂ aqueous solution for narcotizing fauna, fixed in a 4% formaldehyde solution in buffered sea-water and stained with Rose Bengal (0.2 g L^–1). In laboratory, the samples were rinsed with a gentle jet of fresh water through a 0.5-mm sieve to separate the macrofauna (Danovaro et al. 2004Danovaro R., Gambi C., Mirto S., et al. 2004. Meiofauna. Biol. Mar. Mediterr. 11: 55-97.). They were then decanted, sieved 10 times through a 42-μm mesh and centrifuged three times with Ludox HS30 (specific density 1.18) (Pfannkuche and Thiel 1988Pfannkuche O., Thiel H. 1988. Sample processing. In: Higgins R.P., Thiel H. (eds), Introduction to the study of meiofauna. Smithsonian Inst., Washington, pp. 134-145.). According to Moreno et al. (2011)Moreno M., Semprucci F., Vezzulli L., et al. 2011. The use of nematodes in assessing ecological quality status in the Mediterranean coastal ecosystems. Ecol. Ind. 11: 328-336., 100 nematodes were randomly withdrawn from each of the three replicates and subsequently mounted on permanent slides according to Seinhorst (1959)Seinhorst J.W. 1959. A rapid method for the transfer of nematodes from fixative to anhydrous glycerine. Nematologica 4: 67-69.. The specimens were identified under a 100× oil immersion objective using Nomarski differential interference contrast illumination (Optiphot-2 Nikon) and the NeMys online identification key (Deprez et al. 2005Deprez T. et al. 2005. NeMys. World Wide Web electronic publication.).

The trophic and life strategies were among the functional traits considered. Following Wieser (1953)Wieser W. 1953. Die Beziehung zwischen Mundhöhlengestalt. Ernährungsweise und Vorkommen bei freilebenden marinen Nematoden. Eine okologisch-morphologische studie. Ark. Zool. 4: 439-484., the nematodes were divided into four trophic groups: selective (1A) and non-selective (1B) deposit feeders, epistrate feeders (2A), and predators/omnivores (2B).

Bongers (1990)Bongers T. 1990. The maturity index: an ecological measure of environmental disturbance based on nematode species composition. Oecologia 83: 14-19. and Bongers et al. (1991)Bongers T., Alkemade R., Yeates G.W. 1991. Interpretation of disturbance-induced maturity decrease in marine nematode assemblages by means of the Maturity Index. Mar. Ecol. Prog. Ser. 76: 135-142. proposed classifying nematodes as ‘colonizers’ (r-strategists) or ‘persisters’ (k-strategists). Therefore, the maturity index (MI) was calculated as the weighted average of the individual colonizer-persistent (c-p) values.

The biodiversity pattern of the assemblage was described by means of the Shannon (H’, using log-base 2), Pielou, Margalef and Simpson indices.

Statistical analysis

All of the statistical analyses were performed on the relative abundances of the nematode species. Following Armynot du Châtelet et al. (2004)Armynot du Châtelet E., Debenay J.P., Soulard R. 2004. Foraminiferal proxies for pollution monitoring in moderately polluted harbours. Environ. Pollut. 127: 27-40., only taxa with a relative abundance exceeding 5% in at least one sample were taken into consideration. However, though rare species often considered the most sensitive ones were down-weighted in this way, their contribution in the TE system was measured by means of MI calculation. Prior to the statistical analyses, all of the biotic and abiotic data were normalized by applying an additive logarithmic transformation log(x+1). In order to define the different assemblages, which included a group of species with a similar spatial distribution pattern, an R-mode cluster analysis (CA) was performed. Alternatively, a Q-mode CA was carried out for the ordination of the samples based on the relative abundances of species. These tests were performed by applying Ward’s linkage method and the findings were given in terms of the Euclidean distance. A principal component analysis (PCA) was used to determine the assemblage relationships to the abiotic parameters and for the ordination of sample locations.

ResultsTop

Environmental parameters

The data on the environmental parameters are summarized in Supplementary material Table S1. The pH was slightly alkaline, ranging from 7.6 (V1, V15) to 8.0 (V29, V39). The salinity varied from 23.0‰ (V50) to 31.9‰ (V13), with a clear increase in the outer part of the lake and, in particular next to the Capoiale Channel. The DO values varied between 5.7 (V15) and 12.8 (V35), with the highest values found in the central part of the lake and the lowest near to the two channels. The ORP values ranged from 164.0 (V39 and V41) to 286.0 (V1), with the highest values in the outer part of the lagoon. The sediments were primarily composed of silt (47% on average) and sand (51.9% on average), with a low amount of clay (1.1% on average). The outer part of the lake, particularly in front of the two channels and the baymouth bar, was dominated by sand, whereas silt was predominant in the southeastern part of the lake. The CaCO₃ varied between 9.5% (V4) and 80.8% (V37), and was particularly abundant on the western side of the lake. The TN ranged between 0.1 (V4, V7, V11, V41, V48, V49 and V50) and 0.4 (V27 and V47), whereas TS was between 0.9% (V7) and 1.8% (V1). The TC values varied from 1.3 (V4) to 10.9 (V37 and V42), while TH ranged from 0.1 (V4) to 1.3 (V39 and V47). The TOC varied from 0.2% (V4) to 4.2% (V27), with higher values being found in the central and southwestern parts of the lake. A clear spatial distribution of most of the trace elements was visible, with the highest concentrations being found in both the central and southeastern parts of the lake. The PLI values ranged from 1.0 (V4) to 11.4 (V41), and revealed a general enrichment of the concentrations of trace elements in both the central and southeastern parts of the lake.

Nematode assemblage

Fifty-five nematode species, belonging to 37 genera and 17 families, were recorded at Lake Varano (Supplementary material Table S2). The richest families were Chromadoridae (10 species), Linhomoeidae (9) and Xyalidae (8). Almost half of the selected specimens belonged to either the Desmodoridae (35%) or Linhomoeidae (16%) families, followed by Leptolaimidae (12%), Microlaimidae (9%) and Xyalidae (9%). Three species had very high abundances and a wide distribution: Molgolaimus allgeni (35%), Terschellingia longicaudata (13%) and Leptolaimus luridus (11%). They were followed by Daptonema normandicum (6%), Sabatieria pulchra (5%), Aponema torosa (4%) and Calomicrolaimus honestus (4%). Aponema torosa was more abundant in the outer (particularly in front of the two channels) and central parts of the lake, whereas D. normandicum and S. pulchra were more abundant in the stations close to the San Francesco output and in the central area.

The richness of the species ranged from 4 (V29) to 27 (V7). The H’ ranged from 0.9 (V29) to 4.2 (V7), the Margalef index from 0.7 (V29) to 5.6 (V7), the Pielou evenness from 0.4 (V29) to 0.9 (V7), and the Simpson evenness from 0.33 (V29) to 0.94 (V7 and V50) (Supplementary material Table S3). These indices revealed that the nematode assemblage in the central part of the lake was less diversified.

The trophic structure of the assemblage was mainly characterized by 1A and 1B (60.1% and 21.8%, respectively, on average), followed by 2A (16.2%) and 2B (2%) (Supplementary material Table S3). The group 1A was widely distributed, whereas the group 1B was particularly associated with freshwater inputs and the Varano Channel. The group 2B and, in particular, 2A were more abundant in the outer part of the lake.

The MI ranged from 2.1 (V4) to 2.9 (V27), with higher values mainly in the central part of the lake. All of the classes of the colonizers-persisters (c-p) were revealed, with the only exception being class c-p5 (Supplementary material Table S3). On average, the c-p3 species were the dominant component of the nematode assemblage (58.5%) followed by c-p2 (38.4%), c-p1 (2.8 %) and c-p4 (0.4%).

Statistical analysis

The results of the R-mode CA enabled two main (A and B) clusters and two different subclusters (A1 and A2) to be recognized (Fig. 2). The most abundant taxa in subcluster A1 were A. torosa, D. normandicum and S. pulchra, whereas M. allgeni, L. luridus and T. longicaudata were the only three representatives of the subcluster A2. Cluster B grouped together several taxa with very low relative abundances, including Daptonema fistulatus, O. wieseri and P. longicaudata as the main representatives. The Q-mode CA separated two different clusters (A and B) and two subclusters (B1 and B2) (Fig. 3). Cluster A was represented by station 4, which was in front of the baymouth bar, where the main taxa characterizing the nematode assemblages were O. wieseri, P. longicaudatus, P. pellucida and A. caudostriatus. This cluster was characterized by the lowest richness and MI, and was dominated by c-p2 taxa and 1B feeders. This subcluster was dominated by sand and characterized by the lowest values of CaCO₃, TOC, TN, TC, TH and TS. The subcluster B1 included the stations of the southeastern part of the lake: it was mainly dominated by M. allgeni, L. luridus, T. longicaudata, D. normandicum and S. pulchra. The same subcluster was characterized by higher levels of richness and the highest H’ value. It was dominated by c-p2 and c-p3 taxa and by 1A and 1B feeders. The subcluster B1 also showed the lowest salinity and ORP values and was mainly characterized by silty and, secondarily, sandy substrates. It also had higher values of TC and TOC and the highest PLI value.

Subcluster B2 included all the other stations, and can be further subdivided into two subclusters, B2a and B2b. The subcluster B2a grouped a few stations located in front of the Varano Channel, in the southern margin and in the central part of the lake. It was dominated by M. allgeni, C. honestus and D. normandicum, followed by S. pulchra. This subcluster showed the lowest values of Pielou and richness, if cluster A was excluded. It was also dominated by 1A and secondarily by 1B and 2A feeders and by c-p2 and c-p3 taxa. It was also characterized by sandy substrates with a high abundance of silt. It had higher values of OM and the highest value of TS. B2b included plenty of stations located in both the outer margins and the central parts of the lake, which are mainly represented by M. allgeni, along with T. longicaudata, L. luridus and A. torosa. This subcluster showed the highest values of sand and MI, and was also dominated by 1A and, to a lesser extent, by 1B and 2A feeders. It was also dominated by c-p3 taxa. This subcluster was characterized by the highest salinity values and sand-silt dominated sediments with the highest percentages of CaCO₃. It also had the highest values of TOC, TN, TC, and TH.

The eigenvalues of the PCA revealed that the first two factors explained ca. 49.3% of the total variance, whereas the factor variable plans showed how the species (those used for the CA and representing the main variables) and secondary variables (abiotic) were related to the different factors (Fig. 4). Factor 1, which explained 33.3% of the data variability, appeared to be related to the sediment characteristics in terms of grain size. Accordingly, silt and sand were placed on the opposite sides of the factor plan.

Strong relationships among the fine sediment fractions of both silt and clay were evidenced with T. longicaudata, L. luridus and the PLI. On the other hand, an opposite behaviour was highlighted among the finer fractions and the relative abundances of some taxa, including A. caudostriatus, P. pellucida, P. longicaudatus and O. wieseri. Factor 2, which justified 16.0% of the data variability, was mainly related to some of the physicochemical parameters of the bottom water, namely salinity, and the OM quality and quantity within the sediments (TOC, TC, TN, TS and TH). In particular, this factor seemed to be positively related to Microlaimus sp. 2, C. honestus, M. allgeni, and the MI, and negatively with D. normandicum, Thalassomonhystera parva, S. pulchra, N. papillosa and the 1B feeders.

DiscussionTop

Taxonomic biodiversity: comparison with other TEs

The description of the distribution patterns of species is one of the fundamental starting blocks in the ecology of biological assemblages. However, there is very little information available on nematode taxonomic composition in Italian TEs and in the entire southern part of the Adriatic Sea (Semprucci 2013Semprucci F., Frontalini F., Covazzi-Harriague A., et al. 2013. Meio- and Macrofauna in the marine area of the Monte St. Bartolo Natural Park (Central Adriatic Sea, Italy). Sci. Mar. 77: 189-199.).

The nematode assemblage of the study area appeared to be mainly comprised of representative species of the fine sediments (Heip et al. 1985Heip C., Vincx M., Vranken G. 1985. The ecology of marine nematodes. Oceanogr. Mar. Biol. Annu. Rev. 23: 399-489., Travizi and Vidakovic 1997Travizi A., Vidakovic J. 1997. Nematofauna in the Adriatic Sea: Review and check-list of free-living nematode species. Helgol. Meeresunters. 51: 503-519.) or of TEs (Villano and Warwick 1995Villano N., Warwick R.M. 1995. Meiobenthic communities associated with the seasonal cycle of growth and decay of Ulva rigida Arardh in the Palude Della Rosa, Lagoon of Venice. Estuar. Coast. Shelf. Sci. 4: 181-194., Pallo et al. 1998Pallo P., Widbom B., Olafsson E. 1998. A quantitative survey of the benthic meiofauna in the Gulf of Riga (eastern Baltic Sea), with special reference to the structure of nematode assemblages. Ophelia 49: 117-139., Fabbrocini et al. 2005Fabbrocini A., Guarino A., Scirocco T., et al. 2005. Integrated biomonitoring assessment of the Lesina Lagoon (Southern Adriatic Coast, Italy): preliminary results. Chem. Ecol. 21: 479-489., Barnes et al. 2008Barnes N., Bamber R.N., Moncrieff C.B., et al. 2008. Meiofauna in closed coastal saline lagoons in the United Kingdom: Structure and biodiversity of the nematode assemblage. Estuar. Coast. Shelf. Sci. 79: 328-340.). The overall richness revealed values that were highly comparable with those reported for other Italian TEs (Guerrini et al. 1998Guerrini A., Colangelo M.A., Ceccherelli V.U. 1998. Recolonization patterns of meiobenthic communities in brackish vegetated and unvegetated habitats after induced hypoxia/anoxia. Hydrobiologia 375/376: 73-87., Fabbrocini et al. 2005Fabbrocini A., Guarino A., Scirocco T., et al. 2005. Integrated biomonitoring assessment of the Lesina Lagoon (Southern Adriatic Coast, Italy): preliminary results. Chem. Ecol. 21: 479-489.), but were generally lower than those reported for European brackish-water systems (e.g. Hendelberg and Jensen 1993Hendelberg M., Jensen P. 1993. Vertical distribution of the nematode fauna in coastal sediment influenced by seasonal hypoxia in the bottom water. Ophelia 37: 83-94., Barnes et al. 2008Barnes N., Bamber R.N., Moncrieff C.B., et al. 2008. Meiofauna in closed coastal saline lagoons in the United Kingdom: Structure and biodiversity of the nematode assemblage. Estuar. Coast. Shelf. Sci. 79: 328-340., Ferrero et al. 2008Ferrero T.J., Debenham N.J., Lambshead P.J.D. 2008. The nematodes of the Thames estuary: Assemblage structure and biodiversity, with a test of Attrill’s linear model. Estuar. Coast. Shelf. Sci. 79: 409-418.). This is probably due to the low salinity gradient of Lake Varano (Frontalini et al. 2013Frontalini F., Margaritelli G., Francescangeli F., et al. 2013. Benthic foraminiferal assemblages and biotopes in a coastal lake: the case study of Lake Varano (southern Italy). Acta. Protozool. 52: 147-161.), which does not determine a transition from a marine to typical freshwater fauna, leading to a lower level of richness than in the European TEs.

Effect of environmental variables on taxonomic structure and functional traits

Studies carried out in both marine and coastal lagoon ecosystems have revealed that the quantity and quality of the OM and salinity are among the most important factors influencing meiofaunal and nematode assemblages (Gambi et al. 2003Gambi C., Totti C., Manini E. 2003. Impact of organic loads and environmental gradients on microphytobenthos and meiofaunal distribution in a coastal lagoon. Chem. Ecol. 19: 207-223., Barnes et al. 2008Barnes N., Bamber R.N., Moncrieff C.B., et al. 2008. Meiofauna in closed coastal saline lagoons in the United Kingdom: Structure and biodiversity of the nematode assemblage. Estuar. Coast. Shelf. Sci. 79: 328-340.). However, the distribution of the assemblages is often more complex, reflecting multiple interactions between the distance to the sea, water depth, grain size, hydrodynamic turnover time, sediment oxygen and nutrient availability (Frontalini et al. 2013Frontalini F., Margaritelli G., Francescangeli F., et al. 2013. Benthic foraminiferal assemblages and biotopes in a coastal lake: the case study of Lake Varano (southern Italy). Acta. Protozool. 52: 147-161.).

As shown in the PCA, the sediment grain size appeared to be one of the most important environmental variables affecting the nematode species distribution of Lake Varano. This parameter may, in fact, exert significant influence on the structure of the nematode assemblage, which appears to be very sensitive, even to slight sedimentological variations (Steyaert et al. 2003Steyaert M., Vanaverbeke J., Vanreusel A., et al. 2003. The importance of fine-scale, vertical profiles in characterizing nematode community structure. Estuar. Coast. Shelf. Sci. 58: 353-366., Vanaverbeke et al. 2002Vanaverbeke J., Gheskiere T., Steyaert M., et al. 2002. Nematode assemblages from subtidal sandbanks in the southern bight of the north sea: effect of small sedimentological differences. J. Sea. Res. 48: 197-207., Semprucci et al. 2010aSemprucci F., Colantoni P., Baldelli G., et al. 2010a. The distribution of meiofauna on back-reef sandy platforms in the Maldives (Indian Ocean). Mar. Ecol. 31: 592-607.). In particular, the increasing silt content appeared to have the greatest effect on the lake’s nematodes. The presence of silt, even in small amounts, can actually reduce the permeability of sediment. Furthermore, a greater retention of OM is often observed in silty sediments (Albertelli et al. 1999Albertelli G., Covazzi-Harriague A., Danovaro R., et al. 1999. Differential responses of bacteria, meiofauna and macrofauna in a shelf area (Ligurian Sea, NW Mediterranean): role of food availability. J. Sea. Res. 42: 11-26., Frontalini et al. 2011Frontalini F., Semprucci F., Coccioni R., et al. 2011. On the quantitative distribution and community structure of the meio and macrofaunal communities in the coastal area of the Central Adriatic Sea (Italy). Environ. Monit. Assess. 180: 325-344., Semprucci et al. 2013Frontalini F., Margaritelli G., Francescangeli F., et al. 2013. Benthic foraminiferal assemblages and biotopes in a coastal lake: the case study of Lake Varano (southern Italy). Acta. Protozool. 52: 147-161.), leading to OM accumulation and oxygen deficiency (Graf 1992Graf G. 1992. Benthic-pelagic coupling: A benthic view. Oceanogr. Mar. Biol. Annu. Rev. 30: 149-190.). Accordingly, total OM and, in particular, its refractory component, increased in the study area in parallel with the silt fraction, while the highest ORP values detected were only associated with sand (Frontalini et al. 2013Frontalini F., Margaritelli G., Francescangeli F., et al. 2013. Benthic foraminiferal assemblages and biotopes in a coastal lake: the case study of Lake Varano (southern Italy). Acta. Protozool. 52: 147-161.). As suggested by Vanaverbeke et al. (2011)Vanaverbeke J., Merckx B., Degraer S., et al. 2011. Sediment-related distribution patterns of nematodes and macrofauna: Two sides of the benthic coin? Mar. Environ. Res. 71: 31-40., this may prove that grain size only indirectly shapes nematode assemblages, and the biogeochemical environment resulting from the interplay between hydrodynamic features and grain size is a key factor structuring the assemblages. As shown by the PCA, the positive correlation between silt and the PLI is noteworthy. Spagnoli et al. (2002)Spagnoli F., Specchiulli A., Scirocco T., et al. 2002. The Lago di Varano: Hydrologic characteristics and sediment composition. Mar. Ecol. 23: 384-394. and Fabbrocini et al. (2010)Fabbrocini A., Di Stasio M., D’Adamo R. 2010. Computerized sperm motility analysis in toxicity bioassays: A new approach to pore water quality assessment. Ecotoxicol. Environ. Saf. 73: 1588-1595. have already evidenced a correlation between trace elements and the fine sediment in Lake Varano. This greatly depends on the adsorptive properties of clay and silt particles, and supports the notion that grain size is also important in influencing pollutant accumulation in sediment (Bernardello et al. 2006Bernardello M., Secco T., Pellizzato F., et al. 2006. The changing state of contamination in the lagoon of Venice. Part 2: heavy metals. Chemosphere 64: 1334-1345.).

As suggested by the CA and the PCA, the species closely associated with the silt were mainly M. allgeni, T. longicaudata and L. luridus. The high abundance of Desmodoridae in the study area was mainly due to M. allgeni. The Molgolaimus genus is present in all of the marine environments, ranging from shallow waters to the deep sea. It is generally associated with muddy sediment (e.g. Lambshead et al. 2000Lambshead P.J.D., Tietjen J., Ferrero T.J., et al. 2000. Latitudinal diversity gradients in the deep sea with special reference to North Atlantic nematodes. Mar. Ecol. Prog. Ser. 194: 159-167., Muthumbi et al. 2004Muthumbi A.W., Vanreusel A., Duineveld G., et al. 2004. Nematode community structure along the continental slope off the Kenyan Coast, Western Indian Ocean. Int. Rev. Hydrobiol. 89: 188-205.) and its opportunistic life style has been documented (Schratzberger et al. 2003Schratzberger M., Bolam S.G., Whomersley P., et al. 2003. Development of a meiobenthic nematode community following the intertidal placement of various types of sediment. J. Exp. Mar. Biol. Ecol. 303: 79-96.). In the Adriatic Sea, the presence of M. allgeni has been recorded in the Venice Lagoon, as well as in the offshore area of the Rovinj-Po River mouth and the Pored-Venezia transect (Travizi and Vidakovic 1997Travizi A., Vidakovic J. 1997. Nematofauna in the Adriatic Sea: Review and check-list of free-living nematode species. Helgol. Meeresunters. 51: 503-519.).

Terschellingia longicaudata is a typical inhabitant of estuarine, intertidal and subtidal areas, with its highest abundances being in muddy sediment (Warwick and Gee 1984Warwick R.M., Gee J.M. 1984. Community structure of estuarine meiobenthos. Mar. Ecol. Prog. Ser. 18: 97-111., Barnes et al. 2008Barnes N., Bamber R.N., Moncrieff C.B., et al. 2008. Meiofauna in closed coastal saline lagoons in the United Kingdom: Structure and biodiversity of the nematode assemblage. Estuar. Coast. Shelf. Sci. 79: 328-340., Semprucci et al. 2010bSemprucci F., Boi P., Manti A., et al. 2010b. Benthic communities along a littoral of the Central Adriatic Sea (Italy). Helgol. Mar. Res. 64: 101-115.). It is also well known as tolerant to a variety of natural and anthropogenic stressors (e.g. Armenteros et al. 2009Armenteros M., Ruiz-Abierno A., Fernández-Garce R., et al. 2009. Biodiversity patterns of free-living marine nematodes in a tropical bay: Cienfuegos, Caribbean Sea. Estuar. Coast. Shelf. Sci. 85: 179-189., Gollner et al. 2010Gollner S., Riemer B., Martinez Arbizu P., et al. 2010. Diversity of Meiofauna from the East Pacific Rise across a Gradient of Hydrothermal Fluid Emissions. PLoS ONE 5: e12321., Guilini et al. 2012Guilini K., Levin L.A., Vanreusel A. 2012. Cold seep and oxygen minimum zone associated sources of margin heterogeneity affect benthic assemblages, diversity and nutrition at the Cascadian margin (NE Pacific Ocean). Prog. Oceanogr. 96: 77-92.). The physiological and behavioural adaptations of T. longicaudata to poorly oxygenated environments include a low respiratory rate and slow movement (Warwick and Gee 1984Warwick R.M., Gee J.M. 1984. Community structure of estuarine meiobenthos. Mar. Ecol. Prog. Ser. 18: 97-111., Guerrini et al. 1998Guerrini A., Colangelo M.A., Ceccherelli V.U. 1998. Recolonization patterns of meiobenthic communities in brackish vegetated and unvegetated habitats after induced hypoxia/anoxia. Hydrobiologia 375/376: 73-87.). The deposition of insoluble metal sulphides in intracellular inclusions in T. longicaudata has been suggested as a mechanism of detoxification of sulphide (Nicholas et al. 1987Nicholas W.L., Goodchild D.J., Steward A. 1987. The mineral composition of intracellular inclusions in nematodes from thiobiotic mangrove mud-flats. Nematologica 33: 167-179.). Leptolaimus is another genus that is frequently and widely distributed in TE systems (Heip et al. 1985Heip C., Vincx M., Vranken G. 1985. The ecology of marine nematodes. Oceanogr. Mar. Biol. Annu. Rev. 23: 399-489., Pallo et al. 1998Pallo P., Widbom B., Olafsson E. 1998. A quantitative survey of the benthic meiofauna in the Gulf of Riga (eastern Baltic Sea), with special reference to the structure of nematode assemblages. Ophelia 49: 117-139., Ferrero et al. 2008Ferrero T.J., Debenham N.J., Lambshead P.J.D. 2008. The nematodes of the Thames estuary: Assemblage structure and biodiversity, with a test of Attrill’s linear model. Estuar. Coast. Shelf. Sci. 79: 409-418.). In particular, L. luridus has been found in muds characterized by the presence of organic detritus (Timm 1963Timm R.W. 1963. Marine nematodes of the suborder Monhysterina from the Arabian Sea at Karachi. Proc. Helminthol. Soc. Wash. 30: 34-49.).

Sabatieria pulchra inhabits all types of sediment, but its highest densities are in mud (Hendelberg and Jensen 1993Hendelberg M., Jensen P. 1993. Vertical distribution of the nematode fauna in coastal sediment influenced by seasonal hypoxia in the bottom water. Ophelia 37: 83-94., Steyaert et al. 2007Steyaert M., Moodley L., Nadong T., et al. 2007. Responses of intertidal nematodes to short-term anoxic events. J. Exp. Mar. Biol. Ecol. 345: 175-184.). It is generally present in very shallow down to deeper waters. In anoxic and strongly reduced sediment, it is often the only remaining species. It dwells deep in the sediment, and is often found in high densities near the redox potential discontinuity layer (Wetzel et al. 2002Wetzel M.A., Weber A., Giere O. 2002. Re-colonization of anoxic/ sulfidic sediments by marine nematodes after experimental removal of macroalgal cover. Mar. Biol. 141: 679-689., Steyaert et al. 2003Steyaert M., Vanaverbeke J., Vanreusel A., et al. 2003. The importance of fine-scale, vertical profiles in characterizing nematode community structure. Estuar. Coast. Shelf. Sci. 58: 353-366., 2005Steyaert M., Moodley L., Vanaverbeke J., et al. 2005. Laboratory experiments on the infaunal activity of intertidal nematodes. Hydrobiologia 540: 217-223.). Data from a heavily polluted, tropical semi-closed bay (Havana Bay) revealed an enhancement of S. pulchra over T. longicaudata, probably due to the increase in grain size (Armenteros et al. 2009Armenteros M., Ruiz-Abierno A., Fernández-Garce R., et al. 2009. Biodiversity patterns of free-living marine nematodes in a tropical bay: Cienfuegos, Caribbean Sea. Estuar. Coast. Shelf. Sci. 85: 179-189.). This could explain the lack of a close correlation between S. pulchra and silt. However, the same authors also referred to possible interspecific competition for food or space between S. pulchra and T. longicaudata.

As highlighted by the statistical analyses, the enhancement of ORP and the sand percentage values led to a higher abundance of A. caudostriatus, O. wieseri, P. pellucida and P. longicaudatus, especially in the area in front of the baymouth bar. With a progressive increase in grain size (and, above all, a decrease in the silt-clay content), Axonolaimidae and Cyatholaimidae, and some representative species of Xyalidae, become increasingly abundant (Heip et al. 1985Heip C., Vincx M., Vranken G. 1985. The ecology of marine nematodes. Oceanogr. Mar. Biol. Annu. Rev. 23: 399-489.). Indeed, the presence of these taxa has previously been documented, particularly in intertidal or sublittoral sandy substrata (Warwick 1971Warwick R.M. 1971. Nematode associations in the Exe Estuary. J. Mar. Biol. Assoc. U.K. 51: 430-454., Semprucci 2013Semprucci F., Frontalini F., Covazzi-Harriague A., et al. 2013. Meio- and Macrofauna in the marine area of the Monte St. Bartolo Natural Park (Central Adriatic Sea, Italy). Sci. Mar. 77: 189-199.). Diversity indices, unlike the trends observed in the literature (Vanaverbeke et al. 2002Vanaverbeke J., Gheskiere T., Steyaert M., et al. 2002. Nematode assemblages from subtidal sandbanks in the southern bight of the north sea: effect of small sedimentological differences. J. Sea. Res. 48: 197-207., Semprucci et al. 2010aSemprucci F., Colantoni P., Baldelli G., et al. 2010a. The distribution of meiofauna on back-reef sandy platforms in the Maldives (Indian Ocean). Mar. Ecol. 31: 592-607.,bSemprucci F., Boi P., Manti A., et al. 2010b. Benthic communities along a littoral of the Central Adriatic Sea (Italy). Helgol. Mar. Res. 64: 101-115.), were not significantly affected by grain size, maintaining levels that were comparable to those recorded for other TEs. On the other hand, the relative proportions of the feeding guilds appeared to be more related to the sediment features. This depends on the nature of the available food, which in turn depends on the nature of the sediment. In particular, the 1A and 2A feeders were significantly affected by silt and sand increases, respectively. The increase of 1A, mainly bacteriovorus, in the silty sediment may be related to the higher abundances of microorganisms in the finest sediments (Albertelli et al. 1999Albertelli G., Covazzi-Harriague A., Danovaro R., et al. 1999. Differential responses of bacteria, meiofauna and macrofauna in a shelf area (Ligurian Sea, NW Mediterranean): role of food availability. J. Sea. Res. 42: 11-26., Semprucci et al. 2010bSemprucci F., Boi P., Manti A., et al. 2010b. Benthic communities along a littoral of the Central Adriatic Sea (Italy). Helgol. Mar. Res. 64: 101-115.). In sand, a high benthic primary production, a great abundance of diatoms, and wide surfaces suitable for scraping off the algal and bacterial biofilms can be found (Semprucci et al. 2010aSemprucci F., Colantoni P., Baldelli G., et al. 2010a. The distribution of meiofauna on back-reef sandy platforms in the Maldives (Indian Ocean). Mar. Ecol. 31: 592-607.,bSemprucci F., Boi P., Manti A., et al. 2010b. Benthic communities along a littoral of the Central Adriatic Sea (Italy). Helgol. Mar. Res. 64: 101-115.), explaining the association between the 2A feeders and well-sorted fine sand. The only c-p class correlated with sediment type was c-p3 (intermediate colonizers), whose numbers increased with the silt amount, along with the PLI and OM enrichment (Bongers et al. 1991Bongers T., Alkemade R., Yeates G.W. 1991. Interpretation of disturbance-induced maturity decrease in marine nematode assemblages by means of the Maturity Index. Mar. Ecol. Prog. Ser. 76: 135-142.).

In the second factor of the PCA, TN, TOC, TH, TS and salinity were the most important environmental variables. Salinity was not strongly correlated with the nematode assemblages. As stated above, this could be due to the small variations in salinity detected in the area, and also confirms the macrofaunal spatial distribution observed by Specchiulli et al. (2010)Specchiulli A., Renzi M., Scirocco T., et al. 2010. Comparative study based on sediment characteristics and macrobenthic communities in two Italian lagoons. Environ. Monit. Assess. 160: 237-256.. However, Barnes et al. (2008)Barnes N., Bamber R.N., Moncrieff C.B., et al. 2008. Meiofauna in closed coastal saline lagoons in the United Kingdom: Structure and biodiversity of the nematode assemblage. Estuar. Coast. Shelf. Sci. 79: 328-340. suggested that to document the controlling action of salinity on nematode assemblages, longer-term salinity data are fundamental. As highlighted by the CA and the PCA, C. honestus, which was reported by Ferrero et al. (2008)Ferrero T.J., Debenham N.J., Lambshead P.J.D. 2008. The nematodes of the Thames estuary: Assemblage structure and biodiversity, with a test of Attrill’s linear model. Estuar. Coast. Shelf. Sci. 79: 409-418. as being abundant in mixoeuhaline zones, was one of the species more positively affected by salinity. In contrast, T. parva, D. normandicum and S. pulchra appeared to be more negatively affected by this parameter. This is in line with the results of Barnes et al. (2008)Barnes N., Bamber R.N., Moncrieff C.B., et al. 2008. Meiofauna in closed coastal saline lagoons in the United Kingdom: Structure and biodiversity of the nematode assemblage. Estuar. Coast. Shelf. Sci. 79: 328-340., who reported the first two species in an oligohaline salinity regime and the third in a polyhaline regime. Daptonema normandicum was also reported to be widely distributed in the Thames estuary by Ferrero et al. (2008)Ferrero T.J., Debenham N.J., Lambshead P.J.D. 2008. The nematodes of the Thames estuary: Assemblage structure and biodiversity, with a test of Attrill’s linear model. Estuar. Coast. Shelf. Sci. 79: 409-418. and in the ‘Valli di Comacchio’ complex (Northern Adriatic Sea) by Guerrini et al. (1998)Guerrini A., Colangelo M.A., Ceccherelli V.U. 1998. Recolonization patterns of meiobenthic communities in brackish vegetated and unvegetated habitats after induced hypoxia/anoxia. Hydrobiologia 375/376: 73-87..

Pielou and Simpson indices seemed to be more affected by the second factor plan of the PCA, increasing with the reduction in the OM supply and salinity. This could be related to the discharge from the San Antonino and Francesco Canals, as documented by Specchiulli et al. (2008)Specchiulli A., Focardi S., Renzi M., et al. 2008. Environmental heterogeneity patterns and assessment of trophic levels in two Mediterranean lagoons: Orbetello and Varano, Italy. Sci. Tot. Environ. 402: 285-298.. As also suggested by the CA, the 1B feeders were widely distributed in the lake. However, the PCA showed that they were, as an opportunist trophic group, more associated with the freshwater inputs and secondly to sand (Semprucci and Balsamo 2012Semprucci F., Balsamo M. 2012. Key role of free-living nematodes in the marine ecosystem. In: Boeri F., Jordan A.C. (eds). Nematodes: Morphology, Functions and Management Strategies. NOVA Science Publishers, Inc. Hauppauge, New York, pp. 109-134.). Moreover, the MI, along with the other c-p classes, was more affected by the second factor plan, exhibiting a notable correlation with the quality and quantity of OM. Indeed, it appeared to be positively correlated with OM content and salinity enhancements. The lower values of the MI and the increases in abundance of the extreme and general colonizers (c-p1 and 2) in relation to the lower salinity values could corroborate the negative effects of freshwater inputs on the Varano system.

Notes on the EQS of Lake Varano

The nematode assemblage in Lake Varano appeared to be mainly represented by colonizer species. This group has been reported to increase in abundance under stress conditions (Bongers et al. 1991Bongers T., Alkemade R., Yeates G.W. 1991. Interpretation of disturbance-induced maturity decrease in marine nematode assemblages by means of the Maturity Index. Mar. Ecol. Prog. Ser. 76: 135-142.). However, we must highlight that the study area did not show higher values of trace elements when compared with the Effect Range-Low and Effect Range-Median values (Long et al. 1995Long E.R., MacDonald D.D., Smith S.L., et al. 1995. Incidence of adverse biological effects within ranges of chemical concentrations in marine and estuarine sediments. Environ. Manage. 19: 81-97.), as well as high OM enrichment.

Although the MI may generally give good results in ecological assessments (Balsamo et al. 2012Balsamo M., Semprucci F., Frontalini F., et al. 2012. Meiofauna as a tool for marine ecosystem biomonitoring. In: Cruzado A. (ed.), Marine Ecosystems. In Tech Publisher: pp. 77-104., Semprucci et al. 2013Semprucci F., Frontalini F., Covazzi-Harriague A., et al. 2013. Meio- and Macrofauna in the marine area of the Monte St. Bartolo Natural Park (Central Adriatic Sea, Italy). Sci. Mar. 77: 189-199.), Armenteros et al. (2009)Armenteros M., Ruiz-Abierno A., Fernández-Garce R., et al. 2009. Biodiversity patterns of free-living marine nematodes in a tropical bay: Cienfuegos, Caribbean Sea. Estuar. Coast. Shelf. Sci. 85: 179-189. reported unclear results in Cienfuegos Bay (Caribbean Sea), especially for the interaction of anthropogenic and natural factors. A clear relationship between anthropogenic activity (PLI) and the MI was not found in Lake Varano. This could be explained by the environmental complexity of transitional habitats as well as the necessity to update the life style of some species.

In order to classify the EQS of Lake Varano in accordance with the Water Framework Directive (WFD, Directive 2000/60/EC), we used the thresholds for nematode descriptors proposed by Moreno et al. (2011)Moreno M., Semprucci F., Vezzulli L., et al. 2011. The use of nematodes in assessing ecological quality status in the Mediterranean coastal ecosystems. Ecol. Ind. 11: 328-336.. This is the only tool available to classify the EQS by means of marine nematodes and it was integrated considering the percentage of all the c-p classes (Semprucci et al. 2014Semprucci F., Colantoni P., Sbrocca C., et al. 2014. Spatial patterns of distribution of meiofaunal and nematode assemblages in the Huvadhoo lagoon (Maldives, Indian Ocean). J. Mar. Biol. Ass. UK. ) (Fig. 5). On the basis of these thresholds, the EQS of most of the stations can be mainly classified as moderate, followed by good and poor, which is in agreement with the level of organic enrichment and the trace element data (Fig. 5). In detail, moderate conditions of the assemblage were generally revealed in the central area of the lake, while in front of the San Antonino and Francesco Canals, where Specchiulli et al. (2010)Specchiulli A., Renzi M., Scirocco T., et al. 2010. Comparative study based on sediment characteristics and macrobenthic communities in two Italian lagoons. Environ. Monit. Assess. 160: 237-256. also documented the negative influence of the water discharges of the canals, a poor EQS was detected. The notable difference revealed between the western and eastern part of the lake could be related to the stronger influence of the Capoiale Channel and the poor mixing action of the Varano Lake, and this fits perfectly with the findings of Specchiulli et al. (2008)Specchiulli A., Focardi S., Renzi M., et al. 2008. Environmental heterogeneity patterns and assessment of trophic levels in two Mediterranean lagoons: Orbetello and Varano, Italy. Sci. Tot. Environ. 402: 285-298. and Frontalini et al. (2013)Frontalini F., Margaritelli G., Francescangeli F., et al. 2013. Benthic foraminiferal assemblages and biotopes in a coastal lake: the case study of Lake Varano (southern Italy). Acta. Protozool. 52: 147-161.. Although temporal replications are lacking, this research enables there to be a preliminary evaluation of the current ecological conditions of the lake. It may also provide a baseline for the future monitoring of this TE and for assessing long-term changes therein.

AcknowledgEmentsTop

We warmly thank the two anonymous referees for their constructive and critical reviews of the manuscript and Dr. Anabella Harriague-Covazzi for the translation of the abstract into Spanish.

ReferencesTop

Airoldi L., Beck M.W. 2007. Loss, status and trends for coastal marine habitats of Europe. Oceanogr. Mar. Biol. Annu. Rev. 45: 345-405.

Albertelli G., Covazzi-Harriague A., Danovaro R., et al. 1999. Differential responses of bacteria, meiofauna and macrofauna in a shelf area (Ligurian Sea, NW Mediterranean): role of food availability. J. Sea. Res. 42: 11-26.
http://dx.doi.org/10.1016/S1385-1101(99)00012-X

Appeltans W. et al. 2012. The Magnitude of Global Marine Species Diversity. Curr. Biol. 22: 2189-2202.
http://dx.doi.org/10.1016/j.cub.2012.09.036

Armenteros M., Ruiz-Abierno A., Fernández-Garce R., et al. 2009. Biodiversity patterns of free-living marine nematodes in a tropical bay: Cienfuegos, Caribbean Sea. Estuar. Coast. Shelf. Sci. 85: 179-189.
http://dx.doi.org/10.1016/j.ecss.2009.08.002

Armynot du Châtelet E., Debenay J.P., Soulard R. 2004. Foraminiferal proxies for pollution monitoring in moderately polluted harbours. Environ. Pollut. 127: 27-40.
http://dx.doi.org/10.1016/S0269-7491(03)00256-2

Armynot du Châtelet E., Bout-Roumazeilles V., Coccioni R., et al. 2013. Environmental control on shell structure and composition of agglutinated foraminifera along a proximal-distal transect in the Marmara Sea. Mar. Geol. 335: 114-128.
http://dx.doi.org/10.1016/j.margeo.2012.10.013

Balsamo M., Albertelli G., Ceccherelli V.U., et al. 2010. Meiofauna of the Adriatic Sea: current state of knowledge and future perspective. Chem. Ecol. 26: 45-63.
http://dx.doi.org/10.1080/02757541003705492

Balsamo M., Semprucci F., Frontalini F., et al. 2012. Meiofauna as a tool for marine ecosystem biomonitoring. In: Cruzado A. (ed.), Marine Ecosystems. In Tech Publisher: pp. 77-104.
http://dx.doi.org/10.5772/34423

Barnes N., Bamber R.N., Moncrieff C.B., et al. 2008. Meiofauna in closed coastal saline lagoons in the United Kingdom: Structure and biodiversity of the nematode assemblage. Estuar. Coast. Shelf. Sci. 79: 328-340.
http://dx.doi.org/10.1016/j.ecss.2008.03.017

Belmonte G., Scirocco T., Denitto F. 2011. Zooplankton composition in Lake Varano (Adriatic Sea coast, Italy). Ital. J. Zool. 78: 370-378.
http://dx.doi.org/10.1080/11250003.2011.561261

Beneduce L., Vernile A., Spano G., et al. 2010. Occurrence of Vibrio vulnificus in mussel farms from the Varano lagoon environment. Lett. Appl. Microbiol. 51: 443-449.
http://dx.doi.org/10.1111/j.1472-765X.2010.02917.x

Bernardello M., Secco T., Pellizzato F., et al. 2006. The changing state of contamination in the lagoon of Venice. Part 2: heavy metals. Chemosphere 64: 1334-1345.
http://dx.doi.org/10.1016/j.chemosphere.2005.12.033

Bongers T. 1990. The maturity index: an ecological measure of environmental disturbance based on nematode species composition. Oecologia 83: 14-19.
http://dx.doi.org/10.1007/BF00324627

Bongers T., Alkemade R., Yeates G.W. 1991. Interpretation of disturbance-induced maturity decrease in marine nematode assemblages by means of the Maturity Index. Mar. Ecol. Prog. Ser. 76: 135-142.
http://dx.doi.org/10.3354/meps076135

Danovaro R., Gambi C., Mirto S., et al. 2004. Meiofauna. Biol. Mar. Mediterr. 11: 55-97.

Deprez T. et al. 2005. NeMys. World Wide Web electronic publication. www.nemys.ugent.be

Fabbrocini A., Guarino A., Scirocco T., et al. 2005. Integrated biomonitoring assessment of the Lesina Lagoon (Southern Adriatic Coast, Italy): preliminary results. Chem. Ecol. 21: 479-489.
http://dx.doi.org/10.1080/02757540500438623

Fabbrocini A., Di Stasio M., D’Adamo R. 2010. Computerized sperm motility analysis in toxicity bioassays: A new approach to pore water quality assessment. Ecotoxicol. Environ. Saf. 73: 1588-1595.
http://dx.doi.org/10.1016/j.ecoenv.2010.05.003

Ferrero T.J., Debenham N.J., Lambshead P.J.D. 2008. The nematodes of the Thames estuary: Assemblage structure and biodiversity, with a test of Attrill’s linear model. Estuar. Coast. Shelf. Sci. 79: 409-418.
http://dx.doi.org/10.1016/j.ecss.2008.04.014

Forster S.J. 1998. Osmotic stress tolerance and osmoregulation of intertidal and subtidal nematodes. J. Exp. Mar. Biol. Ecol. 224: 109-125.
http://dx.doi.org/10.1016/S0022-0981(97)00192-5

Frontalini F., Coccioni R. 2011. Benthic foraminifera as bioindicators of pollution: a review of Italian research over the last three decades. Rev. Micropal. 54: 115-127.
http://dx.doi.org/10.1016/j.revmic.2011.03.001

Frontalini F., Coccioni R., Bucci C. 2010. Benthic foraminiferal assemblages and trace element contents from the lagoons of Orbetello and Lesina. Environ. Monit. Assess. 170: 245-260.
http://dx.doi.org/10.1007/s10661-009-1229-6

Frontalini F., Semprucci F., Coccioni R., et al. 2011. On the quantitative distribution and community structure of the meio and macrofaunal communities in the coastal area of the Central Adriatic Sea (Italy). Environ. Monit. Assess. 180: 325-344.
http://dx.doi.org/10.1007/s10661-010-1791-y

Frontalini F., Margaritelli G., Francescangeli F., et al. 2013. Benthic foraminiferal assemblages and biotopes in a coastal lake: the case study of Lake Varano (southern Italy). Acta. Protozool. 52: 147-161.

Frontalini F., Semprucci F., Armynot du Châtelet E., et al. 2014. Biodiversity trends of the meiofauna and foraminifera assemblages of Lake Varano (southern Italy). Proc. Biol. Soc. Wash. 127.

Gambi C., Totti C., Manini E. 2003. Impact of organic loads and environmental gradients on microphytobenthos and meiofaunal distribution in a coastal lagoon. Chem. Ecol. 19: 207-223.

Gollner S., Riemer B., Martinez Arbizu P., et al. 2010. Diversity of Meiofauna from the East Pacific Rise across a Gradient of Hydrothermal Fluid Emissions. PLoS ONE 5: e12321.
http://dx.doi.org/10.1371/journal.pone.0012321

Graf G. 1992. Benthic-pelagic coupling: A benthic view. Oceanogr. Mar. Biol. Annu. Rev. 30: 149-190.

Guerrini A., Colangelo M.A., Ceccherelli V.U. 1998. Recolonization patterns of meiobenthic communities in brackish vegetated and unvegetated habitats after induced hypoxia/anoxia. Hydrobiologia 375/376: 73-87.
http://dx.doi.org/10.1023/A:1017096603808

Guilini K., Levin L.A., Vanreusel A. 2012. Cold seep and oxygen minimum zone associated sources of margin heterogeneity affect benthic assemblages, diversity and nutrition at the Cascadian margin (NE Pacific Ocean). Prog. Oceanogr. 96: 77-92.
http://dx.doi.org/10.1016/j.pocean.2011.10.003

Heip C., Vincx M., Vranken G. 1985. The ecology of marine nematodes. Oceanogr. Mar. Biol. Annu. Rev. 23: 399-489.

Hendelberg M., Jensen P. 1993. Vertical distribution of the nematode fauna in coastal sediment influenced by seasonal hypoxia in the bottom water. Ophelia 37: 83-94.
http://dx.doi.org/10.1080/00785326.1993.10429909

Lambshead P.J.D., Tietjen J., Ferrero T.J., et al. 2000. Latitudinal diversity gradients in the deep sea with special reference to North Atlantic nematodes. Mar. Ecol. Prog. Ser. 194: 159-167.
http://dx.doi.org/10.3354/meps194159

Long E.R., MacDonald D.D., Smith S.L., et al. 1995. Incidence of adverse biological effects within ranges of chemical concentrations in marine and estuarine sediments. Environ. Manage. 19: 81-97.
http://dx.doi.org/10.1007/BF02472006

Martins V.A, Frontalini F., Figueira R.C.L., et al. 2013. Assessment of the health quality of Ria de Aveiro (Portugal): heavy metals and benthic foraminifera. Mar. Pollut. Bull. 70: 18-33.
http://dx.doi.org/10.1016/j.marpolbul.2013.02.003

Moreno M., Semprucci F., Vezzulli L., et al. 2011. The use of nematodes in assessing ecological quality status in the Mediterranean coastal ecosystems. Ecol. Ind. 11: 328-336.
http://dx.doi.org/10.1016/j.ecolind.2010.05.011

Muthumbi A.W., Vanreusel A., Duineveld G., et al. 2004. Nematode community structure along the continental slope off the Kenyan Coast, Western Indian Ocean. Int. Rev. Hydrobiol. 89: 188-205.
http://dx.doi.org/10.1002/iroh.200310689

Nicholas W.L., Goodchild D.J., Steward A. 1987. The mineral composition of intracellular inclusions in nematodes from thiobiotic mangrove mud-flats. Nematologica 33: 167-179.
http://dx.doi.org/10.1163/187529287X00308

Pallo P., Widbom B., Olafsson E. 1998. A quantitative survey of the benthic meiofauna in the Gulf of Riga (eastern Baltic Sea), with special reference to the structure of nematode assemblages. Ophelia 49: 117-139.
http://dx.doi.org/10.1080/00785326.1998.10409377

Pfannkuche O., Thiel H. 1988. Sample processing. In: Higgins R.P., Thiel H. (eds), Introduction to the study of meiofauna. Smithsonian Inst., Washington, pp. 134-145.

Pusceddu A., Gambi C., Manini E., et al. 2007. Trophic state, eco system efficiency and biodiversity of transitional aquatic ecosystems: analysis of environmental quality based on different benthic indicators. Chem. Ecol. 23: 505-515.
http://dx.doi.org/10.1080/02757540701760494

Schratzberger M., Warwick R.M. 1998. Effects of the intensity and frequency of organic enrichment on two estuarine nematode communities. Mar. Ecol. Prog. Ser. 164: 83-94.
http://dx.doi.org/10.3354/meps164083

Schratzberger M., Bolam S.G., Whomersley P., et al. 2003. Development of a meiobenthic nematode community following the intertidal placement of various types of sediment. J. Exp. Mar. Biol. Ecol. 303: 79-96.
http://dx.doi.org/10.1016/j.jembe.2003.11.003

Seinhorst J.W. 1959. A rapid method for the transfer of nematodes from fixative to anhydrous glycerine. Nematologica 4: 67-69.
http://dx.doi.org/10.1163/187529259X00381

Semprucci F. 2013. Marine nematodes from the shallow subtidal coast of the Adriatic Sea: species list and distribution. Int. J. Biodiv. 1: 1-9.
http://dx.doi.org/10.1155/2013/187659

Semprucci F., Balsamo M. 2012. Key role of free-living nematodes in the marine ecosystem. In: Boeri F., Jordan A.C. (eds). Nematodes: Morphology, Functions and Management Strategies. NOVA Science Publishers, Inc. Hauppauge, New York, pp. 109-134.

Semprucci F., Colantoni P., Baldelli G., et al. 2010a. The distribution of meiofauna on back-reef sandy platforms in the Maldives (Indian Ocean). Mar. Ecol. 31: 592-607.
http://dx.doi.org/10.1111/j.1439-0485.2010.00383.x

Semprucci F., Boi P., Manti A., et al. 2010b. Benthic communities along a littoral of the Central Adriatic Sea (Italy). Helgol. Mar. Res. 64: 101-115.
http://dx.doi.org/10.1007/s10152-009-0171-x

Semprucci F., Frontalini F., Covazzi-Harriague A., et al. 2013. Meio- and Macrofauna in the marine area of the Monte St. Bartolo Natural Park (Central Adriatic Sea, Italy). Sci. Mar. 77: 189-199.
http://dx.doi.org/10.3989/scimar.03647.26A

Semprucci F., Colantoni P., Sbrocca C., et al. 2014. Spatial patterns of distribution of meiofaunal and nematode assemblages in the Huvadhoo lagoon (Maldives, Indian Ocean). J. Mar. Biol. Ass. UK.
http://dx.doi.org/10.1017/S002531541400068X

Spagnoli F., Specchiulli A., Scirocco T., et al. 2002. The Lago di Varano: Hydrologic characteristics and sediment composition. Mar. Ecol. 23: 384-394.
http://dx.doi.org/10.1111/j.1439-0485.2002.tb00036.x

Specchiulli A., Focardi S., Renzi M., et al. 2008. Environmental heterogeneity patterns and assessment of trophic levels in two Mediterranean lagoons: Orbetello and Varano, Italy. Sci. Tot. Environ. 402: 285-298.
http://dx.doi.org/10.1016/j.scitotenv.2008.04.052

Specchiulli A., Renzi M., Scirocco T., et al. 2010. Comparative study based on sediment characteristics and macrobenthic communities in two Italian lagoons. Environ. Monit. Assess. 160: 237-256.
http://dx.doi.org/10.1007/s10661-008-0691-x

Steyaert M., Vanaverbeke J., Vanreusel A., et al. 2003. The importance of fine-scale, vertical profiles in characterizing nematode community structure. Estuar. Coast. Shelf. Sci. 58: 353-366.
http://dx.doi.org/10.1016/S0272-7714(03)00086-6

Steyaert M., Moodley L., Vanaverbeke J., et al. 2005. Laboratory experiments on the infaunal activity of intertidal nematodes. Hydrobiologia 540: 217-223.
http://dx.doi.org/10.1007/s10750-004-7145-4

Steyaert M., Moodley L., Nadong T., et al. 2007. Responses of intertidal nematodes to short-term anoxic events. J. Exp. Mar. Biol. Ecol. 345: 175-184.
http://dx.doi.org/10.1016/j.jembe.2007.03.001

Timm R.W. 1963. Marine nematodes of the suborder Monhysterina from the Arabian Sea at Karachi. Proc. Helminthol. Soc. Wash. 30: 34-49.

Travizi A., Vidakovic J. 1997. Nematofauna in the Adriatic Sea: Review and check-list of free-living nematode species. Helgol. Meeresunters. 51: 503-519.
http://dx.doi.org/10.1007/BF02908730

Trentesaux A., Recourt P., Bout-Roumazeilles V., et al. 2001. Carbonate grain-size distribution in hemipelagic sediments from a laser particle sizer. J. Sediment. Res. 71: 858-862.
http://dx.doi.org/10.1306/2DC4096E-0E47-11D7-8643000102C1865D

Vanaverbeke J., Gheskiere T., Steyaert M., et al. 2002. Nematode assemblages from subtidal sandbanks in the southern bight of the north sea: effect of small sedimentological differences. J. Sea. Res. 48: 197-207.
http://dx.doi.org/10.1016/S1385-1101(02)00165-X

Vanaverbeke J., Merckx B., Degraer S., et al. 2011. Sediment-related distribution patterns of nematodes and macrofauna: Two sides of the benthic coin? Mar. Environ. Res. 71: 31-40.
http://dx.doi.org/10.1016/j.marenvres.2010.09.006

Villano N., Warwick R.M. 1995. Meiobenthic communities associated with the seasonal cycle of growth and decay of Ulva rigida Arardh in the Palude Della Rosa, Lagoon of Venice. Estuar. Coast. Shelf. Sci. 4: 181-194.
http://dx.doi.org/10.1006/ecss.1995.0060

Warwick R.M. 1971. Nematode associations in the Exe Estuary. J. Mar. Biol. Assoc. U.K. 51: 430-454.
http://dx.doi.org/10.1017/S0025315400031908

Warwick R.M., Gee J.M. 1984. Community structure of estuarine meiobenthos. Mar. Ecol. Prog. Ser. 18: 97-111.
http://dx.doi.org/10.3354/meps018097

Wetzel M.A., Weber A., Giere O. 2002. Re-colonization of anoxic/ sulfidic sediments by marine nematodes after experimental removal of macroalgal cover. Mar. Biol. 141: 679-689.
http://dx.doi.org/10.1007/s00227-002-0863-0

Wieser W. 1953. Die Beziehung zwischen Mundhöhlengestalt. Ernährungsweise und Vorkommen bei freilebenden marinen Nematoden. Eine okologisch-morphologische studie. Ark. Zool. 4: 439-484.

SUPPLEMENTARY MATERIAL

The following material is available through the online version of this article and at the following link:
http://www.icm.csic.es/scimar/supplm/sm04018esm.pdf

Table S1. – Geographic coordinates of the sampling stations, water depth and the physicochemical parameters of the bottom water and sediment characteristics.

Table S2. – Percentages of the nematode species found in the study area.

Table S3. – Data matrix used for the statistical analysis.