Comparing Freshwater Benthic Macroinvertebrate Communities in Forest and Urban Streams of the Coastal Ecological Region of Cameroon

In this study conducted in the coastal zone of Cameroon, biological indices and functional feeding groups of benthic macroinvertebrates were used to assess the health status of two urban streams. For a better diagnosis, two streams located in coastal forest zone were used as a reference. Benthic macroinvertebrates were sampled monthly over a 3-month period (from May to July 2017) in six urban stations and six forest stations. Measurements of the physicochemical variables were done simultaneously. Physicochemical analysis revealed that urban streams are strongly polluted with high content of decaying organic matters, while forest streams are slightly polluted as indicated by the Principal Component Analysis. Concerning benthic macroinvertebrates, urban streams are poorly diversified with the proliferation of taxa tolerant to water pollution and belonging to the functional feeding groups of collectors-gatherers. Inversely, forest streams are more diversified and dominated by sensitive taxa, most belonging to the functional feeding groups of predators and shredders. These marked differences between biological indices and feeding


Introduction
Assessment of the impacts of urbanization on rivers and streams involves the research and use of a wide range of physical, chemical, and biological indicators [1]. However, [2] recommend the use of biological communities at the expense of physicochemical variables which provide a punctual approximation of the quality of the environment. The benthic macroinvertebrates are metazoans that are widely represented in rivers and streams and constitute a good tool for this diagnosis [3]. Their preference for biological monitoring of rivers is justified for to the following characteristics [2]: 1) they occupy a wide range of trophic levels and pollution tolerance, providing solid information for the interpretation of cumulative effects due to exposure to pollutants; 2) their sedentary life style make it possible to reflect the local conditions of the environment, and to give a significant image of the quality of the habitat over time.
In Cameroon, the use of benthic macroinvertebrates as bioindicators of the health status of rivers has expanded considerably in the last 15 years. However, most of the works on this group of organisms has been conducted in urban and forest streams in the South-Central Forest Ecological region with bimodal rainfall [4] [5] [6]. The data available for the coastal eco-zone with unimodal rainfall relate to only a few urban and peri-urban rivers in the industrial-port city of Douala [7] [8] [9]. To date and to our knowledge, no study has been conducted on the ecology of benthic macroinvertebrates of coastal forest streams in the unimodal rainfall ecological region. However, in order to understand the functioning of a biological system, it must be understood and interpreted according to its particular ecological situation and implicitly implies a reference model adapted to the ecological region considered [10].
The main objective of this study is to assess the impact of urban pollution during the rainy season on the diversity of benthic macroinvertebrates communities of two streams of the industrial-port city of Douala, as compared to the referential status of two other streams chosen in a natural forested river basin at Yabassi. The specific objectives of this study are: 1) to evaluate the physicochemical quality of the waters of the forest and urban streams; 2) to determine the variation in the structure of benthic macroinvertebrates communities of naturally vegetated watersheds and urban watersheds using the biological indices and functional feeding groups; 3) to determine the main physicochemical variables that influence the distribution and the dynamic of benthic macroinvertebrates in the studied streams.

Study Area and Sampling Stations
The present study was conducted in the littoral region of Cameroon, which is subject to a unimodal equatorial climate of Guinean type and Cameroonian coastal sub-type, characterized by two seasons: a long rainy season of 9 months (March-November) and a shot dry season of 3 months (December to February) [11]. Rainfalls are abundant and regular y (2596 -5328 mm). The air temperature is relatively high with a monthly average of 28˚C. The streams chosen for study are situated in the localities of Douala and Yabassi.    (Figure 1(b)). The Mboppi river basin is characterized by

Physicochemical Analyzes
The measurements of the physicochemical parameters at each of the twelve sampling stations were done monthly, between 7 and 10 am, from May to July 2017.
These measurements were conducted in accordance with standard methods described by [12]. Thus, a total of 36 measurements were carried out for each pa-  Table 2). The OPI corresponds to the arithmetic mean value of the class number of each parameter and the values obtained are apportioned into 5 levels of pollution (Table 3).

Sampling and Identification of Benthic Macroinvertebrates
Benthic macroinvertebrates samples were collected at each station using a longhandled kick net (30 cm × 30 cm side, 400-µm mesh size and 45 cm of depth).
All the available habitats within the study station including leaf packs, branches, macrophyte beds, sand, mud and rafts were sampled. Samplings were done in a 100-m stretch for each station, following protocols described by [14]. At each sampling station, 20 drags of the kick net were done in different micro-habitats, each corresponding to a surface area of 0.3 m 2 (30 cm × 100 cm), for a total sampled surface area of 6 m 2 per station. Each time, the benthic macroinvertebrates retained by the meshes of the net were sorted and collected using a pair of fine forceps, and then introduced into plastic bottles containing 70% alcohol. In the laboratory, all of the benthic macroinvertebrates at each sampling station were identified to the family taxonomic level and counted, using a Wild M5 stereomicroscope and appropriate taxonomic identification keys proposed by [15] [16] [17]. Recent research indicates that extensive identification up to the family provides a sufficient level of taxonomic resolution to detect responses of the communities studied to anthropogenic disturbances [18].

Benthic Macroinvertebrates Metrics and Indices
The diversity index of Shannon and Weaver (H') and the Pielou's evenness index (J) were calculated from the taxonomic richness (S) and abundance (TA) of macroinvertebrates. These indices informed on the distribution of individuals within the taxa so as to compare the diversity of the communities of the different sampling stations [19]. The assumption is that undisturbed environments are characterized by high diversity or richness and an even distribution of individuals among the Taxa [20].
The adaptation of benthic macroinvertebrates to organic pollution was assessed through six metrics: relative abundance of taxa tolerant to pollution (%-To; tolerance score > 6), relative abundance of taxa sensitive to pollution Open Journal of Ecology taxa [21]. The tolerance score assigned to each organism in this study was taken from [22] [23].
Lastly, we assigned each taxon to one of the five functional feeding groups (FFGs): shredders (shr), scrapers (scr), collectors-gatherers (c-g), collectors-filterers (c-f) and predators (prd); based on the food habit descriptions in literature [16]- [24]. In fact, functional approach which the goal is to characterize ecosystem condition uses FFGs to describe ecosystem attributes in streams [24].

Physicochemical Quality of Water
Except for orthophosphates, the values of pH, temperature, electrical conductivity, color, Suspended Solids, turbidity, ammonia and nitrite were significantly lower (P < 0.01) at the forest stations compared those recorded in the urban stations. Inversely, the highest values of dissolved oxygen were recorded at the forest stations, while the urban stations showed very high oxygen content (Table   4). Moreover, the organic pollution level of the water was moderate at the forest sites (OPI = 3.41 ± 0.54) whereas in the urban area it was high (OPI = 2.11 ± 0.58) ( Table 4).
The PCA carried out with the physicochemical parameters data set shows that the first two axes F1 (58.36%) and F2 (18.21%) accounted for 76.57% of the total variance ( Figure 2). The distribution of the 36 water samples on factorial plane of the PCA allowed to discriminate the sampling sites into two groups ( Figure   2   Mann Whitney tests were used to evaluate differences among the two groups of streams (forest streams and urban streams). The values followed by the superscript (a, b) for a given physicochemical parameter indicate that this parameter is significantly different between forest and urban streams (P < 0.01).
from the Simbi and Mboppi streams, located in Douala. These urban sampling sites have very poor water quality with higher values of water color, turbid, ions content, and important organic matter load.

Benthic Macroinvertebrates Communities' Structure
Considering spatial distribution of the taxonomic richness, 36 families divided into 4 classes (Crustacea Gastropoda, Insecta, Oligochaeta) and 12 orders were   Mann Whitney tests were used to evaluate differences among the two groups of streams (forest streams and urban streams). The values followed by the superscript (a, b) for a given metric indicate that this metric is significantly different between forest and urban streams (P < 0.01). account for 82.53% of the total abundance in urban streams. In the forest streams,  (Table 5). Although indicating, a better distribution of individuals within different taxa in forest (0.75 ± 0.08) and urban (0.69 ± 0.19) stations, the Pielou evenness index (J) did not vary significantly (P > 0.05) ( Table 5). The %-Into and FMWP indices measuring sensitivity of taxa to pollution were higher in forest streams than in urban ones. For %-To and FBI, the values recorded in urban streams are significantly high (P < 0.05) than those registered in forested sites (Table 5). In addition, the FBI revealed that forest streams waters are of medium quality with fairly substantial organic pollution (FBI = 5.43 ± 0.71) whereas waters from urban streams are of very poor quality with severe organic pollution (FBI = 7.54 ± 0.37).
Regarding FFGs, in the forest streams, predators and shredders are the most represented with 38.04% and 30.02% of relative abundances respectively ( Figure   3). The collectors-gatherers and the scrapers represent 16.18% and 15.40% of the total abundance respectively. The collectors-filterers are very poorly represented with 0.36% of individuals. In urban streams, macroinvertebrates are dominated by collectors-gatherers, which constitute 82.15% of the population. However, we also find collectors-filterers, predators and shredders which respectively account for 12.85%, 4.80% and 0.20% of macroinvertebrates.

Relationship between Benthic Macroinvertebrates Metrics and Physicochemical Water Quality
The PCA based on the physicochemical variables and the benthic macroinvertebrates metrics showed that the first two axes (F1, F2) account for 67.60% of the total variance ( Figure 4). The distribution of the observations according to the factorial plane 1 -2 discriminates two groups. The group I is consisted of physicochemical parameters such as nitrites, ammonium, orthophosphates, color, SS, turbidity, pH and temperature, which are positively and significantly correlated to %-Dip, %-Oligo, %-To, FBI and ASPT. It may be thought that benthic  macroinvertebrates individuals belonging to the orders of Diptera, Oligochaeta and all other individuals belonging to pollution-tolerant families (tolerance score > 6) proliferate in very turbid waters, with high temperature and heavy organic matter load, characteristics of urban streams located in Douala. This is confirmed by the positive and significant correlations between metrics such as %-Dip, %-Oligo, %-To, FBI and ASPT and environmental variables such as temperature, electrical conductivity, color, turbidity, SS, nitrites, ammonium and orthophosphates (Table 6). Inversely, the group II, gathered metrics such as %-Ins-N-Dip, %-Into, %-EPT, S, S-EPT and S-POET which are positively and significantly correlated with dissolved oxygen and OPI (Figure 4), indicating low organic matter load in the forest streams of Yabassi. Similarly, the Spearman correlation matrix revealed that the installation and the proliferation of non-Dipterans insects, EPT, POET and all other sensitive individuals (tolerance score < 4) are favored in well oxygenated water with low organic matter concentration (Table 6). Mboppi market, Sandaga market), households and traditional toilets of the spontaneous habitats. The resulting aerobic degradation of organic matter is responsible of the hypoxic condition of these waters [25]. As regards forest streams, they showed a better water quality due to the presence of riparian vegetation. In fact, the important riparian forest limits the exposure of the streams to solar rays, which considerably reduces the temperature of the water. It also reduces soils erosion and indirectly the amount of nutrients and organic matter loaded in the streams [26]. The values of parameters indicating organic matter load (nitrites, ammonium and phosphates) are lower in forest streams; however, these parameters revealed through the OPI, a moderate organic pollution of water. [27] has attributed such a situation in the forest zone, principally due to the degradation of the litter coming from the riparian vegetation.

Discussion
Concerning forest streams of Yabassi, there was a numerical preponderance and taxonomic richness of sensitive taxa which includes Decapoda (Atiydae, Desmocarididae) and non-Dipterans Insects (Leptophlebiidae, Baetidae, Potamanthidae, Ephemerellidae, Perlidae, Hydropsychidae, Veliidae, Gomphidae, Gyrinidae, Noteidae, Hydrophilidae, Chrysomelidae). Their preponderance would be related to the good oxygenation of the forest streams. The positive and significant correlations between S-EPT, S-POET, %-EPT, %-POET, %-Ins-N-Dip and the high concentrations of dissolved oxygen in forest streams confirm this observation. This numerical and taxonomic preponderance of Decapoda and non-Dipterans Insects found in the forest streams as compared to their low richness and abundance in the urban streams could reveal a state of advanced degradation of the Mboppi and Simbi stations. [28] showed that taxonomic richness particularly that of aquatic insects can be a good indicator of the influence of anthropogenic disturbances on rivers. Indeed, [16] indicate that aquatic In a study conducted in the United States, [29] indicate that EPT taxa are among the most sensitive benthic macroinvertebrates and therefore their species richness decreases drastically with the urbanization of watersheds. Our results are in line with those of [30] who showed that a decrease in the number of EPT taxonomic richness in human impacted streams is clearly related to changes in water quality and habitat suitability. For freshwater Decapoda, it is recognized that these organisms proliferate in well-oxygenated waters [16] [31]. In addition, the great prominence of the Decapoda order would also be related to their belonging to the functional feeding groups of shredders and to the presence of large amounts of litter in forest streams. They find an important part of their food resources in the litter. [32] showed that macroinvertebrates associated with litter in Guinea's rivers were largely dominated by shredder Decapoda.
The proliferation of Physidae, Lymnaeidae and Chironomidae in the urban streams of the city of Douala is explained by the fact that these waters systems are subjected to heavy organic matter load. To live in such disturbed environments, these organisms have developed special adaptations. Indeed, [33] showed that Physidae and Lymnaeidae have a rudimentary lung through which they breathe atmospheric air and thus can live in hypoxic aquatic environments.
Moreover, their rapid maturation and their ability to self-fertilize in difficult conditions are also at the origin of their important development [34]. As for Chironomidae, studies by [35] [36] have shown that these organisms have adapted to hypoxic media with the hemoglobin in their hemolymph which allows them to store oxygen in the cuticle when it is abundant in the medium and to release it when its concentration becomes low. Furthermore, Chironomidae have what appears to be finger-like "gills" at the posterior ends of their bodies, which are actually involved in osmoregulation and permitted them to live in polluted waters with high ionic concentrations [8]. In addition, the high abundance of Physisae, Lymnaeidae and Chironomidae is related to their feeding mode; indeed, they are collectors-gatherers and proliferate in turbid aquatic milieu with important decaying organic matters load [24]. This is primarily due to various anthropogenic pollution sources found in the rivers basins of Douala township.

Conclusion
In this study, the set of benthic macroinvertebrates metrics and FFGs used to