# Variability in the Structure of Phytoplankton Assemblages in relation to Human Disturbance in Southern Coast of Tunisia

Authors: Lotfi Mabrouk; Lamia Dammak; Asma Hamza; Mabrouka Mahfoudhi; Med-Najmeddine Bradai
Journal: Journal of Marine Biology
 (2014)

---

## Abstract

We investigated the impact of industrial effluents on phytoplankton assemblage in southern Tunisia (Skhira). We specifically addressed changes in microalgae composition caused by this anthropogenic interference. A hierarchical sampling design was used to compare planktonic microalgae structure between one disturbed station and one control station. Samples were collected by scuba diving at 5 m depth in August 2012. A total of 76 microalgae taxa were identified. Dinoflagellates abundance was low in the disturbed station, especially Gonyaulacales and Prorocentrales due to P-limitation, whereas diatoms and cyanobacteria abundance were low in control station which is characterized by Si-limitation.


---

## Body

## 1. Introduction

Marine ecosystems are particularly vulnerable to environmental change and many are, at present, severely degraded [1]. The availability of good quality water is an indispensable feature for preserving biodiversity. In recent years, increase in human population and industrial activity have led to faster degradation of many coastal ecosystems. The discharge of urban, industrial, and agricultural wastes has added the quantum of various harmful chemicals to the water body considerably altering their inherent physicochemical characteristics [2].Phytoplankton constitutes the basic components of the aquatic food chain. They act as primary producers and represent themselves as a direct food source for other aquatic animals [3]. Quantifying the phytoplankton community composition is essential for understanding the structure and the dynamics of marine ecosystems [4]. It is also important to evaluate the role of the physical and chemical processes on the biological responses of the ecosystem. The interplay of physical, chemical, and biological properties of water most often leads to the production of phytoplankton, while their assemblage (composition, distribution, diversity, and abundance) is also structured by these factors.It is, therefore, of importance to undertake a study of the effects of urban and industrial pressure on phytoplankton dynamics under the impact of human activities versus control conditions. In southern Tunisia, we analyzed the phytoplankton diversity and abundance and asked if anthropogenic disturbance such as industrial effluents causes their variations.We attempted to answer the following question: how nutrient enrichment induced by anthropogenic disturbance affects abundance and species composition of phytoplankton community? We thus compared phytoplankton assemblages between a station exposed to industrial effluents and a control station and examined their variability at spatial scales using a hierarchical sampling design [5].


## 2. Materials and Methods

### 2.1. Study Area

The study was carried out in two sampling stations (Figure1), 100 km apart, were chosen in August 2012; unpolluted station in the east of Tunisia, Mahdia (M) (35° 35′43′′N 11° 05′42′′E), was selected as control station, and Skhira (K) (34°19′.179′N 10°11′.706′′E) was chosen as disturbed station, localised in the Gulf of Gabes. This station is exposed, since 1988, to industrial effluents from phosphoric acid and fertilizer production of the “Groupe Chimique Tunisien” (GCT). These industrial activities affect this marine ecosystem [6] as in deterioration of seagrass bed [7, 8] and decline in fishing [9]. The sampling depth in all stations was 5 m.Figure 1
Map of the study area, showing the sampling stations with M, Mahdia (control station), and K, Skhira (disturbed station).


### 2.2. Sampling and Data Collection

In accordance with the hierarchical sampling design, three sites (500 m apart) were chosen randomly at each station. At each site, three replicate water samples were collected from each site for nutrient analyses. Samples were taken using a 125 mL plastic bottle, having received prior treatment with hydrochloric acid. On returning to the surface, samples were shaken and then filtered with a 0.45μm filter (cellulose acetate, 17 mm). Samples were frozen in liquid nitrogen for transportation to the laboratory, where concentrations of NO2-, NO3-, NH4+, PO4-, Si(OH)4, total dissolved nitrogen (TN), and total dissolved phosphorus (TP) were measured following standard colorimetric techniques [10].In each site, temperature, salinity, and pH were measured immediately after sampling using a multiparameter kit (Multi 340i/SET). For determination of suspended particulate matter (SM) concentration, water samples were filtered through preweighed GF/C Whatman filters (pore size 0.45 pm), which were subsequently dried at 80°C for about 24 h and reweighed [11].For phytoplankton enumeration, three replicate water samples (about 100 m apart at 5 m depth) were selected from the same sites. Water column samples, using 1-litre glass bottle sampler lowered from the surface to the near bottom, were conducted by scuba diving. All samples were collected at noon.Samples were fixed with Lugol’s solution and finally preserved in 5% formalin. All sampling waters were kept in the dark at ambient temperature until their microscopic observation. Settling long glass tubes used for sedimentation procedure were 2 cm wide by 21 cm long and have a base plate that contains a coverslip on which the algae settle. To mix the sample, the bottle was gently tilted back and forth 10 times before pouring [12]. A 50 mL subsample was poured into the settling chamber and left to settle for 24 h. subsamples were examined in an inverted microscope at medium (× 200) magnification by scanning the entire surface of the settling chamber to enumerate epiphytic microalgae [12, 13]. The total number of microalgae individuals (N) contained in 1 litre expressed as number of individuals per litre is obtained by the following conversion: N = (n  × 1000)/(v), with n = the number of individuals counted and v = the volume of the sedimentation chamber (50 mL) [13]. The identified taxa were divided into groups (diatoms, dinoflagellates, and cyanobacteria).


### 2.3. Data Analysis

Data were tested for normality using the Kolmogorov-Smirnov test [14] and for heteroscedasticity using Cochran’s C test and transformed if necessary [5].Relationships between species abundance and abiotic parameters were examined using theRELATE procedure in PRIMER. RELATE is the equivalent of a nonparametric Mantel test [15]; it assesses the degree of correspondence between matrices, and, via a randomization test, it provides a measure of statistical significance of the relationship [16]: the matrix of similarities between phytoplankton species abundances (based on Bray-Curtis coefficient from Log(x+1)-transformed data) was compared with a matrix of the similarity between abiotic parameters (based on Euclidean distance from Log(x+1)-transformed data). The significance of any correlation between matrices is assessed with a randomization test.Analyses of similarity (ANOSIM) randomization tests were used to test for differences in community composition (with presence/absence transformed data) and for differences in species abundance (with Log(x+1)-transformed data) between control and disturbed stations [17]. Differences found using ANOSIM were followed up using theSIMPER analysis to identify which species primarily accounted for the observed differences between sites. SIMPER generates a ranking of the species responsible for the significant differences. These analyses used a matrix composed of Bray-Curtis similarity coefficient generated with Log(x+1)-transformed species abundance data [13].Analysis of variance (ANOVA) was used to test the hypothesis that the abundance of each group of taxa differed between disturbed and control stations [13]. Factor station was fixed. The same design was employed for nutrient concentrations, pH, temperature, and suspended matter. Tukey's HSD test was applied for multiple comparisons of means.Canonical correspondence analysis (CCA), a direct gradient analysis technique [18], was used to investigate the relationship between microalgae species and environmental parameters. Microalgae abundance data were Log(x+1)-transformed prior to the analysis in order to stabilise the variance and to optimise the signal-to-noise ratio in the data set. Downweighting for rare species was performed. Forward selection and associated Monte Carlo permutation tests (499 unrestricted permutations) were used to identify a subset of environmental variables that contributed mostly to the species abundances in the data set. The environmental parameters which better described the distribution of the species data were a priori identified by forward selection [18]. Only significant environmental variables are included in the model. The same procedure was repeated using species abundance and nutrient ratios (Si : N, N : P, and Si : P). CANOCO 4.5 (Scientia Software) software was used.


## 2.1. Study Area

The study was carried out in two sampling stations (Figure1), 100 km apart, were chosen in August 2012; unpolluted station in the east of Tunisia, Mahdia (M) (35° 35′43′′N 11° 05′42′′E), was selected as control station, and Skhira (K) (34°19′.179′N 10°11′.706′′E) was chosen as disturbed station, localised in the Gulf of Gabes. This station is exposed, since 1988, to industrial effluents from phosphoric acid and fertilizer production of the “Groupe Chimique Tunisien” (GCT). These industrial activities affect this marine ecosystem [6] as in deterioration of seagrass bed [7, 8] and decline in fishing [9]. The sampling depth in all stations was 5 m.Figure 1
Map of the study area, showing the sampling stations with M, Mahdia (control station), and K, Skhira (disturbed station).


## 2.2. Sampling and Data Collection

In accordance with the hierarchical sampling design, three sites (500 m apart) were chosen randomly at each station. At each site, three replicate water samples were collected from each site for nutrient analyses. Samples were taken using a 125 mL plastic bottle, having received prior treatment with hydrochloric acid. On returning to the surface, samples were shaken and then filtered with a 0.45μm filter (cellulose acetate, 17 mm). Samples were frozen in liquid nitrogen for transportation to the laboratory, where concentrations of NO2-, NO3-, NH4+, PO4-, Si(OH)4, total dissolved nitrogen (TN), and total dissolved phosphorus (TP) were measured following standard colorimetric techniques [10].In each site, temperature, salinity, and pH were measured immediately after sampling using a multiparameter kit (Multi 340i/SET). For determination of suspended particulate matter (SM) concentration, water samples were filtered through preweighed GF/C Whatman filters (pore size 0.45 pm), which were subsequently dried at 80°C for about 24 h and reweighed [11].For phytoplankton enumeration, three replicate water samples (about 100 m apart at 5 m depth) were selected from the same sites. Water column samples, using 1-litre glass bottle sampler lowered from the surface to the near bottom, were conducted by scuba diving. All samples were collected at noon.Samples were fixed with Lugol’s solution and finally preserved in 5% formalin. All sampling waters were kept in the dark at ambient temperature until their microscopic observation. Settling long glass tubes used for sedimentation procedure were 2 cm wide by 21 cm long and have a base plate that contains a coverslip on which the algae settle. To mix the sample, the bottle was gently tilted back and forth 10 times before pouring [12]. A 50 mL subsample was poured into the settling chamber and left to settle for 24 h. subsamples were examined in an inverted microscope at medium (× 200) magnification by scanning the entire surface of the settling chamber to enumerate epiphytic microalgae [12, 13]. The total number of microalgae individuals (N) contained in 1 litre expressed as number of individuals per litre is obtained by the following conversion: N = (n  × 1000)/(v), with n = the number of individuals counted and v = the volume of the sedimentation chamber (50 mL) [13]. The identified taxa were divided into groups (diatoms, dinoflagellates, and cyanobacteria).


## 2.3. Data Analysis

Data were tested for normality using the Kolmogorov-Smirnov test [14] and for heteroscedasticity using Cochran’s C test and transformed if necessary [5].Relationships between species abundance and abiotic parameters were examined using theRELATE procedure in PRIMER. RELATE is the equivalent of a nonparametric Mantel test [15]; it assesses the degree of correspondence between matrices, and, via a randomization test, it provides a measure of statistical significance of the relationship [16]: the matrix of similarities between phytoplankton species abundances (based on Bray-Curtis coefficient from Log(x+1)-transformed data) was compared with a matrix of the similarity between abiotic parameters (based on Euclidean distance from Log(x+1)-transformed data). The significance of any correlation between matrices is assessed with a randomization test.Analyses of similarity (ANOSIM) randomization tests were used to test for differences in community composition (with presence/absence transformed data) and for differences in species abundance (with Log(x+1)-transformed data) between control and disturbed stations [17]. Differences found using ANOSIM were followed up using theSIMPER analysis to identify which species primarily accounted for the observed differences between sites. SIMPER generates a ranking of the species responsible for the significant differences. These analyses used a matrix composed of Bray-Curtis similarity coefficient generated with Log(x+1)-transformed species abundance data [13].Analysis of variance (ANOVA) was used to test the hypothesis that the abundance of each group of taxa differed between disturbed and control stations [13]. Factor station was fixed. The same design was employed for nutrient concentrations, pH, temperature, and suspended matter. Tukey's HSD test was applied for multiple comparisons of means.Canonical correspondence analysis (CCA), a direct gradient analysis technique [18], was used to investigate the relationship between microalgae species and environmental parameters. Microalgae abundance data were Log(x+1)-transformed prior to the analysis in order to stabilise the variance and to optimise the signal-to-noise ratio in the data set. Downweighting for rare species was performed. Forward selection and associated Monte Carlo permutation tests (499 unrestricted permutations) were used to identify a subset of environmental variables that contributed mostly to the species abundances in the data set. The environmental parameters which better described the distribution of the species data were a priori identified by forward selection [18]. Only significant environmental variables are included in the model. The same procedure was repeated using species abundance and nutrient ratios (Si : N, N : P, and Si : P). CANOCO 4.5 (Scientia Software) software was used.


## 3. Results

The disturbed station (K) exhibited quite higher concentrations of ammonium, Si(OH)4, TP, and TN and higher pH and temperature. Phosphate was higher in control station (M) (Table 1). We also calculate nutrient ratios; N/P ratio [DIN (DIN = NO2-+NO3-+NH4+) to DIP (DIP = PO43-)] in disturbed station (K) was 75.73±35.43 (mean ± sd) and 15.24±1.22 in control station (M). This average in disturbed station was higher than the Redfield ratio (16), which suggests potential P-limitation in this station. Si/N ratio was 0.41±0.20 and 0.09±0.01 in disturbed and control stations, respectively. Si/P ratio was 26.38±9.99 and 1.39±0.14 in disturbed station and control station, respectively, which suggests potential Si-limitation in control station.Table 1
The mean ± SD values of physical and chemical variables measured at all stations during the sampling period. Difference between control station (M) and disturbed station (K) was tested by ANOVA.
Variables
Unit
N
Disturbed station (K)
Control station (M)
F
(
1,16
)
Tukey’s test
SM
g·L−1
9
12.00 ± 6.76
7.00 ± 0.70
4.32
n.s.
NO
2
-
μmol·L−1
9
0.301 ± 0.12
0.283 ± 0.09
0.119
n.s.
NO
3
-
μmol·L−1
9
4.91 ± 2.85
5.91 ± 0.04
0.98
n.s.
NH
4
+
μmol·L−1
9
23.36 ± 16.80
4.17 ± 0.19
10.42*
M < K
PO
4
3
-
μmol·L−1
9
0.307 ± 0.17
0.684 ± 0.05
25.33***
K < M
Si(OH)4
μmol·L−1
9
8.28 ± 1.93
0.94 ± 0.05
115.41***
M < K
TP
μmol·L−1
9
3.26 ± 0.84
3.34 ± 0.34
9.86*
M < K
TN
μmol·L−1
9
45.29 ± 26.08
16.32 ± 0.47
9.86*
M < K
pH
9
7.97 ± 0.15
7.32 ± 0.34
25.70***
M < K
Salinity
g·L−1
9
38.91 ± 0.44
38.24 ± 0.88
3.96
n.s.
T
°
°C
9
31.03 ± 1.05
27.26 ± 0.96
59.02*
M < K
*
P
<
0.05; ***P<0.001; SM: suspended matter; n.s.: not significant.The results of RELATE tests indicated that there is a correlation between physicochemical parameters and species abundance (Spearman rank correlation statistic, Rho = 0.484). None of the 999 random permutations resulted in a correlation equal to or greater than the measured value of 0.484, indicating that the correlation was significant at theP<0.001 level.For phytoplankton, 76 species were counted including 54 in the control station (M) and 74 species in the disturbed station (K) (Table2). Analysis of similarity (ANOSIM) of phytoplankton microalgae species abundances log⁡(x+1)-transformed showed significant differences (R=0.954; P=0.01) between the control station and the disturbed station. Similar results are obtained (R=0.913; P=0.01) using presence/absence transformed data showing that the dissimilarities between stations are due to both species composition and their abundances.Table 2
Presence-absence species list of planktonic microalgae identified in prospected stations.
Stations
Skhira (K)(disturbed station)
Mahdia (M)(control station)
Dinoflagellates
Alexandrium minutum Halim
1
1
Alexandriumsp.
1
1
Amphidinium carterae Hulburt
1
1
Amphidiniumsp.
1
1
Ceratium candelabrum (Ehrenberg) Stein
1
0
Ceratium furca (Ehrenberg) Claparède & Lachmann
1
1
Ceratium fusus (Ehrenberg) Dujardin
1
0
Ceratium lineatum (Ehrenberg) Cleve
1
0
Coolia monotis Meunier
1
1
Dinophysis rotundata Claparède & Lachmann
1
0
Diplopsalopsis sp.
1
1
Goniodoma polyedricum (Pouchet) Jørgensen
1
1
Gonyaulax polyedra F.Stein
1
1
Gonyaulaxsp.
1
1
Gonyaulax spinifera (Claparède & Lachmann) Diesing
1
1
Gymnodinium veneficum
1
0
Gyrodiniumsp.
1
1
Hermesinumsp.
1
1
Karenia selliformis A.J.Haywood, K.A.Steidinger & L.MacKenzie in Haywood et al.
1
0
karlodinium veneficum
1
1
Ostreopsis ovata Fukuyo
1
1
Peridiniumsp.
1
1
Podolampas palmipes Stein
1
1
Polykrikossp.
1
1
Prorocentrum gracile Schütt
1
1
Prorocentrum lima (Ehrenberg) Dodge
1
1
Prorocentrum micansEhrenberg
0
1
Prorocentrum minimum(Pavillard) Schiller
1
1
Prorocentrum rathymum(Loeblich) Shirley and Schmidt
1
1
Prorocentrum triestinum J.Schiller
1
1
Protoperidinium divergens(Ehrenberg) Balech
1
1
Protoperidinium granii (Ostenfeld) Balech
1
1
Protoperidinium ovatum Pouchet
1
1
Protoperidinium pellucidum Bergh ex Loeblich Jr. & Loeblich III
1
0
Protoperidinium quinquecorne (Abé) Balech
1
1
Protoperidiniumsp.
1
1
Protoperidinium steinii (Jørgensen) Balech
1
1
Protoperidinium  diabolus(Cleve) Balech
1
0
Pyrophacussp.
1
0
Scrippsiellasp.
1
1
Diatoms
Achnanthes sp.
1
0
Amphiprorasp.
1
0
Amphora sp.
1
1
Asterionellopsis glacialis (Castracane) Round
1
0
Bacteriastrum sp.
1
1
Biddulphia sp.
1
1
Bacteriosira fragilis (Gran) Gran
0
0
Coscinodiscussp.
1
0
Chaetocerossp.
1
1
Climacosphenia sp.
1
1
Guinardiasp.
1
1
Grammatophorasp.
1
1
Licmophorasp.
1
1
Leptocylindrussp.
1
1
Melosirasp.
1
1
Naviculasp.
1
1
Nitzschiasp.
1
1
Nitzshoidessp.
1
0
Plagiotropissp.
1
1
Pseudo-nitzschiasp.
1
0
Pinnulariasp.
1
1
Pleurosigmasp.
1
1
Rhizosoleniasp.
1
1
Rhizosolenia fragilissima Bergon
1
1
Rhabdonemasp.
1
1
Striatella unipunctata (Lyngbye) C.Agardh
1
0
Skeletonema costatum (Greville) Cleve
1
0
Thalassionema nitzschioides (Grunow) Mereschkowsky
1
1
Thalassiosirasp.
1
0
Synedra sp.
1
1
Triceratiumsp.
1
0
Cyanobacteria
Anabaenasp.
1
1
Pseudo-anabaenasp.
1
0
Oscillatoriasp.
1
1
Spirulina subsalsa Oersted ex Gomont
1
0
Croccocus sp.
1
0Dinoflagellates, diatoms, and cyanobacteria were common in water column (Figure2). Abundances of those groups were included in the univariate analyses of variance. The number of species (S) was higher in polluted station (F(1;17) = 6.95, P=0.017). Dinoflagellates differed significantly between stations (F(1;17) = 5.66, P=0.029) with high abundance in control station (Tukey’s test). When dinoflagellates were grouped in order, significant differences were detected for Prorocentrales and Gonyaulacales (F(1;17) = 6.17, P=0.024 andF(1; 17) = 12.49, P=0.003, resp.) which were more abundant in control station. Diatoms abundance was higher (F(1;17) = 16.39, P<0.001) in disturbed station (K) and specially centric diatoms (F(1;17) = 7.16, P=0.016). Cyanobacteria abundance was also higher (F(1;17) = 12.61, P=0.002) in disturbed station (K). No significant differences were recorded for total abundance and forH′ index.Figure 2
Average abundance of planktonic microalgae in prospected stations, K: Skhira (disturbed station), M: Mahdia (control station).Analyses of similarity percentage (SIMPER) showed that the average dissimilarity between the control and disturbed groups is high (86.31%). This procedure also allowed determining the species that contribute to this dissimilarity; they areCoscinodiscussp.,Nitzshoidessp.,Chaetocerossp.,Anabaenasp.,Gyrodiniumsp.,Biddulphiasp.,Striatella unipunctata,Thalassiosirasp.,Leptocylindrussp., andScrippsiellasp. Species which were more abundant in disturbed station (K) areBiddulphiasp.,Polykrikos sp.,Pleurosigma sp.,Coolia monotis, andAmphidinium carterae (Table 4). Species which had higher abundance in control station areOstreopsis ovata,Coolia monotis,Prorocentrum lima,Pleurosigma sp.,Peridinium sp.,Pinnularia sp., andPolykrikos sp.CCA analysis using all measured physicochemical parameters (Table5) indicated that the axis I (eigenvalue λ1 = 0.493) and axis II (eigenvalue λ2 = 0.154) expressed 47.5% of the cumulative variance species-environmental variable. Eigenvalues ​​(P value) calculated with the CCA associated with the analysis in Monte Carlo test were used to select four environmental variables statistically significant (P<0.05) that explain better variations (81.88%) of species composition: they are Si(OH)4, NH4+, TN, and suspended matter (33.72%, 9.41%, 10.18%, and 11.86% of total variance, resp.). In the triplot diagram (Figure 3) disturbed station samples (K) were grouped in the right and are correlated with the axis I, which is defined by the environmental variables Si(OH)4, NH4+, TN, and suspended matter (intraset correlation −0.955, −0.613, −0.588, and −0.497, resp.).Figure 3
Diagram of canonical correspondence analysis (CCA) showing the effects of environmental variables on microalgae species ordination according to the first and second axes. Species associated with K (disturbed station) samples are as follows: 4:Amphidinium sp., 24:Prorocentrum minimum, 48:Achnanthes sp., 20:Karenia selliformis, 56:Chaetoceros sp., 60:Licmophora sp., 76:Thalassiosira sp., 55:Coscinodiscus sp., 53:Biddulphiasp., 37:Protoperidinium pellucidum, 61:Leptocylindrus sp., and 66:Plagiotropis sp. Species associated with M (control station) samples are as follows: 9:Coolia monotis, 11:Diplopsalopsis sp., 21:Ostreopsis ovata, 22:Peridinium sp., 28:Prorocentrum lima, 29:Prorocentrum triestinum, 50:Amphorasp., and 68:Pinnulariasp. Species whose coverage and frequency are less than 40% were eliminated.CCA analysis using nutrients' ratios (Table5) indicated that axis I (eigenvalue λ1 = 0.153) and axis II (eigenvalue λ2 = 0.104) expressed 45.20% of the cumulative variance species-environmental variable. Monte Carlo test showed that all nutrient ratios were statistically significant (P=0.002) and explain 45.22% of variation of species ordination detailed as Si : P (55.04%), N : P (24.65%), and Si : P (20.3%) of total explained variance. In the triplot diagram (Figure 4) disturbed station samples (K) were grouped in the right and are correlated with axis I, which is defined by the nutrient ratios N : P, Si : N, and Si : P (intraset correlation 0.720, −0.761, and 0.877, resp.). Most of diatoms species were associated with K samples (disturbed station), while M samples (control station) were associated with dinoflagellates species (Figure 4).Figure 4
Diagram of canonical correspondence analysis (CCA) showing the effects of nutrient ratios on microalgae species ordination according to the first and second axes. Species associated with K (disturbed station) were the dinoflagellates 3:Amphidinium carterae, 19:Karlodinium veneficum, 8:Ceratium candelabrum, 6:Ceratium fusus, 38:Protoperidinium diabolum, 34:Protoperidinium divergens, 20:Karenia selliformis, and 41:Pyrophacus sp., the diatoms 51:Climacosphenia sp., 47:Biddulphia sp., 54:Licmophora sp., 49:Coscinodiscus sp., 50:Chaetoceros sp., and 45:Amphora sp., and the cyanobacteria 71:Thalassiosira sp. Species associated with M (control station) were the dinoflagellates 9:Coolia monotis, 11:Diplopsalopsis sp., 21:Ostreopsis ovata, 22:Peridinium sp., 28:Prorocentrum lima, 31:Protoperidinium sp., 29:Prorocentrum triestinum and the diatoms 45:Asterionellopsis glacialis, 58:Navicula sp. Species whose coverage and frequency were less than 40% were eliminated.


## 4. Discussion

Our results show that the abundance and composition of phytoplankton varied between the disturbed and control stations. Dinoflagellates abundance, especially the orders Prorocentrales and Gonyaulacales, was low in disturbed station contrary to diatoms that were more abundant. Phytoplankton variation between control and disturbed conditions is similar to those found by Ben Brahim et al. [19] in the south of Tunisia (Gulf of Gabes), by Wang et al. [20] in Baiyangdian Lake (China), and by Davies and Ugwumba [21] in Okpoka Creek (Nigeria). In addition, some fertilization experiments have shown shifts in relative dominance among algal species, with a shift towards cyanobacteria [2] and diatoms [22] under nutrient enrichment.The high species number in disturbed station reflects an enhanced nutrient enrichment, as reported in other studies with a high abundance of centric over pennate diatoms [23, 24]. The same observation was made in Lebanese coastal waters [25]. Indeed, Egge and Aksnes [26] and Fouillaron et al. [27] showed that diatoms always numerically dominated the phytoplankton community when silicic acid (DSi) concentrations were above 2 μM. Previous studies have shown that nutrient enrichment was associated with changes in species composition [1, 21, 28]. However, each microalgal group displayed a unique spatial pattern in response to N and P enrichment [2]. Our data showed that there is an increase of diatoms and cyanobacteria in disturbed station compared to control. For insistence,Anabaenasp. that was abundant in Skhira station (Table 3) is pollution indicator species usually associated with eutrophic water bodies [29] and has been recorded in Okpoka Creek (Niger) due to the increased influx of nutrients from anthropogenic inputs [21].Table 3
ANOVA on microalgae abundances between controls and disturbed stations. M: Mahdia (control station); K: Skhira (disturbed station); n.s.: not significant.
DF
MS
F
P
Tukey’s test
Total abundance
Station
1
781250
0.002
0.963
n.s.
Residuals
17
356395368
Species number
Station
1
338
6.953
0.017943
K > M
Residuals
17
48.61
H
′
Station
1
0.25515
1.552
0.230796
n.s.
Residuals
17
0.16442
Dinoflagellates
Station
1
504560556
5.662
0.029
K < M
Residuals
17
89109085
Gonyaulacales
Station
1
297273472
12.488
0.003
K < M
Residuals
17
23802884
Gymnodiniales
Station
1
6805.6
0.036
0.85
n.s.
Residuals
17
185040.8
Total
18
Dinophysiales
Station
1
80000
4
0.062
n.s.
Residuals
17
20000
Prorocentrales
Station
1
40951250
6.177
0.024
K < M
Residuals
17
6628897
Peridiniales
Station
1
1250
0.001
0.974
n.s.
Residuals
17
1148603
Diatoms
Station
1
4.81
E
+
08
16.398
>0.001
K > M
Residuals
17
2.93
E
+
07
Pennate diatoms
Station
1
48840139
0.924
0.35
n.s.
Residuals
17
52850433
Centric diatoms
Station
1
196020000
7.158
0.016
M < K
Residuals
17
27384118
Cyanobacteria
Station
1
533889
12.613
0.002
K > M
Residuals
17
42326Table 4
SIMPER results.
Species
Group K (disturbed station)
Group M (control station)
Av. Diss.
Diss./SD
Contrib.%
Cum.%
Av. abund.
Av. abund.
Coscinodiscus sp.
5566.67
0
15.47
1.14
17.93
17.93
Ostreopsis ovata
33.33
4566.67
13.93
3.89
16.14
34.07
Coolia monotis
0
3844.44
11.4
1.98
13.21
47.28
Prorocentrum lima
77.78
3666.67
10.88
2.11
12.61
59.88
Nitzshoides sp.
2416.67
0
5.35
0.46
6.19
66.08
Nitzschia sp.
933.33
122.22
2.58
1.26
2.99
69.07
Naviculasp.
1633.33
1444.44
2.48
1.38
2.88
71.95
Pleurosigma sp.
300
838.89
1.97
1.15
2.28
74.24
Peridinium sp.
133.33
538.89
1.49
1.29
1.72
77.72
Pinnularia sp.
55.56
483.33
1.38
1.61
1.6
79.32
Prorocentrum rathymum
477.78
0
1.28
1.34
1.48
80.8
Karenia selliformis
400
0
0.91
0.52
1.06
81.86
Polykrikossp.
255.56
338.89
0.9
1.45
1.05
82.91
Licmophora sp.
311.11
0
0.89
1.49
1.03
83.93
Chaetoceros sp.
244.44
0
0.75
0.62
0.87
84.8
Anabaena sp.
288.89
116.67
0.64
1.43
0.75
85.55
Protoperidiniumsp.
88.89
150
0.57
1.28
0.66
86.93
Gyrodinium sp.
111.11
61.11
0.54
0.56
0.62
87.56
Biddulphia sp.
166.67
0
0.51
0.58
0.59
88.15
Striatella unipunctata
200
0
0.49
0.88
0.57
88.72
Thalassiosira sp.
166.67
0
0.45
0.88
0.52
89.23
Leptocylindrus sp.
133.33
0
0.42
0.57
0.49
89.72
Scrippsiella sp.
133.33
0
0.41
0.71
0.48
90.2Table 5
CCA results.
Using environmental variables
Using nutrient ratios
Axis 1
Axis 2
Axis 1
Axis 2
Eigenvalues
0.493
0.154
0.153
0.104
Species-environment correlations
0.994
0.943
0.975
0.935
Cumulative percentage variance of species data
29.7
38.9
38.8
45.2
Cumulative percentage variance of species-environment relation
36.2
47.5
85.8
100.0
Sum of all eigenvalues
1.662
1.623
Sum of all canonical eigenvalues
1.361
0.734
Explained variance by CCA
81.88%
45.22%
Intraset correlations of variables with axes
NO2-
−0.135
−0.450
Si(OH)4
−0.955
−0.026
NO3-
0.281
−0.339
NH4+
−0.613
−0.445
PO43-
0.832
−0.060
NT
−0.588
−0.384
PT
0.145
0.210
T
−0.879
−0.072
Salinity
−0.386
0.143
pH
−0.821
0.033
SM
−0.497
−0.411
N : P
0.7204
0.6059
Si : N
0.7606
−0.3145
Si : P
0.8765
0.2950Phytoplankton species respond in different ways to nutrient enrichment, probably based on their life cycle characteristics, such as growth rate and their absorptive capacity [30]. According to Lepoint et al. [31], it is difficult to predict which group would be favored by the moderate intake of nutrients. For example, Lin et al. [32] found that nutrient enrichment (NO3-, NH4+, and PO4-, either alone or in combination) did not result in an increase in the abundance of phytoplankton. These communities respond in complex ways to the addition of nutrients. The seasonal influence on the response of microalgae to nutrient enrichment seems important [33].Our results showed that Si(OH)4 and nitrogen (NH4+, NT) are factors that contribute to differences of species abundance between stations (CCA; Figure 3). Indeed, silicon (Si) played the most important role in the growth and development of diatoms while dinoflagellates were mostly controlled by phosphorus (P) availability [34, this study]. When phosphorus loading increased, a shift from diatoms to dinoflagellates was observed [35]. Krumme and Liang [36] and Wear et al. [28] have found that diatoms were high in terms of abundance and species composition in nutrient enrichment station. Pennate and centric diatoms respond differently to disturbance. Indeed, our results showed that centric diatoms were higher in disturbed station, result found also by Harrison et al. [37]. Indeed, diatoms have the capacity of rapid cell division with growth rates (up to 5.9 day−1) generally much higher than those for dinoflagellates (up to 2.7 day −1) based on equivalent body mass [38].The canonical molar Redfield ratio of 16 : 1 (N : P) in the water column has traditionally been considered the determining factor for whether there is potential for N- or P-limitation in phytoplankton for a specific ecosystem [39, 40]. If the ratio of dissolved inorganic N (DIN) compared to dissolved inorganic P (DIP) is greater than 16 : 1 and nutrients are limiting, the system is considered P-limited [41]. Our data allow us to highlight that disturbed station was characterized by a slight P-limitation, which is known to induce the synthesis of alkaline phosphatases, AP, for numerous phytoplankton species such asAlexandrium catenella [42],Ptychodiscus brevis [43],Karenia mikimotoi [44],Gymnodinium catenatum,andAlexandrium tamarense [45]. AP is an enzyme used to convert DOP (dissolved organic P) to bioavailable DIP (dissolved inorganic P) by cleaving P from the DOP molecule [42, 46]. Some diatoms also do produce AP, includingPseudo-nitzschia,Chaetoceros, andSkeletonema [41, 47]. But previous studies [41, 47] have indicated that some phytoplankton groups (dinoflagellates and coccolithophores) may have a greater P-requirement and utilize DOP more readily than diatoms through the production of AP.We also found a strong Si-limitation in the control station which seems to be the main driver of differences in diatoms’ assemblages between control and disturbed conditions. This station showed low abundance of diatoms (Figure4, Table 4) compared to disturbed station. This result was also found by previous study [1, 21, 34, 35] that showed that increase in the N : Si ratio was proportional with the increase in flagellates and the decrease of diatoms abundance [1].Differences in species abundance and composition between the two stations cannot be assigned only to nutrients enrichment; many other factors could influence their distributions, such as water motion [1], light intensity [2], and pollution by heavy metals. For instance, in Skhira station (K), there is phosphate treatments factory that contaminates water and sediments by cadmium, mercury, arsenic, and so forth, as reported by Smaoui-Damak et al. [6].The phytoplankton-grazer interaction also plays an important role in controlling abundance and diversity of microalgae. Indeed, plankton is a food source for a range of grazers which influence their diversity and their abundance [48].The present study underlines the use of phytoplankton assemblages as general indicators of anthropogenic disturbance in marine ecosystems. It was clear that the effect of nutrient loading and turbidity could be depicted in the community structure and diversities of phytoplankton. During period of study (August), nutrient enrichment (mainly Si(OH)4 and NH4+) and turbidity were associated with increase of diatoms and cyanobacteria and decrease of dinoflagellates abundance, especially Gonyaulacales and Prorocentrales. The existence of Si- and P-limitation in study sites seems to be the main driver of differences in phytoplankton assemblages between control and disturbed conditions.



---

*Source: 241410-2014-06-25.xml*