| - Biomarkers are an important tool in laboratory assays that link exposure or effect of specific toxicants to key molecular and cellular events, but they have not been widely used in invertebrate populations exposed to complex mixtures of environmental contaminants in their natural habitats. The present study focused on a battery of biomarkers and their comparative analysis in natural populations of the benthic larvae of Chironomus riparius (Diptera), sampled in three differentially polluted rivers (the Con, Sar, and Louro in Galicia, Spain). In our study, some parameters were identified, such as hsp70 gene activity, GST enzymatic activity, total glycogen content and mouthpart deformities, which showed significant differences among populations from the three rivers that differed in the levels and types of sedimentary contaminants analyzed (metals, organic-chlorine pesticides, alkylphenols, pharmaceutical, and personal care products). In contrast to these sensitive biomarkers, other parameters showed no significant differences (hsc70 gene, EcR gene, P450 gene, RNA:DNA ratio, total protein content), and were stable even when comparing field and nonexposed laboratory populations. The hsp70 gene seems to be particularly sensitive to conditions of pollutant exposure, while its constitutive counterpart hsc70 showed invariable expression, suggesting that the hsc70/hsp70 ratio may be a potential indicator of polluted environments. Although further studies are required to understand the correlation between molecular responses and the ecological effects of pollutants on natural populations, the results provide new data about the biological responses to multiple-stressor environments. This field study adds new molecular endpoints, including gene expression, as suitable tools that, complementing other ecotoxicological parameters, may help to improve the methodologies of freshwater monitoring under the increasing burden of xenobiotics.