First evidence of SARS-CoV-2 genome detection in zebra mussel (Dreissena polymorpha).

The uses of bivalve molluscs in environmental biomonitoring have recently gained momentum due to their ability to indicate and concentrate human pathogenic microorganisms. In the context of the health crisis caused by the COVID-19 epidemic, the objective of this study was to determine if the SARS-CoV-2 ribonucleic acid genome can be detected in zebra mussels (Dreissena polymorpha) exposed to raw and treated urban wastewaters from two separate plants to support its interest as bioindicator of the SARS-CoV-2 genome contamination in water. The zebra mussels were exposed to treated wastewater through caging at the outlet of two plants located in France, as well as to raw wastewater in controlled conditions. Within their digestive tissues, our results showed that SARS-CoV-2 genome was detected in zebra mussels, whether in raw and treated wastewaters. Moreover, the detection of the SARS-CoV-2 genome in such bivalve molluscans appeared even with low concentrations in raw wastewaters. This is the first detection of the SARS-CoV-2 genome in the tissues of a sentinel species exposed to raw and treated urban wastewaters. Despite the need for development for quantitative approaches, these results support the importance of such invertebrate organisms, especially zebra mussel, for the active surveillance of pathogenic microorganisms and their indicators in environmental waters.

The genome of the zebra mussel, Dreissena polymorpha: a resource for comparative genomics, invasion genetics, and biocontrol.

The zebra mussel, Dreissena polymorpha, continues to spread from its native range in Eurasia to Europe and North America, causing billions of dollars in damage and dramatically altering invaded aquatic ecosystems. Despite these impacts, there are few genomic resources for Dreissena or related bivalves. Although the D. polymorpha genome is highly repetitive, we have used a combination of long-read sequencing and Hi-C-based scaffolding to generate a high-quality chromosome-scale genome assembly. Through comparative analysis and transcriptomics experiments, we have gained insights into processes that likely control the invasive success of zebra mussels, including shell formation, synthesis of byssal threads, and thermal tolerance. We identified multiple intact steamer-like elements, a retrotransposon that has been linked to transmissible cancer in marine clams. We also found that D. polymorpha have an unusual 67 kb mitochondrial genome containing numerous tandem repeats, making it the largest observed in Eumetazoa. Together these findings create a rich resource for invasive species research and control efforts.

Faster monitoring of the invasive alien species (IAS) Dreissena polymorpha in river basins through isothermal amplification.

Zebra mussel (Dreissena polymorpha) is considered as one of the 100 most harmful IAS in the world. Traditional detection methods have limitations, and PCR based environmental DNA detection has provided interesting results for early warning. However, in the last years, the development of isothermal amplification methods has received increasing attention. Among them, loop-mediated isothermal amplification (LAMP) has several advantages, including its higher tolerance to the presence of inhibitors and the possibility of naked-eye detection, which enables and simplifies its potential use in decentralized settings. In the current study, a real-time LAMP (qLAMP) method for the detection of Dreissena polymorpha was developed and tested with samples from the Guadalquivir River basin, together with two real-time PCR (qPCR) methods using different detection chemistries, targeting a specific region of the mitochondrial gene cytochrome C oxidase subunit I. All three developed approaches were evaluated regarding specificity, sensitivity and time required for detection. Regarding sensitivity, both qPCR approaches were more sensitive than qLAMP by one order of magnitude, however the qLAMP method proved to be as specific and much faster being performed in just 9 min versus 23 and 29 min for the qPCR methods based on hydrolysis probe and intercalating dye respectively.

Structure, function and parallel evolution of the bivalve byssus, with insights from proteomes and the zebra mussel genome.

The byssus is a structure unique to bivalves. Byssal threads composed of many proteins extend like tendons from muscle cells, ending in adhesive pads that attach underwater. Crucial to settlement and metamorphosis, larvae of virtually all species are byssate. By contrast, in adults, the byssus is scattered throughout bivalves, where it has had profound effects on morphological evolution and been key to adaptive radiations of epifaunal species. I compare byssus structure and proteins in blue mussels (Mytilus), by far the best characterized, to zebra mussels (Dreissena polymorpha), in which several byssal proteins have been isolated and sequenced. By mapping the adult byssus onto a recent phylogenomic tree, I confirm its independent evolution in these and other lineages, likely parallelisms with common origins in development. While the byssus is superficially similar in Dreissena and Mytilus, in finer detail it is not, and byssal proteins are dramatically different. I used the chromosome-scale D. polymorpha genome we recently assembled to search for byssal genes and found 37 byssal loci on 10 of the 16 chromosomes. Most byssal genes are in small families, with several amino acid substitutions between paralogs. Byssal proteins of zebra mussels and related quagga mussels (D. rostriformis) are divergent, suggesting rapid evolution typical of proteins with repetitive low complexity domains. Opportunities abound for proteomic and genomic work to further our understanding of this textbook example of a marine natural material. A priority should be invasive bivalves, given the role of byssal attachment in the spread of, and ecological and economic damage caused by zebra mussels, quagga mussels and others. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

Environmental (e)RNA advances the reliability of eDNA by predicting its age.

Environmental DNA (eDNA) analysis has advanced conservation biology and biodiversity management. However, accurate estimation of age and origin of eDNA is complicated by particle transport and the presence of legacy genetic material, which can obscure accurate interpretation of eDNA detection and quantification. To understand the state of genomic material within the environment, we investigated the degradation relationships between (a) size of fragments (long vs short), (b) genomic origins (mitochondrial vs nuclear), (c) nucleic acids (eDNA vs eRNA), and (d) RNA types (messenger (m)RNA vs ribosomal (r)RNA) from non-indigenous Dreissena mussels. Initial concentrations of eRNA followed expected transcriptional trends, with rRNAs found at > 1000 × that of eDNA, and a mitosis-associated mRNA falling below detection limits within 24 h. Furthermore, the ratio of eRNA:eDNA significantly decreased throughout degradation, potentially providing an estimate for the age of genomic material. Thus, eRNA quantification can increase detection due to the high concentrations of rRNAs. Furthermore, it may improve interpretation of positive detections through the eRNA:eDNA ratio and/or by detecting low abundant mitosis-associated mRNAs that degrade within ~ 24 h.

Environmental transcriptomes of invasive dreissena, a model species in ecotoxicology and invasion biology.

Dreissenids are established model species for ecological and ecotoxicological studies, since they are sessile and filter feeder organisms and reflect in situ freshwater quality. Despite this strong interest for hydrosystem biomonitoring, omics data are still scarce. In the present study, we achieved full de novo assembly transcriptomes of digestive glands to gain insight into Dreissena polymorpha and D. rostriformis bugensis molecular knowledge. Transcriptomes were obtained by Illumina RNA sequencing of seventy-nine organisms issued from fifteen populations inhabiting sites that exhibits multiple freshwater contamination levels and different hydrosystem topographies (open or closed systems). Based on a recent de novo assembly algorithm, we carried out a complete, quality-checked and annotated transcriptomes. The power of the present study lies in the completeness of transcriptomes gathering multipopulational organisms sequencing and its full availability through an open access interface that gives a friendly and ready-to-use access to data. The use of such data for proteogenomic and targeted biological pathway investigations purpose is promising as they are first full transcriptomes for this two Dreissena species.

The quagga mussel genome and the evolution of freshwater tolerance.

Freshwater dreissenid mussels evolved from marine ancestors during the Miocene ∼30 million years ago and today include some of the most successful and destructive invasive species of freshwater environments. Here, we sequenced the genome of the quagga mussel Dreissena rostriformis to identify adaptations involved in embryonic osmoregulation. We provide evidence that a lophotrochozoan-specific aquaporin water channel, a vacuolar ATPase subunit and a sodium/hydrogen exchanger are involved in osmoregulation throughout early cleavage, during which time large intercellular fluid-filled 'cleavage cavities' repeatedly form, coalesce and collapse, expelling excess water to the exterior. Independent expansions of aquaporins coinciding with at least five freshwater colonization events confirm their role in freshwater adaptation. Repeated aquaporin expansions and the evolution of membrane-bound fluid-filled osmoregulatory structures in diverse freshwater taxa point to a fundamental principle guiding the evolution of freshwater tolerance and provide a framework for future species control efforts.

Shedding light on eDNA: neither natural levels of UV radiation nor the presence of a filter feeder affect eDNA-based detection of aquatic organisms.

The use of environmental DNA (eDNA) as a species detection tool is attracting attention from both scientific and applied fields, especially for detecting invasive or rare species. In order to use eDNA as an efficient and reliable tool, however, we need to understand its origin and state as well as factors affecting its degradation. Various biotic and abiotic environmental factors have been proposed to affect degradation of eDNA in aquatic environments and thus to influence detection rates of species. Here, we were interested in two of them, namely UV light, which can break down DNA, and the presence of filter feeders, which can remove DNA and DNA-bound particles. A few, mostly laboratory-based studies have found minor effects of UVB on the degradation of eDNA. Ultraviolet A radiation (UVA), however, has been neglected although it also causes DNA lesions and is 10- to 100-fold more prevalent than UVB when reaching the earth's surface. Filter feeders are common in aquatic ecosystem, but their effects on eDNA has hitherto been ignored. We conducted a full-factorial aquatic mesocosm experiment under near-natural outdoor conditions manipulating UV radiation as well as the presence of Dreissena polymorpha, a strong filter feeder capable of filtering cells or organelles containing DNA. Surprisingly, we found that neither UV radiation nor the presence of the filter feeder affected eDNA-based detection rates of macroinvertebrates, even though the experiment took place in summer when UV radiation intensity and filtration activity is high for the chosen experimental site and conditions. These results, in combination with studies from marine or laboratory settings finding no effect of sunlight and its UV components on the detectability of eDNA, suggest that eDNA based species assessments could be relatively robust with respect to our two factors studied.

Differentiation of sympatric zebra and quagga mussels in ecotoxicological studies: A comparison of morphometric data, gene expression, and body metal concentrations.

The zebra mussel is among the best studied freshwater molluscs in ecotoxicology, but information on the quagga mussel is lacking. Considering its potential spread, we selected a river in France in which zebra and quagga mussels coexisted, and then we used genetic markers to differentiate the two species and compared morphological parameters. cDNA sequencing assays of ten genes already used in zebra mussels were performed on quagga mussels to obtain functional specific primers. Then we analyzed the expression of genes involved in cellular metabolic activities (Cytochrome-c-oxidase - cox, and ATP synthase - atp), detoxification processes (Glutathione-S-Transferase - gst), oxidative stress (Catalase - cat), and digestive functions (Amylase - amy) on the two species. Whereas morphometric analysis underlined similarities in shape between the two species, relative gene expression profiles and metal concentrations evidenced strong differences. Quagga mussels notably presented half as high concentrations in Cd and Pb, two particularly toxic elements, as zebra mussels. These results imply that i) particular attention should be paid to properly distinguish the two species considering their similar external appearance, and ii) zebra mussels cannot be replaced by quagga mussels in ecotoxicological studies without preliminary investigations on biomarker response patterns. To our knowledge, this study is the first to have undertaken such an approach in gene expression analysis in quagga mussels, and more generally to have compared such biomarker responses of zebra and quagga mussels in the field.

Isothermal amplification of environmental DNA (eDNA) for direct field-based monitoring and laboratory confirmation of Dreissena sp.

Loop-mediated isothermal amplification (LAMP) of aquatic invasive species environmental DNA (AIS eDNA) was used for rapid, sensitive, and specific detection of Dreissena sp. relevant to the Great Lakes (USA) basin. The method was validated for two uses including i) direct amplification of eDNA using a hand filtration system and ii) confirmation of the results after DNA extraction using a conventional thermal cycler run at isothermal temperatures. Direct amplification eliminated the need for DNA extraction and purification and allowed detection of target invasive species in grab or concentrated surface water samples, containing both free DNA as well as larger cells and particulates, such as veligers, eggs, or seeds. The direct amplification method validation was conducted using Dreissena polymorpha and Dreissena bugensis and uses up to 1 L grab water samples for high target abundance (e.g., greater than 10 veligers (larval mussels) per L for Dreissena sp.) or 20 L samples concentrated through 35 µm nylon screens for low target abundance, at less than 10 veligers per liter water. Surface water concentrate samples were collected over a period of three years, mostly from inland lakes in Michigan with the help of a network of volunteers. Field samples collected from 318 surface water locations included i) filtered concentrate for direct amplification validation and ii) 1 L grab water sample for eDNA extraction and confirmation. Though the extraction-based protocol was more sensitive (resulting in more positive detections than direct amplification), direct amplification could be used for rapid screening, allowing for quicker action times. For samples collected between May and August, results of eDNA direct amplification were consistent with known presence/absence of selected invasive species. A cross-platform smartphone application was also developed to disseminate the analyzed results to volunteers. Field tests of the direct amplification protocol using a portable device (Gene-Z) showed the method could be used in the field to obtain results within one hr (from sample to result). Overall, the direct amplification has the potential to simplify the eDNA-based monitoring of multiple aquatic invasive species. Additional studies are warranted to establish quantitative correlation between eDNA copy number, veliger, biomass or organismal abundance in the field.

Validated methodology for quantifying infestation levels of dreissenid mussels in environmental DNA (eDNA) samples.

The zebra mussel (Dreissena polymorpha Pallas, 1771) and the quagga mussel (D. rostriformis Deshayes, 1838) are successful invasive bivalves with substantial ecological and economic impacts in freshwater systems once they become established. Since their eradication is extremely difficult, their detection at an early stage is crucial to prevent spread. In this study, we optimized and validated a qPCR detection method based on the histone H2B gene to quantify combined infestation levels of zebra and quagga mussels in environmental DNA samples. Our results show specific dreissenid DNA present in filtered water samples for which microscopic diagnostic identification for larvae failed. Monitoring a large number of locations for invasive dreissenid species based on a highly specific environmental DNA qPCR assay may prove to be an essential tool for management and control plans focused on prevention of establishment of dreissenid mussels in new locations.

Seasonal and size-related variation of subcellular biomarkers in quagga mussels (Dreissena bugensis) inhabiting sites affected by moderate contamination with complex mixtures of pollutants.

The size-related differences in subcellular biomarker responses were assessed in Dreissena bugensis mussels inhabiting harbours moderately affected by pollution with complex mixtures of heavy metals and polycyclic aromatic hydrocarbons (PAHs). Adult D. bugensis samples were collected from three harbours of Lake Balaton (Hungary) characterized by moderate shipping activity, and as reference site, from a highly protected remote area of the lake. Biomarkers of exposure (metallothioneins (MTs), ethoxyresorufin-o-deethylase (EROD)), oxidative stress (lipid peroxidation (LPO), DNA strand breaks (DNAsb)) and possible endocrine disruption (vitellogenin-like proteins (VTG)) were analysed in whole-tissue homogenates of differently sized groups of mussels in relation to environmental parameters and priority pollutants (heavy metals and polycyclic aromatic hydrocarbons). Integrated biomarker response (IBR) indices were calculated for biomarker responses gained through in situ measurements to signalize critical sites and to better distinguish natural tendencies from biological effects of contaminants. Biomarker responses showed close positive correlation in case of MT, EROD, LPO, and DNAsb and negative correlation with VTG levels with mussel shell length in autumn, when higher levels of biomarkers appeared, possibly due to natural lifecycle changes of animals.

Use of a combined effect model approach for discriminating between ABCB1- and ABCC1-type efflux activities in native bivalve gill tissue.

Aquatic organisms, such as bivalves, employ ATP binding cassette (ABC) transporters for efflux of potentially toxic chemicals. Anthropogenic water contaminants can, as chemosensitizers, disrupt efflux transporter function enabling other, putatively toxic compounds to enter the organism. Applying rapid amplification of cDNA ends (RACE) PCR we identified complete cDNAs encoding ABCB1- and ABCC1-type transporter homologs from zebra mussel providing the molecular basis for expression of both transporter types in zebra mussel gills. Further, efflux activities of both transporter types in gills were indicated with dye accumulation assays where efflux of the dye calcein-am was sensitive to both ABCB1- (reversin 205, verapamil) and ABCC1- (MK571) type specific inhibitors. The assumption that different inhibitors targeted different efflux pump types was confirmed when comparing measured effects of binary inhibitor compound mixtures in dye accumulation assays with predictions from mixture effect models. Effects by the MK571/reversin 205 mixture corresponded better with independent action, whereas reversin 205/verapamil joint effects were better predicted by the concentration addition model indicating different and equal targets, respectively. The binary mixture approach was further applied to identify the efflux pump type targeted by environmentally relevant chemosensitizing compounds. Pentachlorophenol and musk ketone, which were selected after a pre-screen of twelve compounds that previously had been identified as chemosensitizers, showed mixture effects that corresponded better with concentration addition when combined with reversine 205 but with independent action predictions when combined with MK571 indicating targeting of an ABCB1-type efflux pump by these compounds.

Amphetamine exposure imbalanced antioxidant activity in the bivalve Dreissena polymorpha causing oxidative and genetic damage.

Illicit drugs have been recognized as emerging aquatic pollutants due to their presence in aquatic ecosystems up to µg/L level. Among these, the synthetic psycho-stimulant drug amphetamine (AMPH) is commonly found in both surface and wastewaters worldwide. Even though the environmental occurrence of AMPH is well-known, the information on its toxicity towards non-target freshwater organisms is completely lacking. This study investigated the imbalance of the oxidative status and both oxidative and genetic damage induced by a 14-day exposure to two concentrations (500 ng/L and 5000 ng/L) of AMPH on the freshwater bivalve Dreissena polymorpha by the application of a biomarker suite. We investigated the activity of antioxidant enzymes (SOD, CAT and GPx), the phase II detoxifying enzyme GST, the lipid peroxidation level (LPO) and protein carbonyl content (PCC), as well as primary (Single Cell Gel Electrophoresis assay) and fixed (DNA diffusion assay and Micronucleus test) genetic damage. Our results showed that a current realistic AMPH concentration (500 ng/L) did neither cause notable imbalances in enzymatic activities, nor oxidative and genetic damage to cellular macromolecules. In contrast, the bell-shaped trend of antioxidants showed at the highest tested concentration (5000 ng/L) suggested an overproduction of reactive oxygen species, leading to oxidative damage, as confirmed by the significant increase of protein carbonylation and DNA fragmentation.

DNA oxidation and DNA repair in gills of zebra mussels exposed to cadmium and benzo(a)pyrene.

Freshwater bivalve molluscs are considered as effective indicators of environmental pollution. The comet assay allows the detection of DNA damage such as DNA strand breaks and alkali-labile sites. The main oxidative lesion, 8-oxo-7,8-dihydro-2'-deoxyguanosine (8-oxodG), which is a pre-mutagenic lesion, can be detected by the comet assay coupled with the hOGG1 DNA repair enzyme. With this modified assay we recently observed that BaP induced 8-oxodG lesions and with the modified comet-Fpg assay we observed that Cd induced oxidative DNA damage. The aim of this study was to determine the stability of DNA lesions in Cd and BaP exposed zebra mussels using the comet-hOGG1 assay. Mussels were exposed for 24 h to these two chemicals and then placed in clean water for 6 days. We observed that BaP (7, 12 and 18 µg/L) induced an increase of DNA strand break levels as soon as 6 h of exposure and that the two highest concentrations of BaP induced a low level of hOGG1-sensitive sites. After 2 days of depuration, BaP induced DNA lesions returned to the basal level, indicating an effective DNA repair. Cd (3, 32 and 81 µg/L) induced an increase of the DNA strand break levels and a low level of hOGG1-sensitive sites. This study revealed that BaP-induced DNA lesions are repaired more efficiently than Cd-induced DNA lesions. As the level of hOGG1 sensitive sites was increased in Cd and BaP exposed mussels, it seems that these chemicals induce 8-oxo-dG.

qRT-PCR evaluation of the transcriptional response of zebra mussel to heavy metals.

BACKGROUND: The transcriptional response of adult zebra mussels (Dreissena polymorpha) to heavy metals (mercury, copper, and cadmium) was analyzed by quantitative Real-Time Polymerase Chain Reaction (qRT-PCR) to study the coordinated regulation of different metal-, oxidative stress- and xenobiotic defence-related genes in gills and digestive gland. Regulatory network analyses allowed the comparison of this response between different species and taxa. RESULTS: Chemometric analyses allowed identifying the effects of these metals clearly separating control and treated samples of both tissues. Interactions between the different genes, either in the same or between both tissues, were analysed to identify correlations and to propose stress-related genes' regulatory networks. These networks were finally compared with existing data from human, mouse, zebrafish, Drosophila and the roundworm to evaluate their mechanistically-known response to metals (and to stressors in general) with the correlations observed in the still poorly-known, invasive zebra mussel. CONCLUSIONS: Our analyses found a general conservation of regulation genes and of their interactions among the different considered species, and may serve as a guide to extrapolate regulatory data from model species to lesser-known environmentally (or medically) relevant species.

Using Massive Parallel Sequencing for the development, validation, and application of population genetics markers in the invasive bivalve zebra mussel (Dreissena polymorpha).

The zebra mussel (Dreissena polymorpha, Pallas, 1771) is one of the most invasive species of freshwater bivalves, due to a combination of biological and anthropogenic factors. Once this species has been introduced to a new area, individuals form dense aggregations that are very difficult to remove, leading to many adverse socioeconomic and ecological consequences. In this study, we identified, tested, and validated a new set of polymorphic microsatellite loci (also known as SSRs, Single Sequence Repeats) using a Massive Parallel Sequencing (MPS) platform. After several pruning steps, 93 SSRs could potentially be amplified. Out of these SSRs, 14 were polymorphic, producing a polymorphic yield of 15.05%. These 14 polymorphic microsatellites were fully validated in a first approximation of the genetic population structure of D. polymorpha in the Iberian Peninsula. Based on this polymorphic yield, we propose a criterion for establishing the number of SSRs that require validation in similar species, depending on the final use of the markers. These results could be used to optimize MPS approaches in the development of microsatellites as genetic markers, which would reduce the cost of this process.

Genotoxicity and activation of cellular defenses in transplanted zebra mussels Dreissena polymorpha along the Seine river.

The aim of the present study was to confirm the relevance of studying DNA adduct formation in a field study. In that context, freshwater mussels Dreissena polymorpha, collected in a reference station, were transplanted in different sites with a pollution gradient. After one and two months, mussels were collected and DNA adduct formation was analyzed using the (32)P post labelling technique on both gills and digestive glands. In addition, the expression of genes involved in the detoxification system (catalase (CAT), superoxide dismutase (SOD), glutathione S-transferase (GST), HSP70, aryl hydrocarbon receptor (AHR), P glycoprotein (PgP), metallothionein (MT)) was assessed by RT-PCR. DNA adducts were observed at amount comparable to data from literature. Increase of DNA adducts after two months of transplantation could be correlated with strong modulation of gene expression implicated in detoxification processes. Indeed, PgP and HSP70 gene expressions were similarly induced in gills and digestive glands while SOD and CAT expressions were down regulated in both tissues. AHR, GST and MT genes were differently regulated depending upon the tissue studied and the level of contamination in the different sites. We demonstrated that mussels transplanted in the different stations with pollution gradient were able to biotransform PAHs, assessed by DNA adduct formation and the high decrease of detoxification genes. Specific DNA adducts pattern obtained after one and two month mussel transplantations demonstrated the relevance of DNA adduct as biomarker of environmental pollution.

DNA adduct formation and induction of detoxification mechanisms in Dreissena polymorpha exposed to nitro-PAHs.

Derived polycyclic aromatic hydrocarbons (PAHs) such as nitro-PAHs are present in the environment and are known to be much more toxic than PAHs compounds. However, very few studies have analysed their effects on the aquatic environment and none have investigated the freshwater environment. In the present study, we determined whether 1-nitropyrene (1-NP), a model of nitro-PAHs, can induce DNA adducts in gills and digestive glands of the freshwater mussel Dreissena polymorpha. Two concentrations of 1-NP (50 and 500 µM) were tested. In addition, in order to understand the metabolic pathways involved in 1-NP genotoxicity, mRNA expression of genes implicated in biotransformation mechanisms was assessed by quantitative reverse transcription-PCR. Results showed the presence of DNA adducts in both gills and digestive glands, with highest levels obtained after 5 days of exposure to 500 µM. Metallothionein mRNA levels were enhanced in digestive glands exposed to 50 µM. Surprisingly, at the higher concentration (500 µM), aryl hydrocarbon receptor and HSP70 genes were only up-regulated in digestive glands while PgP mRNA levels were increased in both tissues. Results suggested a cytotoxic and genotoxic effect of 1-NP. Mussels seemed to be able to partially detoxify this compound, in view of the low amount of DNA adducts observed after 5 days exposure to 50 µM. For the first time, 1-NP biotransformation and detoxification systems have been characterised in D. polymorpha.

Effects of freshwater pollution on the genetics of zebra mussels (Dreissena polymorpha) at the molecular and population level.

Revealing long-term effects of contaminants on the genetic structure of organisms inhabiting polluted environments should encompass analyses at the population, molecular, and cellular level. Following this concept, we studied the genetic constitution of zebra mussel populations from a polluted (Dp) and reference sites (Cl) at the river Drava, Croatia, and applied microsatellite and DNA damage analyses (Comet assay, micronucleus test (MNT)). Additionally, mussels from both populations were exposed to polluted wastewater in the laboratory for three days, and DNA damage was analyzed to evaluate acclimatization and genetic adaptation of the investigated populations to the polluted environment. The two populations differed in their genetic constitution. Microsatellite analysis suggested that Dp had undergone a genetic bottleneck. Comet assay did not indicate any difference in DNA damage between the two populations, but MNT revealed that Dp had an increased percentage of micronuclei in hemocytes in comparison to Cl. The laboratory experiment revealed that Dp had a lower percentage of tail DNA and a higher percentage of micronuclei than Cl. These differences between populations were possibly caused by an overall decreased fitness of Dp due to genetic drift and by an enhanced DNA repair mechanism due to acclimatization to pollution in the source habitat.