Biodiversity of Talitridae family (Crustacea, Amphipoda) in some Tunisian coastal lagoons

Abstract

Background

Although wetlands were remarkable habitats with their fauna and flora diversity, few studies have been devoted to the study of amphipod distribution in this type of environment.

To study both qualitatively and quantitatively amphipod community, surveys were conducted during the spring season in ten coastal lagoons ranging from subhumid to arid bioclimatic stage. At each station, eight quadrats of 50 × 50 cm were randomly placed. Amphipods were preserved in alcohol 70°C. In the laboratory, the specimens collected were identified and counted. Meanwhile, analyses of organic matter, particle size, and heavy metals from the soil taken from each station were made.

Results

A total of 1,340 specimens of amphipods were collected, and eight species belonging to Talitridae family were identified.

Species richness ranges from one species collected in the supralittoral zone of El Bcherliya (Ghar El Melh lagoon) and eight species in the supralittoral zone of Bizerte lagoon. In this last station, the relative abundance of amphipods was significantly higher (36.04%, N = 483). In addition, the diversity indices of Simpson, Shannon-Weaver, and equitability shows that the highest species diversity characterizes this same station while the community was more balanced in opposite El Boughaz (Ghar El Melh lagoon) (J″ = 0.996).

Conclusions

The spatial distribution of different amphipod species depends on edaphic (heavy metals, granulometry, organic matter) and climatic (temperature, humidity) factors.

Background

Coastal areas, the natural interface between water and land, were valuable ecosystems since they host high biodiversity levels (Defeo and McLachlan 2005; McLachlan and Brown 2006). Many human activities, such as fishing, land reclamation, engineering, shipping, and recreational activities, affect coastal ecosystems (McLachlan and Brown 2006; Schlacher et al. 2007). The impact of pollution in coastal areas was generally strong and the contaminants of main concern include persistent organic pollutants, oil, radionuclides, fertilizers, trace metals, and pathogens (Islam and Tanaka 2004). Trace metals were toxic and nondegradable elements that have adverse effects on living organisms (e.g., immunodeficiency, negative effect on metabolic processes and cell membrane permeability) (Ikem and Egiebor 2005). Many of these elements occur naturally, but their input has been enhanced through the millennia by various human activities (mining, metallurgical and tanning industries, chemical plants, paper mills) (Islam and Tanaka 2004).

In the Mediterranean, there was a high diversity of wetlands (lagoon, lake, sebkhas, oued, hill reservoir and dam) that were of great importance in conservation biology, and they were considered among the most biologically diverse and productive ecosystems (Medail and Quezel 1999). They offer a wide variety of natural habitats for plants and aquatic animals as well as semiterrestrial and terrestrial. The interactions of biological (plants, animals, microorganisms, etc.) and physicochemical components (granulometry, temperature, humidity, etc.) of wetlands enable them to perform many ecological functions such as shoreline stabilization and water purification.

In Tunisia, semi-closed shallow lagoons were among the most sensitive areas to environmental stresses (Benrejeb-Jenhani and Romdhane 2002). In general, lagoon sediments were considered reservoirs of many chemical pollutants, especially heavy metals, which represent the most prominent marine pollution agents (Phillips and Rainbow 1994) affecting local communities and human health (Förstner and Wittmann 1981; Boucheseiche et al. 2002). Hence, quality management of marine coastal environments becomes a priority for many countries.

Peracarid crustaceans have received special attention, because many of them were important components of soft sediment faunas (Dauvin et al. 1994) and were considered good indicators of water and sediment quality (Corbera and Cardell 1995; Alfonso et al. 1998). Further, detritivorous peracarids play an important role in the degradation of organic matter, in both aquatic and terrestrial habitats. Talitrid amphipods were key species in the energy flow of sandy shore ecosystems (Griffiths et al. 1983). Feeding on terrestrial and marine material, these species integrate the different sources of contamination and constitute an important source of food for many species of invertebrates, fishes, and birds (Griffiths et al. 1983; Wildish 1988; Bergerard 1989; Koch 1989). This group constitutes one of the dominant macrofaunal groups in sandy beaches (Dahl 1946; McLachlan and Jaramillo 1995). Their ecological relevance has justified worldwide studies, for instance with respect to their behavioral plasticity (Scapini and Fasinella 1990; Scapini et al. 1993; Scapini et al. 1995), locomotor activity rhythms (Nasri-Ammar and Morgan 2005, 2006; Ayari and Nasri-Ammar 2012a, b; Jelassi and Nasri-Ammar 2013; Jelassi et al. 2013a, b), the factors influencing their spatial distribution and oriented movements in sandy beach (Scapini and Quochi 1992; Borgioli et al. 1999; Scapini et al. 1999; Ayari and Nasri-Ammar 2011; Jelassi et al. 2012; Jelassi et al. 2013c), their behavioral strategies (Fallaci et al. 1999), and genetic structure of different populations (De Matthaeis et al. 1995; Bulnheim and Schwenzer 1999). With regard to biodiversity, along the European coasts, talitrid populations have been compared genetically to assess inter- and intraspecific variations (De Matthaeis et al. 1995). In terms of applied research, a number of articles have been published on trace metal (Cu, Zn, Fe, Cd, Pb, Mn, and Ni) concentrations and bioaccumulation by talitrids and on their role in biomonitoring (Fialkowski et al. 2000; Rainbow et al. 1989; Weeks 1992). Studying the amphipod diversity in three different complexes, lagoons in Northern Tunisia, Jelassi et al. (2013c) showed that the abundance of different amphipod species at the three lagoon complexes in northern Tunisia could be best explained by the soil contents of several heavy metals.

In Tunisia, amphipod communities inhabiting the supralittoral zone of wetlands, other than sandy beaches (Jelassi and Nasri-Ammar 2013; Jelassi et al. 2012; Jelassi et al. 2013a), have not received much attention. Thus, aims of this study were to produce a list of amphipod species collected in the supralittoral zones of different lagoons, to estimate species richness, abundance, and density, and we focus on environmental factors that may control their spatial distribution.

Methods

Study stations

In total, ten stations belonging to the supralittoral zone of eight coastal lagoons distributed along the Tunisian coasts were studied (Figure 1). The lagoon is a stretch of salt water partially or completely separated from the open ocean by barriers of sand. These stations were S1: Bizerte lagoon (S1: 37°13′8″N009°55′1″E), three stations belonging to Ghar E Melh lagoon (S2: El Bcherliya (37°10′03″N 010°09′57″E), S3: the old harbor (37°10′04″N 010°11′40″E), S4: opposite El Boughaz (37°10′09″N 010°13′12″E), S5: Sidi Ali Mekki lagoon (37°09′50″N 010°14′45″E)), S6: North (36°48′01″N 010°12′27″E) and S7: South (36°47′59″N 010°12′26″E) lagoons of Tunis, S8: Korba lagoon (36°38′12″N 010°54′11″E), S9: Tazarka lagoon (36°32′20″N 010°50′38″E), and S10: El Bibane lagoon (33°15′57″N 011°08′28″E).

Figure 1

Localization of sampling lagoon stations. S1: Bizerte, S2: El Bcherliya, S3: Old harbor of Ghar El Melh, S4: Opposite El Boughaz, S5: Sidi Ali Mekki, S6: North lagoon of Tunis, S7: South lagoon of Tunis, S8: Korba, S9: Tazarka, S10: El Bibane.

Sampling method and soil analysis

Quantitative samples of amphipods were manually collected in spring (April 2010) during ten consecutive days in the early morning hours. In the supralittoral zone of each site, eight quadrats of 50 × 50 cm were randomly placed and 20 min were devoted to collect amphipods inside the quadrat. The content of each quadrat (7 cm depth) was placed in individual boxes. Humidity and temperature of soil were measured in situ at each site. In the laboratory, specimens were preserved in 70% ethanol. Afterwards, they were identified and counted. The identification of these species was carried out under a Leica MS 5 binocular microscope (Leica Microsystems, Wetzlar, Germany), using the key of Ruffo (1993) and Ruffo et al. (2014) (for the new name of Orchestia cavimana).

At each station, a soil sample was taken from a depth of 0 to 10 cm. Grain size distribution of these composite samples was analyzed using different sieves in descending order (from 2 mm to 25 μm). A subsample was brought to the ICP-MS laboratory at University of Kiel and sieved to obtain the <250-μm grain size fraction which was then dried and milled (Bat and Raffaelli 1999). Heavy metals were extracted from a 250-mg sample of powder with 10 mL 7 N nitric acid on a hot plate at 80°C (2.5 h). The solution was made up to 20 mL, centrifuged at 3,500 rpm for 15 min, and the supernatant transferred to a 20-mL sample vial. The metals vanadium (V), chromium (Cr), manganese (Mn), cobalt (Co), nickel (Ni), copper (Cu), zinc (Zn), arsenic (As), cadmium (Cd), tin (Sn), thallium (Tl), lead (Pb), lithium (Li), rubidium (Rb), and strontium (Sr) were analyzed by inductively coupled plasma-mass spectrometry (ICP-MS). Average analytical reproducibility was estimated from replicate analyses of some samples and was found to be better than 2% RSD (1 sigma relative standard deviation) for all elements. The accuracy of analytical results was monitored by analyzing certified reference materials (CRM): GSMS-2 (marine sediment; Chinese Academy of Geological Sciences, PR China) and PACS-1 (coastal sediment; NRCC Canada) as unknowns along with the samples. Organic matter content was determined by weighing before and after ashing at 450°C for 3 h at the University of Salzburg.

Soil sodium content was determined by atomic absorption technique based on the methodology of the IC2MP (Institute of Chemistry of Poitiers: materials and natural resources). To obtain a usable sample, we put 5 g of soil in 50 ml of distilled water. After 1 h, the solution was filtered through a 50-lm mesh. To obtain reliable measurements, it was necessary to carry out successive dilutions.

Data analysis

In order to characterize the amphipod communities, the following ecological parameters were calculated: species richness (S) expressed by the number of species at each station, occurrence frequency (F = (n/N) * 100) where ‘n’ was the number of times that species appears in the sample and ‘N’ was the total number of samples, and relative species abundance (A _r = (n _i/N) * 100) where ‘n _i’ was the number of individuals of each species and ‘N’ was the total number of amphipods at each station. Mean density of the amphipod community at each station and the mean density of each species at each station were expressed as number of individuals per m².

Species diversity and evenness were calculated by the Simpson index (Is = 1/∑ p _i ²) with p _i = (n _i/N) (N = total number of amphipods; n _i = the number of individuals for each species) (Simpson 1949), the Shannon-Weaver index (H′ = −Σ p _i log₂ p _i) (Frontier 1983), and Pielou's evenness index (J′ = H′/ log₂S) (Pielou 1966).

The degree of similarity among sampling stations was evaluated using similarity cluster dendrograms. For this analysis, the data matrix consisting of total abundances of species at each site was converted into a symmetric matrix using the Bray-Curtis similarity index. The similarity matrix was agglomerately clustered using complete linkage based on presence/absence of species. The analysis above was performed with the PRIMER software package (Clarke and Warwick 1994).

Differences in abundance, species richness, and densities of amphipod species among lagoons were tested using ANOVA test. Principal correspondence analysis (PCA) of amphipod distribution was performed using the free version of XLSTAT.

Results

Heavy metals, sodium, organic matter, and granulometry of different stations

The highest concentration of the majority of heavy metals, namely vanadium, nickel, zinc, arsenic, cadmium, thallium, and lead, was observed in the supralittoral zone of Tunis North lagoon. For the chromium and manganese, the maximum of content characterizes the supralittoral zone of Bizerte lagoon with respectively 26.393 ppm and 281.748 ppm. The maximum content in copper (39.098 ppm) was observed in the supralittoral zone of El Bcherliya. The supralittoral zone of Korba lagoon revealed the highest concentration in cobalt (8.311 ppm) and rubidium (15.814 ppm). Whereas, the supralittoral zone of El Bibane lagoon was characterized by the highest content in lithium (29.087 ppm), in strontium (2,101.549 ppm), and in tin (7.340 ppm). Furthermore, the lowest concentration for all heavy metals studied was observed in the supralittoral zone of Sidi Ali Mekki.

The different heavy metals analyzed do not exceed the maximum value tolerated except for the lead that exceeded 100 ppm in the supralittoral zone of Tunis North lagoon (Henin 1983).

In these different stations, soil sodium content ranged between 1.26 mg/g and 7.2 mg/g of soil in the supralittoral zone of Bizerte lagoon and the old harbor of Ghar El Melh, respectively. For the organic matter, the percentage varied between 0.6% in the supralittoral zone of Sidi Ali Mekki and 9.46% in the supralittoral zone of Bizerte lagoon.

Concerning the granulometry, the supralittoral zones of Bizerte lagoon, El Bcherliya, old harbor, Tunis North, and Korba were characterized by loamy sand substrates. The supralittoral zones of Tunis South and El Bibane lagoons were characterized by sandy substrates. The sandy loam, sandy clay loam, and sandy silt loam substrates characterized respectively the supralittoral zones of opposite El Boughaz, Tazarka, and Sidi Ali Mekki lagoons.

Species richness and occurrence frequency

A total of eight species were collected in lagoons namely Orchestia montagui Audouin, 1826, Orchestia mediterranea Costa, 1853, Orchestia gammarellus (Pallas, 1766), Orchestia stephenseni Cecchini, 1928, Cryptorchestia cavimana Heller, 1865, Platorchestia platensis (Kroyer, 1845), Deshayesorchestia deshayesii (Audouin, 1826) and Talitrus saltator (Montagu, 1808). Species richness varied both between quadrats of the same station and between lagoons. It ranged from one species in the supralittoral zone of El Bcherliya and eight species in the supralittoral zone of Bizerte lagoon (Table 1). The differences observed between lagoons were highly significant (F = 3.495, DF = 7, p = 0.004).

Table 1 Species composition, mean density and diversity indices of amphipod community in different lagoons

The study of the occurrence frequency of different species of amphipods on lagoon shores showed that O. gammarellus was qualified as common species (F = 80%) and, O. mediterranea (F = 70%) and O. stephenseni (F = 60%) were classified as constant species. As for P. platensis (F = 50%), O. montagui (F = 40%), and D. deshayesii (F = 30%), they were called accessory species. Finally, T. saltator (F = 20%) and C. cavimana (F = 10%) were qualified as rare species.

Relative abundance and density

A count of 1,340 amphipod individuals was collected in different lagoons. In this type of wetlands, the supralittoral zone of Bizerte lagoon revealed the highest abundance (36.04%) (Table 1). Moreover, the difference between lagoons was statistically significant (Anova test: F = 4.371; DF = 7; p = 0.001).

The study of relative abundance of the different species showed that O. mediterranea was the most abundant species in the supralittoral zone of Bizerte lagoon (25.7%). Whereas, in El Bcherliya, the old harbor of Ghar El Melh, the North and South lagoons of Tunis, O. gammarellus dominated the community. In the supralittoral zone of El Bibane lagoon, O. montagui was the most important species (28.3%). Anova test revealed that differences in species observed between lagoons were highly significant (F = 4.371; DF = 7; p = 0.001).

Concerning the mean density of amphipod community, it varied between lagoons but it was the supralittoral zone of Bizerte lagoon that presented the important density with 241.5 ind. m⁻². The lowest density was observed in the supralittoral zone of El Bcherliya (0.5 ind. m⁻²) (Figure 2).

Figure 2

Variation of the mean density (ind. m − 2) of amphipod communities in the supralittoral zone of lagoons.

Furthermore, results showed that O. mediterranea presented the most important density in the supralittoral zone of Bizerte lagoon (62 ind. m⁻²). Whereas, O. gammarellus showed the most important density in the supralittoral zones of El Bcherliya (0.5 ind. m⁻²), the old harbor of Ghar El Melh (19 ind. m⁻²), the North (34 ind. m⁻²) and South (34.5 ind. m⁻²) lagoons of Tunis.

Diversity and community similarity

According to Simpson index, the highest diversity was observed in the supralittoral zone of Bizerte lagoon (Is = 6.059) (Table 1). Unlike to the previous index, the Shannon-Weaver index takes into account much more the rare species, it confirmed this result and it ranged between 1.287 in the supralittoral zone of Sidi Ali Mekki lagoon and 2.771 in Bizerte lagoon. Concerning the equitability, it was important in opposite El Boughaz where we noted that the species were equitably distributed.

A dendrogram of similarity between lagoons, obtained by calculating the Bray-Curtis similarity index based on presence/absence of species, showed the presence of two clusters: the first associate Sidi Ali Mekki and El Bcherliya where we noted the lowest count of amphipods, the second associate the six others stations and it was divided into two sub-groups; the first combines the opposite El Boughaz and the old harbor of Ghar El Melh, and the other comprises North and South lagoons of Tunis, Bizerte lagoon, and El Bibane lagoon. The combination of these four stations could be explained on the one hand by the presence of the species O. montagui and on the other hand by a considerable concentration of lead exceeding 50 ppm (Figure 3).

Figure 3

Clustering of sampling data, based on total number of amphipods collected (Bray-Curtis index, presence/absence, complete linkage).

Amphipod distribution according to environmental factors

Principal component analysis (PCA) was performed to better understand which of the environmental factors influenced the distribution of different species in lagoons (Figure 4).

Figure 4

Principal Component Analysis (PCA) performed on the environmental factors. V: Vanadium, Cr: Chromium, Mn: Manganese, Co: Cobalt, Ni: Nickel, Cu: Copper, Zn: Zinc, As: Arsenic, Cd: Cadmium, Sn: Tin, Tl: Thallium, Pb: Lead, Li: Lithium, Rb: Rubidium, Sr: Strontium, Tsol: soil temperature, Hsol: soil humidity, Gra: granulometry, Na+: sodium, MO: organic matter.

The first three factorial axes (F1, F2, and F3) extract 54.39%, 19.57%, and 9.32% of the variance respectively; it was a cumulative percentage of 83.28% (Figure 4).

On the first axis (F1) extracting almost half of the total inertia, there were Bizerte lagoon, North and south lagoons of Tunis which were correlated positively with the majority of heavy metals namely vanadium, chromium, manganese, cobalt, nickel, copper, zinc, arsenic, rubidium, cadmium, thallium, and lead. These stations were characterized negatively by the soil temperature.

On the second axis which extracts only 19.57% of the total variance, there were El Bcherliya, El Bibane and Korba lagoons which were related to soil humidity, particle size, and the content of strontium, lithium, and tin.

On the third axis, were projected negatively the organic matter and positively the concentration on sodium. These two parameters do not characterize these different lagoon stations.

Discussion

The study of biodiversity of amphipods communities highlighted differences in species richness, relative abundance and diversity between and within lagoons.

Globally the species richness was higher in the supralittoral zone of Bizerte lagoon than in the other stations. In this station, eight amphipod species were identified. This would be related to the presence of important vegetation in spring as well as the Cymodocea nodosa leaf litter and a high percentage of organic matter. In this same station, the highest diversity, abundance, and density were observed. Moreover, it was O. mediterranea that presented the most important density in this supralittoral zone. In other type of wetlands such as oueds, species richness varied between one species in oueds Laakarit, Khniss, and El Fared and six species in Tinja (Korsi station) whereas in sebkhas, only two species were collected in sebkhas Gargour and Moknine (Jelassi 2014; Jelassi et al. 2013c). In supralittoral zones of dams and hill lakes, no species has been found (Jelassi 2014). Using the same method of collection, Ayari (2012) showed the presence of five species namely T. saltator, D. deshayesii, O. gammarellus, O. montagui, and O. mediterranea in the supralittoral zone of Bizerte Corniche beach and only one species, T. saltator, in the supralittoral zone of Gabès beach. This author attributed the lowest species richness in the last station (Gabès) to the highest tidal amplitude reported in this beach and the highest species richness observed in Bizerte Corniche to the presence of Posidonia oceanica benches associated with other algae that provide shelter and food for different amphipod species. In Zouaraa beach, only two species were identified: T. saltator and Talorchestia brito living in sympatry with a total density generally low and never exceeded 860 ind.m⁻² (Charfi-Cheikhrouha et al. 2000). At Ouderef, Gabes, and Zarrat beaches (Gulf of Gabes, Tunisia), the amphipod T. saltator was the most abundant species, present in the three beaches studied; with the highest densities in Ouderef beach (Pérez-Domingo et al. 2008). O. montagui and D. deshayesii in the Bay of Bou Ismail reached more than 45.000 ind. m⁻² (Louis 1980), and O. mediterranea in the estuary of Bou Regreg, Morocco, attained 7,000 ind. m⁻² (Elkaïm et al. 1985). On the Isle of Man, densities of T. saltator were estimated at 80 to 400 ind. m⁻² (Williams, 1995). Marsden (1991) and Cardoso and Veloso (1996) showed that fluctuations in population density were frequent in talitrids and indicated periods of intense reproduction. They observed, in fact, similar patterns of variation for Talorchestia quoyana and Pseudorchestoidea brasiliensis, respectively, with highest densities in summer and late winter.

The distribution of amphipod species in the different lagoons was investigated according to environmental factors using PCA analysis. Our statistical analyses indicated that the variation in the spatial distribution of amphipod species depended on some climatic (temperature, humidity) as well as edaphic factors, which was influenced by soil quality. In the different lagoons, T. saltator, found only in the supralittoral zone of Bizerte and El Bibane lagoons, was correlated with temperature, humidity, particle size, and some heavy metals. In sebkhas and oueds, this species was correlated only with edaphic factors (unpublished data). Several studies were focused on the role of environmental factors and showed the influence of factors than others (Borgioli et al. 1999; Scapini et al. 1995; Scapini and Fasinella 1990; Scapini et al. 1999; Scapini and Quochi 1992). Jelassi et al. (2012) showed that Talitridae abundance was tightly controlled by air temperature. The highest abundance observed in this site could be also explained by the important recruitment in spring (Jelassi 2014). Studying the amphipod diversity at three Tunisian lagoon complexes, Jelassi et al. (2013c) showed that according to the constrained correspondence analysis, the abundance of O. mediterranea and O. gammarellus at the three lagoon complexes in northern Tunisia could be best explained by the soil contents of several heavy metals, namely zinc, thallium, and cadmium and the proportion of the coarse sand fraction. T. saltator abundance, by contrast, negatively corresponded to these station characteristics. Air and soil temperature were the best predictors for O. stephenseni abundance that negatively corresponded with the proportion of the fine sand fraction and the organic matter content of the soil. O. montagui and C. cavimana abundances corresponded positively with air humidity and the soil lithium and rubidium contents but negatively with the soil tin content and the proportion of the silt and clay fraction. D. deshayesii and P. platensis did not exhibit any clear correspondence with station characteristics. Bouslama et al. (2009) showed that the temperature was the most important factor influencing the zonation; indeed, the increase of temperature induced the migration of T. saltator population from the top to the bottom of beach. This result was similar to that found by Fallaci et al. (2003) that reported that the mean zonation of T. saltator was exclusively influenced by the temperature during the activity period. Marques et al. (2003) have highlighted a positive correlation between temperature and the density or biomass of Talitrus saltator population. These authors showed that the positive correlation may be interpreted as a cause-and-effect relation, with temperature favoring recruitment and, consequently, the increase in density and biomass. On the other hand, Colombini et al. (2002) noted the importance of the parameters of the sediment in the choice of specific distribution area especially for young individuals.

Conclusions

A total of 1,340 talitrids individuals was collected in different coastal lagoons. Abundance, density and diversity were more pronounced in the supralittoral zone of Bizerte lagoon characterized by the presence of eight amphipod species: O. montagui, O. gammarellus, O. mediterranea, O. stephenseni, O. cavimana, P. platensis, D. deshayesii and T. saltator. Moreover, according to principal component analysis, the spatial distribution of species in the different lagoons in Tunisia depends on edaphic (heavy metals, granulometry, organic matter) and climatic (temperature, humidity) factors.