Study Area and sampling design
The study took place in the Nyanza Gulf on the Kenyan side of Lake Victoria (Fig. 1). This catchment area contributes 30% of the riverine inflow into the lake, sourced from five major rivers: Nzoia, Yala, Kuja, Nyando, and Sondu, all originating in the highlands and passing through agricultural areas (Gikuma-Njiru & Hecky, 2005; Gikuma-Njiru et al., 2013). Additionally, smaller crucial rivers such as R. Kisat, R. Kisian, and R. Nyamasaria (a highly polluted river receiving influents from the City’s sewage treatment plant) flow from agricultural and suburban areas, respectively. Sondu and Miriu rivers are utilized for hydropower generation. These rivers traverse sugarcane agricultural and industrial zones, as well as rice paddies where rice cultivation is integrated with aquaculture, albeit on a small scale. The region, however, faces challenges due to the impact of industries like breweries, tanning, fish processing, agro-processing, and abattoirs, which discharge significant volumes of pollutants into the lake. Moreover, the area is known for activities such as livestock farming, gold mining, sand harvesting, and oil depots.
Nyanza gulf, is densely populated at the rate of about 440 people per km2 (KPHC, 2009), putting pressure on available land resulting to deforestation and increasing levels of pollutants in the gulf. Major towns along the lake shores, which include: Kisumu City, Homa Bay and Kendu Bay are particularly under greater pressure in terms of sewage facilities and refuse deposition which translates into nonpoint sources of pollution on lake water quality. The effluents through storm waters find their way to the lake, compromising on the water quality. The pollutants impact on the gulf ecosystem, causing destruction of habitats for aquatic flora and fauna. The zonation of lake ecosystem forming growth of emergent and submergent macrophytes in the littoral zone and the water column, create insect habitats, which form part of the aquatic food chains and food web, allowing energy flow up the trophic levels.
Sampling was done in September 8th -11th, 2020. Three sampling stations, Kisumu bay, Homabay and Kendu-bay located in urban environs and with an approximate maximum depth of 3.5m were identified to capture the effects of urban pollution. This was to include effects of domestic effluents, raw sewage, industrial wastes and the nonpoint sources like storm water from the towns. The water contains suspended and dissolved solids like oils, detergents from carwash, wasted batteries etc. Kisumu bay has a highly polluted river in-let Kisat, which is mainly sewage effluent. Ndere island was identified to represent the inshore stations in the gulf and presumably observed to be relatively clean, with limited human activities due to the presence of Ndere island National Reserve Park under the management of Kenya wildlife service, KWS.
A systematic random design was employed in the sampling. Approximately 50 m belt along the lake shores was estimated and a first point randomly located at the center. A transect was developed across and three points identified for random sampling of water and sediments in triplicates. The samples were pooled to obtain composite homogeneous sample for heavy metal analysis. Sediment samples were collected from the three sampling points per sampling station for heavy metal analysis. Insect samples were collected from three sampling points on the transects, pooled together to form a representative sample per stations.
Water samples for heavy metal analyses were collected using a polycarbonate sampling bar of polytetrafluoroethylene container which was thoroughly pre-washed with acid and deionized water. The samples were preserved in liquid HNO3 at a pH of < 2, placed in cooler boxes at 4ºC and transported to the laboratory for analysis.
The U.S. EPA (2020) standard operating procedure (ID: LSASDROC-200-R4) and (Ohio EPA, 2001) was employed in collection of sediment samples. A 400cm2 (20 x 20 cm) size AISI 316 stainless steel fabricated Eckman bottom corer grab was used. The spring-tensioned, scoop jaw-like part was mounted on pivot points and set with a trigger assembly activated from the surface by a messenger. A total of three grabs was made per sample, dumped on a tray and pooled to form a composite sample. The sample was then placed in wide bottom 1000 ml Polyvinyl chloride, PVC bottles and covered with aluminum foil paper, placed in a cooler box at 4ºC and transported to the laboratory for heavy metal analysis.
Analysis of insect samples was done using a profundal lake sampling procedure as outlined in the standard SFS 5076. A boat was used to reach a 50 m distance inshore for littoral zone sampling. At the anchoring site, Eckman grab–Birge dredge sampler was used to make random triplicate grabs of submerged insect larvae placed into the plastic bucket through a bucket sieve with um mesh and pooled to obtain a composite sample. The contents emptied for sorting aided by washing bottles to flash the remaining content using alcohol, then placed in paper slips and labeled (location, date, time, collector, sampling method, habitat, habitat description, weather and photographs of every site taken, sample number). These were then filled with 80% alcohol as in ISO-EN 5667-3, closed and packed in readiness for transportation in cooler boxes at 20ºC for further identification (Moller, 2013; Mehmet et al., 2015; McMurtrie et al., 2014). Sampling was carried out in the morning hours between 7 am-11.30 am.
Sample Analyses (water, sediment and insect)
Heavy metal analyses in water were done by Vacuum filtration with high particulate matter, high turbidity, suspended solids and with low levels in concentration. Pre-treatment of the water sample was carried out by adding a mixture of HNO3:H2O2, in a ratio 1:3 for acid digestion. The contents were placed in a hollow solid fiber micro-extractant with graphene oxide silica coating. The extract was raised to the 100 ml mark in readiness for analysis by Inductively Coupled Plasma-Mass Spectrometry, ICP-MS technique (EPA Methods, 2008). Mercury analysis was done using the Cold Vapor Atomic Fluorescence Spectroscopy, CV-AFS (EPA Methods, 245.1).
The sediment samples were dried at room temperature, ground and sieved using 2 mm sieve. About 2.00 g obtained from the dried sample was placed in a 100 ml beaker, 15 ml HNO3 added and heated at 130 ºC till boiling point for 5 hrs. to obtain 3 ml. The heated content was filtered, washed in 0.1M HNO3, made up to a volume of 100 ml using distilled water in readiness for analysis using Graphite Furnace Atomic Absorption Spectrophotometer - GF-AAS (EPA Methods, 200.9) for metals except Mercury and Arsenic, analyzed by Cold Vapor Absorption technique-CV- AFS, (EPA methods, 245.1)
Heavy metal analyses on preserved insect samples were done by oven drying at 120ºC and weighed on a Mettler Toledo microbalance, a highly precise weighing instrument used for determining the weight of extremely small samples in the microgram weighing range. Microbalances offer six decimal place readability (1 µg). Triplicate samples were digested in a flask by adding HNO3 and H2O2, in the ratio 1:1, and heated at 130ºC till dissolution. Dilution was done by adding distilled water to obtain 100ml sub-samples. Then followed by metal analysis using Graphite furnace- Atomic Absorption spectrophotometer-GF-AAS (EPA Method 200.9) for metals except Mercury and Arsenic, analyzed by Cold Vapor Absorption technique-CV- AFS, (EPA method 245.1)
Molecular analysis of chironomids
Genomic DNA was extracted from insects’ isolates using the DNeasy Blood and Tissue Kit Qiagen™ kit according to the manufacturer’s specifications. The concentration and purity of DNA was estimated using a DNA Nanodrop™ Lite Spectrophotometer (Thermo Scientific Inc., USA) at 260–280 nm and by horizontal gel electrophoresis (Thistle Scientific Ltd, USA) on a 0.8% (w/v) agarose gel at 100 V for 30 min and visualized under UV after staining with GelRed™ (Thermo Scientific, USA).
Amplification was performed in a programmable master cycler thermocycler (C1000-BioRad, USA) using specific primers that target COI gene of mitochondrial DNA, LCO1490 5GGTCAACAAATCATAAAGATATTGG and HC02198 TAAACTTCAGGGTGACCAAAAAATCA. The thermal cycler conditions were as follows: initial denaturation at 98°C for 2 minutes followed by 35 cycles at 98°C for 30 seconds, annealing at 47.3°C for30 seconds, elongation at 75°C for 30 seconds and final elongation at 75°C for 10 minutes. 50µL PCR mix constituted DNA polymerase enzyme 1U/50µL reaction, 5X Buffer 10µL, 10pm dNTP, 50Mm Mg Cl₂ 1µL,10pm primers 1µL each and distilled water.
PCR products were separated by horizontal gel electrophoresis on 1.5% (w/v) agarose gel at 100 V for 45 mins and visualized under UV after staining with 2 µl Gel Red™ (Thermo Scientific). PCR amplicons purified using the Thermo Scientific®GeneJET Purification Kit (EU Lithuania). A ratio of 1:1 volume of binding buffer added to the completed PCR mixture and vortexed to mix properly. When the color of the mixture remained orange or- violet, 10 µl of 3M Sodium acetate (pH 5.2) was added to alter the color to yellow. Eight hundred microliters of the solution transferred to the Gene-JET purification column and centrifuged at 10,000rpm for 30 sec. and the flow-through was discarded. Seven hundred microliters of the buffer (diluted with ethanol) added and centrifuged at 10,000 rpm on a rotor for 30 sec. and the flow-through discarded. Additional centrifugation done to completely remove any residual buffer. The purification column was transferred to clean 1.5 ml microcentrifuge tubes and 50 µl of elution buffer added followed by centrifugation at 10,000 rpm for 1 min to obtain pure DNA amplicons.
Purified PCR products were sequenced by capillary sequencing on ABI 3730xl DNA Analyzer, (Applied Biosystems) using the same forward and reverse primers and ABI BigDye® Terminator v3.1 Cycle Sequencing reaction kit (Applied Biosystems, USA). Forward and reverse sequences were obtained and used to produce consensus sequences after editing using Chromas Pro v3.1. and further editing and alignment was achieved using Bio-edit Sequence Alignment Editor ver. 7.0. The sequences were converted into FASTA file format in preparation for analyses. Sequences obtained from the present study were submitted to the Gen Bank accession number ON455096-ON455103. Additional CO1 sequences for Chironomus from other parts of the world were obtained from the Gen Bank (www.ncbi.nlm.nih.gov) for phylogenetic analyses (Table 1).
Descriptive statistics were employed to evaluate data on heavy metals across stations. The data was analyzed by ANOVA and a Tukey’s pairwise post hoc test was used to confirm significant variations. The composition of aquatic insects was analyzed independently based on morphological and molecular approaches. The non-biting midge, chironomids a species representative sample was identified across stations for genetic analysis. Molecular analysis, assembled sequences were transferred to MEGA v.11 software and pairwise sequence alignment of the nucleotides done using CLUSTAL W according to Tamura et al., (2011). Sequences were submitted to the NCBI BLAST portal (www.ncbi.blast.nlm.nih.gov) for a sequence homology search, and sequences with greater than 97% similarity retrieved for phylogenetic analysis. Evolutionary history was inferred using the Neighbor-Joining algorithm and distances computed using the Maximum Composite Likelihood (Tamura et al., 2004; Tamura et al., 2011). Bootstrap tests (1000 replicates) were used to cluster associated taxa, with replicate trees showing likelihoods above 50% indicated on the branches.