Next Article in Journal
Study on the Horizontal Distribution Law of Flood Water and Sediment Factors under the Effect of Vegetation on a Curved Beach
Next Article in Special Issue
Zooplankton Community Dynamics in Temporary Mediterranean Wetlands: Which Drivers Are Controlling the Seasonal Species Replacement?
Previous Article in Journal
Assessing Water Withdrawals in Scarce-Data Transboundary Areas by Use of Dynamic Precipitation–Flow Relationships: The Case of the Hasbani River Basin
Previous Article in Special Issue
Epiphytic Diatom-Based Biomonitoring in Mediterranean Ponds: Traditional Microscopy versus Metabarcoding Approaches
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Promotion of the Development of Sentinel Species in the Water Column: Example Using Body Size and Fecundity of the Egg-Bearing Calanoid Copepod Eurytemora affinis

1
Université de Lille, CNRS, Université du Littoral Côte d’Opale, UMR 8187 LOG, Laboratoire d’Océanologie et de Géosciences, Station Marine de Wimereux, F-59000 Lille, France
2
Université Littoral Côte d’Opale, UMR 1158 BioEcoAgro, TERRA Viollette, USC Anses, INRAe, Université Lille, Université Artois, Université Picardie Jules Verne, Université Liège, Yncréa, F-62200 Boulogne-sur-Mer, France
*
Author to whom correspondence should be addressed.
Water 2021, 13(11), 1442; https://doi.org/10.3390/w13111442
Submission received: 31 March 2021 / Revised: 18 May 2021 / Accepted: 19 May 2021 / Published: 21 May 2021
(This article belongs to the Special Issue Advances in Mountain and Mediterranean Wetlands Conservation)

Abstract

:
The development of sentinel species in aquatic ecosystems is mostly based on benthic organisms; however, organisms living in water column such as zooplankton have received less attention, except for some cladocerans. In this paper, a new ecological indicator based on simple measurements of the size and fecundity of egg-bearing calanoid copepods is developed. The well-studied estuarine copepod Eurytemora affinis is used to illustrate this new framework. A large database obtained from laboratory experiments developed under different conditions is used to define a reference regression line between clutch size (CS) and prosome length (PL). The same database allowed one to confirm that the coefficient of variation (CV) of CS is an adequate estimator of the accumulated stress at population level. The CV of PL shows very little variability in all experimental and field conditions. The values of CS and PL obtained from the Seine, Loire, Gironde, Scheldt and Elbe estuaries in Europe are compared to the reference regression line. A quality index (QI) is calculated as a percentage of difference between the observed and the predicted CS. The QI classified 19 samples collected in the Seine estuary between 2004 and 2010 into four classes according to the physiological condition of the copepod female. A single sampling from June 2004 (5.26%) showed a very good condition, whereas 57.9% of the sampling dates confirmed good conditions. On the other hand, four sampling dates were associated to very bad conditions and three sampling dates indicated bad conditions. Seven additional samples obtained from other European estuaries between 2006 and 2009 were also used. Females showing poor conditions were observed in the early spring of 2005 and 2008 as well as during the month of November. These years were characterized by very strong climatic anomalies with a very cold late winter in 2005 and a warm winter in 2008. Therefore, it seems that the QI perfectly reflected the strong stress caused by the sudden change in hydro-climatic conditions that have certainly affected the physiology of copepod females and probably the availability of food. The new indicator is very simple to calculate and can be generalised to several aquatic ecosystems (fresh water and brackish water) by targeting the dominating egg-bearing calanoid copepods. As in the case of E. affinis, the development of sentinel species based on copepods or cladocerans can enrich ecological and ecotoxicological studies given their capacity to integrate the variability of their habitats’ quality at the individual and population levels.

1. Introduction

The majority of aquatic ecosystems are facing multiple anthropogenic pressures (e.g., eutrophication, pollution, habitats fragmentation, freshwater demand, riverine and maritime transportation, etc.) combined to climatic fluctuations (e.g., global warming, extreme events, drought seasons, storms, etc.) [1]. Consequently, it is important to develop a link between ecological studies and future management strategies to better understand the functioning of most aquatic ecosystems, evaluate their global services and to adopt the most appropriate management and conservation strategies [2]. To reach this goal, a comprehensive monitoring strategy of the whole ecosystem is required, but unfortunately, this was never achieved for most aquatic ecosystems. The studies comparing different aquatic systems and using adequate modelling tools are scarce [3]. Most often, the heterogeneity of the measured variables as well as the absence of an adequate modelling framework capable of linking individual, population and community levels, limit the development of a holistic approach. Some global initiatives, such as the European Water Framework Directive (WFD; 2000/60/EC) aiming to protect the quality of the water by using an integrated river basin management approach, tried to promote biological indicators at large scale. Although the WFD is important in enhancing a global debate on the ways of developing and harmonizing ecological indicators in the context of water quality, some key components of all aquatic ecosystems, such as zooplankton, were missing [4]. Furthermore, the extension of the WFD to the transition zones, such as estuaries, was extremely limited. The limitations of the WFD for the transition zones including estuaries and the coastal zones were partially completed by the European strategy for marine environment directive (2008/56/EC). In that directive, the zooplankton is retained among the other biological groups (phytoplankton, benthos, fish). From a conceptual point of view, zooplankton is a suitable group to develop different ecological indicators at different levels of organisation (individual, population and community). In fact, zooplankton is a very diverse group of small-sized organisms, with a fast growth rate and short generation length. All these characteristics allowed zooplankton to be a good indicator of natural and/or anthropogenic variabilities [4,5]. Particularly, in the freshwater habitats, cladocerans and mainly Daphnia have been intensively used as a good ecological indicator [6]. In addition to the advantages of studying cladocerans, other crustacean groups of zooplankton, such as copepods, are also present in all aquatic ecosystems and their study at community and/or population levels can also be useful to enrich the battery of ecological indicators.
In this paper, we focus on freshwater and brackish water habitats at a basin level and their extension to the estuarine systems. Although there is a low number of ecological studies targeting zooplankton communities from the river source to the estuarine habitat, Le Coz et al. [7] tried to use different ecological concepts originated from either riverine or estuarine ecology to the continuum of the Scheldt river and estuary. Zooplankton communities are also very good indicators to characterize the trophic status of Mediterranean wetlands and to suggest some conservation strategies [8]. The ecological status of some large lakes can be assessed by using several ratios, such as cladoceran/calanoid and cyclopoid/calanoid [9]. The ratios can be based on the size of species; for example, the ratio between large and small rotifers in Lake Geneva can be simultaneously affected by eutrophication and global warming [10]. The body size of zooplankton groups was more useful than taxonomic resolution to rapidly assess the effect of environmental factors on temporal wetlands [11]. Furthermore, the body size of the Daphnia cucullata was suggested to be a good ecological indicator of the trophic conditions of lakes [6]. Consequently, the body size of zooplankton species seems to be an important trait to be considered when developing ecological indicators.
However, most freshwater ecosystems experience different sources of anthropogenic pressure with different intensities of climatic variability that challenge the development of ecological indicators. Whereas the studies of zooplankton communities are relatively very common [7,8], ecological indicators can be conceived at either individual or population levels. The link between individual and population levels can be achieved by studying the variability of the life traits of sentinel species [2]. For example, the use of condition indices based on measuring size (or morphometry) and weight of individuals was used in studies targeting fish juveniles [12] or benthic invertebrates [13]. The allometric relationship between weight and size of individuals is also used as a good ecological indicator of the health status of intertidal mudflats [14]. The same study explored the reproductive status by dissecting mature individuals and observing oocyte development. Although reproductive traits are important, this approach needs specific skills and cannot be easily applied by nonexperts.
Nevertheless, the importance of considering life traits and size of individuals seems to be a good option to enrich the development of ecological indicators, but such an approach received less attention in zooplankton species. In fact, the use of body size alone in a single species of a calanoid copepod [15] or caladoceran [6] allowed one to compare different water bodies with contrasting ecological statuses. However, these studies did not combine body size with other life traits and particularly the fecundity of females. We hypothesise that the simultaneous consideration of body size and fecundity of some zooplankters can enrich the potential of developing new ecological indicators.
In this paper, we target calanoid copepods and especially egg-bearing species to suggest a general framework that can be applied to almost all freshwater and brackish water ecosystems. In order to illustrate this approach, we select the widely distributed species Eurytemora affinis and focus on the well-studied population [16,17,18] from the low salinity zone of the Seine estuary. In several freshwater and estuarine ecosystems, the diversity of calanoid copepods is very low and most often a single species is dominating even the whole zooplankton community [19,20]. Moreover, when adult females carry eggs, the clutch size can be estimated in addition to other morphometric parameters [15]. However, the interpretation of these data remains complicated given their high variability. The first objective of this paper is to demonstrate the way that the use of this relationship between fecundity and female size can help to build ecological indicators. The second objective is to show that inter-individual variability that can be measured simply by the coefficient of variation (CV) provides a reliable indicator of the magnitude of stress experienced by the population.

2. Materials and Methods

2.1. The Selected Species

In this study, we used the calanoid copepod Eurytemora affinis, an euryhaline and eurytherm species, that has a large geographical distribution and can be found in the low salinity zone of most estuarine ecosystems in the northern Atlantic as well as several lakes [21]. In Europe, three populations (Baltic, English Channel–North Sea and Atlantic) were identified using a phylogenetic approach [22,23]. Other studies in the Baltic confirmed the co-existence of the native E. affinis with other close species recently described as Eurytemora carolleeae [24]. Although there are close morphological similarities between these two species, we only focus on the European populations of E. affinis and particularly the one from the Seine estuary. In fact, the life cycle properties of this population were studied in the laboratory under controlled natural conditions [25,26,27] or in the presence of different classes of pollutants [28,29,30,31,32] as well as in the field [16,17,18]. All these studies allowed for the development of individual-based models [33,34]. Very recently, the whole genome of E. affinis from the Seine estuary maintained in laboratory culture was sequenced and compared with other available genomes of three copepods, including E. carolleeae [35]. The complete barcoding of the species using the mitochondrial genome [36] as well as the annotation of the transcriptome [37] were realised. This gives a great background on the selected species to test the hypotheses of this study and promote the development of sentinel species.

2.2. Experimental Data

Eurytemora affinis was sampled in the low salinity zone of the Seine estuary and transferred to the laboratory. Then, after manually sorting a large number of late developmental stages, the culture was initiated in a large 25 L aquarium at a temperature of 15 °C and salinity of 15. After two generations of acclimation (~1 month), the batch culture was transferred to 20 °C and a salinity of 15 and maintained during 2 additional generations (~20 days) for acclimation. Copepods were fed every two days using Rhodomonas marina as described in [27]. Then, the experiment started by allowing 40 ovigerous females to lay their nauplii for the first generation (F0) at a temperature of 20 °C and salinity of 15 (T20S15). Then, adult females were removed in order to avoid the formation of a new egg-sac. The produced nauplii were followed in a 2 L beaker until adulthood. The water was changed after one week when copepodite stages appeared. After approximatively two weeks of development, a high number of ovigerous females appeared in the beaker, allowing one to stop running the generation experiment. Therefore, 40 live ovigerous females were randomly sorted and used to initiate the following generation F1. At the same time, 40 additional ovigerous females were sorted and used in another beaker containing a freshly prepared culture medium at the same temperature of 20 °C but a higher salinity of 25 (T20S25). The rest of the population was fixed in buffered formalin with final concentration 4% for subsequent morphometric analyses and clutch size counting. Then, the same procedure was repeated for each generation and for both beakers independently. After 18 generations at T20S15 and 17 generations at T20S25, the follow up of the multigeneration protocol was continued in a new incubator with a higher temperature of 24 °C, leading to two final conditions: T24S15 and T24S25. We should notice that in each generation, between 20 and 40 ovigerous females were sorted randomly to measure prosome length (PL) and estimate the clutch size (CS).
Although in the previous study, the first 11 generations were not used in the analysis (see Figure 1 in [27]), they were included in this analysis to obtain more data on the size and fecundity of E. affinis. The objective of using this very long-term experiment is to obtain a reference data on 4 different experimental conditions (T20S15, T20S25, T24S15 and T24S25) to test our hypotheses under controlled conditions. For each generation, the mean values of prosome length and clutch size were calculated. The variability in PL and CS between females in each generation was assessed by computing the coefficient of variation (CV = 100 × SD/mean), where SD is the standard deviation. The CV can be used as an additional indicator of stress within each generation.

2.3. Field Data

E. affinis samples collected between June 2004 and October 2010 in the same tidal conditions (during the ebb) near the Normandy bridge (mouth of the Seine estuary, as in [16]) were combined in this study to develop the ecological indicator. The dates, numbers of ovigerous females, mean values and CV(%) of PL and CS are provided.

2.4. Ecological Indicator Based on Prosome Length–Clutch Size Relationship

The experimental data allowed one to define a reference regression line that is supposed to represent good conditions that copepods can encounter in the field.
CSpred = a × PLobsb
where CSpred is the predicted CS and PLobs is the observed prosome length.
Then, a simple ecological indicator named the quality index (QI) was computed for all data obtained from the field. This QI corresponds to the percentage of deviation of the observed clutch size compared to the predicted clutch size.
QI = 100 × (CSobs − CSpred)/CSpred
where CSobs is the observed mean clutch size.
To illustrate the principle of this method, 4 classes were suggested (Table 1):
All data obtained from the Seine estuary (19 dates) including the year 2005 (7 dates) characterised by a late winter/early spring unfavourable conditions to E. affinis development [17] are used to compute this index. Furthermore, additional samples collected in the Loire estuary near Paimboeuf (the same site as in the study of [22]) between 2006 and 2009 (one sampling per year), as well as three other samples obtained between March and April 2006 (one sample per estuary) in the Gironde, Scheldt and Elbe estuaries were added as a supplementary data.

2.5. Data Analysis

The mean values of prosome length (PL) and clutch size (CS) obtained from each generation were gathered to perform a linear regression (CS = a × PL − b). In order to assess the robustness of this model, the input data (CS and PL) were resampled using a bootstrap technique available on Matlab Softaware (Mathworks, Inc., Portola Valley, CA, USA). This allows one to obtain 5000 resampled values of the slope (a), correlation coefficient (r) and the intercept (b/a corresponding to CS = 0).
The values of CV(%) for prosome length and clutch size were aggregated in the following categories (Table 2):
Then, the nonparametric Wilcoxon rank sum test was applied to compare the three paired data samples (In situ good vs. In situ bad), (T20S15 vs. T24S15) and (T20S25 vs. T24S25).

3. Results

Figure 1A shows that the different generations from all experimental conditions follow a clear linear trend between clutch size and prosome length with r2 = 0.83. Excepting two synchronous generations from T20S15 and T25S25, that produced significantly smaller females and clutches due to a unique episode of low food quality (i.e., Rhodomonas marina used not at its exponential phase of growth), the increase in temperature reduced both prosome length and clutch size. Figure 1B shows the distribution of the coefficient of correlation (r) with a mean value of r = 0.91. All reconstructed r values by the bootstrapping technique are highly significant (p < 0.001), confirming the robustness of the linear relationship between CS and PL. Figure 1C,D shows the distributions of the values of the slope and intercept, respectively. The mean value of the intercept is around 730 µm, confirming that logically, egg-carrying females of E. affinis from the Seine population cannot be smaller than that size.
The data obtained from the Seine estuary and calculated in the last column of Table 3 are shown in Figure 2 (circle symbols). In addition, seven data obtained from other European estuaries (Loire, n = 4; Gironde, n = 1; Scheldt, n = 1 and Elbe, n = 1) are also used. We notice that the highest value of the QI was observed in June 2004 (24.85%), suggesting very good conditions (Figure 2). Additionally, the largest females with good conditions were observed in February 2009 in the Seine estuary (PL = 1127.42 µm, Table 3). This contrasted with the smallest females observed in the Seine estuary in November 2007 (PL = 846.29 µm). The size differences between these extreme situations is around 33.22% but the differences in clutch size is much higher (458.76%). In the Seine estuary, 11 sampling dates corresponded to good conditions. The QI index classified 63.16% of the dates sampled in the Seine estuary in good or very good conditions. The other sampling dates belong to bad (n = 3, 15.79%) and very bad (n = 4, 21.05%) conditions. Figure 2 shows that excepting one sampling date of November 2007, with very bad conditions, all other sampling dates with low conditions are observed in the years 2005 and 2008. In particular, the critical months are March, April, November and December (Figure 2). We notice that among the seven samples from the other European estuaries, only the one obtained in Loire in April 2008 showed bad condition. It seems that the early spring and late autumn of the years 2005 and 2008 corresponded to unfavourable conditions for E. affinis development. It is also interesting to see that the seasonal pattern of the QI was similar in both years 2005 and 2008. In 2005, the very bad conditions were observed in March and April, followed by four good conditions in May, June, July and September, and ended by a bad condition in December (late November too). In 2008, a very bad condition was observed in March followed by a bad condition in April and two consecutive good conditions in June and July, before again experiencing a bad condition in November.
Figure 3A shows that the variability of measures of prosome length (PL) are very low (CV < 5%) with slightly higher values in the field samples compared to the experimental conditions. This confirms that PL is a very robust morphometric measurement that can be used as an independent variable in morphological or trait analyses. Although the CV of PL in bad conditions (in situ), and in the most stressed condition T24S25, are slightly higher than the less stressed conditions, the observed differences are not statistically significant. On the contrary, the values of CV of CS are very high and particularly in field conditions. Figure 3B shows that the variability of CS is on average higher in the field (36.12%) than in the laboratory conditions (20.35%, average of all conditions). In field conditions, the CV of CS when QI indicates good or very good conditions is equal to 27.52%, whereas the CV when QI indicated bad or very bad conditions is equal to 43.05% (Figure 3B). Although this is a big difference in CV between favourable and unfavourable field conditions, it is not statistically significant at p = 0.05, and this is most likely due to the low number of samples in our study. Under controlled conditions, the increase in water temperature from 20 °C to 24 °C at each salinity (15 and 25) is associated with a statistically significant increase in the CV of CS (Figure 3B). The value of CV in the condition T24S25 with a double thermal and saline stresses (30.66%) is higher than the value of CV obtained in only thermal stress condition T24S15 (25.86%).

4. Discussion

4.1. Assessing the Environmental Quality Using Morphometric Measurements

In aquatic ecology, the individual growth parameters, such as body mass and size are often used to assess the condition of benthic organisms [12,14]. More generally, measuring several morphometric traits can be useful in comparing different sites and populations. This approach was applied to the cladoceran Daphnia cucullata [6] and the calanoid copepod Arctodiaptomus salinus [15]. In addition to morphometric measurements that can be easily achieved, some reproductive traits (e.g., sexual maturity, fecundity) can be good indicators of water quality and ecosystem health in general. In fact, the reproductive traits of benthic organisms (i.e., sexual maturity and fecundity) are good physiological indicators, very useful to compare between different sites with contrasted anthropogenic pressure such as pollution [13]. More generally, most sentinel species and the use of biometric measurements are developed for benthic organisms and fishes. On the contrary, the organisms living in the water column such as zooplankters received less attention in the framework of developing sentinel species.
One difficulty in promoting sentinel species based on crustacean zooplankters is certainly related to the diversity of this group, favouring the use of the community level to assess the quality of habitats [8]. Even if several species are models in aquatic ecotoxicology [38], most toxicity studies realised in the laboratory under controlled conditions cannot be easily extrapolated to in situ conditions. Particularly, cladocerans and copepods are fast-growing organisms which allow one to study the acute toxicity of several pollutants [28,30], but also the intergenerational effects of natural factors [27] or pollutants [32,39]. In field conditions, several factors can simultaneously affect body size and life traits. In freshwater ecosystems, other benthic crustaceans such as the Gammarus fossarum are used as sentinel species. This status has been favoured by an intensive research program aiming to better understand the sensitivity of the species to different environmental conditions while maintaining a balance between in vivo and in situ studies [40]. In addition to Daphnia, this study confirmed that some species of calanoid copepods can be promoted as good candidates for developing sentinel species of freshwater and estuarine ecosystems. Indeed, the diversity within this group is often low, unlike cyclopoids, which favours the work at the level of a single population [16,19,20]. Moreover, the majority of these species carry their eggs, which offers a considerable advantage to easily measure morphological traits (e.g., body size) and fecundity (e.g., clutch size). It is sufficient to focus on ovigerous females to easily and operationally make reliable measurements of prosome length and estimate fecundity by counting eggs in the sac [27]. This will offer a great advantage to assess the fecundity of copepods compared to other benthic organisms that require dissection and histological techniques [14]. For copepods, several ecological studies conducted in estuarine environments [41] and in lakes [42] already measured female prosome length and clutch size. However, the potential of using an allometric relationship between female size and clutch size has not been yet explored to develop indices of individual condition. Nonetheless, an integral indicator based on some morphometric measurements applied to the calanoid copepod Arctodiaptomus salinus was developed and used to assess the population state at larger scale [15]. In our study, the potential of using simple allometric relationships between size and fecundity to identify the individual condition was demonstrated by using the well-studied estuarine European species E. affinis.

4.2. Promoting Certain Calanoid Copepods as Good Sentinels: Case of the Species Eurytemora affinis

The copepod E. affinis has a large geographical distribution [22] and can dominate zooplankton in the low salinity zone of several estuaries [16]. However, given the geographical discontinuity of the aquatic habitats of this cryptic species, its strong genetic and morphological heterogeneity [43], the description of a new species E. carolleeae [44] which is very close morphologically to E. affinis, requires great caution when using data from different in situ populations. In this study, we focused on the population of E. affinis from several European estuaries. For example, in the Seine estuary the effects of natural factors, such a temperature and salinity on the dynamics of this population are well described at tidal [16], annual [17] and inter-annual [18] scales. Generally, the populations of E. affinis are distributed in the low salinity zones of estuaries near the maximum turbidity zone and can bioaccumulate different pollutants [45,46]. Due to the complexity of aquatic ecosystems, the separation between natural factors and pollution remains difficult, although this challenge is crucial in order to better understand the possible combined effects of global warming and pollution on aquatic ecosystems [1]. Parallel to in situ studies of E. affinis, an intensive experimental work has been developed allowing the switch from the acquisition of demographic or ecotoxicological parameters using individual observations [25,26,30] to the use of multi-generation observations at the population level [27,32]. The use of E. affinis cultured in the laboratory offers great potential to obtain individuals raised in a standardized way to conduct interdisciplinary research such as the effect of turbulence [47] or pollution [29,48] on swimming behaviour. Furthermore, the kinetics of interaction between heavy metals and E. affinis exposed in large volumes confirmed the hypothesis of competition between ions of cadmium in presence of nickel and copper [31]. Recently, the same population of E. affinis originated from the Seine estuary was promoted as a good model for molecular ecotoxicological studies using genomic approaches [36,37].

4.3. Development of a Quality Index (QI) Based on the Relationship Body Size-Fecundity

Due to the large plasticity of copepods, and particularly E. affinis living in a highly variable environment, it is difficult to obtain accurate reference values for biomarkers used in ecotoxicology and life traits (i.e., fecundity) measured in the field. The use of experimentation is important to understand the effects of natural parameters such as salinity and temperature on enzymatic activity [49], as well as body size and fecundity changes [27]. In field conditions, a high variability of the prosome length and clutch size of ovigerous females of E. affinis was observed. The species is known to have a high plasticity, mainly regarding size and fecundity [27]. In addition, the age of females that may affect the fecundity [25] can be heterogeneous in the field. Consequently, it might be difficult to build a reference relationship between female prosome length and clutch size when all data are gathered together. An average regression line may not be the best solution to obtain the reference physiological relationship between body size and fecundity. As in our study, a database gathering more than 40 generations was available, a first reference line from lab conditions was suggested. These experiments that aimed to test two sources of stress related to an experimentally simulated global warming of +4 °C and salinity increase helped to validate the importance of using the coefficient of variation (CV) as an indicator of stress.
Contrary to the average relationship between CS and PL published in [27], based only on four data points, in this study, all data obtained from each single generation were used. Moreover, the early generations of the acclimation phase at a temperature of 20 °C were also added. This provided a very robust linear regression that defined the reference line to be compared to data from field conditions. In order to illustrate this approach, a simple quality index (QI), based on the normalised difference between the observed and predicted values of clutch size, was calculated. Then, the QI values were classified arbitrarily into four classes (Table 3). We assumed that our experimental conditions based on an optimised feeding protocol [27,32] reflected the good conditions in the field. Generally, the best mono-algal diet for copepods is composed by Rhodomonas sp. ([50], our personal observations). This confirms that E. affinis in our long-term experiment (>1 year) received good diet excepting during one generation that produced smaller females with a lower clutch size (see Figure 1A). However, this low-quality food episode significantly reduced PL and CS but did not affect the general trend of the CS–PL relationship.
In the field conditions, the lowest values of the QI were observed in March 2005, April 2005 and March 2008. Devreker et al. [17] studied the annual variability of E. affinis population in the Seine estuary and all environmental parameters during 2005. They showed that the negative anomaly of temperature that was observed in late winter was responsible for the very bad conditions of the population of E. affinis. This extreme climatic event observed at a larger scale in southern Europe [51] negatively affected the E. affinis population in the Seine estuary [18], as well as the seasonal succession of phytoplankton communities in the English Channel [52]. The QI clearly identified the critical period of poor conditions in March and April 2005, and the good conditions observed from May to September 2005. In fact, Devreker et al. [18] observed a high mortality in the beginning of the study (March–April 2005) that did not allow one to obtain a high density of the population of E. affinis during May and June, as it was expected from previous studies [17]. During 2005, the surviving individuals developed to adult stages in May and encountered favourable conditions. Furthermore, the QI confirmed that all ovigerous females observed between May and September were in good conditions. Unfortunately, no sampling was realised in October and November 2005, but it seems that the females sampled on 1 December 2005 showed bad conditions, but nevertheless with large females compared to those observed during autumn (October–November) in the years 2006, 2007 and 2008 (see Figure 2). This confirms that the annual cycle of 2005 was extremely disturbed at both population [18] and individual (this study) scales. It also seems that November is a critical month for E. affinis development in the Seine estuary. The year 2008 appeared to be unfavourable for the spring development of E. affinis in the Seine estuary (March–April) and also in the Loire estuary (April). According to the earth observatory of the NASA, the year 2008 was much warmer in Europe, Russia and the Arctic compared to the Pacific Ocean, which was cooler than usual (https://earthobservatory.nasa.gov/images/36699/2008-global-temperature, accessed on 15 March 2021). This climatic situation, with an unusual warming of winter in Europe observed in 2008, could also disturb the regular development of E. affinis population, due to the direct effect on metabolism and/or the indirect effect on food quality. However, this hypothesis suggested by the new QI developed in this study should be explored deeply in the future. Even if the number of data was limited in years 2007 and 2008, it seems that the condition of E. affinis using the QI was bad or very bad during spring and autumn (see Figure 2). Both years were characterised by a significant warming during winter. However, other factors related to food quality, as well as other anthropogenic factors, could be involved in such pattern.
The river discharge is an easily accessible global parameter that plays a key role in the dynamics of the estuary [17]. Some studies already confirmed that higher river discharges usually enhance the lower trophic level production in the estuarine environment such as Eurytemora species [53] or the benthic sentinel Hediste (=Nereis) diversicolor [54]. During the period 2004–2010, the hydrological regime in the Seine river was characterised by low annual discharges (419.97 m3/s; SD = 225.85) compared to the previous decade 1994–2003 (mean = 596.97, m3/s; SD = 424.01). This situation of low river discharges could exacerbate the observed effects of extreme climatic events encountered during this period. However, these hypotheses on the link between hydro-climatic situations and the production of the lower trophic components need to be confirmed in the future with more dedicated experimental and/or sampling designs.

4.4. Generalisation of the QI to Other Freshwater and Estuarine Ecosystems

Several studies in marine ecology tried to establish a relationship between the reproductive output of copepods, generally expressed in terms of carbon content (e.g., for a clutch of eggs), and the female body weight. By aggregating several published data from different species, a global trend between fecundity and female body weight was observed [55]. This power law relationship (linear in log–log plot) could be useful in biogeochemical modelling but it is not accurate when a single species is considered. However, a better taxonomic resolution, particularly at species level, is required to increase the ecological relevance of the study and provide species-specific relationship.
In this study, the European populations of E. affinis, and particularly the population from the Seine estuary, were used to develop our approach. The additional data points from the Elbe, Scheldt, Loire and Gironde estuaries perfectly fitted to the regression line, excepting a single data point of April 2008 showing a bad condition (see Figure 2). This means that the approach suggested in this study can be applied to all E. affinis populations at least at European level, although phylogenetically, three lineages were described [22]. However, for other species, even those very close genetically, a new calibration of the reference line, that could be obtained from in situ data, is necessary. For example, Lloyd et al. [56] used field data obtained in the Chesapeake Bay (USA) in 2002 and 2003 and fitted a good relationship between CS and PL of the dominating copepod which was named before E. affinis but whose name was corrected later to E. carolleeae [57]. This robust relationship (r2 = 0.85) obtained during two favourable years, 2002 and 2003, could be used as a reference line for the species E. carolleeae. By assuming that 2006 was also a favourable year for E affinis development in Europe, we performed a simple regression line using five estuaries sampled in the same year. This provided an excellent relationship, explaining 99.35% of the total variability (PL = 0.1931 × CS − 139.32, r2 = 0.9935), with an intercept value of 721.49 mm. This later value fell within the range of possible values generated by the bootstrap technique (see Figure 1D). Therefore, this relationship obtained from the field data for a single season in 2006, was not statistically different from the relationship established in the lab and can also be used as reference line. The linear regression line obtained in [56] for E. acrolleeae in the Chesapeake Bay (CS = 0.229 × PL − 126) was statistically different from the relationship obtained in the present study (CS = 0.2 × PL − 145.3). The value of the slope (0.229) projected in the histogram showing the distribution of the slope values produced by bootstrap (Figure 1C) was significantly higher than the mean value of slope (0.2) obtained in this study (p = 0.0244). The intercept value of E. carolleeae was 550.22 µm, which was much smaller than the lowest value generated by bootstrap in this study 677.07 µm (p = 0). It seems that the conditions of E. carolleeae in the study of [56] were good, because the in situ data did not show a huge dispersion, allowing one to suggest a linear decreasing relationship between PL and temperature with r2 = 0.59. Pierson et al. [57] measured the PL and carbon contents of females and males of E. carolleeae during winter in 2007 and 2008. They compared their results with the predicted value of PL vs. temperature of [56], and confirmed that their values obtained in 2007 were slightly above the regression line, whereas the values of 2008 were below the regression line (see Figure 4b in [56]). This means that females of E. carolleeae in winter 2008 indicated bad conditions. Unfortunately, such observations during winter in 2007 and 2008 were not available in the studied European estuaries. However, the available samples realised in March and April 2008 in the Seine and Loire indicated bad conditions in April (both estuaries) and very bad conditions in March (Seine estuary). This means that a meta-analysis can be easily conducted by using the QI in several aquatic ecosystems even if the key sentinel species are different. The inter-comparison between different sites and ecosystems could be conducted in the future with much larger database. We can speculate that large-scale climatic events (e.g., global warming, extreme climatic anomalies) can lead to synchronous observations, such as the examples observed in early 2008 in three transatlantic estuaries (the Seine, Loire and Chesapeake Bay). However, a larger database has to be built in the future to accurately test such hypothesis.
The reference lines obtained from field observations for E. affinis in Europe during spring 2006 (this study) and in the Chesapeake Bay during 2002 and 2003 [56] showed very low dispersion of the data, which contrasted with the regular situation when field data are used. For a scatter plot with a high dispersion of the data, it is still possible to identify the reference line by modelling the maximum clutch size (e.g., 95th percentile) instead of fitting the average (or median) response. This approach was applied in Daphnia [58], where a large dataset allowed one to use 50µm width of size classes and adjust the maximum of brood size (see Figure 1 in [58]).
The QI suggested in this paper can easily be generalised to other freshwater and brackish water ecosystems. In several aquatic ecosystems, the calanoid copepods, and sometimes the whole zooplankton community, are dominated by a single or very low number of species. For example, the egg-bearing calanoid copepod Eudiaptomus gracilis is present along the year in several European lakes [20,59]. Gilbert et al. [60] studied the crustacean zooplankton in 36 Mediterranean wetlands in the south of Spain and confirmed that the number of calanoid copepods in each wetland is very low (1 or 2). For that reason, Jiménez-Melero et al. [61] targeted the dominant egg-bearing copepod Arctodiaptomus salinus to study its population dynamics and secondary production in a shallow saline pond. Moreover, the morphometric measurements based on the same species were used to compare 12 different populations [15]. In order to apply the QI, it is important to target a single species present almost along the year and avoid a mixture of species that may have different morphologies and reproductive strategies. This will allow one to have different reference regression lines between CS and PL, as we confirmed above for the very close species E. affinis from Europe and E. carolleeae from the Chesapeake Bay. When a single species is used, several sites can be compared. For example, Chow-Fraser and Maly [62] studied the clutch size variability of two Diaptomus species in 11 lakes. They found that the CS is positively correlated with chlorophyll a and negatively correlated with temperature, which is a very common pattern observed in other marine studies [55]. For the species Diaptomus oregonensis, they showed a positive correlation between log-transformed data of both CS and PL, with a high dispersion of the data (see Figure 2a in [62]). Without establishing a reference line of regression between CS and PL, it is difficult to explain the sources of this dispersion and variability in CS within and between lakes.
The QI is calculated by using only the mean values of CS and PL. In order to evaluate the variability of the mean values of CS and PL, the coefficient of variation (CV) was also estimated. The inverse of the CV (1/CV) was used at the community level to estimate ecosystem stability [3]. In this study, the CV can be considered as a standardized estimator of the individual variability between females, and could be an indicator of stress intensity. In fact, in optimal and non-stressful conditions, the individual variability in most life traits was often lower than the one observed during stressful conditions (personal observations). This study confirmed that in controlled conditions, the simulated warming (+4 °C) induced a stressful response by significantly increasing the CV at two different salinities (see Figure 3). In field conditions, the CV was generally higher than that observed in laboratory conditions, and particularly in the sampling dates associated to bad or very bad conditions according the QI classification. This means that the variability associated to the mean values should be considered when developing simple ecological indicators such as the QI. However, this conclusion remains valid for all indices of conditions often based on average values of biometric measurements. The variability around the average values should not be neglected for a better ecological interpretation.
This study mainly focussed on adult females of egg-bearing copepods. For free spawning copepods and adult males, only morphometric measurements are possible. In that context, it is possible to use the CV of different morphometric parameters to develop indicators of variability [15]. In fact, stochastic variation is an important component of phenotypic variance, at the same level as genotypic variation and phenotypic plasticity, and usually increases under stress [63]. Without using the fluctuating asymmetry (FA) of bilateral traits, that may need very accurate morphometric data [63], the most common traits such as prosome length, prosome width, and urosome length can be easily measured and used in copepod [15] or in cladoceran studies [6]. This approach can complete the QI developed in this study. In addition, for free spawning copepods in freshwater ecosystems, it is possible to obtain a relationship between daily egg production rate (equivalent to CS) and PL by using the incubation techniques commonly deployed in marine systems [64].

Author Contributions

S.S. and A.S. conceived the study and wrote the MS. A.S. completed all morphological measurements and the associated database. S.S. performed statistical analyses. Both authors have read and agreed to the published version of the manuscript.

Funding

This research promoting the copepod Eurytemora affinis as a sentinel species received funding from the GIP Seine-Aval under several research projects: ZOOSEINE, ZOOGLOBAL and SENTINELLES. This work is a contribution to the project CPER 2014–2020 MARCO (funded by Europe FEDER, French government the region Hauts-de-France and IFREMER) and the International Associated Laboratory between Université de Lille and National Taiwan Ocean University (IAL MULTIFAQUA). This work benefitted from the French GDR “Aquatic Ecotoxicology” framework which aims at fostering stimulating scientific discussions and collaborations for more integrative approaches.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to some rules of the funding agencies.

Acknowledgments

We thank all past and present members of S. Souissi’s group who helped in keeping micro-algae and copepod cultures for several years and particularly those who helped in initiating the multigeneration experiment. We are grateful to all team members and the research vessel crew that contributed to all sampling campaigns in the Seine and Loire estuaries. We are grateful to Gesche Winkler for providing the samples collected in Spring 2006 in the Elbe, Scheldt, Seine, Loire and Gironde.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Rodríguez-Romero, A.; Viguri, J.R.; Calosi, P. Acquiring an evolutionary perspective in marine ecotoxicology to tackle emerging concerns in a rapidly changing ocean. Sci. Total Environ. 2021, 764, 142816. [Google Scholar] [CrossRef]
  2. Tlili, S.; Mouneyrac, C. New challenges of marine ecotoxicology in a global change context. Mar. Pollut. Bull. 2021, 166, 112242. [Google Scholar] [CrossRef]
  3. Chang, C.W.; Ye, H.; Miki, T.; Deyle, E.; Souissi, S.; Anneville, O.; Adrian, R.; Chiang, Y.R.; Ichise, S.; Kumagai, M.; et al. Long-term warming destabilizes aquatic ecosystems through weakening biodiversity-mediated causal networks. Glob. Chang. Biol. 2020, 26, 6413–6423. [Google Scholar] [CrossRef]
  4. Jeppesen, E.; Nõges, P.; Davidson, T.A.; Haberman, J.; Nõges, T.; Blank, K.; Lauridsen, T.L.; Søndergaard, M.; Sayer, C.; Laugaste, R.; et al. Zooplankton as indicators in lakes: A scientific-based plea for including zooplankton in the ecological quality assessment of lakes according to the European Water Framework Directive (WFD). Hydrobiologia 2011, 676, 279–297. [Google Scholar] [CrossRef]
  5. Beaugrand, G.; Brander, K.M.; Lindley, J.A.; Souissi, S.; Reid, P.C. Plankton effect on cod recruitment in the North Sea. Nature 2003, 426, 661–664. [Google Scholar] [CrossRef]
  6. Karpowicz, M.; Sługocki, L.; Kozłowska, J.; Ochocka, A.; López, C. Body size of Daphnia cucullata as an indicator of the ecological status of temperate lakes. Ecol. Indic. 2020, 117, 106585. [Google Scholar] [CrossRef]
  7. Le Coz, M.; Chambord, S.; Meire, P.; Maris, T.; Azémar, F.; Ovaert, J.; Buffan-Dubau, E.; Kromkamp, J.C.; Sossou, A.C.; Prygiel, J.; et al. Test of some ecological concepts on the longitudinal distribution of zooplankton along a lowland water course. Hydrobiologia 2017, 802, 175–198. [Google Scholar] [CrossRef]
  8. Gilbert, J.D.; de Vicente, I.; Jiménez-Melero, R.; Parra, G.; Guerrero, F. Selecting priority conservation areas based on zooplankton diversity: The case of Mediterranean wetlands. Mar. Freshw. Res. 2014, 65, 857–871. [Google Scholar] [CrossRef]
  9. Hsieh, C.H.; Sakai, Y.; Ban, S.; Ishikawa, K.; Ishikawa, T.; Ichise, S.; Yamamura, N.; Kumagai, M. Eutrophication and warming effects on long-term variation of zooplankton in Lake Biwa. Biogeosciences 2011, 8, 1383–1399. [Google Scholar] [CrossRef] [Green Version]
  10. Molinero, J.C.; Anneville, O.; Souissi, S.; Balvay, G.; Gerdeaux, D. Anthropogenic and climate forcing on the long-term changes of planktonic rotifers in Lake Geneva. J. Plankton Res. 2006, 28, 287–296. [Google Scholar] [CrossRef]
  11. Gilbert, J.D.; de Vicente, I.; Jiménez-Melero, R.; Parra, G.; Guerrero, F. Zooplankton body size versus taxonomy in Mediterranean wetlands: Implication for aquatic ecosystem evaluation. Freshw. Sci. 2017, 36, 774–783. [Google Scholar] [CrossRef]
  12. Schloesser, R.W.; Fabrizio, M.C. Temporal dynamics of condition for estuarine fishes in their nursery habitats. Mar. Ecol. Prog. Ser. 2016, 557, 207–2019. [Google Scholar] [CrossRef] [Green Version]
  13. Tlili, S.; Métais, I.; Ayache, N.; Boussetta, H.; Mouneyrac, C. Is the reproduction of Donax trunculus affected by their sites of origin contrasted by their level of contamination? Chemosphere 2011, 84, 1362–1370. [Google Scholar] [CrossRef]
  14. Durou, C.; Mouneyrac, C.; Amiard-Triquet, C. Environmental quality assessment in estuarine ecosystems: Use of biometric measurements and fecundity of the ragworm Nereis diversicolor (Polychaeta, Nereididae). Water Res. 2008, 42, 2157–2165. [Google Scholar] [CrossRef]
  15. Shadrina, N.V.; Anufriieva, E.V. Integral Indicators of Variability of Arctodiaptomus salinus (Daday, 1885) (Copepoda, Diaptomidae) and Their Possible Use in Assessing the Population State. Inland Water Biol. 2018, 11, 456–464. [Google Scholar] [CrossRef]
  16. Devreker, D.; Souissi, S.; Molinero, J.C.; Nkibuto, F. Trade-offs of the copepod Eurytemora affinis in mega-tidal estuaries. Insights from high frequency sampling in the Seine Estuary. J. Plankton Res. 2008, 30, 1329–1342. [Google Scholar] [CrossRef]
  17. Devreker, D.; Souissi, S.; Molinero, J.C.; Beyrend-Dur, D.; Gomez, F.; Forget-Leray, J. Tidal and annual variability of the population structure of Eurytemora affinis in the middle part of the Seine Estuary during 2005. Estuar. Coast. Shelf Sci. 2010, 89, 245–255. [Google Scholar] [CrossRef]
  18. Dur, G.; Souissi, S. Ontogenetic optimal temperature and salinity envelops of the copepod Eurytemora affinis in the Seine estuary (France). Estuar. Coast. Shelf Sci. 2018, 200, 311–323. [Google Scholar] [CrossRef]
  19. Yoshida, T.; Kagami, M.; Gurung, T.B.; Urabe, J. Seasonal succession of zooplankton in the north basin of Lake Biwa. Aquat. Ecol. 2001, 35, 19–29. [Google Scholar] [CrossRef]
  20. Anneville, O.; Molinero, J.C.; Souissi, S.; Balvay, G.; Gerdeaux, D. Long-term changes in the copepod community of Lake Geneva. J. Plankton Res. 2007, 29, i49–i59. [Google Scholar] [CrossRef] [Green Version]
  21. Lee, C.E. Rapid and Repeated Invasions of Fresh Water by the Copepod Eurytemora affinis. Evolution 1999, 53, 1423–1434. [Google Scholar] [CrossRef] [PubMed]
  22. Winkler, G.; Souissi, S.; Poux, C.; Castric, V. Genetic heterogeneity among several Eurytemora affinis populations in Western Europe. Mar. Biol. 2011, 158, 1841–1856. [Google Scholar] [CrossRef]
  23. Sukhikh, N.; Souissi, A.; Souissi, S.; Winkler, G.; Castric, V.; Holl, A.C.; Alekseev, V. Genetic and morphological heterogeneity among populations of Eurytemora affinis (Crustacea: Copepoda: Temoridae) in European waters. Comptes Rendus Biol. 2016, 239, 197–206. [Google Scholar] [CrossRef]
  24. Sukhikh, N.; Souissi, A.; Souissi, S.; Holl, A.C.; Schizas, N.V.; Alekseev, V. Life in sympatry: Coexistence of native Eurytemora affinis and invasive Eurytemora carolleeae in the Gulf of Finland (Baltic Sea). Oceanologia 2019, 61, 227–238. [Google Scholar] [CrossRef]
  25. Devreker, D.; Souissi, S.; Forget-Leray, J.; Leboulenger, F. Effects of salinity and temperature on the post embryonic development of Eurytemora affinis (Copepoda; Calanoida) of the Seine estuary: A laboratory study. J. Plankton Res. 2007, 29, i117–i133. [Google Scholar] [CrossRef] [Green Version]
  26. Devreker, D.; Souissi, S.; Winkler, G.; Forget-Leray, J.; Leboulenger, F. Effects of salinity and temperature on the reproduction of Eurytemora affinis (Copepoda; Calanoida) from the Seine estuary: A laboratory study. J. Exp. Mar. Biol. Ecol. 2009, 368, 113–123. [Google Scholar] [CrossRef]
  27. Souissi, A.; Souissi, S.; Hwang, J.S. Evaluation of the copepod Eurytemora affinis life history response to temperature and salinity increases. Zool. Stud. 2016, 55, e4. [Google Scholar]
  28. Cailleaud, K.; Budzinski, H.; Le Menach, K.; Souissi, S.; Forget-Leray, J. Uptake and elimination of hydrophobic organic contaminants in estuarine copepods: An experimental study. Environ. Toxicol. Chem. 2009, 28, 239–246. [Google Scholar] [CrossRef]
  29. Michalec, F.G.; Holzner, M.; Menu, D.; Hwang, J.S.; Souissi, S. Behavioral responses of the estuarine calanoid copepod Eurytemora affinis to sub-lethal concentrations of waterborne pollutants. Aquat. Toxicol. 2013, 138, 129–138. [Google Scholar] [CrossRef]
  30. Kadiene, E.U.; Bialais, C.; Ouddane, B.; Hwang, J.S.; Souissi, S. Differences in lethal response between male and female calanoid copepods and life cycle traits to cadmium toxicity. Ecotoxicology 2017, 26, 1227–1239. [Google Scholar] [CrossRef]
  31. Zidour, M.; Boubechiche, Z.; Pan, Y.J.; Bialais, C.; Cudennec, B.; Grard, T.; Drider, D.; Flahaut, C.; Ouddane, B.; Souissi, S. Population response of the estuarine copepod Eurytemora affinis to its bioaccumulation of trace metals. Chemosphere 2019, 220, 505–513. [Google Scholar] [CrossRef]
  32. Das, S.; Ouddane, B.; Hwang, J.S.; Souissi, S. Intergenerational effects of resuspended sediment and trace metal mixtures on life cycle traits of a pelagic copepod. Environ. Pollut. 2020, 267, 115460. [Google Scholar] [CrossRef]
  33. Dur, G.; Souissi, S.; Devreker, D.; Ginot, V.; Schmitt, F.G.; Hwang, J.S. Effect of temperature and individual variability on the reproduction of the estuarine copepod Eurytemora affinis Individual-based experimental and modelling approach. Ecol. Model. 2009, 220, 1073–1089. [Google Scholar] [CrossRef]
  34. Dur, G.; Jiménez-Melero, R.; Beyrend-Dur, D.; Hwang, J.S.; Souissi, S. Individual-based model of the phenology of egg-bearing copepods: Application to Eurytemora affinis from the Seine estuary, France. Ecol. Model. 2013, 269, 21–36. [Google Scholar] [CrossRef]
  35. Choi, B.S.; Kim, D.H.; Kim, M.S.; Park, J.C.; Lee, Y.H.; Kim, H.J.; Jeong, C.B.; Hagiwara, A.; Souissi, S.; Lee, J.S. The genome of the European estuarine calanoid copepod Eurytemora affinis: Potential use in molecular ecotoxicology. Mar. Pollut. Bull. 2021, 166, 112190. [Google Scholar] [CrossRef]
  36. Choi, B.S.; Han, J.; Hwang, D.S.; Souissi, S.; Hagiwara, A.; Lee, J.S. Complete mitochondrial genome of the calanoid copepod Eurytemora affinis (Calanoida, Temoridae). Mitochondrial DNA Part B 2019, 4, 2731–2733. [Google Scholar] [CrossRef] [Green Version]
  37. Lee, B.Y.; Lee, M.C.; Jeong, C.B.; Kim, H.J.; Hagiwara, A.; Souissi, S.; Han, J.; Lee, J.S. RNA-Seq-based transcriptome profiling and expression of 16 cytochrome P450 genes in the benzo[α]pyrene-exposed estuarine copepod Eurytemora affinis. Comp. Biochem. Physiol. Part D 2018, 28, 142–150. [Google Scholar] [CrossRef] [PubMed]
  38. Kwok, K.W.H.; Souissi, S.; Dur, G.; Won, E.J.; Lee, J.S. Copepods as references species in estuarine and marine waters. In Aquatic Ecotoxicology: Advancing Tools for Dealing with Emerging Risks; Amiard-Triquet, C., Amaird, J.C., Mouneyrac, C., Eds.; Elsevier: Amsterdam, The Netherlands, 2015; pp. 281–308. [Google Scholar]
  39. Zhang, C.; Jeong, C.B.; Lee, J.S.; Wang, D.; Wang, M. Transgenerational Proteome Plasticity in Resilience of a Marine Copepod in Response to Environmentally Relevant Concentrations of Microplastics. Environ. Sci. Technol. 2019, 53, 8426–8436. [Google Scholar] [CrossRef]
  40. Chaumot, A.; Coulaud, R.; Adam, O.; Quéau, H.; Lopes, C.; Geffard, O. In Situ Reproductive Bioassay with Caged Gammarus fossarum (Crustacea): Part 1—Gauging the Confounding Influence of Temperature and Water Hardness. Environ. Toxicol. Chem. 2020, 39, 667–677. [Google Scholar] [CrossRef]
  41. Crawford, P.; Daborn, G.R. Seasonal Variations in Body Size and Fecundity in a Copepod of Turbid Estuaries. Estuaries 1986, 9, 133–141. [Google Scholar] [CrossRef]
  42. Ali, A.K.; Primicerio, R.; Folstad, I.; Liljedal, S.; Berge, J. Morphological correlates of mating frequency and clutch size in wild caught female Eudiaptomus graciloides (Copepoda: Calanoida). J. Plankton Res. 2009, 31, 389–397. [Google Scholar] [CrossRef] [Green Version]
  43. Sukhikh, N.; Abramova, E.; Holl, A.C.; Souissi, S.; Alekseev, V. A comparative analysis of genetic differentiation of the E. affinis species complex and some other Eurytemora species, using CO1, ITSn and 18SrRNA genes (Copepoda, Calanoida). Crustaceana 2020, 93, 931–955. [Google Scholar] [CrossRef]
  44. Alekseev, V.; Souissi, A. A new species within the Eurytemora affinis complex (Copepoda: Calanoida) from the Atlantic Coast of USA, with observations on eight morphologically different European populations. Zootaxa 2011, 2767, 41–56. [Google Scholar] [CrossRef]
  45. Cailleaud, K.; Forget-Leray, J.; Pehulet, L.; Le Menach, K.; Souissi, S.; Budzinski, H. Tidal influence on the distribution of hydrophobic organic contaminants in the Seine Estuary and on their biological effects on the copepod Eurytemora affinis. Environ. Pollut. 2009, 157, 64–71. [Google Scholar] [CrossRef]
  46. Cailleaud, K.; Forget-Leray, J.; Souissi, S.; Hilde, D.; Le Menach, K.; Budzinski, H. Seasonal variations of hydrophobic organic contaminant concentrations in the water-column of the Seine Estuary and their transfer on a planktonic species Eurytemora affinis (Calanoïda, copepoda). Part1: PCBSs and PAHs. Chemosphere 2007, 70, 270–280. [Google Scholar] [CrossRef]
  47. Michalec, F.G.; Fouxon, I.; Souissi, S.; Holzner, M. Zooplankton can actively adjust their motility to turbulent flow. Proc. Natl. Acad. Sci. USA 2017, 114, E11199–E11207. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Michalec, F.G.; Holzner, M.; Souissi, A.; Stancheva, S.; Barras, A.; Boukherroub, R.; Souissi, S. Lipid nanocapsules for behavioural testing in aquatic toxicology: Time response of Eurytemora affinis to environmental concentrations of PAHs and PCB. Aquat. Toxicol. 2016, 170, 310–322. [Google Scholar] [CrossRef] [Green Version]
  49. Cailleaud, K.; Maillet, G.; Budzinski, H.; Souissi, S.; Forget-Leray, J. Effects of salinity and temperature on the expression of enzymatic biomarkers in Eurytemora affinis. Comp. Biochem. Physiol. A 2007, 147, 841–849. [Google Scholar] [CrossRef]
  50. Dayras, P.; Bialais, C.; Lee, J.S.; Souissi, S. Effects of microalgal diet on the population growth and fecundity of the cyclopoid copepod Paracyclopina nana. J. World Aquac. Soc. 2020, 51, 1386–1401. [Google Scholar] [CrossRef]
  51. Somavilla, R.; González-Pola, C.; Rodriguez, C.; Josey, S.A.; Sánchez, R.F.; Lavín, A. Large changes in the hydrographic structure of the Bay of Biscay after the extreme mixing of winter 2005. J. Geophys. Res. 2009, 114, C01001. [Google Scholar] [CrossRef] [Green Version]
  52. Gómez, F.; Souissi, S. The impact of the 2003 summer heat wave and the 2005 late cold wave on the phytoplankton in the north-eastern English Channel. Comptes Rendus Biol. 2008, 331, 678–685. [Google Scholar] [CrossRef] [PubMed]
  53. Kimmel, D.G.; Miller, W.D.; Roman, M.R. Regional Scale Climate Forcing of Mesozooplankton Dynamics in Chesapeake Bay. Estuaries Coasts 2006, 29, 375–387. [Google Scholar] [CrossRef]
  54. Bessineton, C. Historical records of the Nereis diversicolor population in the Seine estuary. In Environmental Assessment of Estuarine Ecosystems; Amiard-Triquet, C., Rainbow, P.S., Eds.; A Case Study; CRC Press—Taylor & Francis Group: Boca Raton, FL, USA, 2009; pp. 183–197. [Google Scholar]
  55. Bunker, A.J.; Hirst, A.G. Fecundity of marine planktonic copepods: Global rates and patterns in relation to chlorophyll a, temperature and body weight. Mar. Ecol. Prog. Ser. 2004, 279, 161–181. [Google Scholar] [CrossRef]
  56. Lloyd, S.S.; Elliott, D.T.; Roman, M.R. Egg production by the copepod, Eurytemora affinis, in Chesapeake Bay turbidity maximum regions. J. Plankton Res. 2013, 35, 299–308. [Google Scholar] [CrossRef] [Green Version]
  57. Pierson, J.J.; Kimmel, D.G.; Roman, M.R. Temperature Impacts on Eurytemora carolleeae Size and Vital Rates in the Upper Chesapeake Bay in Winter. Estuaries Coasts 2016, 39, 1122–1132. [Google Scholar] [CrossRef]
  58. Hülsmann, S. Reproductive potential of Daphnia galeata in relation to food conditions: Implications of a changing size-structure of the population. Hydrobiologia 2001, 442, 241–252. [Google Scholar] [CrossRef]
  59. Seebens, H.; Straile, D.; Hoegg, R.; Stich, H.B.; Einsle, U. Population dynamics of a freshwater calanoid copepod: Complex responses to changes in trophic status and climate variability. Limnol. Oceanogr. 2007, 52, 2364–2372. [Google Scholar] [CrossRef]
  60. Gilbert, J.D.; de Vicente, I.; Ortega, F.; Jiménez-Melero, R.; Parra, G.; Guerrero, F. A comprehensive evaluation of the crustacean assemblages in southern Iberian Mediterranean wetlands. J. Limnol. 2015, 74, 169–181. [Google Scholar] [CrossRef] [Green Version]
  61. Jiménez-Melero, R.; Gilbert, J.D.; Guerrero, F. Secondary production of Arctodiaptomus salinus in a shallow saline pond: Comparison of methods. Mar. Ecol. Prog. Ser. 2013, 483, 103–116. [Google Scholar] [CrossRef]
  62. Chow-Fraser, P.; Maly, E.J. Factors Governing Clutch Size in Two Species of Diaptomus (Copepoda: Calanoida). Can. J. Fish. Aquat. Sci. 1991, 48, 364–370. [Google Scholar] [CrossRef]
  63. Lajus, D.; Sukhikh, N.; Alekseev, V. Stochastic phenotypic variation: Empirical results and potential use in Eurytemora research (Copepoda, Calanoida). Crustaceana 2020, 93, 317–336. [Google Scholar] [CrossRef]
  64. Head, E.J.H.; Harris, L.R.; Ringuette, M.; Campbell, R.W. Characteristics of egg production of the planktonic copepod, Calanus finmarchicus, in the Labrador Sea: 1997–2010. J. Plankton Res. 2013, 35, 281–298. [Google Scholar] [CrossRef] [Green Version]
Figure 1. (A) Experimental relationship between clutch size (CS) and prosome length (PL) of Eurytemora affinis from the Seine estuary. The different experimental conditions T20S15, T24S15, T20S25, T24S25 giving in total 43 generations are represented by different symbols. The continuous line corresponds to the linear regression CS = 0.2PL − 145.3. (B). Histogram of the distribution of the 5000 values of coefficient of correlation (r) obtained by sampling the original data using a bootstrap technique. (C). Similar to B. but showing the distribution of values of the slope parameter (a). (D). Similar to A. but showing the distribution of values of the intercept (value corresponding to predicted CS = 0).
Figure 1. (A) Experimental relationship between clutch size (CS) and prosome length (PL) of Eurytemora affinis from the Seine estuary. The different experimental conditions T20S15, T24S15, T20S25, T24S25 giving in total 43 generations are represented by different symbols. The continuous line corresponds to the linear regression CS = 0.2PL − 145.3. (B). Histogram of the distribution of the 5000 values of coefficient of correlation (r) obtained by sampling the original data using a bootstrap technique. (C). Similar to B. but showing the distribution of values of the slope parameter (a). (D). Similar to A. but showing the distribution of values of the intercept (value corresponding to predicted CS = 0).
Water 13 01442 g001
Figure 2. Projection of the mean values of clutch size (CS) and prosome length (PL) of ovigerous females sampled in the Seine estuary (circle symbols), Loire estuary (square symbols) and other (Elbe, Scheldt and Gironde) European estuaries (triangle symbols). The continuous line corresponds to the reference regression line obtained from Figure 1. The month and year of each sampling date is shown in the figure (black text for the Seine estuary and grey colour and italic for Loire, Gironde, Scheldt and Elbe estuaries). The different colours correspond to the values of the QI computed for each sampling date and indicated in Table 3. Blue colour: very good conditions (QI ≥ +20%), green colour: good conditions (−20% ≤ QI < −20%), yellow colour: bad conditions (−40% ≤ QI < −20%) and red colour: very bad conditions (QI ≤ −40%).
Figure 2. Projection of the mean values of clutch size (CS) and prosome length (PL) of ovigerous females sampled in the Seine estuary (circle symbols), Loire estuary (square symbols) and other (Elbe, Scheldt and Gironde) European estuaries (triangle symbols). The continuous line corresponds to the reference regression line obtained from Figure 1. The month and year of each sampling date is shown in the figure (black text for the Seine estuary and grey colour and italic for Loire, Gironde, Scheldt and Elbe estuaries). The different colours correspond to the values of the QI computed for each sampling date and indicated in Table 3. Blue colour: very good conditions (QI ≥ +20%), green colour: good conditions (−20% ≤ QI < −20%), yellow colour: bad conditions (−40% ≤ QI < −20%) and red colour: very bad conditions (QI ≤ −40%).
Water 13 01442 g002
Figure 3. Mean values of the coefficient of variation (CV) of prosome length (A) and clutch size (B) obtained in field and experimental conditions. Error bars correspond to the standard error. Stars indicate significant difference using Wilcoxon rank sum nonparametric test (p < 0.05).
Figure 3. Mean values of the coefficient of variation (CV) of prosome length (A) and clutch size (B) obtained in field and experimental conditions. Error bars correspond to the standard error. Stars indicate significant difference using Wilcoxon rank sum nonparametric test (p < 0.05).
Water 13 01442 g003
Table 1. Definition of 4 classes of the quality index (QI) values and their associated score/colour.
Table 1. Definition of 4 classes of the quality index (QI) values and their associated score/colour.
Class NumberQI ValuesEcological StatusColour
1QI ≥ +20%very goodblue
2−20% ≤ QI < +20%goodgreen
3−40% ≤ QI < −20%badyellow
4QI < −40%very badred
Table 2. Description of the different categories of data obtained from the field and in the laboratory.
Table 2. Description of the different categories of data obtained from the field and in the laboratory.
CategoryDescriptionSample Size (n)
In situ allAll sampling dates19
In situ goodSampling dates with good or very good QI12
In situ badsampling dates with bad or very bad QI7
T20S15generations followed at T = 20 °C and S = 1517
T24S15generations transferred from T20S15 to T = 24 °C and S = 155
T20S25generations followed at T = 20 °C and S = 2516
T24S25generations transferred from T20S25 to T = 24 °C and S = 255
Table 3. Summary of the sampling dates realised in the lower Seine estuary (salinity gradient zone) and the number of ovigerous females used to calculate the quality index (QI) between 2004 and 2010. Mean values of prosome length (PL; in µm) and clutch size (CS; number of eggs per female) for each date were associated with the coefficient of variation (CV = 100 × SD/mean, where SD is the standard deviation). The predicted clutch size (CSpred) for each date is obtained using the regression line of Figure 1 (CS = 0.1988 × PL − 145.3). Colours in the last column of the QI are the same as in Figure 2, and they represent arbitrary 4 classes: very good (blue), QI is greater than 20%; good (green), QI ranged between −20% and 20%; bad (yellow), QI ranged between −40% and −20%, and very bad (red), QI is lower than −40%.
Table 3. Summary of the sampling dates realised in the lower Seine estuary (salinity gradient zone) and the number of ovigerous females used to calculate the quality index (QI) between 2004 and 2010. Mean values of prosome length (PL; in µm) and clutch size (CS; number of eggs per female) for each date were associated with the coefficient of variation (CV = 100 × SD/mean, where SD is the standard deviation). The predicted clutch size (CSpred) for each date is obtained using the regression line of Figure 1 (CS = 0.1988 × PL − 145.3). Colours in the last column of the QI are the same as in Figure 2, and they represent arbitrary 4 classes: very good (blue), QI is greater than 20%; good (green), QI ranged between −20% and 20%; bad (yellow), QI ranged between −40% and −20%, and very bad (red), QI is lower than −40%.
Sampling Date (DD/MM/YYYY)Number of Females (n)Mean PLCV (%) of PLMean CSCV (%) of CSPredicted CS (CSpred)Quality Index (QI)
16/06/200450886.45.8338.624.1330.9224.85
08/03/2005701006.766.2825.7338.9854.84−53.09
12/04/2005641072.786.2132.644.7667.97−52.04
21/05/2005103974.294.0541.1631.3948.39−14.94
10/06/200560988.33.5557.0218.2051.1711.42
21/07/200558937.73.1743.1229.2641.114.88
16/09/200528940.13.9638.6428.4741.59−7.10
01/12/200529981.753.9532.5225.9549.87−34.79
23/05/200620949.854.0044.812.8943.532.92
20/11/200650840.865.4718.0621.1521.86−17.39
05/10/200740846.334.3525.551.0022.9511.11
16/11/200780846.294.4313.736.0722.94−40.29
19/03/200840978.963.6619.654.2949.32−60.26
10/04/200859988.165.3439.348.4551.15−23.16
18/06/200850899.843.8634.637.9833.593.01
08/07/200886879.175.1126.867.5729.48−9.09
04/11/200850831.794.0414.952.8320.06−25.72
25/02/2009201127.423.1476.5531.6478.83−2.89
02/10/2010181042.414.9952.2231.2061.93−15.68
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Souissi, S.; Souissi, A. Promotion of the Development of Sentinel Species in the Water Column: Example Using Body Size and Fecundity of the Egg-Bearing Calanoid Copepod Eurytemora affinis. Water 2021, 13, 1442. https://doi.org/10.3390/w13111442

AMA Style

Souissi S, Souissi A. Promotion of the Development of Sentinel Species in the Water Column: Example Using Body Size and Fecundity of the Egg-Bearing Calanoid Copepod Eurytemora affinis. Water. 2021; 13(11):1442. https://doi.org/10.3390/w13111442

Chicago/Turabian Style

Souissi, Sami, and Anissa Souissi. 2021. "Promotion of the Development of Sentinel Species in the Water Column: Example Using Body Size and Fecundity of the Egg-Bearing Calanoid Copepod Eurytemora affinis" Water 13, no. 11: 1442. https://doi.org/10.3390/w13111442

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop