Benthic Macroinvertebrates as Ecological Indicators: Their Sensitivity to the Water Quality and Human Disturbances in a Tropical River

Macroinvertebrate metrics are helpful tools for the assessment of water quality and overall aquatic ecosystem health. However, their degree of sensitivity and the most reliable metrics for the bioassessment program development are very poorly studied in Togo. This study aimed to test the sensitivity of metrics calculated at the family and genus levels. A total of 21 water quality parameters and macroinvertebrates' data were collected during three periods at 20 sampling sites within the Zio River. The canonical correspondence analysis (CCA), factor analysis (FA), and Spearman's correlation analysis were conducted on water quality parameters and macroinvertebrates' data. The results reveal that macroinvertebrate structure and composition were affected by water quality parameters related to human disturbances. In this study, three groups of macroinvertebrate communities were identified including sensitive taxa such as Ephemeroptera, Plecoptera, Trichoptera, and Odonata (EPTO) taxa; the resistant or resilient taxa such as Oligochaeta, Hirudinea, Diptera, and Pulmonates (OHDP) taxa; and tolerant taxa such as Prosobranchia, Bivalvia, Lepidoptera, Heteroptera, and Coleoptera (PBLHC). All the 13 macroinvertebrate-based metrics were found to be sensitive in the detection of water quality and human disturbance gradients. However, metrics related to EPTO and the tolerance measure [multimetric index of the Zio River basin (MMIZB), Average Score per Taxon (ASPT), and Biological Monitoring Working Party (BMWP)] are the most robust in discrimination of pressure gradients. This study reveals that macroinvertebrates are sensitive and can be used for the bioassessment program development at the order, family, or genera taxonomic level.


INTRODUCTION
Rivers are important ecosystems with high ecological value (Nguyen et al., 2018), and the health of these ecosystems is important for the human societies that depend on them (Dickens et al., 2018). Rivers are an important source of renewable water supply for human beings and freshwater ecosystems (Vorosmarty et al., 2010) and provide many goods and services such as domestic uses, navigation, recreational activities, and nursing grounds and food for many organisms (Berger et al., 2016). Until the late 1960s, the overriding interest in water has been on the available amount for consumption. Except when undesirable water quality conditions persist, the available water was considered acceptable for consumption. Only during the last three decades of the 20th century has water quality been deteriorated to the point where it is considered as important as managing for the availability of water (Abbasi and Abbasi, 2012).
Globally, human impact is changing the availability of freshwater (Rodell et al., 2018). Among freshwater ecosystems, streams and rivers are the most influenced or threatened by a range of anthropogenic stresses (Allan, 2004;Best and Darby, 2020). In the same line, Dudgeon et al. (2006) grouped the major threats to freshwater species under five interacting categories: overexploitation; water pollution; flow modification; destruction or degradation of habitat; and invasion by exotic species. The degree of these stressors led to a geographical variation in the threat to freshwater species with almost one in three freshwater species threatened with extinction worldwide (Collen et al., 2013), while estimates suggest that at least 10,000-20,000 freshwater species are extinct or at risk of extinction worldwide (Balian et al., 2008;Vorosmarty et al., 2010). For example, 1,116 freshwater migratory fish species are threatened and with 102 additional extinctions (Hogan, 2011). In tropical regions, particularly in Sub-Saharan countries, rivers are under pressure due to human activities that deteriorate water quality, limiting water availability for drinking and other uses (Traoré et al., 2016;Kaboré et al., 2018;Tampo, 2018;Chetty and Pillay, 2019;Agboola et al., 2020a). Furthermore, the river pollution problem, related to human disturbances through anthropogenic activities and urbanization (Azrina et al., 2006;Faridah et al., 2012;Edegbene et al., 2021), has many adverse impacts on river ecosystems that have required managers to increase water assessment efforts (Clifford and Tariro, 2005). Many African countries have concerns about the ecological status of their rivers and are increasing their investments on the restoration of degraded rivers Fayiga et al., 2018;Wantzen et al., 2019).
In many studies, biological methods are valuable for determining natural and anthropogenic influences on water resources and habitats because biota responds to stressors from multiple spatial or temporal scales (Weigel and Robertson, 2007;Resende et al., 2010;He et al., 2020;Kurthen et al., 2020). In addition, the use of aquatic organisms in ecological studies has proven more effective than using environmental variables alone, because the aquatic community integrates structural and functional characteristics and reflects the health of the studied streams (Bonada et al., 2006;He et al., 2020). Accordingly, there is an increasing interest in the application of ecological thresholds for natural resource management (King and Baker, 2010;Forio et al., 2018). Several metrics and biotic indices are developed and used across the world (Barbour et al., 1996;Hering et al., 2006;Agboola et al., 2019;He et al., 2020;Ko et al., 2020. Gonçalves and argued that the use of biotic indices as a tool for river quality assessment was more useful in evaluating river health than the conventional national water quality assessment standard practices in many countries. Metcalfe (1989) and Alba-Tercedor (1996) on the other hand pointed at the disadvantage of the physicochemical assessment, which measures the water quality only at the time of sampling.
The freshwater macroinvertebrate taxa vary in response to organic pollution; and thus, their diversity and composition have been used to make inferences about pollution loads (Kenney et al., 2009;Wan Abdul Ghani et al., 2018;Tampo et al., 2020). Generally, in natural pristine rivers, high abundances and richness of species could be found (Barbour et al., 1996(Barbour et al., , 1999. However, high impact due to human activities caused many changes of the assemblages and biodiversity of the river fauna (Wright et al., 1993;Pinel-Alloul et al., 1996). Despite the development and application of a variety of biotic indices, scores, and metrics based on macroinvertebrates for water quality and ecosystem health assessment in America and Europe, the literature provides little information on biological assessment and monitoring tools of freshwater ecosystems in Sub-Saharan Africa especially in West Africa (Adams, 1993;Dallas, 1997;Kaboré et al., 2016). In Togo, there is paucity of information on bioassessment of riverine ecosystems health, except for some works on hydrochemistry, macroinvertebrates as indicators, and multimetric index of the Zio River basin (MMIZB) carried out by Tampo et al. (2015Tampo et al. ( , 2020 and Tampo (2018). The measure of macroinvertebrate sensitivity at different taxonomic scales is however very important for bioassessment (Bailey et al., 2001;Schmidt-Kloiber and Nijboer, 2004;Marshall et al., 2006;Metzeling et al., 2006;Jones, 2008;Costas et al., 2018) and essential for developing a biomonitoring program in Togo. Therefore, the aim of this study was to test the effectiveness and the robustness of a set of macroinvertebrate-based metrics in the Zio River.

Study Area and Sampling Sites
This study was undertaken in the main stem of the Zio River. It is located in a tropical climate between the latitude 6 • 5 ′ and 7 • 18 ′ N ( Figure 1) and includes three of the five ecoregions of Togo. The agriculture and other human activities described previously (Tampo et al., 2015(Tampo et al., , 2020Tampo, 2018) are more intensive in the lower reaches of the Zio River. The river crosses some villages and the outskirts of the city of Lomé; and its water is the major source for domestic and agricultural purposes in rural and semi-urban areas. Moreover, the main activity in the basin is agriculture (maize, cassava, bean, yams, rice, etc.), fishing, and few industrial activities such as gravel washing and sand extraction (Bawa et al., 2018). According to Tampo (2018), the river can be subdivided into three sections, such as the upper section, the middle section, and the lower section, which is influenced by seawater during the low flow. Twenty sampling sites were selected across the study area, which included a range of impaired and unimpaired ecological states within the Zio River following anthropogenic disturbances (Tampo et al., 2020).

Sampling and Water Parameter Analyses
Twenty sites were sampled three times between 2013 and 2014 during different seasons (the rainy season, the dry season, and a transition period between the rainy season and dry season). Therefore, 60 samples were collected and 21 water quality parameters measured. At each site, the water quality parameters such as temperature, electrical conductivity (EC), dissolved oxygen (DO), and pH were measured in situ using Wissenschaftliche Technische Werkstätten (WTW) multipurpose water quality probe. After the in situ measurement, 1.5 L of water was taken in a plastic bottle and stored in a cool environment for analysis of chemical parameters in the , and NH + 4 ), while K + and Na + were determined by flame emission spectrophotometer. Total suspended solid (TSS) was determined by gravimetric method (dried at 105 • C). All these parameters were measured in the Laboratory of Applied Hydrology and Environment (LAHE) of the Université de Lomé (Togo) with an accuracy ranking from 1 to 2% according to the standard methods as prescribed by AFNOR (1997) and Rodier et al. (2009).

Macroinvertebrate Collection and Identification
Benthic macroinvertebrates were sampled using a dip net (circular opening, 33 cm of diameter; mesh size, 320 µm) in lentic habitats and a Surber Sampler (rectangular opening, 20 cm × 25 cm; mesh size, 320 µm) for lotic habitats. At each site, substrate samples were taken and combined to one composite sample following protocol described by AFNOR (1997) and Rodier et al. (2009). All animals were identified at the family level and at the genus level for mollusk gastropod, Annelida, and insects using taxonomic manuals and keys (Durand and Levêque, 1981;Merritt and Cummins, 1996;Tachet et al., 2010).

Metric Selection
In this study, 13 metrics classified in three categories were used to assess the ecological status of the Zio River and to test their sensitivities (Table 1). These metrics were calculated using macroinvertebrate features (number of taxa, diversity indices, abundance, and tolerance score). The selection of these metrics was based on their simplicity and reliability for assessing the water quality of the river as well as their suitability to detect anthropogenic disturbances (Raburu et al., 2009;Jun et al., 2012;Nguyen et al., 2014;Kaboré et al., 2016;Agboola et al., 2019;Tampo et al., 2020). Table 1 indicates the metrics used and their definition.

Statistical Analysis
The canonical correspondence analysis (CCA) was applied in order to establish the relationship between water quality parameters and macroinvertebrate abundance and to identify water quality parameters affecting macroinvertebrate community. The analysis was performed under PAST software (version 3.0) based on dataset of water quality parameters and a dataset of abundance of macroinvertebrate community at phylum-class, subclass-order, and family and genus levels. In addition, we assessed the potential of macroinvertebrate taxa and metrics detected in this study to serve as bioindicators for the river's environmental condition investigated. Therefore, factor analysis (FA) and Spearman's correlation analysis between macroinvertebrate data and water quality variables were used, and with expert's consensus following Kaboré et al. (2016) and Tampo et al. (2015). The FA was performed using principal components as factor extraction method without any rotation. The factor loadings were considered for the explanation of correlations among variables and the detection of reliable metrics as indicators of water quality. Spearman's correlation and FA were computed using STATISTICA (version 7.0) for Windows.

Status of Water Quality in the Zio River
The descriptive analysis gives an overview on the variation of water quality parameters in the Zio River during the three sampling periods ( Table 2). EC expresses the degree of water mineralization and salinity. It varies in this study from 10,506 to 38.10 µS/cm, with high standard deviations (±1,755.07 µS/cm) indicating the large variability of EC from upstream to downstream of the Zio River. The pH ranged from 6.4 to 7.8, with the mean value varying from neutral toward a state of alkaline. The variation trends of most of major ions such as Ca 2+ , Mg 2+ , Na + , K + , HCO − 3 , Cl − , and SO 2− 4 are similar to EC and in the range of natural water quality. The concentration of minor ions such as NH + 4 , NO − 2 , Fe, PO 3− 4 , and Mn 2+ was below the detection limits in some sites and very low in other sites. The COD and the BOD are indicators of organic pollution in the water. COD values ranged from 4 to 30.7 mg/L with a mean of 8.68 mg/L, while BOD ranged from 2 to 20 mg/L with a mean value of 4.88 mg/L. The ratio COD/BOD is also an indicator of river pollution in terms of domestic or industrial effluents and biodegradability. This ratio ranges from 1.16 to 4.30 with a mean value of 1.92. The DO is one of the common parameters used to assess water quality and aquatic ecosystems' health. In the present study, the DO value varied from 0.6 to 12.8 mg/L. Many values recorded are high or close to 7 mg/L, showing that water in the Zio River can be qualified as good-to-excellent DO quality according to surface water standards except for a few polluted sites at the lower section of the river. The TSSs are defined as solids in the water including organic and inorganic materials that can be trapped by a filter. TSS ranged from 12 to 320 mg/L with a mean value and 75th percentile value of 108.6 and 158 mg/L, respectively.
Among microbiological parameters, total coliforms (TC) and fecal coliforms (FC) are often used as indicators of bacterial contamination in the water. Their presence is an indication of fecal contamination in water. In the present study, TC and FC range from 8 to 710 Cfu/100 ml and from 0 to 177 Cfu/100 ml, respectively. These values can express a risk of a bacterial contamination of water in the Zio River.

Sensitivity at Phylum-Class Level
The CCA indicates the environmental parameters that affect macroinvertebrate community at phylum-class level (Figure 2A). Table 3 shows the key indicator taxa and the total taxa identified at different taxonomic levels. The Supplementary Data in Appendix A indicate Spearman's correlation between water quality parameters and phylum-class abundance. From Figure 2A, phylum Annelida is affected by BOD, COD, and NH + 4 . Appendix A indicates that Annelida abundance is significantly and positively correlated with these parameters (r > 0.60; p < 0.05) and negatively correlated with DO (r < −70; p < 0.05). The significant correlation (r > 0.60; p < 0.05) of these parameters with TC and FC suggests that they can affect microbiological quality in the Zio River water. From Figure 2A, the class Crustacea is affected by EC and Cl and Na contents, which are core salinity parameters. Moreover, their abundance is significantly correlated with EC, Na, and Cl (r > 0.60; p < 0.05) and can suggest the degree of crustacean tolerance to the water salinity. From the results of CCA (Figure 2A) and the relationship (Spearman's correlation between water quality parameters and taxon abundance), this study identified at phylum-class level five taxa as indicators including three taxa as tolerant, one taxon as resistant, and one taxon as an indicator of water salinity (Table 3).

Response at Subclass-Order Level
The CCA indicates the environmental parameters that affect macroinvertebrate community at subclass-order level ( Figure 2B). The Supplementary Data in Appendix B indicate Spearman's correlation between water quality parameters and subclass-order composition. The orders such as Ephemeroptera, Plecoptera, Trichoptera, and Odonata (EPTO) are significantly influenced by DO with a significant and positive correlation (r > 0.60; p < 0.05). Their abundances are negatively influenced by BOD, COD, TSS, and NH4 (r < −0.60; p < 0.05). The subclasses such as Hirudinea (Leeches) and Oligochaeta are positively affected by BOD, TSS, and NH4 (r > 0.60; p < 0.05)  but negatively correlated with DO (r < −0.65; p < 0.05). The Pulmonates mollusks and Diptera are significantly correlated with FC and TC (r > 0.60; p < 0.05). From Figure 2B, EPTO are in an opposite trend in comparison with Oligochaeta, Hirudinea, Diptera, and Pulmonates (OHDP). Decapoda is significantly and positively correlated with EC, Na, K, and Cl and then positively correlated with the salinity. The other subclass-order such as Prosobranchia, Bivalvia, Lepidoptera, Heteroptera, and Coleoptera (PBLHC) do not reveal a significant association with water quality parameters. From Figure 2B and Spearman's correlation, at subclass-order resolution, four groups of macroinvertebrates can be distinguished as indicator taxa with different sensitivity to water quality parameters. The first group is composed of EPTO, which are sensitive to the decrease of DO, and to the increase of BOD, COD, TSS, and NH4. The second group is composed of OHDP, which has the opposite response in comparison with the first group. The third group is represented by Decapoda, which seems to be influenced by the salinity parameters (EC, Na, Cl, and K). The fourth group is composed of PBLHC, which does not show significant correlation with water quality parameters. The key indicator taxa and total number of taxa in each group at subclass-order resolution are mentioned in Table 3.

Response at Family Level
The CCA indicates the environmental parameters that affect macroinvertebrate community at the family level ( Figure 3A). The Supplementary Data in Appendix C indicate Spearman's correlation between water quality parameters and families composition. In Figure 3A, except some few families (Baetidae, Caenidae, Libellulidae, and Hydropsychidae), most of the families are from EPTO orders, and other families such as Atyidae from Decapoda, and Elmidae and Gyrinidae from Coleoptera are also strongly affected by water quality parameters. The abundance of these families increases with the increasing DO but decreases with the increasing BOD, COD, NH 4 , TC, and FC concentrations. Spearman's correlation confirmed this trend by a significant correlation between the abundance of these families and DO (r > 0.50; p < 0.05), BOD, COD, NH 4 , FC, and TC (r < −0.50; p < 0.05). Therefore, the families such as Atyidae, Elmidae, and Gyrinidae and those from EPTO can be classified as sensitive taxa (Table 3). In contrast, the abundance of Syrphidae, Psychodidae, Culicidae, and Chironomidae from Diptera order; Hirudinidae from Hirudinea; Lymnaeidae and Bulinidae from Pulmonates; and Tubificidae and Naididae from Oligochaeta increases when BOD, COD, NH 4 , TC, and FC concentrations increase and decreases when DO concentration increases. The abundances of these taxa are significantly correlated with DO (r < −0.50; p < 0.05) and with BOD, COD, NH 4 , TC, and FC (r > 0.50; p < 0.05). These relationships suggest that the families such as Syrphidae, Psychodidae, Culicidae, Chironomidae, Hirudinidae, Lymnaeidae, Bulinidae, Tubificidae, and Naididae can be classified as resistant or resilient taxa (Table 3). According to Figure 3A, the abundance of Palaemonidae is influenced positively by EC, Na, and Cl and is significantly and positively correlated with these parameters (r > 0.55; p < 0.05). This family from Decapoda can be considered as an indicator taxon of water salinity. The remaining families' taxa mainly from PBLHC; some families' taxa from Diptera and Decapoda; and also the few families from EPTO did not reveal a significant relationship with water quality variables (r < |0.50|; p < 0.05). These taxa were classified as tolerant or indifferent taxa as indicated in Table 3.

Response at Genus Level
The CCA indicates the environmental parameters that affect macroinvertebrate community at the genus level ( Figure 3B). From Figure 3B and Table 3, it is seen that at the genus level, there is increase of sensitive taxa number for each group. Most of the genera from EPTO; genera Atya (Atyidae) from Decapoda, Limnius and Elmis (Elmidae) and Orectogyrus and Aulonogyrus (Gyrinidae) from Coleoptera; and some genera from Heteroptera such as Ranatra (Nepidae) and Hydrometra (Hydrometridae) were identified as sensitive taxa (Table 3). However, some genera from Ephemeroptera such as Pseudocloeon, Cloeon (Baetidae), Caenodes (Caenidae), Eatonica, and Ephemera (Ephemeridae) and others from Odonata such as Urothemis, Zygonyx, Olpogastra, Pantala (Libellulidae), and Pseudagrion (Coenagriinidae) are found to be tolerant taxa in this study. Some genera from Diptera are also identified as tolerant taxa (Chironomus, Culicoides, Dixa, Simulium, Stenochironomus, Cryptochironomus, Atherix, etc.), while other genera from Diptera are identified as resistant or resilient taxa (Erystalis, Orthocladius, Tanytarsus, Tanypus, Cricotopus, Ablabesmyia, Procladius, Aedes, Culex, etc.). The genera from Pulmonate (Limnaea, Ceratophallus, Biomphalaria, Bulinus, and Ceratophallus) were observed as resilient or resistant taxa, while the genera from Prosobranchia were distributed among tolerant taxa (Potadoma, Cleopatra, and Melania) and resilient taxa (Lanistes, Pila, and Theodoxus). Most of the genera from Coleoptera and Heteroptera identified here belong to tolerant and resistant taxa, with more genera from Heteroptera found in the tolerant group ( Table 3). The genus Macrobrachium from Palaemonidae is positively influenced by salinity parameters and positively correlated with these parameters. In this study, it can be considered as an indicator taxon of water salinity ( Table 3).
The tolerant taxa were the most represented in terms of total number and at the different levels of identification. At the family level, 30 tolerant taxa, 24 sensitive taxa, 11 resistant taxa, and two taxa influenced by water salinity (salinity sensitive taxa) were identified. At the genus level, the same trend was observed, with 68 tolerant taxa, 53 sensitive taxa, 20 resistant taxa, and four salinity indicator taxa recorded ( Table 3).

Response of Macroinvertebrate Metrics
The CCA indicates the association between water quality parameters and metrics (Figure 4). Table 4 presents the correlation coefficients (factor loadings) between observed variables (raw water quality parameters and macroinvertebrate metrics) and the latent or underlying variables (factors). The Supplementary Data in Appendix D indicate Spearman's correlation between water quality parameters and metrics. From Figure 4, three groups of metrics can be distinguished according to their relationship with water quality parameters. First, the metrics for taxonomic richness are composed of ETOF, ETOG, EPTF, and EPTG. These metrics are positively and significantly correlated with DO (r > 0.65, p < 0.05), but they are negatively and significantly correlated with EC, BOD, COD, TC, FC, Cl, Na, and NH 4 (r < −0.60, p < 0.05). Then, the tolerant metrics are composed of ASPT and Biological Monitoring Working Party (BMWP) including MMIZB. These metrics are positively correlated with DO (r > 0.65, p < 0.05) and negatively correlated with EC, BOD, COD, TC, FC, and NH 4 (r < −0.55, p < 0.05). Finally, the diversity indices are represented by Marg_G, Marg_F, Sha_H_F, and Sha_H_G, as well as the total number of family and genera that reveal a strong positive association with DO (r ≤ 0.5, p < 0.05) and significant negative correlation with EC, TC, Cl, Na, and NH 4 (r ≥ −0.50, p < 0.05). From Table 4, the first factor is significantly and positively correlated with TSS, COD, BOD, TC, and NH 4 , but a significant negative association was observed with pH and DO. Therefore, the first factor can represent the water quality deterioration due to anthropogenic pressures. There is a positive association between the second factor and EC, Na, Cl, and K. This factor can represent the water salinity degree due mainly to mineralization processes or seawater. All the metrics tested in this study displayed a significant negative association with the first factor, but they did not reveal a significant association with the second factor. This result means that all metrics tested in this study are sensitive to water quality deterioration or anthropogenic disturbances affecting water quality.

DISCUSSIONS Water Quality Parameters and Their Interpretation
Knowledge about water quality parameters is an important part of environmental monitoring and determining the condition of habitats. When water quality is poor, it affects not only aquatic life but also the aquatic ecosystem health. This section details with some parameters that affect water quality and aquatic ecosystem health. The values of EC show that more than 75% of water samples were in the range of natural freshwater, which varies from 0.5 to 1,500 µS/cm (Rodier et al., 2009).
The high values of EC can affect freshwater organisms such as macroinvertebrate communities (Environmental Protection Agency of Ireland, 2001) as found in the present study; it is revealed that EC affects Palaemonidae taxon particularly the genus Macrobrachium. As a result, this taxon can detect the salinity or EC in the water of the Zio River. The pH has been considered as an important parameter in the ecology of aquatic macroinvertebrates (Thomsen and Friberg, 2002;Yuan, 2004). Benthic macroinvertebrates are sensitive to pH variation, and values below 5 or >9 are considered harmful (Yuan, 2004). The pH value found in the present study is in the range of natural water pH 6.5-8.5 according to some standards [U.S. Environmental Protection Agency (USEPA), 1980; Environmental Protection Agency of Ireland, 2001;Rodier et al., 2009;WHO, 2011]. However, according to Thomsen and Friberg (2002), the low pH values are associated with lower diversity of benthic macroinvertebrates and cause decreasing emergence rates of this community (Hall et al., 1980). The range of pH found in the present study is associated with high diversity of some macroinvertebrates such as EPTO taxa. The DO is an excellent indicator of water quality and one of the most sensitive to water quality degradation and therefore indicates the degree of water body self-purification (WHO, 1996;Sanchez et al., 2006;Makhoukh et al., 2011), whereas BOD indicates the quantity of biodegradable substances that mainly originate from organic matters into the water body due to human activities. In this study, the mean DO value observed in the range of saturated water (6-7 mg/L) with the 75th percentile value around 9.75 mg/L can indicate that most of the water samples can be classified from good to excellent quality according to international standards (Environmental Protection Agency of Ireland, 2001;Rodier et al., 2009). In the same vein, our study showed that DO and pH affect EPTO taxa found in the present study as sensitive taxa. Therefore, EPTO taxa or related metrics are able to detect variables such as DO and pH in the water of the Zio River. The BOD value increases when DO value decreases and often indicates an organic pollution. We found that in the pristine rivers, this value is below 2 mg/L; and in moderately polluted rivers, it can range from 2 to 8 mg/L, although the rivers above 8 mg/L may be considered as severely polluted (Rejsek, 2002;Rodier et al., 2009). The COD is an indicative measurement of oxygen consumed by reducing substances (organic, nitrite, sulfide, ferrous salts, and others) in the water. It is also an important parameter of water quality assessment (Rejsek, 2002).
The 75th percentile values and the mean values of BOD and COD indicate that the waters of this study can be considered as moderately polluted. The ratio COD/BOD is a meaningful parameter used to assess biodegradation in natural water. It is also used in environmental monitoring to assess the selfpurification power of a river, which is the suitability for a river containing organic matter to transform this organic matter into mineral matter by natural processes (Mara and Pearson, 1999;Mara, 2004). A ratio of COD/BOD ranging from 1.5 to 2.5 indicates a good biodegradability and self-purification power (Mara et al., 2007;Rodier et al., 2009). Accordingly, the Zio River can be considered as a suitable river for selfpurification process.
The TSS in water may cause pathogen problems, as the growth of anaerobic bacteria attached to suspended material can increase the risk of disease outbreaks and a significant portion of toxic organics (Rosewarne et al., 2014;Rono, 2017). In addition, the trace elements can be adsorbed onto organic matter present in water as suspended materials (Rono, 2017;Helmecke et al., 2020). TSS can also cause ecosystem disturbance through reduced water clarity, limiting photosynthesis and asphyxiation of some gill-breathing organism by clogging gills (Rodier et al., 2009;Rosewarne et al., 2014;Swinkels et al., 2014). The high levels of suspended solids can be considered as a form of pollution, which will have the effect of reducing the quality of the habitat for cold-water organisms (Rejsek, 2002;Rodier et al., 2009). All sampling sites within the Zio River contain acceptable concentration of TSS for aquatic ecosystem functioning life. Ammonium usually reflects a process of anaerobic degradation of nitrogenous organic matter and can be a good indicator of river pollution by urban effluents (Chapman and Kimstach, 1996). The range of NH + 4 values is near to natural water, but the standard deviation indicates a slight punctual source pollution of the study river. The present study revealed that COD, BOD, TSS, and NH + 4 affect OHDP taxa found here as resistant taxa. Thus, the findings from this study suggest that OHDP taxa or related metrics are able to detect COD, BOD, TSS, and NH + 4 in the water of the Zio River. Chloride is frequently associated with wastewater; and it is often included in assessments as an indicator of possible fecal contamination or as a measure of the amount of dispersion of wastewater discharges into the natural environment (WHO, 1996). In pristine freshwater, chloride concentrations are generally less than 10 mg/L and sometimes less than 2 mg/L (WHO, 1996;Rodier et al., 2009). Sodium concentration is often related to the nature of contact rock, evaporation, and seawater intrusion phenomenon. In this latter case, sodium content is often correlated with chloride content. The 75th percentile and mean values of sodium and chloride indicate that about 75% of water samples were not affected by sea intrusion, while 25% were affected as reported by Tampo (2018) at the lower reaches. TC and FC are parameters that indicate the possible fecal contamination of water. The recommended values in freshwater should be below 10 Cfu/100 ml. The values recorded here reveal a risk of fecal contamination at some sites. The Zio River water quality trend follows a gradient of pristine-to-poor water quality from the upper reaches to the lower reaches with a mosaic of good and poor water quality in the lower reaches based on proximity to anthropogenic pressure and saltwater intrusion (Tampo, 2018).

Response of Macroinvertebrates at Different Taxonomic Levels
In the present study, some phylum-class taxa could effectively indicate the ecological integrity associated with water quality. This could be explained by the fact that these phylum-classes were dominated by taxa with the same or very few different lifestyles that could point out a single category of pollution (Tachet et al., 2010). From this study, crustaceans were found to be good indicators of the freshwater salinity. Indeed, crustaceans were dominated by two families (Palaemonidae and Atyidae) and well known to be influenced by EC or water salinity (Attrill et al., 1999;Cuesta et al., 2006;N'Zi et al., 2008;Collocott et al., 2014). At the genus level, it was found that Macrobrachium from Palaemonidae was the most affected by salinity gradients and can be considered as the key indicator of water salinity. This corroborates the findings of other studies, which indicate that the life cycle or style of Macrobrachium is mainly related to estuary or water salinity (N'Zi et al., 2008;Collocott et al., 2014;Gangbe et al., 2016;Chen et al., 2018;Adam et al., 2019). Meanwhile, Mollusca taxon is associated with water eutrophication due to organic pollution. The findings revealed that phylum Annelida could be considered as generalists with high tolerance to stream degradation. At subclass-order level, the identified group of EPTO was found to be a sensitive group due to high number of sensitive taxa belonging to this group. This sensitivity of EPTO taxa is well known and described by earlier authors (Hilsenhof, 1988;Rodier et al., 2009;Tachet et al., 2010;Ko et al., 2020). The high abundance and diversity of EPTO taxa indicate good or excellent water quality with good ecological conditions (Tampo et al., 2015(Tampo et al., , 2020Kaboré, 2016). Many metrics and multimetric indices calculated through EPTO taxa are widely used in biomonitoring programs as reported by many authors (Rodier et al., 2009;Mereta et al., 2012Mereta et al., , 2013Tampo et al., 2015Tampo et al., , 2020Kaboré et al., 2016). The results support other studies that found that the high abundance of Oligochaeta, Hirudinea, Pulmonates, and some Diptera often indicates organic pollution as seen in the present study. Globally, this has been seen with possible fecal contamination and substantial anthropogenic pressures (Hilsenhof, 1987;Rodier et al., 2009;Tachet et al., 2010;Tampo et al., 2020). This resistance to pollution of some macroinvertebrate groups can be explained by their lifestyle and the feeding strategy. For example, the Oligochaeta feed on detritus, bacteria, and dead remains of other aquatic organisms. Pulmonates are not only detritivores but also the intermediate host of many parasites (Tachet et al., 2010). According to Rodier et al. (2009), Oligochaeta and Leeches are indicators of polluted water and deteriorated habitat. Their number increase with the increase of eutrophication, which affects the overall ecosystem health (Kazanci et al., 2015). Within the insect fauna, Hilsenhof (1987) has demonstrated that Libellulidae, Lestidae, and Coenagrionidae from Odonata order are tolerant taxa. Furthermore, Rodier et al. (2009) showed that in a biomonitoring program based on a standardized global biotic index, Baetidae (algae grazer) and Caenidae from Ephemeroptera are resistant taxa, while Mary (1999) has reported some sensitive taxa in Heteroptera families. Here, we found within the EPTO order a notable diversity, with several families observed in good habitat conditions. This could be explained by the presence of families less tolerant or even sensitive to pollution. At the genus level, more taxa emerged as sensitive with specific indications. Some studies reported this increase of the sensitivity degree and the specific indication of macroinvertebrates when moving from phylum-class to species taxonomic level (Bailey et al., 2001;Jones, 2008;Menezes et al., 2011). However, according to the review paper of Jones (2008), many experts showed that the family level is sufficient or may be better in bioassessments even if genus and species taxonomic levels are required. The present study agrees with this consensus and shows that the family level is sufficient and better in the detection of nutrient pollution, organic pollution, and human disturbances and which pressures may drive shifts when using indices or metrics in macroinvertebrate communities. In addition to this, the present study showed that subclass and order levels are less sufficient of specific pollution but better for rapid evaluation of pollution at the global scale. These findings also agreed with those of other studies carried out in Africa (Dickens and Graham, 2002;Mereta et al., 2013;Elias et al., 2014;Kaboré et al., 2016;Dalu and Chauke, 2020;Ochieng et al., 2020;Tampo et al., 2020;Edegbene et al., 2021) in Asia (Blakely et al., 2014;Nguyen et al., 2014) and Europe (Poquet et al., 2009;Rodier et al., 2009;Costas et al., 2018) where these studies suggested the use of macroinvertebrates at the family level for the development of biomonitoring programs. The findings of the present study highlight the importance of using the family level for biomonitoring program in Togo and even elsewhere in Africa because of cost-effectiveness and the lack of systematic knowledge on macroinvertebrates.

Macroinvertebrate-Based Metrics
According to the results of the ordination (CCA and FA) and Spearman's correlation analysis, the metrics calculated at the family level and those calculated at the genus level were strongly correlated (r ≥ 0.80) and with the first factor (r ≥ |0.75|). This strong correlation is a redundant association between metrics calculated at family and genus levels. Thus, these two taxonomic levels can translate the same or common information. Our results also indicate that the taxonomic detail (family or genus level) does not substantially affect metrics used in bioassessment. These results are near the consensus that the lowest taxonomic level (species) is not always required for bioassessments (Jones, 2008). For example, in the case of differentiating impacted from unimpacted sites, some researchers showed that species-level identification was not always the best, but total richness and EPT richness performed better at family than at species, and family-level diversity measures also performed well (Schmidt-Kloiber and Nijboer, 2004;Schmidt-Kloiber et al., 2006). The negative correlation observed in the present study with all metrics confirmed the fact that they appear to be highly sensitive in the variation of water quality variables linked to the degree of human pressures. This sensitivity of macroinvertebrate metrics to the water quality and anthropogenic pressures is widely reported (Hering et al., 2006;Jun et al., 2012;Helson and Williams, 2013;Lakew and Moog, 2015;Kaboré et al., 2016;Aura et al., 2017;Tampo et al., 2020). The degree and the type of sensitivity of each one can explain the three groups of metrics identified by the CCA. In this way, reporting the sensitivity of EPTO group, many authors found that the diversity and abundance of this group increase with the increase of DO and with the decrease of nutrients, BOD, COD, and TSS (Mereta et al., 2013;Tampo et al., 2015;Kaboré et al., 2016). Thus, the EPTO group is the most suitable in the evaluation of water quality and aquatic ecosystems health and is used worldwide in bioassessments and monitoring programs (Armitage et al., 1983;Hilsenhof, 1987;Barbour et al., 1996Barbour et al., , 1999Schmidt-Kloiber et al., 2006;Rodier et al., 2009). The metrics related to tolerance measure (MMIZB; ASPT and BMWP) are also highly correlated with the first factor and some water quality parameters. This result revealed their sensitivity in water quality variation and mostly in the detection of anthropogenic pressures. Furthermore, in the previous studies, MMIZB revealed its sensitivity by discriminating impaired sites and unimpaired sites (Tampo et al., 2020). The two metrics ASPT and BMWP were also used in many studies as a single metric or integrated in a multimetric index for the monitoring of watersheds (Ferreira et al., 2011) for assessing water quality (Nguyen et al., 2014) and for assessing ecological integrity (Solimini et al., 2007;Lakew and Moog, 2015). That is likely because of their robust sensitivity and high discrimination power; the tolerance measure metrics are mostly recommended for ecological status evaluation under many conditions (Barbour et al., 1996;Hering et al., 2006). The total taxonomic number and the diversity indices were also calculated at family and genus levels. The results show that for the same metric of this group calculated, using familylevel and genera-level data did not reveal a big difference about their relationship with water quality parameters and the first factor. Indeed, in bioassessments, richness is used to summarize biological condition, not to generate exhaustive taxon checklists, so coarse taxonomic resolution is often sufficient. For example, Marshall et al. (2006) reported only a 6% information loss when benthos data were rolled up to family from species. These results agree with those of several authors who found that correlations between ordination site-scores and environmental variables are slightly affected by taxonomic detail (Hewlett, 2000;Waite et al., 2000;Metzeling et al., 2006). In the same vein, Metzeling et al. (2006) reported that reducing taxonomic detail from species to family had little effect on ordination plots and gave similar principal axis correlations with environmental variables. Stronger evidence was provided by Hewlett's (2000) classification and ordination-based study of Australian stream sites. Even if the different taxonomic resolutions are known to influence the responses of macroinvertebrates through several mechanisms, our findings suggest that the use of macroinvertebrate metrics approach based on the family level has proven useful for Zio River's aquatic ecosystem monitoring.

CONCLUSION
In many developing countries across Africa such as Togo, rivers are threatened by intense agriculture (using fertilizers and pesticides), urbanization, and severe pollution, leading to habitat loss and water quality and negatively affecting aquatic organisms. With ongoing multiple pressures, the urgent need of biomonitoring tools is crucial for ensuring valuable biodiversity and water resource conservation. Along the Zio River gradient including several types of pressures, it has been demonstrated that macroinvertebrate indices or metrics have proven their sensitivity to water quality variation and human disturbance from order-level identification up to the genus level. Specifically, those related to EPTO taxa and tolerance were suitable in detection of water quality and human disturbances. The findings of this study confirm the importance of maintaining the family level for bioassessment and biomonitoring programs development in developing countries, especially in Sub-Saharan Africa, due to its cost-effectiveness and the lack of systematic knowledge on macroinvertebrates. However, the genus and species taxonomic levels are needed to improve the understanding of responses on the family level and the detection of specific pollution. We, therefore, recommend the use of the family level in metric and index formulation to monitor the Zio River.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Materials, further inquiries can be directed to the corresponding author.

ETHICS STATEMENT
The animal study was reviewed and approved by Wendengoudi Guenda, University of Joseph KI-ZERBO, Laboratory of Animals Biology and Ecology, 03 BP 7021 Ouagadougou 03, Burkina Faso.

AUTHOR CONTRIBUTIONS
LT conceived the manuscript and wrote the paper with IK and EHA. LT and IK collected and analyzed the data. AO assisted LT in the laboratory analyses. LMB and GD-B supervised the laboratory analyses. EHA corrected the English language. LMB and GD-B supervised the study and contributed valuable comments to the manuscript. All authors contributed to the article and approved the submitted version.