Next Article in Journal
Inhibition of Rice Stripe Virus Accumulation by Polyubiquitin-C in Laodelphax striatellus
Previous Article in Journal
Cyantraniliprole and Thiamethoxam Exposure Changes Expression of Transcripts Associated with Small Non-Coding RNA Processing in the Colorado Potato Beetle
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Response of Chironomidae (Diptera) to DDT, Mercury, and Arsenic Legacy Pollution in Sediments of the Toce River (Northern Italy)

1
National Research Council-Water Research Institute (CNR-IRSA), Via del Mulino 19, 20861 Brugherio, MB, Italy
2
Department of Biology, Chemistry, and Pharmacy, Freie Universität Berlin, Arnimallee 22, 14195 Berlin, Germany
3
Department of Agricultural and Environmental Sciences (DISAA), University of Milan, Via Celoria 2, 20133 Milan, MI, Italy
*
Author to whom correspondence should be addressed.
Current address: Department of Ecohydrology, IGB Leibniz-Institute of Freshwater Ecology and Inland Fisheries, Müggelseedamm 310, 12587 Berlin, Germany.
Submission received: 15 January 2024 / Revised: 18 February 2024 / Accepted: 20 February 2024 / Published: 22 February 2024
(This article belongs to the Section Insect Ecology, Diversity and Conservation)

Abstract

:

Simple Summary

Historical contamination may still pose a toxicity risk to living organisms. The insecticide DDT and the heavy metals mercury (Hg) and arsenic (As) were discharged in the last century (ca. 1915–1996) into the Toce River (Northern Italy) by a factory, and nowadays residual contamination is still detectable in bottom sediments. Various aquatic invertebrates, including non-biting midges (Diptera, Chironomidae), inhabit these sediments and can serve as bioindicators of potential adverse effects caused by the toxicants. We collected and analyzed sediments and Chironomidae at riverine sites located upstream and downstream of the industrial area. Contamination downstream reaches levels potentially toxic for aquatic organisms. A total of 32 chironomid species were identified. Most species were distributed among different sites based on natural factors such as water temperature, water flow, and substrate type. Few species proved to be sensitive to contamination, such as Diamesa spp., Sympotthastia spinifera, and Prodiamesa olivacea to DDT; Potthastia longimanus to Hg; Odontomesa fulva and Microtendipes pedellus to As. These species may serve as bioindicators of contamination in other freshwater ecosystems.

Abstract

The Toce River (Northern Italy) is characterized by legacy contamination of dichloro-diphenyl-trichloroethane (DDT), mercury, and arsenic deriving from an industrial plant active between ca. 1915 and 1996. Chironomidae taxa assemblages and sediments were collected in 2014 and 2019 upstream and downstream of the industrial area to analyze species responses to toxic substances in a river stretch with relatively uniform natural (i.e., hydro-morphological) characteristics. A total of 32 chironomid taxa were identified. Sediment concentrations reached levels potentially toxic for benthic invertebrates: 15.7 µg kg−1 1% organic carbon for DDT, 197 µg kg−1 dry weight (d.w.) for Hg, and 55.7 mg kg−1 d.w. for As. Canonical Correspondence Analysis (CCA) revealed a predominant seasonal gradient, followed by an upstream-downstream gradient. Partial CCA indicated that 5.2% of the total variation was associated with sediment contamination. Self-Organizing Maps (SOMs) were used to represent species responses to toxicants. Most species appeared to be tolerant, e.g., Chironomus riparius, Micropsectra atrofasciata, Conchapelopia pallidula, and Polypedilum spp. Sensitivity to contaminants was observed in only a few species: Diamesa spp., Sympotthastia spinifera, and Prodiamesa olivacea to DDT; Potthastia longimanus to Hg; Odontomesa fulva and Microtendipes pedellus to As. The chironomid community was characterized in presence of contamination levels commonly observed in freshwater ecosystems.

Graphical Abstract

1. Introduction

The response of chironomids from different lentic and lotic habitats to environmental factors has been extensively analyzed in previous studies [1]. Different species have been shown to be sensitive to natural factors, such as water temperature, oxygen level, river flow, substrate particle size, and nutrient content [2]. Thus, the analysis of chironomid taxa assemblages may reveal anthropogenic alterations of these parameters. For example, the relationship with water quality, organic pollution, nutrient enrichment, and oxygen shortage has been observed [3,4,5].
While these organisms are good indicators of overall anthropogenic pressure, their relationship to specific stressors has often proven difficult to disentangle. Regarding the effects of toxic substances, such as trace metals, pesticides, or other pollutants, responses have primarily been described through analysis of specific biomarkers, mouth part deformities [6,7], gene expression [8,9], or alterations in life history parameters [10] in laboratory [11] or in microcosm studies [12]. Conversely, field studies at the community level have rarely shown clear relationships between species and toxicants, as other variables may act as confounding factors. For example, long-term shifts in taxa assemblages were described along sediment cores collected from historically contaminated sites, yet variations in contaminant concentrations have often been associated with trends of other significant drivers, such as nutrient loads or water temperature [13,14]. Studies conducted in contaminated rivers have frequently revealed community structures strongly correlated with variables expressing the longitudinal axis of the river, such as distance from the source, altitude, and water temperature [15]. However, chironomids demonstrate high potential as bioindicators of the presence of toxicants, being one of the most abundant and species-rich macroinvertebrate groups at sites with elevated levels of contamination [16].
The aim of this study is to find potential relations between the chironomid community and toxic contaminants in the Toce River (Northern Italy), focusing on a river stretch characterized by a gradient of dichloro-diphenyl-trichloroethane (DDT), mercury (Hg), and arsenic (As) industrial legacy contamination in conditions of relative uniformity of natural factors such as distance from the source, width, depth, and discharge. A DDT factory equipped with a mercury-cell chlor-alkali plant, located close to the riverbank, caused heavy DDT and Hg contamination in water and sediments during the last century (ca. 1915–1996). Additionally, documentation indicates arsenopyrite roasting to produce sulfuric acid in the industry, leading to As enrichment into the Toce [17]. Pollutant concentrations in different environmental compartments (sediments, biota) of this river are constantly monitored by the International Commission for the Protection of the Italian-Swiss Waters (CIPAIS; https://www.cipais.org/web/, accessed on 10 January 2024). Data indicate that contamination is slowly decreasing due to the natural recovery determined by the transport of uncontaminated sediments from upstream [18]. At the river mouth, total DDT concentrations decreased from 142 µg kg−1 normalized to 1% organic carbon (1% OC) in 2001 to 3.2 µg kg−1 1% OC in 2018; Hg shifted from 949 µg kg−1 dry weight (d.w.) in 2008 to 68 µg kg−1 d.w. in 2018; As decreased from 48 mg kg−1 d.w. in 2008 to 17 mg kg−1 d.w. in 2018 [18], but concentrations in sediments and in the aquatic biota still remain above background levels [19,20].
The response of the entire benthic macroinvertebrate community in this river section was analyzed in 2014 [17]. Partial Redundancy Analysis (RDA) revealed that 5% of the total variance was associated with sediment contamination. Different biotic metrics were calculated, including non-stressor-specific metrics such as the Multimetric Intercalibration Index STAR_ICMi [21] and stressor-specific metrics such as Species at Risk for pesticides (SPEARpesticide [22]) and mean Sensitivity to Hg (SHg [23]), but significant differences between communities upstream and downstream of the factory were not emphasized.
In this study, we have furthered the taxonomic identification of chironomids to the most detailed level possible, aiming to understand if this species-rich group may exhibit relationships with residual contamination. According to published literature, these organisms may be sensitive to DDT, Hg, and As. DDT exposure is known to affect Chironomus emergence, mortality, reproduction, growth, and sex ratio [24,25]. Exposure to Chironomus riparius to Hg has been reported to reduce growth and emergence, delay development time, and decrease activity of the larvae [10], and increase morphological deformities [26]. The latter effect is commonly observed also with exposure to other trace elements, such as As [7,27]. However, community-level responses have rarely been reported. For example, Diggins and Steward [28] showed reduced richness and densities at increasing metal contamination levels, with Procladius, Chironomus, and Cricotopus sp. as the most tolerant taxa. Therefore, the ultimate objectives of this study are to assess the potential toxic effect of residual contamination on chironomid taxa assemblages and identify potential indicators of DDT, Hg, and As pollution.

2. Materials and Methods

2.1. Study Area

The Toce River is 84 km in length and is one of the primary tributaries of Lake Maggiore, situated along the boundary between Italy and Switzerland. Flowing through the Ossola Valley in the Central-Western Alps of the Piedmont Region, Northern Italy, the river exhibits an average annual flow of approximately 65 m3 s−1 (Figure 1). Within its catchment area of 1600 km2, natural land use predominates (92%), while urban and industrial areas only account for 2% and 0.4%, respectively [18].
The river section selected for this study is approximately 23 km long, with the most upstream station (Domo) located 8.6 km from the industrial site and the furthest downstream site (Mergozzo) situated 14.4 km away. This particular river stretch was selected due to its relatively uniform hydro-morphological characteristics, classified as WFD intercalibration type R-A2 and Italian river type 01SS4G, i.e., it belongs to the Western Alps hydroecoregion, situated at elevations ranging between 500 and 1000 m above sea level, characterized by runoff origins, a nival-glacial flow regime, siliceous substrate, and large size (distance from the source ranging between 75 and 150 km). Based on river discharge, this stretch comprises two distinct water bodies: type G1 upstream from Pieve Vergonte (average annual flow at Domo site: 32 m3 s−1) and type G2 downstream (average annual flow at Ornavasso site: 55 m3 s−1), influenced by the inflow of two large water channels (data provided by the Environmental Protection Agency of Piedmont Region, ARPA Piemonte, www.arpa.piemonte.it, accessed on 10 January 2024). The riverbed width ranges from 40 to 60 m, with a maximum depth exceeding 1.5 m.
Background values were estimated as 0.044 ± 0. 026 mg kg−1 d.w. for Hg and 34.3 ± 3.0 mg kg−1 d.w. for As in river sediments [29]. The relatively high geological background for As is due to the presence of arsenopyrite formations [30].
Six sampling sites were selected (Figure 1): two stations located upstream from the industrial site, namely Domo and Prata (8.6 km and 3.4 km from the industrial site, respectively); and four downstream stations, namely Bosco Tenso, Premosello, Ornavasso, and Mergozzo (3.7, 8.7, 13.1, and 14.4 km from the industrial site, respectively).

2.2. Insect Collection

Chironomids were collected in April, July, and October 2014 at all six stations and in February 2019 at all stations except Mergozzo. This latter sampling campaign was aimed at increasing the number of samples for statistical analyses and to account for potential interannual variations in the chironomid community and/or contaminant concentrations.
Sampling of chironomids followed the Italian standard protocol for benthic macroinvertebrate community [21], employing a multihabitat proportional method. Within each site (a 50 m long river stretch), a depositional area (pool mesohabitat) and an erosional area (riffle mesohabitat) were visually identified, where possible. Pools are characterized by finer substrates, greater depth, and lower flow velocity compared to riffles. Based on dominant substrates and flow types, different microhabitat types were identified within each of the two areas: silt (<63 µm grain size), sand (63 µm–2 mm grain size), microlithal (stones with 2–6 cm length), mesolithal (6–20 cm), macrolithal (20–40 cm), coarse particulate organic matter (CPOM), fine particulate organic matter (FPOM), algae, aquatic macrophytes, parts of terrestrial plants (Tp), and submerged dead wood (xylal) (Table S1). Ten replicate units were collected within both the pool and riffle, based on the relative occurrence of each microhabitat type [21]. Each replicate unit consisted of a 0.1 m2 area and was sampled using a Surber net (300 µm mesh size). In total, 114 quantitative samples were collected.
Additionally, drift nets were positioned at each site for at least one hour to collect pupal exuviae and floating adults. These qualitative samples were used solely to achieve a more precise taxonomic identification of taxa and were not included in statistical analyses.
Samples were preserved in 70% ethanol and sorted in the laboratory under a stereomicroscope. Chironomids were mounted onto slides and identified under microscope magnification to the most detailed taxonomic level possible, i.e., genus, species group, or species, following taxonomic keys [31,32,33]. Mature larvae of Chironomus were identified to species level by examining the karyotype of polytene chromosomes extracted from salivary glands, following the method described by Bettinetti et al. [34].

2.3. Water Quality and Sediment Chemical Analyses

Sediment sampling was carried out at the same sites and dates as the chironomid samplings. Different sub-samples were collected at each site using a metal spoon and mixed to obtain a representative sample. The sediments were preserved in acid-washed dark glass bottles at 4 °C until freeze-drying (72 h at 0.2 mbar and −45 °C; Telstar LyoQuest system, Telstar, Élancour, France). Sediments were sieved to separate the finest fraction (<63 µm grain size) for chemical analyses.
Concentrations of DDT (sum) and trace elements (Hg, As, Cd, Cu, Ni, and Pb) in sediments collected in 2014 were analyzed and published in Marziali et al. [17], and concentrations of DDT and Hg in 2019 samples were reported in CIPAIS [35], thus these data were used for statistical analyses. Concentrations of the other trace elements (As, Cd, Cu, Ni, and Pb) in 2019 samples are new original data presented in this work. For this analysis, the following procedure was followed: aliquots of 150–200 mg of freeze-dried sediments (<63 µm grain size fraction) were homogenized with a ball mill (Retsch MM2000, Retsch Technology GmbH, Haan, Germany) and mineralized in closed pressurized Teflon vessels with 6 mL of concentrated HNO3 and 2 mL of ultrapure water using a microwave system (Preekem Excel, Fulltech Instruments, Arenato, Rome, Italy). Solutions were diluted to 50 mL with ultrapure water. After centrifugation, quantification was carried out by Inductively Coupled Plasma-Optical Emission Spectrometry (ICP-OES) (iCap7200 Duo, Thermo Fisher Scientific, Rodano, Milan, Italy). For quality assurance, multielement standard solutions at 20 µg L−1 and 100 µg L−1 were analyzed every 10 samples, with recoveries ranging between 93 and 110% for all elements. Blank values were always below the Limit of Detection (LOD). The certified reference material GBW07305 sediment powder from the National Standard Centre of China was analyzed, with recoveries ranging between 85% and 99% for all elements. Analyses were run in triplicate, obtaining relative standard deviations ≤ 10%.
Total organic carbon (OC) in sediments was determined in 0.5 g d.w. sample aliquots by back-titration after oxidation with potassium dichromate in the presence of sulfuric acid, following the Walkley-Black procedure [36].
During each sampling, water depth was measured with a graduated rod, water velocity was estimated with a flow probe (PCE, PCE Italia s.r.l., Lucca, Italy), and water parameters were measured using field multiprobes (Radiometer Analytical, Hach, Milan, Italy): pH, water temperature, water conductivity, and percent oxygen saturation (Table 1). The dominant flow type was visually estimated as smooth, rippled, unbroken standing waves, broken standing waves, or chute.
Regarding water quality analysis, data on total alkalinity, nitrates (N-NO3), ammonia (N-NH4), total phosphorous and dissolved DDT (sum of DDT, DDE, and DDE), and trace elements were provided by ARPA Piemonte [14], with site and sampling season used as criteria to match chemical data with chironomid samplings (Table 1). Those data were included to test the potential influence of water parameters on chironomids. Dissolved DDT, Hg, and As were mostly below the LOD (Table 1), therefore they were not considered in the data analysis.

2.4. Data Analysis

A multivariate analysis was conducted to examine the relationship between chironomid community composition and environmental factors. The environmental variables considered are outlined in Table 1 and include altitude, distance from the source, mean water depth, current velocity, water temperature, alkalinity, water conductivity, pH, oxygen saturation, ammonia, nitrates, total phosphorous, dominant flow type (coded as follows: 1 = smooth, 2 = rippled, 3 = unbroken standing waves, 4 = broken standing waves, 5 = chute), percent fine fraction (<63 μm grain size) in whole sediments, organic carbon in sediments, concentration of micropollutants and trace metals in sediments, including DDT, As, Hg, Cu, Cd, Ni, and Pb. Additionally, two mesohabitats (pool and riffle), each with different microhabitat types (silt, sand, microlithal, mesolithal, macrolithal, CPOM, FPOM, algae, aquatic macrophytes, Tp, and xylal), were included, totaling 114 samples. In terms of species, 32 taxa identified to genus or species level were included.
The gradient length calculated using Detrended Correspondence Analysis (DECORANA) suggested a bimodal response as more appropriate than a linear one. Therefore, Canonical Correspondence Analysis (CCA) was preferred over Redundancy Analysis (RDA) [37]. The “formula interface” method was chosen for analysis [38,39]. A nXp matrix, comprising n (=114) rows (i.e., microhabitat samplings) and p (=32) columns (chironomid taxa), and a nZs matrix, including the same rows and s (=19) columns (environmental variables), were analyzed. The nXp values were log(x + 1) transformed before calculation.
A forward selection of environmental variables was run to test reduced models containing only the variables contributing most to total variance. With this method, collinear variables were removed. Subsequently, partial Canonical Correspondence Analysis (pCCA) [38,39] was carried out, running a CCA analysis with a species data matrix nXp explained by an environmental data matrix nZs, which constitutes the constrained variables in the presence of conditioning variables nWt (covariates). This pCCA aimed to assess the influence of sediment contamination (DDT, Hg, and As concentrations) on overall community variability. Subsequently, a variation partitioning test was carried out. With this method, the variation bound to sediment contamination was calculated by removing the influence of other variables.
The significance of CCA, pCCA, and variation partitioning was tested by running a Monte Carlo simulation with 999 permutations.
A Self-Organizing Map (SOM) was trained based on the nXp matrix [40]. SOM approaches k-mean clustering in aggregating sites in units and has the advantage that it is not affected by outliers [41,42]. SOM creates a map in which sites are plotted according to similarities in species composition: sites with similar species assemblages are aggregated in clusters of cells. SOM allows a simple representation of sites and species in two dimensions and can be used to easily represent the response of single taxa to target parameters. In fact, SOM analysis was trained as a Supervised Self-Organizing Map (SSOM) [43]. In SSOM, an external variable is included to guide the clustering. Different SSOMs were trained, each including a different column of the nZs matrix, representing an environmental variable or factor. In the present analysis, the site maps were created by dividing the values of each environmental variable into classes to create a factor with different levels, represented by different colors. Cells with the same environmental factor level and similar species composition were clustered together. The empty cells emphasize the distance between clusters (filled cells). Species maps were produced based on codebook values. The relative abundance of each species was figured out on the map. SOM map size was selected considering quantization (QE) and topographic error (TE). QE is a measure of the average distance between the data points and the map nodes to which they are mapped, with smaller values indicating a better fit. TE is a measure of how well the structure of the input space is modeled by the map [44].
All data analyses were carried out in the R environment (R version 4.3.2, 2023-10-31 ucrt) [38,45] supplemented by the following additional packages: vegan [46], adespatial, gclus, cluster, Hmisc, labdsv, pvclust, kohonen [47], awesome [48], plotrix, and RColorBrewer.

3. Results

The ranges of the environmental parameters measured at the sampling sites are presented in Table 1. The target parameters included DDT, showing values in sediments comprised between 0.4 and 15.7 µg kg−1 1% OC; Hg, comprised between 20 and 197 µg kg−1 d.w.; and As, ranging between 3.6 and 55.7 mg kg−1 d.w. (Table 1). Concentrations of contaminants in sediments collected at the downstream sites were significantly higher than those at the upstream sites (t-test, p < 0.05; Table 2).
Table 1. Environmental parameters used in analyses with range and mean values. Abbreviations used in Figures and Tables are also reported.
Table 1. Environmental parameters used in analyses with range and mean values. Abbreviations used in Figures and Tables are also reported.
CompartmentParameterUnitAbbreviationMinMaxMean
Altitudem a.s.l.alt196236217
Distance from the sourcekmdist497259
SedimentsFine sediments%sed0.2277
Organic carbon% d.w.sostorg0.473.71.4
DDTµg kg−1 1% OCDDT0.415.76.2
Hgmg kg−1 d.w.Hg0.020.1970.081
Asmg kg−1 d.w.As3.655.716.6
Cdmg kg−1 d.w.Cd0.10.30.2
Cumg kg−1 d.w.Cu13.755.928.7
Nimg kg−1 d.w.Ni12.347.231.5
Pbmg kg−1 d.w.Pb8.717.713.4
WaterSampling depthcmdepth55024
Current velocitycm s−1velcorr110020
Water temperature°Ctemp715.710.9
ConductivityµS cm−1cond122269181
Alkalinitymg L−1 CaCO3alcal3610268
pHpH unitpH7.187.5
N-NO3mg L−1NO3<0.010.50.4
N-NH4mg L−1NH4<0.010.030.015
dissolved DDTng L−1-<2<2<2
dissolved Hgng L−1-<20<20<20
dissolved Asµg L−1-<34<3
Total Phosphorusmg L−1TP<0.005<0.005<0.005
Oxygen saturation%O296124102
A total of 32 taxa were collected using the multihabitat proportional method and were included in the multivariate analyses (Table 3). The most frequent species were Micropsectra atrofasciata, Microtendipes pedellus, Polypedilum spp., Tvetenia spp., Odontomesa fulva, and Macropelopia spp.
A comprehensive list of species collected in the Toce River, considering also samples obtained using drift nets, is presented in Table S2.
Correlations between environmental parameters and species are presented in Table S3.
A preliminary DECORANA analysis yielded axis lengths of 5.060, 4.762, 3.826, and 4.010 in the first four axes, respectively, suggesting a unimodal species-environment relationship. According to the CCA analysis, the chironomid community exhibited a strong correlation with the 19 environmental factors included in the analysis, explaining 32% of the total inertia (Table S4a). The scores of species and environmental variables for the first two canonical axes are depicted in Figure 2a,b (Tables S5 and S6a).
The CCA highlighted a spatio-temporal gradient (Figure 2c,d). The temporal gradient, evident along the first axis, reflected seasonality, with February samples distinct from July samples, consistent with variations in water temperature values (Figure 2c, Table S5). The spatial gradient manifested as an upstream-downstream ordination along the second axis, separating sites based on increasing depth (Figure 2d, Table S5). However, this trend was confounded by the predominant seasonal gradient. The species driving ordination were Orthocladius frigidus and Eukiefferiella spp., which clearly prevailed in the upstream stations, except in February, when Diamesa spp. and Orthocladius (Euorthocladius) spp. dominated (Table S6a). C. riparius predominated in October. Potthastia longimanus, Conchapelopia pallidula, Cricotopus bicinctus, and C. riparius showed high scores on the second axis with no clear spatial distribution (Table S6a).
The contaminant DDT, along with the trace elements Cu, Cd, and Pb, were associated with percent fine sediments, while As and Ni correlated with organic carbon (Figure 2b) (Tables S3 and S5). Conversely, Hg followed a different pattern, likely influenced by the seasonal gradient.
The stepwise forward analysis conducted to reduce the number of environmental variables identified water temperature, Pb, and flow type as the most significant variables (i.e., with the highest Variance Inflation Factor) to be included in the model, with an adjusted R squared of 0.066, 0.061, and 0.025, respectively. These variables confirmed the spatio-temporal gradient, with a seasonal gradient (water temperature) prevailing on the spatial one (Pb). Temperature was associated with temporal factors (sampling month), Pb was correlated with spatial factors (altitude and distance from the source), and flow type was associated with current velocity (Figure 2, Table S3). This trace element, as well as Cd, Cu, and Ni, were proved to be mostly bound to geogenic origin [18], with slight differences between concentrations at upstream and downstream sites (Table 2). For this reason, they were considered natural variables.
Partial CCA was subsequently performed (Table S4b). Variables measuring sediment contamination (DDT, Hg, and As concentrations) were included as constrained variables in the nZs matrix, while temperature, Pb, and flow type were included as conditioning variables in the nWt matrix (Figure 3). The analysis aimed to determine the extent to which the constrained variables explained the chironomid species pattern when the effect of the natural variables was removed. The results indicated that constrained variables explained 5.2% of total inertia, while conditioning ones accounted for 3.7% (Table S4b).
Partial CCA allowed to separate the upstream sites Domo and Prata, which clustered near the center of the plot, from the contaminated stations, namely Premosello, Bosco Tenso, Ornavasso, and Mergozzo, which were positioned at some distance from the center. DDT and As exerted opposite influences on the first axis, while Hg was correlated with the second axis (Figure 3, Table S7).
Some species were positioned far from the center of the graph and exhibited an inverse correlation with DDT, Hg, and As, respectively, indicating sensitivity to the presence of these toxicants (Figure 3, Tables S6b and S7). Cricotopus (Paratrichocladius) rufiventris, Thienemannimyia sp., Sympotthastia spinifera, and Diamesa spp. showed an inverse correlation with DDT; Cardiocladius sp. and P. longimanus with Hg; C. riparius, Prodiamesa olivacea, Chaetocladius spp., Orthocladius (Euorthocladius) spp., and Paratrissocladius excerptus with As. Conversely, factor score values indicated tolerance of C. riparius to DDT and of P. olivacea to Hg (Figure 3, Tables S6b and S7).
To aid in the interpretation of the response of each species to the toxicants, SSOM maps were used. The selected map size was 10 × 7 cells, which optimized quantization (QE) and topographic (TE) errors for the 70-cell map size carried out with 1000 iterations (Table S8). The species responses to contamination are reported in Figure 4, Figure 5 and Figure 6.
Species maps showed that almost all taxa seem to avoid the sites with the highest DDT concentrations, while they could tolerate moderate DDT contamination (Figure 4). Only C. riparius and P. longimanus showed the highest abundance at sites with the highest DDT levels. Some taxa, such as Diamesa spp., S. spinifera, and P. olivacea, were mostly abundant at sites with the lowest DDT levels, suggesting sensitivity to this toxicant.
Species tolerant to high Hg concentrations included M. atrofasciata, Polypedilum spp., O. fulva, Rheocricotopus spp., C. bicinctus, and P. olivacea, while only P. longimanus avoided sites with the highest Hg concentrations (Figure 5).
Arsenic had a limited influence on species distribution. Only M. pedellus and Orthocladius (Euorthocladius) spp. avoided sites with the highest As concentrations. Conversely, S. spinifera demonstrated tolerance to the toxicant, being most abundant at sites with the highest As concentrations (Figure 6).
Species responses to other variables are reported in Supplementary material: responses to the longitudinal gradient of the river (distance from the source) are reported in Figure S1; to seasonality (sampling month and water temperature) in Figures S2 and S3; to habitat characteristics (mesohabitat and microhabitat type) in Figures S4 and S5; to hydrological factors (flow type, current velocity, and depth) in Figures S10–S12; to physical-chemical parameters of water (conductivity, oxygen saturation) in Figures S13 and S14; to organic carbon in sediments in Figure S15. These responses may be used to derive the ecological traits of species.

4. Discussion

The values of contaminants (DDT, Hg, and As) in sediments downstream of the industrial plant are typical of legacy contamination [19,20,49,50]. Despite covering a relatively narrow range, concentrations exceeded the consensus-based Threshold Effect Concentrations (cb-TECs), indicating thresholds below which toxic effects for benthic invertebrates are not expected (5.28 µg kg−1 1%OC for total DDT, 180 µg kg−1 d.w. for Hg, 9.79 mg kg−1 d.w. for As) [51], showing maxima up to 15.7 µg kg−1 1%OC for DDT, 197 µg kg−1 d.w. for Hg, and 55.7 mg kg−1 d.w. for As. Similar values of DDT and Hg are often observed in freshwater ecosystems, even in the absence of local anthropogenic sources, due to the atmospheric long-distance transport of these contaminants [52,53]. For As, the main source in freshwater is generally geogenic [30,54]. However, in the case of the Toce River, contributions from past anthropogenic activities are reported [17,55] and downstream values exceed the background value of 34.3 ± 3.0 mg kg−1 d.w. [29], reaching concentrations above the consensus-based Probable Effect Concentration (cb-PEC, 33 mg kg−1 d.w.), which is the threshold above which effects on benthic organisms are expected [51].
The Toce River hosts species-rich benthic invertebrate communities [17], as confirmed by chironomid taxa assemblages, which comprised a total of 56 taxa resulting from the analysis of both Surber net and drift net samples (Table S2). This number aligns with taxa richness values obtained in other Italian lowland rivers, such as 56 chironomid taxa in the Adda River (Northern Italy, [56,57]), 68 taxa in the Adige River (Northern Italy [58]), and 79 taxa in Southern rivers (Cilento National Park [59]). The presence of many species, including some sensitive to anthropogenic alterations such as O. fulva and S. spinifera, offers the potential to identify indicators of residual toxic contamination.
The CCA results emphasized that chironomid taxa assemblages were primarily structured by temporal (seasonal) factors rather than the spatial (upstream-downstream) gradient. In fact, species separation n was more evident when grouping sites by sampling month rather than by distance from the source.
Temporal gradients are critical in structuring chironomid communities due to their association with the adult emergence seasonality [60,61]. Most species show their optimum for emerge in winter, spring, summer, or, more rarely, autumn. In this study, we considered different sampling seasons with the aim of collecting a representative set of species. Additionally, samples were collected in different years (2014 and 2019) to account for potential interannual variability related to climatic and hydrological factors [62]. This may explain the strong temporal-seasonal gradient observed.
The spatial gradient within a watershed is generally related to altitude, distance from the source, and it may also be influenced by water depth, water temperature, dissolved oxygen, current velocity, substrate grain size, substrate type, etc. These variables significantly influence chironomid species composition [1,2,15]. In our study, the river stretch considered was selected to minimize the influence of these parameters, with any spatial gradient expected to be primarily bound to the contamination gradient.
To analyze the response of chironomids to contamination, it was necessary to remove the effects of natural spatio-temporal factors. Partial CCA (pCCA) showed that only a small proportion of the total variance could be attributed to pollutants (5.2%). Chironomids were generally tolerant to a wide range of conditions, and only in some cases was it possible to detect a differential response among species.
With the aim of SSOM maps, the distribution of species according to toxicant concentrations was analyzed. The prevalence of C. riparius at sites with higher DDT concentrations agrees with the identification of this species as one of the most tolerant [63]. Similarly, the prevalence of M. atrofasciata at the sites with the highest Hg values confirms the tolerance of this species [3]. C. pallidula and Polypedilum spp. (including Polypedilum nubeculosum and Polypedilum laetum) were confirmed as tolerant species [3,5]. Conversely, O. fulva, S. spinifera, and P. longimanus appeared less tolerant than previously thought [4]. In fact, O. fulva, as well as M. pedellus, showed maxima at sites with the lowest As concentrations. P. longimanus was sensitive to Hg and to DDT. S. spinifera avoided DDT-contaminated sites, as did Diamesa spp. (including Diamesa zernyi and Diamesa tonsa), confirming them to be sensitive taxa [4].
The results also allowed for a preliminary outline of some ecological traits of the most frequent and/or peculiar species, as reported in Table 4. The statistical significance of these traits cannot always be tested due to the need for a larger dataset. However, these descriptions may serve as a basis for future studies aimed at using chironomids as indicators of specific habitat alterations [64].

5. Conclusions

Chironomids are primarily responsive to various natural factors, including season, water temperature, flow type, and substrate type. However, certain taxa may exhibit sensitivity to contaminants that accumulate in river sediments, although this finding requires further investigation. This study represents a preliminary effort to analyze the response of chironomid communities to contamination levels commonly observed in freshwater ecosystems.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/insects15030148/s1, Table S1: List of microhabitats; Table S2: List of chironomid taxa collected in the Toce River; Table S3: Correlations between taxa and environmental parameters; Table S4: Inertia and eigenvalues of CCA and partial CCA axes; Table S5: Results of CCA: scores of environmental factors; Table S6: Scores of species according to CCA and partial CCA; Table S7: Results of pCCA: constraining variables; Table S8: Results of SOM: Quantization and Topographic errors; Figure S1: SSOM map according to distance from the source; Figure S2: SSOM map according to sampling month; Figure S3: SSOM map according to water temperature; Figure S4: SSOM map according to mesohabitat; Figure S5: SSOM map according to microhabitats; Figure S6: SSOM map according to Pb; Figure S7: SSOM map according to Ni; Figure S8: SSOM map according to Cd; Figure S9: SSOM map according to Cu; Figure S10: SSOM map according to flow type; Figure S11: SSOM map according to current velocity; Figure S12: SSOM map according to water depth; Figure S13: SSOM map according to water conductivity; Figure S14: SSOM map according to oxygen saturation; Figure S15: SSOM map according to organic carbon.

Author Contributions

Conceptualization, L.M. and B.R.; methodology, L.M. and B.R.; software, B.R.; formal analysis, B.R. and L.M.; investigation, L.M., N.P. and A.S.; resources, L.M.; data curation, L.M.; writing, review, and editing, L.M. and B.R.; project administration, L.M.; funding acquisition, L.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partially funded by the International Commission for the Protection of Italian-Swiss Waters (CIPAIS, www.cipais.org/web, accessed on 10 January 2024), Research Programs 2013–2015 and 2019–2021.

Data Availability Statement

Data regarding sediment contamination analysis and benthic invertebrates of the Toce River are published in the project reports regarding the analysis of contamination in Lake Maggiore, available at the website www.cipais.org/web (accessed on 10 January 2024).

Acknowledgments

We thank Lucia Valsecchi, Federica Rosignoli, Antonina Drago, and Simona Pascariello for their help in the field and laboratory work. We thank Emiliano Rinaldo Veronesi, Atlantic Technological University, Department of Natural Resources and Environment, Galway, Ireland for reviewing the English text.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Rossaro, B.; Marziali, L.; Montagna, M.; Magoga, G.; Zaupa, S.; Boggero, A. Factors Controlling Morphotaxa Distributions of Diptera Chironomidae in Freshwaters. Water 2022, 14, 1014. [Google Scholar] [CrossRef]
  2. Armitage, P.D.; Cranston, P.S.; Pinder, L.C.V. The Chironomidae; Springer: Dordrecht, The Netherlands, 1995; ISBN 978-94-010-4308-3. [Google Scholar]
  3. Pillot, H.K.M.M. Chironomidae Larvae, Vol. 2: Chironomini: Biology and Ecology of the Chironomini; KNNV Publishing: Leiden, The Netherlands, 2009; ISBN 978-90-04-27804-2. [Google Scholar]
  4. Pillot, H.K.M.M. Chironomidae Larvae, Vol. 3: Orthocladiinae: Biology and Ecology of the Aquatic Orthocladiinae; KNNV Publishing: Leiden, The Netherlands, 2014; ISBN 978-90-04-27805-9. [Google Scholar]
  5. Vallenduuk, H.J.; Pillot, H.K.M.M. Chironomidae Larvae, Vol. 1: Tanypodinae: General Ecology and Tanypodinae; KNNV Publishing: Leiden, The Netherlands, 2014; ISBN 978-90-04-27803-5. [Google Scholar]
  6. Di Veroli, A.; Goretti, E.; Paumen, M.L.; Kraak, M.H.S.; Admiraal, W. Induction of Mouthpart Deformities in Chironomid Larvae Exposed to Contaminated Sediments. Environ. Pollut. 2012, 166, 212–217. [Google Scholar] [CrossRef]
  7. Di Veroli, A.; Santoro, F.; Pallottini, M.; Selvaggi, R.; Scardazza, F.; Cappelletti, D.; Goretti, E. Deformities of Chironomid Larvae and Heavy Metal Pollution: From Laboratory to Field Studies. Chemosphere 2014, 112, 9–17. [Google Scholar] [CrossRef]
  8. Arambourou, H.; Llorente, L.; Moreno-Ocio, I.; Herrero, Ó.; Barata, C.; Fuertes, I.; Delorme, N.; Méndez-Fernández, L.; Planelló, R. Exposure to Heavy Metal-Contaminated Sediments Disrupts Gene Expression, Lipid Profile, and Life History Traits in the Midge Chironomus riparius. Water Res. 2020, 168, 115165. [Google Scholar] [CrossRef]
  9. Stefani, F.; Rusconi, M.; Valsecchi, S.; Marziali, L. Evolutionary Ecotoxicology of Perfluoralkyl Substances (PFASs) Inferred from Multigenerational Exposure: A Case Study with Chironomus riparius (Diptera, Chironomidae). Aquat. Toxicol. 2014, 156, 41–51. [Google Scholar] [CrossRef] [PubMed]
  10. Azevedo-Pereira, H.M.V.S.; Soares, A.M.V.M. Effects of Mercury on Growth, Emergence, and Behavior of Chironomus riparius Meigen (Diptera: Chironomidae). Arch. Environ. Contam. Toxicol. 2010, 59, 216–224. [Google Scholar] [CrossRef]
  11. Planelló, R.; Servia, M.J.; Gómez-Sande, P.; Herrero, Ó.; Cobo, F.; Morcillo, G. Transcriptional Responses, Metabolic Activity and Mouthpart Deformities in Natural Populations of Chironomus riparius Larvae Exposed to Environmental Pollutants. Environ. Toxicol. 2015, 30, 383–395. [Google Scholar] [CrossRef] [PubMed]
  12. Choung, C.B.; Hyne, R.V.; Stevens, M.M.; Hose, G.C. The Ecological Effects of a Herbicide–Insecticide Mixture on an Experimental Freshwater Ecosystem. Environ. Pollut. 2013, 172, 264–274. [Google Scholar] [CrossRef] [PubMed]
  13. Marchetto, A.; Lami, A.; Musazzi, S.; Massaferro, J.; Langone, L.; Guilizzoni, P. Lake Maggiore (N. Italy) Trophic History: Fossil Diatom, Plant Pigments, and Chironomids, and Comparison with Long-Term Limnological Data. Quat. Int. 2004, 113, 97–110. [Google Scholar] [CrossRef]
  14. Armstrong, I.; Moir, K.E.; Ridal, J.J.; Cumming, B.F. Subfossil Chironomid Assemblages as Indicators of Remedial Efficacy in the Historically Contaminated St. Lawrence River at Cornwall, Ontario. Arch. Environ. Contam. Toxicol. 2023, 85, 191–207. [Google Scholar] [CrossRef] [PubMed]
  15. Ruse, L.P.; Ruse, L.R.; Herrmann, S.J.; Sublette, J.E. Chironomidae (Diptera) Species Distribution Related to Environmental Characteristics of the Metal-Polluted Arkansas River, Colorado. West. N. Am. Nat. 2000, 60, 34–56. [Google Scholar]
  16. Wright, I.A.; Burgin, S. Effects of Organic and Heavy Metal Pollution on Chironomids within a Pristine Upland Catchment. Hydrobiologia 2009, 635, 15–25. [Google Scholar] [CrossRef]
  17. Marziali, L.; Rosignoli, F.; Drago, A.; Pascariello, S.; Valsecchi, L.; Rossaro, B.; Guzzella, L. Toxicity Risk Assessment of Mercury, DDT and Arsenic Legacy Pollution in Sediments: A Triad Approach under Low Concentration Conditions. Sci. Total Environ. 2017, 593–594, 809–821. [Google Scholar] [CrossRef] [PubMed]
  18. Marziali, L.; Guzzella, L.; Salerno, F.; Marchetto, A.; Valsecchi, L.; Tasselli, S.; Roscioli, C.; Schiavon, A. Twenty-Year Sediment Contamination Trends in Some Tributaries of Lake Maggiore (Northern Italy): Relation with Anthropogenic Factors. Environ. Sci. Pollut. Res. 2021, 28, 38193–38208. [Google Scholar] [CrossRef] [PubMed]
  19. Marziali, L.; Roscioli, C.; Valsecchi, L. Mercury Bioaccumulation in Benthic Invertebrates: From Riverine Sediments to Higher Trophic Levels. Toxics 2021, 9, 197. [Google Scholar] [CrossRef] [PubMed]
  20. Tasselli, S.; Marziali, L.; Roscioli, C.; Guzzella, L. Legacy Dichlorodiphenyltrichloroethane (DDT) Pollution in a River Ecosystem: Sediment Contamination and Bioaccumulation in Benthic Invertebrates. Sustainability 2023, 15, 6493. [Google Scholar] [CrossRef]
  21. Spitale, D. Performance of the STAR_ICMi Macroinvertebrate Index and Implications for Classification and Biomonitoring of Rivers. Knowl. Manag. Aquat. Ecosyst. 2017, 418, 20. [Google Scholar] [CrossRef]
  22. Liess, M.; Liebmann, L.; Vormeier, P.; Weisner, O.; Altenburger, R.; Borchardt, D.; Brack, W.; Chatzinotas, A.; Escher, B.; Foit, K.; et al. Pesticides Are the Dominant Stressors for Vulnerable Insects in Lowland Streams. Water Res. 2021, 201, 117262. [Google Scholar] [CrossRef]
  23. Malaj, E.; Grote, M.; Schäfer, R.B.; Brack, W.; Von Der Ohe, P.C. Physiological Sensitivity of Freshwater Macroinvertebrates to Heavy Metals. Environ. Toxicol. Chem. 2012, 31, 1754–1764. [Google Scholar] [CrossRef]
  24. Ma, P.; Li, H.; You, J. Full–Life Cycle Toxicity Assessment of Sediment-Bound DDT and Its Degradation Products on Chironomus dilutus. Environ. Toxicol. Chem. 2019, 38, 2698–2707. [Google Scholar] [CrossRef]
  25. Rakotondravelo, M.L.; Anderson, T.D.; Charlton, R.E.; Zhu, K.Y. Sublethal Effects of Three Pesticides on Larval Survivorship, Growth, and Macromolecule Production in the Aquatic Midge, Chironomus tentans (Diptera: Chironomidae). Arch. Environ. Contam. Toxicol. 2006, 51, 352–359. [Google Scholar] [CrossRef]
  26. Gremyatchikh, V.A.; Tomilina, I.I.; Grebenyuk, L.P. The Effect of Mercury Chloride on Morphofunctional Parameters in Chironomus riparius Meigen (Diptera, Chironomidae) Larvae. Inland Water Biol. 2009, 2, 89–95. [Google Scholar] [CrossRef]
  27. Martinez, E.A.; Moore, B.C.; Schaumloffel, J.; Dasgupta, N. The Potential Association Between Menta Deformities and Trace Elements in Chironomidae (Diptera) Taken from a Heavy Metal Contaminated River. Arch. Environ. Contam. Toxicol. 2002, 42, 286–291. [Google Scholar] [CrossRef] [PubMed]
  28. Diggins, T.P.; Stewart, K.M. Chironomid Deformities, Benthic Community Composition, and Trace Elements in the Buffalo River (New York) Area of Concern. J. N. Am. Benthol. Soc. 1998, 17, 311–323. [Google Scholar] [CrossRef]
  29. Vignati, D.A.L.; Guilizzoni, P. Metalli Nel Lago Maggiore: Livelli Naturali e Antropici. Acqua Aria 2011, 1, 22–27. [Google Scholar]
  30. Pfeifer, H.-R.; Derron, M.-H.; Rey, D.; Schlegel, C.; Atteia, O.; Piazza, R.D.; Dubois, J.-P.; Mandia, Y. Chapter 2 Natural Trace Element Input to the Soil-Sediment-Water-Plant System: Examples of Background and Contaminated Situations in Switzerland, Eastern France and Northern Italy. In Trace Metals in the Environment; Elsevier: Amsterdam, The Netherlands, 2000; Volume 4, pp. 33–86. ISBN 978-0-444-50532-3. [Google Scholar]
  31. Andersen, T.; Cranston, P.S.; Epler, J.H. The Larvae of Chironomidae (Diptera) of the Holarctic Region—Keys and Diagnoses; Cederholm, L., Ed.; Insect Systematics & Evolution: Lund, Sweden, 2013; Volume 66. [Google Scholar]
  32. Wiederholm, T. Chironomidae of the Holarctic Region: Pupae; Chironomidae of the Holarctic Region: Keys and Diagnoses; Entomologica Scandinavica: Lund, Sweden, 1983. [Google Scholar]
  33. Wiederholm, T. Chironomidae of the Holarctic Region: Adult Males; Chironomidae of the Holarctic Region: Keys and Diagnoses; Entomologica Scandinavica: Lund, Sweden, 1983. [Google Scholar]
  34. Bettinetti, R.; Colombo, S.; Rossaro, B. Cytogenetic Characteristics of Different Populations of Chironomus riparius Meigen, 1804. SIL Proc. 1922–2010 2000, 27, 2363–2366. [Google Scholar] [CrossRef]
  35. CIPAIS Indagini Sulle Sostanze Pericolose nell’Ecosistema del Lago Maggiore. Programma 2019–2021. Rapporto Annuale 2019 2020. Available online: www.cipais.org/net (accessed on 10 January 2024).
  36. US-EPA. Methods for the Determination of Total Organic Carbon (TOC) in Soils and Sediments; US Environmental Protection Agency, Ecological Risk Assessment Support Center, Office of Research and Development NCEA-C-1282, EMASC-001: Washington, DC, USA, 2002.
  37. Ter Braak, C.J.F.; Prentice, I.C. A Theory of Gradient Analysis. In Advances in Ecological Research; Begon, M., Fitter, A.H., Ford, E.D., Macfadyen, A., Eds.; Academic Press: Cambridge, MA, USA, 1988; Volume 18, pp. 271–317. ISBN 0065-2504. [Google Scholar]
  38. Borcard, D.; Gillet, F.; Legendre, P. Numerical Ecology with R; Use R! 2nd ed.; Springer International Publishing: Cham, Switzerland, 2018; ISBN 978-3-319-71403-5. [Google Scholar]
  39. Oksanen, J.; Blanchet, F.G.; Kindt, R.; Legendre, P.; Minchin, P.; O’Hara, B.; Simpson, G.; Solymos, P.; Stevens, H.; Wagner, H. Vegan: Community Ecology Package. R Package Version 2.2-1, 2015. Available online: https://cran.r-project.org/web/packages/vegan/vegan.pdf (accessed on 10 January 2024).
  40. Giraudel, J.L.; Lek, S. A Comparison of Self-Organizing Map Algorithm and Some Conventional Statistical Methods for Ecological Community Ordination. Ecol. Model. 2001, 146, 329–339. [Google Scholar] [CrossRef]
  41. Park, Y.-S.; Céréghino, R.; Compin, A.; Lek, S. Applications of Artificial Neural Networks for Patterning and Predicting Aquatic Insect Species Richness in Running Waters. Ecol. Model. 2003, 160, 265–280. [Google Scholar] [CrossRef]
  42. Lek, S.; Guégan, J.-F. Artificial Neuronal Networks; Springer: Berlin/Heidelberg, Germany, 2000; ISBN 978-3-642-63116-0. [Google Scholar]
  43. Wehrens, R.; Kruisselbrink, J. Flexible Self-Organizing Maps in Kohonen 3.0. J. Stat. Soft 2018, 87, 1–18. [Google Scholar] [CrossRef]
  44. Breard, G. Evaluating Self-Organizing Map Quality Measures as Convergence Criteria; University of Rhode Island: Kingston, RI, USA, 2017. [Google Scholar]
  45. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2021; Available online: https://www.R-project.org/ (accessed on 10 January 2024).
  46. Dixon, P. VEGAN, a Package of R Functions for Community Ecology. J. Veg. Sci. 2003, 14, 927–930. [Google Scholar] [CrossRef]
  47. Kohonen, T. Self-Organizing Maps; Springer Series in Information Sciences; Springer: Berlin/Heidelberg, Germany, 2001; Volume 30, ISBN 978-3-540-67921-9. [Google Scholar]
  48. Boelaert, J.; Ollion, E.; Megdoud, M.; Naji, O.; Lemba Kote, A.; Renoud, T.; Hym, S. aweSOM: Interactive Self-Organizing Maps, R. R Package Version 1.3; 2022. Available online: https://CRAN.R-project.org/package=aweSOM (accessed on 10 January 2024).
  49. Ivanova, A.; Wiberg, K.; Ahrens, L.; Zubcov, E.; Dahlberg, A.-K. Spatial Distribution of Legacy Pesticides in River Sediment from the Republic of Moldova. Chemosphere 2021, 279, 130923. [Google Scholar] [CrossRef]
  50. Janssen, S.E.; Tate, M.T.; Krabbenhoft, D.P.; DeWild, J.F.; Ogorek, J.M.; Babiarz, C.L.; Sowers, A.D.; Tuttle, P.L. The Influence of Legacy Contamination on the Transport and Bioaccumulation of Mercury within the Mobile River Basin. J. Hazard. Mater. 2021, 404, 124097. [Google Scholar] [CrossRef]
  51. MacDonald, D.D.; Ingersoll, C.G.; Berger, T.A. Development and Evaluation of Consensus-Based Sediment Quality Guidelines for Freshwater Ecosystems. Arch. Environ. Contam. Toxicol. 2000, 39, 20–31. [Google Scholar] [CrossRef]
  52. Rimondi, V.; Costagliola, P.; Lattanzi, P.; Morelli, G.; Cara, G.; Cencetti, C.; Fagotti, C.; Fredduzzi, A.; Marchetti, G.; Sconocchia, A.; et al. A 200 Km-Long Mercury Contamination of the Paglia and Tiber Floodplain: Monitoring Results and Implications for Environmental Management. Environ. Pollut. 2019, 255, 113191. [Google Scholar] [CrossRef] [PubMed]
  53. Ahmed, G.; Anawar, H.M.; Takuwa, D.T.; Chibua, I.T.; Singh, G.S.; Sichilongo, K. Environmental Assessment of Fate, Transport and Persistent Behavior of Dichlorodiphenyltrichloroethanes and Hexachlorocyclohexanes in Land and Water Ecosystems. Int. J. Environ. Sci. Technol. 2015, 12, 2741–2756. [Google Scholar] [CrossRef]
  54. Shaji, E.; Santosh, M.; Sarath, K.V.; Prakash, P.; Deepchand, V.; Divya, B.V. Arsenic Contamination of Groundwater: A Global Synopsis with Focus on the Indian Peninsula. Geosci. Front. 2021, 12, 101079. [Google Scholar] [CrossRef]
  55. Rosso, J.J.; Schenone, N.F.; Pérez Carrera, A.; Fernández Cirelli, A. Concentration of Arsenic in Water, Sediments and Fish Species from Naturally Contaminated Rivers. Environ. Geochem. Health 2013, 35, 201–214. [Google Scholar] [CrossRef] [PubMed]
  56. Groppali, R. Atlante Della Biodiversità Del Parco Adda Sud. Primo Elenco Delle Specie Viventi Nell’area Protetta; Conoscere il parco; Parco Adda Sud: Lodi, Italy, 2006. [Google Scholar]
  57. Marziali, L.; Casalegno, C.; Rossaro, B. The First Record of the Subfamily Buchonomyiinae (Diptera, Chironomidae) from Italy. Ital. J. Zool. 2004, 71, 341–345. [Google Scholar] [CrossRef]
  58. Lencioni, V.; Moubayed-Breil, J. Urban Chironomid Fauna (Diptera: Chironomidae) Visiting the Streets of Verona along the Adige River Banks and Adjacent Spring. Stor. Nat. Della Città Di Verona-Mem. Del Mus. Civ. Di Stor. Nat. Di Verona-2. Ser.-Monogr. Nat. 2021, 6, 157–170. Available online: https://museodistorianaturale.comune.verona.it/media/_Musei/_StoriaNaturale/_Allegati/Biblioteca/Memorie%202.%20serie/Memorie%20MN/MSN_Storia_naturale_Verona_2021.pdf (accessed on 10 January 2024).
  59. Marziali, L.; Armanini, D.G.; Cazzola, M.; Erba, S.; Toppi, E.; Buffagni, A.; Rossaro, B. Responses of Chironomid Larvae (Insecta, Diptera) to Ecological Quality in Mediterranean River Mesohabitats (South Italy). River Res. Appl. 2010, 26, 1036–1051. [Google Scholar] [CrossRef]
  60. Reiss, F. Ökologische Und Systematische Untersuchungen an Chironomiden (Diptera) Des Bodensees: Ein Beitrag Zur Lakustrischen Chironomidenfauna Des Nördlichen Alpenvorlandes. Arch. Für Hydrobiol. 1968, 64, 176–323. [Google Scholar]
  61. Ringe, F. Chironomiden-Emergenz 1970 in Breitenbach Und Rohrwiesenbach. Schlitzer Produktionsbiologische Studien (10). Arch. Für Hydrobiol./Suppl. 1974, 45, 212–304. [Google Scholar]
  62. Lehmann, J. Die Chironomiden Der Fulda (Systematische, Ökologische Und Faunistische Untersuchungen). Arch. Für Hydrobiol. Suppl. 1971, 37, 466–555. [Google Scholar]
  63. Thienemann, A. Chironomus, Leben, Verbreitung Und Wirtschaftliche-Bedeutung Der Chironomiden; Schweizerbart’s Verlagbuchhandlung: Stuttgart, Germany, 1974. [Google Scholar]
  64. Serra, S.R.Q.; Cobo, F.; Graça, M.A.S.; Dolédec, S.; Feio, M.J. Synthesising the Trait Information of European Chironomidae (Insecta: Diptera): Towards a New Database. Ecol. Indic. 2016, 61, 282–292. [Google Scholar] [CrossRef]
Figure 1. Map of the Toce River (Piedmont Region, Northern Italy) with the location of sampling stations (blue markers) and of the industrial site (red marker).
Figure 1. Map of the Toce River (Piedmont Region, Northern Italy) with the location of sampling stations (blue markers) and of the industrial site (red marker).
Insects 15 00148 g001
Figure 2. Results of CCA, including all variables. Plot of the first two axes with scores of (a) species and (b) environmental variables. Cases are colored according to: (c) sampling month; (d) station. Abbreviations of variable and species names are reported in Table 1 and Table 3.
Figure 2. Results of CCA, including all variables. Plot of the first two axes with scores of (a) species and (b) environmental variables. Cases are colored according to: (c) sampling month; (d) station. Abbreviations of variable and species names are reported in Table 1 and Table 3.
Insects 15 00148 g002
Figure 3. Partial CCA, including DDT, Hg, and As as constrained variables, with water temperature (temp), Pb concentration in sediments, and flow type as conditioning variables. Sampling stations (centroids) are plotted with different colors. Abbreviations of species and variable names are reported in Table 1 and Table 3.
Figure 3. Partial CCA, including DDT, Hg, and As as constrained variables, with water temperature (temp), Pb concentration in sediments, and flow type as conditioning variables. Sampling stations (centroids) are plotted with different colors. Abbreviations of species and variable names are reported in Table 1 and Table 3.
Insects 15 00148 g003
Figure 4. Results of SSOM analysis using DDTsum as an ordering external factor: (a) sites mapped in 5 clusters, with different colors: <3 µg kg−1 1%OC (blue cells), 3–4.5 µg kg−1 1%OC (green cells), 4.5–7.5 µg kg−1 1%OC (yellow cells), 7.5–12 µg kg−1 1%OC (orange cells), >12 µg kg−1 1%OC (red cells); (b) distances between clusters: high distances in yellow, low distances in red; (c) species graphs representing the different abundances of a species in the cells (values on y axis represent log10-tranformed abundance): very abundant (purple cells), abundant (cyan cells), present (green cells), absent (red cells). Abbreviations of species names are reported in Table 3.
Figure 4. Results of SSOM analysis using DDTsum as an ordering external factor: (a) sites mapped in 5 clusters, with different colors: <3 µg kg−1 1%OC (blue cells), 3–4.5 µg kg−1 1%OC (green cells), 4.5–7.5 µg kg−1 1%OC (yellow cells), 7.5–12 µg kg−1 1%OC (orange cells), >12 µg kg−1 1%OC (red cells); (b) distances between clusters: high distances in yellow, low distances in red; (c) species graphs representing the different abundances of a species in the cells (values on y axis represent log10-tranformed abundance): very abundant (purple cells), abundant (cyan cells), present (green cells), absent (red cells). Abbreviations of species names are reported in Table 3.
Insects 15 00148 g004
Figure 5. Results of SSOM analysis using Hg as an ordering external factor: (a) sites mapped in 5 clusters, with different colors: <40 µg kg−1 d.w. (blue cells), 40–60 µg kg−1 d.w. (green cells), 60–100 µg kg−1 d.w. (yellow cells), 100–160 µg kg−1 d.w. (orange cells), >160 µg kg−1 d.w. (red cells); (b) distances between clusters: high distances in yellow, low distances in red; (c) species graphs representing the different abundances of a species in the cells (values on y axis represent log10-tranformed abundance): very abundant (purple cells), abundant (cyan cells), present (green cells), absent (red cells). Abbreviations of species names are reported in Table 3.
Figure 5. Results of SSOM analysis using Hg as an ordering external factor: (a) sites mapped in 5 clusters, with different colors: <40 µg kg−1 d.w. (blue cells), 40–60 µg kg−1 d.w. (green cells), 60–100 µg kg−1 d.w. (yellow cells), 100–160 µg kg−1 d.w. (orange cells), >160 µg kg−1 d.w. (red cells); (b) distances between clusters: high distances in yellow, low distances in red; (c) species graphs representing the different abundances of a species in the cells (values on y axis represent log10-tranformed abundance): very abundant (purple cells), abundant (cyan cells), present (green cells), absent (red cells). Abbreviations of species names are reported in Table 3.
Insects 15 00148 g005
Figure 6. Results of SSOM analysis using As as an ordering external factor: (a) sites mapped in 5 clusters, with different colors: <8 mg kg−1 d.w. (blue cells), 8–10 mg kg−1 d.w. (green cells), 10–16 mg kg−1 d.w. (yellow cells), 16–25 mg kg−1 d.w. (orange cells), >25 mg kg−1 d.w. (red cells); (b) distances between clusters: high distances in yellow, low distances in red; (c) species graphs representing the different abundances of a species in the cells (values on y axis represent log10-tranformed abundance): very abundant (purple cells), abundant (cyan cells), present (green cells), absent (red cells). Abbreviations of species names are reported in Table 3.
Figure 6. Results of SSOM analysis using As as an ordering external factor: (a) sites mapped in 5 clusters, with different colors: <8 mg kg−1 d.w. (blue cells), 8–10 mg kg−1 d.w. (green cells), 10–16 mg kg−1 d.w. (yellow cells), 16–25 mg kg−1 d.w. (orange cells), >25 mg kg−1 d.w. (red cells); (b) distances between clusters: high distances in yellow, low distances in red; (c) species graphs representing the different abundances of a species in the cells (values on y axis represent log10-tranformed abundance): very abundant (purple cells), abundant (cyan cells), present (green cells), absent (red cells). Abbreviations of species names are reported in Table 3.
Insects 15 00148 g006
Table 2. Comparison between concentrations of DDT and trace elements (mean ± standard deviation) in sediments collected at sites upstream and downstream of the industrial site (t-test, p-value).
Table 2. Comparison between concentrations of DDT and trace elements (mean ± standard deviation) in sediments collected at sites upstream and downstream of the industrial site (t-test, p-value).
ContaminantUnitUpstream SitesDownstream Sitesp-Value
DDTµg kg−1 1% OC4.0 ± 3.17.5 ± 3.9<0.0001
Hgmg kg−1 d.w.0.05 ± 0.020.10 ± 0.05<0.0001
Asmg kg−1 d.w.7.2 ± 3.022.1 ± 9.8<0.0001
Cdmg kg−1 d.w.0.15 ± 0.040.17 ± 0.040.111
Cumg kg−1 d.w.31.2 ± 12.927.2 ± 8.10.044
Nimg kg−1 d.w.37.4 ± 5.928.0 ± 6.2<0.0001
Pbmg kg−1 d.w.12.6 ± 2.813.8 ± 2.20.015
Table 3. List of chironomid taxa collected in the Toce River using the multihabitat proportional method, with abbreviations used in Figures, density, and frequency (i.e., percent of samples where the species was observed on a total of 114 samples).
Table 3. List of chironomid taxa collected in the Toce River using the multihabitat proportional method, with abbreviations used in Figures, density, and frequency (i.e., percent of samples where the species was observed on a total of 114 samples).
TaxaAuthorsAbbreviationDensity (ind m−2)Frequency (%)
Micropsectra atrofasciata(Kieffer, 1911)Matro17735
Microtendipes pedellus(De Geer, 1776)Microt16931
Polypedilum spp.Kieffer, 1912Polyp32131
Tvetenia spp.Kieffer, 1922Tvet22231
Odontomesa fulva(Kieffer, 1919)Odont202629
Macropelopia spp.Thienemann, 1916Mpelop13622
Sympotthastia spiniferaSerra-Tosio, 1968Symp36422
Diamesa spp.Meigen, 1835Diam16024
Conchapelopia pallidula(Meigen, 1818)Concha8116
Prodiamesa olivacea(Meigen, 1818)Prod10216
Eukiefferiella spp.Thienemann, 1926Euk22923
Paracladopelma spp.Harnish, 1923Pclad3813
Monodiamesa batyphila(Kieffer, 1918)Monod9112
Orthocladius (Orthocladius) spp.Van der Wulp, 1874Orth3112
Cricotopus (Cricotopus) tremulus(Linnaeus, 1756)Ctrem4211
Cricotopus (Cricotopus) bicinctus(Meigen, 1818)Cbici16810
Rheocricotopus spp.Brundin, 1956Rheoc1310
Orthocladius (Euorthocladius) spp.Thienemann, 1935Euorth4910
Chironomus ripariusMeigen, 1804Cripa597
Synorthocladius semivirens(Kieffer, 1909)Synor296
Cardiocladius sp.Kieffer, 1912Cardio125
Cricotopus (Paratrichocladius) rufiventris(Meigen, 1830)Ptrich355
Potthastia longimanus(Kieffer, 1922)Potth275
Chaetocladius spp.Kieffer, 1911Chaet234
Thienemannimyia sp.Fittkau, 1957Thimmyia454
Corynoneura spp.Winnertz, 1846Cory64
Metriocnemus spp.van der Wulp, 1874Metr33
Brillia spp.Kieffer, 1913Bril93
Orthocladius (Mesorthocladius) frigidus(Zetterstedt, 1838)Ofrig112
Psectrocladius sp.Kieffer, 1906Psectr52
Tanytarsini spp.-Tany22
Paratrissocladius excerptus(Walker, 1856)Ptriss21
Table 4. Ecological traits of the most typical species collected in the Toce River.
Table 4. Ecological traits of the most typical species collected in the Toce River.
TaxaEcological Traits
Macropelopia spp.The genus prefers waters with higher oxygen saturation and moderate water flow. No particular sensitivity to toxicants was observed. It is a predator of other chironomid larvae, or Oligochaetes [5].
Conchapelopia pallidulaThe species prevails at the downstream stations in the summer. It prefers the pool habitat. It shows no flow type preference and no sensitivity to toxicants in the Toce River. It is a common predator of other chironomid larvae and Oligochaetes [5].
Diamesa spp.The taxon includes Diamesa zernyi and Diamesa tonsa. It prefers running waters with high flow velocity and turbulence, such as upstream riffle habitats with lower water temperatures. It tolerates a moderate organic carbon concentration in sediments and high Hg levels, while it is mostly absent at sites with higher DDT concentrations.
Potthastia longimanusThis species prefers pool microhabitat, silty sediments, broken standing wave flow type, moderate organic matter content in sediments, and high DDT concentrations, while it is moslty absent at sites with higher Hg concentrations.
Sympotthastia spiniferaAbundant in February, this species is rare in October and April. It prefers pools rich in submerged macrophytes, but it is present at all microhabitat types and stations. It tolerates high organic carbon content (>3.2%) and is present at sites most contaminated by As and Hg, while it is mostly absent at DDT-contaminated sites.
Odontomesa fulvaAbundant in April, it is rare in October and February. It prefers pool areas with rippled flow type and sand and mesolithal microhabitats, tolerating organic carbon in sediments. It is present at all stations, including sites with higher Hg and DDT concentrations but lower As levels.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Marziali, L.; Pirola, N.; Schiavon, A.; Rossaro, B. Response of Chironomidae (Diptera) to DDT, Mercury, and Arsenic Legacy Pollution in Sediments of the Toce River (Northern Italy). Insects 2024, 15, 148. https://doi.org/10.3390/insects15030148

AMA Style

Marziali L, Pirola N, Schiavon A, Rossaro B. Response of Chironomidae (Diptera) to DDT, Mercury, and Arsenic Legacy Pollution in Sediments of the Toce River (Northern Italy). Insects. 2024; 15(3):148. https://doi.org/10.3390/insects15030148

Chicago/Turabian Style

Marziali, Laura, Niccolò Pirola, Alfredo Schiavon, and Bruno Rossaro. 2024. "Response of Chironomidae (Diptera) to DDT, Mercury, and Arsenic Legacy Pollution in Sediments of the Toce River (Northern Italy)" Insects 15, no. 3: 148. https://doi.org/10.3390/insects15030148

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