Composition, Structure and Functional Feeding of Aquatic Entomofauna in Kodjoboué Lake: Water Quality Assessment

The composition, structure and Functional Feeding Groups (FFGs) of aquatic insects were investigated in Kodjoboué lake (South-East of Côte d’Ivoire) to assess water quality. Four sites (K1 - K4) were sampled during eight campaigns. At each site aquatic insects were collected using a kick-net (25 × 20 cm, 500 µm mesh size). Environmental variables were also recorded during this period. Water was collected from each sampling site for nutrient (phos-phorus, nitrate) analysis in laboratory. A total of 68 aquatic insect taxa be-longing to six orders and 32 families were recorded during the study. Site K4, subject to intense agricultural activity, had the lower diversity and least balanced stand with the lowest values of Shannon-Weaver index, Pielou’s Evenness index and rarefied richness. The functional approach revealed five major FFGs: predators, scrapers, shredders, gatherers-collectors and fil-ter-collectors. Among the different functional feeding groups, the percentage of predators was higher than that of other groups at all the sites. The P/R re-port indicated that our four sites are characterized as heterotrophic with K1 as the least heterotrophic. Focused Principal Component Analysis (FPCA) showed that parameters of water, like pH, dissolved oxygen, nitrate were sta-tistically correlated with species abundance. The diversity and Functional Feeding Groups metrics of aquatic insects used to describe the state of Kodjoboué lake reveal its disturbed state.


Introduction
Aquatic macroinvertebrates are the main group of organisms using to assess the ecological state of continental hydrosystems. They are recognized as the best indicators of the ecological health of hydrosystems through the use of the structure of their communities as an indicator of the effects of human activity on hydrosystems [1] [2] [3] but also, by having within them, bio-indicator species of water quality. Of the biological quality elements that have been used for water quality assessment, aquatic insects are most often recommended [4] due to their varying sensitivity to several abiotic and biotic factors in the environment [5]. They constitute the most diverse taxonomic group and the most abundant in aquatic environments [6]. The use of these organisms in bioindication is commonly done with taxonomic descriptors (diversity indices, taxonomic richness, relative abundances, indicator species) [7] [8] [9]. However, the functional approach is more appropriate and faster than the taxonomic approach to characterize ecosystem conditions [10]. As a result, the functional approach is increasingly used in the assessment of surface water quality [11] [12] [13]. Also, functional feeding groups measures do not require standardized sampling and different sampling methods can be used to assign aquatic invertebrates to functional feeding groups [14].
Despite the growing interest in the study of functional diversity, no attempt to our knowledge has been made to establish a diagnosis of water quality in Côte d'Ivoire.
In addition, while much work has been done on the macrofauna of some lakes in Côte d'Ivoire [15] [16] [17] [18] [19], very little data exists on Kodjoboué lake [20]. Given the importance of the functional feeding groups approach in biomonitoring and conservation, the evaluation of the functional organization of macroinvertebrate communities is essential [21], particularly that of aquatic insects. The aim of this study is to use the taxonomic composition and the functional feeding groups of the aquatic insects to assess the ecological state of Kodjoboué lake (South-East, Côte d'Ivoire).

Study Area and Sampling Sites
Kodjoboué lake is located in south-eastern Côte d'Ivoire in the South-Comoé region about 4 km south of the town of Bonoua between longitudes 3˚35'9"W and latitudes 5˚14'11"N. This area is contiguous to the Grand-Bassam wetland and covers an area of 3.404 Km 2 [20]. This lake communicates with the Comoé river by a thin stream of water. This region has an average altitude of 90 m, its climate is subequatorial with an average annual temperature ranged between 25˚C to 33˚C. The rainfall varies from 1400 to 2500 mm/year and an annual rate of humidity about 80% to 90% [22]. This climate is characterized by four seasons: two dry seasons (mid-July to mid-September and December to March) and two rainy seasons (April to mid-July and mid-September to mid-November) [23]. Four sampling sites were selected in Kodjoboué lake (K1 to K4) ( Figure 1). These sites were sampled during eight sampling campaigns from June 2014 to May 2017. Table 1 summarizes the characteristics of these sampling sites.

Data Collection
Aquatic insects of the native vegetation, covering the banks of the lake were  sampled at the sites during the eight sampling campaigns with a kick-net (25 × 20 cm, 500 µm mesh size). At each site 5 net strokes are sampled. A net stroke consists of pulling the net over a distance approximately 1 m. The initial net is followed by two round trips on the surface to capture the dislodged organisms. The samples were then pooled to make one composite sample per site. Samples were fixed on 70% ethanol and stored in labelled polyethylene flasks.
In the laboratory, specimens were sorted and identified to the lowest possible taxonomic level using a stereomicroscope Olympus SZ (40× magnification) and a series of identification keys [24]- [29].

Data Analysis
Species diversity was calculated using the Shannon-Weaver diversity index and Pielou evenness. Taxonomic richness was rarefied in each site per sampling period. Precisely, rarefied richness was used to avoid any bias related to differences in abundances between samples [31]. Before performing comparison analyses, data normality was checked using Shapiro test. As the biotic and environment data distribution follow non-normal distribution (P < 0. Focused Principal Component Analysis (FPCA) [33] was used to assess relationships between taxa abundance and environmental variables. Focused Prin-cipal Component Analysis (FPCA) is a variation of the traditional principal component analysis. It uses the same types of matrix as the PCA but differs from it in that it is centred or focused on a variable xi. Indeed, it allows a graphical representation of the correlations that exist between this variable xi and the other variables. The graph provides access not only to the nature (positive or negative) but also to the significance (P < 0.05) of the correlations between variable xi and the other variables. It is also possible to observe on the graph the correlations between the other variables.
In this study, the FPCA was used to determine the variables that influence the abundance of macroinvertebrate taxa that contribute at least 5% of the total abundance. This analysis was performed using the psy package on the R software. Table 2

Composition and Distribution of Aquatic Insects
A total of 5562 aquatic insect individuals from six orders, 32 families and 68 taxa have been recorded. The order Coleoptera was the most diversified with 19 taxa and eight families. The order of the Diptera was as follows with 15 taxa and  The qualitative data on sampling at the lake sites considered are summarized in Table 5. Site K1, with 54 taxa, has the highest taxonomic richness. It is followed by the site K3 with 36 taxa. Site K2, with 32 taxa, has the lowest taxonomic richness of all the stations considered. Of the taxa identified, only 13 are common to all four sites.

Frequency of Occurrence
The different sites are compared with each other based on their taxonomic composition (Table 5). Sites K1 and K4 contain more common taxa (26), site pairs (K1 -K3) and (K1 -K3) contain fewer common taxa with 18 taxa each. It should also be noted that site K3 contains more specific taxa (04), while sites K1, K2 and K4 each contain only one specific taxon ( Figure 2).

Diversity Metrics
The lowest values of the three diversity metrics, Shannon-Weaver index, Pielou's Evenness Index and rarefied richness were recorded at site K4. However, the highest values were observed at the site K3 for these three metrics. Of these three metrics, only rarefied richness does not show significant variation between the different sites (P > 0.05; Mann-Whitney) (Figure 3).

Functional Feeding Groups
Five major FFGs were identified in this study: predators (pre), scrapers (scr), shredders (shr), gatherers-collectors (gc), and filter-collectors (fc). The analysis     of the functional composition of the assemblage revealed that the predators were the most abundant ( Figure 4) with a proportion of 53.74% of the insects harvested. Gatherers-collectors ranked second with 42.05% of the total insect population. A proportion of 2.84% was observed for the shredders, while filter-collectors and scrapers had the lowest proportions respectively 0.91% and 0.44% of the total abundance.  (Table 6).

Correlation between Aquatic Insects and Environment Variables
Focused Principal Component Analysis (FPCA) was performed using environmental parameters with a focus on the most abundant taxa (at least 5% of total abundance). Thus, the FPCA states that Hydrovatus sp.is significantly and negatively influenced (P < 0.05) by dissolved oxygen (O 2 ). Also, this insect is positively correlated with parameters such as conductivity (CND), nitrate (NO 3 ), temperature (T) and transparency (Transp) but negatively correlated with pH and phosphate (PO 4 ). Polypedilum deletum is significantly and positively correlated with nitrate (P < 0.05). This taxon is negatively correlated with temperature, pH and dissolved oxygen but positively correlated with other parameters.
Cloeon sp. is significantly correlated with dissolved oxygen and pH (P < 0.05).
This Ephemeroptera is positively correlated with transparency, but negatively correlated with other parameters. As for Esolus sp., it is negatively correlated with pH and temperature but positively correlated with conductivity, nitrate, phosphate and transparency. Esolus sp. is significantly correlated with dissolved oxygen (P < 0.05) ( Figure 6).

Discussion
In this study, the water temperature range of Kodjoboué Lake was between 25.9˚C and 31.9˚C and was considered relatively warm. This trend reflects the temperature in the tropics. This range is close to that observed in Taabo Lake by Kouamé et al. (2010) [17]. Our low conductivity levels (14 to 37.1 μS/cm) would indicate low levels of natural or anthropogenic sediments throughout Kodjoboué Lake. With regard to dissolved oxygen, the range is between 0.10 and 11.3 mg/l. This range is greater than that obtained at Taabo Lake (3.62 -9.2 mg/l) by Kouamé et al., (2010) [17]. The amplitude of temperature and dissolved oxygen variations can be linked to the lack of coverage plant on the lake. The transparency of the water varies between 10 and 218 cm at Lake Kodjoboué. This range is higher than that observed (73 and 112 cm) at Ayamé Lake by Diétoa (2002) [16].  conditions become hostile to reduced interspecific competition in other species [34].
In this study, the lowest values of the three metrics of diversity calculated the Chironomidae, Polypedilum deletum is the only taxon that was most common at all sites. The regular harvest of these chironomids would be linked to their broad spectrum of resilience and morphological, physiological and behavioural adaptations.
Among the different functional feeding groups, the percentage of predators was higher than that of the other groups in our sites, followed by gatherers-collectors.
Predator preponderance would probably be related to prey availability, especially gatherer-collectors such as Chironomidae ranked second in abundance. In fact, predatory organisms reproduce proportionally to the availability of food. The growth rate of predators is therefore an increasing function of prey [35].
Shredders, filter-collectors and scrapers had the lowest proportions of total abundance. The low proportion of shredders in this lake is thought to be related to the fact that in tropical areas high temperatures increase the decomposition of litter and leaves by microbial activity [36].  [42], specialized feeders, such as shredders and scrapers, are presumed to be more sensitive to disturbance. Their low proportions are due to a disturbed environment.
The P/R ratio (<0.75 at all sites) indicates that our four sites are characterized as heterotrophic. The K1 site is the least heterotrophic compared to the other three sites. According to Savić et al. (2017) [43], this would mean that the dominant staple food chain for aquatic insect communities has been found to be non-native detritus. Concerning the Top-down predator control, all values are above 0.2 at all our sites. According to Masese et al. (2014) [44], the value of the predator-prey ratio must be between 0.1 and 0.2 for there to be a balance between predators and prey. In other words, none of our sites have predator-prey balance or top-down control. So, we have an overabundance of predators in Lake Kodjoboué.
Correlations between the most abundant taxa and environmental parameters using the Focused Principal Component Analysis (FPCA) show a strong negative correlation between Hydrovatus sp. (Predator) and dissolved oxygen. This correlation reveals the ability of this Beetle to live in anoxia. Unlike the Ephemeroptera Cloeon sp. is positively correlated with dissolved oxygen. Indeed, Ephemeroptera are known to live in highly oxygenated environments [45].

Conclusion
This work enabled us to identify 68 aquatic insect taxa divided into six orders and 32 families. Beetles were the most abundant in both quality and quantity. Site K4, subject to intense agricultural activity, had the least diversified and least balanced stand with the lowest values of the Shannon-Weaver index, the Pielou Evenness index and rarefied richness. Functionally, five major FFGs were identified in this study: predators (pre), scrapers (scr), shredders (shr), gatherers-collectors (gc) and filter-collectors (fc). Among the different functional feeding groups, the percentage of predators was higher than that of other groups at our sites. The P/R ratio indicated that our four sites, are characterized as heterotrophic with K1 as the least heterotrophic. Top-down predator control ratio revealed an overabundance of predators. The diversity and Functional Feeding Groups metrics of aquatic insects used to describe the state of Kodjoboué lake reveal a disturbed state, therefore of poor quality. It would be important to implement a sustainable management program for this lake in order to conserve its biodiversity.