Phylogeographic Substructuring in the Southernmost Refugium of the European Common Frog Rana temporaria.

Rana temporaria is one of the most widespread Palearctic brown frogs. We aimed to clarify distribution pattern of two main genetic clades in the understudied Balkan peninsula by using 16SrRNA and MT-CYTB sequences, already widely applied in analyses of populations from other parts of Europe, while focusing on the broad area along the Morava river (central Balkans) as a known gap in the species distribution. Additionally, we were interested in revealing the extent of haplotype diversity within the main genetic clades in the Balkans, particularly around the supposed suture zone. The results revealed a suture zone between the Western and Eastern Clades in the central part of the Balkan Peninsula. This indicated the existence of a historical barrier between the Balkan Mountain Belt and geographically close mountains surrounding the Vlasina Plateau (Rhodope/Serbian-Macedonian Massif). The overall observed haplotype diversity in populations of R. temporaria from the Balkan Peninsula seems high. Harboring both main genetic clades of R. temporaria qualifies the Balkan Peninsula as another important center of species' genetic diversity, as well as rich in unique haplotypes. This points out the necessity of applying conservation measures focused on the common European frog populations and habitats in this part of the species' distribution area.

In situ detection of the genotoxic potential as one of the lines of evidence in the weight-of-evidence approach-the Joint Danube Survey 4 Case Study.

Environmental studies which aim to assess the ecological impact of chemical and other types of pollution should employ a complex weight-of-evidence approach with multiple lines of evidence (LoEs). This study focused on in situ genotoxicological methods such as the comet and micronucleus assays and randomly amplified polymorphic DNA analysis as one of the multiple LoEs (LoE3) on the fish species Alburnus alburnus (bleak) as a bioindicator. The study was carried out within the Joint Danube Survey 4 (JDS4) at nine sites in the Danube River Basin in the Republic of Serbia. Out of nine sampling sites, two were situated at the Tisa, Sava, and Velika Morava rivers, and three sites were at the Danube River. The three additionally employed LoEs were: SumTUwater calculated based on the monitoring data in the database of the Serbian Environmental Protection Agency (SEPA) (LoE1); in vitro analyses of JDS4 water extracts employing genotoxicological methods (LoE2); assessment of the ecological status/potential by SEPA and indication of the ecological status for the sites performed within the JDS4 (LoE4). The analyzed biomarker responses in the bleak were integrated into the unique integrated biomarker response index which was used to rank the sites. The highest pollution pressure was recorded at JDS4 39 and JDS4 36, while the lowest was at JDS4 35. The impact of pollution was confirmed at three sites, JDS4 33, 40, and 41, by all four LoEs. At other sampling sites, a difference was observed regarding the pollution depending on the employed LoEs. This indicates the importance of implementing a comprehensive weight-of-evidence approach to ensure the impact of pollution is not overlooked when using only one LoE as is often the case in environmental studies.

Wastewater-based epidemiology in countries with poor wastewater treatment - Epidemiological indicator function of SARS-CoV-2 RNA in surface waters.

Wastewater-based epidemiology (WBE) surveillance of COVID-19 and other future outbreaks is a challenge for developing countries as most households are not connected to a sewerage system. In December 2019, SARS-CoV-2 RNA was detected in the Danube River at a site severely affected by wastewaters from Belgrade. Rivers are much more complex systems than wastewater systems, and efforts are needed to address all the factors influencing the adoption of WBE as an alternative to targeting raw wastewater. Our objective was to provide a more detailed insight into the potential of SARS-CoV-2 surveillance in Serbian surface waters for epidemiological purposes. Water samples were collected at 12 sites along the Sava and Danube rivers in Belgrade during the fourth COVID-19 wave in Serbia that started in late February 2021. RNA was concentrated using Amicon Ultra-15 centrifugal filters and quantified using RT-qPCR with primer sets targeting nucleocapsid (N1 and N2) and envelope (E) protein genes. Microbiological (faecal indicator bacteria and human and animal genetic faecal source tracking markers), epidemiological, physicochemical and hydromorphological parameters were analysed in parallel. From 44 samples, SARS-CoV-2 RNA was detected in 31, but only at 4 concentrations above the level of quantification (ranging from 8.47 × 103 to 2.07 × 104 gc/L). The results indicated that surveillance of SARS-CoV-2 RNA in surface waters as ultimate recipients could be used as an epidemiological early-warning tool in countries lacking wastewater treatment and proper sewerage infrastructure. The performance of the applied approach, including advanced sampling site characterization to trace and identify sites with significant raw sewage influence from human populations, could be further improved by adaptation of the methodology for processing higher volumes of samples and enrichment factors, which should provide the quantitative instead of qualitative data needed for WBE.

Detection of SARS-CoV-2 RNA in the Danube River in Serbia associated with the discharge of untreated wastewaters.

In Serbia less than 13% of collected municipal wastewaters is being treated before their release in the environment. This includes all municipal wastewater discharges from Belgrade (capital city of Serbia; population 1,700,000). Previous research has identified the impacts of raw wastewater discharges from Belgrade on the Danube River, and this study investigated if such discharges also provided a pathway for SARS-CoV-2 RNA material. Samples were collected during the most critical circumstances that occurred so far within the COVID-19 pandemics in Serbia. Grab and composite samples were collected in December 2020, during the peak of the third wave (in terms of reported cases) at the site which receives the wastewater loads in Belgrade. Grab samples collected upstream and downstream of Belgrade were also analyzed. RNA was quantified using RT-qPCR with primer sets targeting nucleocapsid (N1 and N2) and envelope (E) protein genes. SARS-CoV-2 RNA (5.97 × 103 to 1.32 × 104 copies/L) was detected only in samples collected at the site strongly impacted by the wastewaters where all three applied primer sets gave positive signals. Determined concentrations correspond to those reported in wastewater influents sampled at treatment plants in other countries indicating an epidemiological indicator function of used approach for rivers with high pollution loads in countries with poor wastewater treatment.

Unravelling the effects of multiple stressors on diatom and macroinvertebrate communities in European river basins using structural and functional approaches.

Rivers suffer from more severe decreases in species diversity compared to other aquatic and terrestrial ecosystems due to a variety of pressures related to human activities. Species provide different roles in the functioning of the ecosystem, and their loss may reduce the capacity of the ecosystems to respond to multiple stressors. The effects on diversity will differ based on the type, combination and severity of stressors, as well as on the characteristics of the community composition and tolerance. Multiple trait-based approaches (MTBAs) can help to unravel the effects of multiple stressors on communities, providing a mechanistic interpretation, and, thus, complementing traditional biodiversity assessments using community structure. We studied the relationships between diversity indexes and trait composition of macroinvertebrate and diatom communities, as well as environmental variables that described the hydrological and geomorphological alterations and toxic pollution (pesticides and pharmaceuticals) of three different European river basins: the Adige, the Sava, and the Evrotas. These river basins can be considered representative cases of different situations in European freshwater systems. Hydrological variables were the main drivers determining the community structure and function in the rivers, for both diatoms and macroinvertebrates. For diatom communities, pharmaceutical active compound (PhAC) toxic units were also identified as a very important driver of diversity changes, explaining up to 57% of the variance in taxonomic richness. For macroinvertebrates, river geomorphology was an important driver of structural changes, particularly affecting Plecoptera richness. In addition, PhAC and pesticide toxic units were also identified as stressors for macroinvertebrate communities. MTBA provided a detailed picture of the effects of the stressors on the communities and confirmed the importance of hydrological variables in shaping the functional attributes of the communities.

Sensitivity of Bleak (Alburnus alburnus) in Detection of the Wastewater Related Pressure in Large Lowland Rivers.

This study deals with bleak (Alburnus alburnus) sensitivity in detecting of the wastewater related pressure in large lowland rivers. The major objective was to investigate if the response measured in bleak should be linked to a certain stretch of the river and characterised as "stretch specific", or it should be linked to the sampling site and characterised as "site specific". The response was evaluated via condition index, metal pollution index, DNA damage and cell viability using integrated biomarker response approach. The study was conducted at 3 sub-sites characterized by different pollution levels in a relatively short stretch (2 km) of the Sava River (Serbia). Results indicated that the response of the biomarkers in bleak can be interpreted as "site specific". Among the studied biomarkers, DNA damage assessed by comet assay and micronucleus test has shown high sensitivity in differentiation of the sites.

Selection of assay, organism, and approach in biomonitoring significantly affects the evaluation of genotoxic potential in aquatic environments.

In this study, few different evaluation concepts were used for the assessment of genotoxic potential at the stretch of the Danube River identified as a significant hotspot of pollution originated through the untreated wastewaters. Three sites were chosen: one site upstream of the wastewater outlet in Novi Sad (Serbia), one at the outlet of wastewaters, and one site few kilometer downstream. Ex situ approach comprised prokaryotic SOS/umuC test on Salmonella typhimurium TA1535/pSK1005 and comet assay on human hepatoma cell line (HepG2). In situ approach was based on the active monitoring (cage approach) using freshwater mussels Sinanodonta woodiana and fish Cyprinus carpio. The comet and micronucleus assays were selected for evaluation of DNA damage in mussel haemocytes and fish blood cells. Within the ex situ part of the study, our results indicated that the eukaryotic model system is more sensitive compared to the prokaryotic one. In situ bioassays are recommended for obtaining a better insight into ecosystem status and in the case of our study the complete insight of genotoxic pressure. However, the choice of animals as bioindicators also has a significant impact on the quality of the obtained information. Differential response between fish and mussels was observed at the highly polluted site suggesting possible involvement of additional protective mechanism such as valve closure in mussels.

Simulium reptans (Linnaeus, 1758) and Simulium reptantoides Carlsson, 1962 from the Balkan Peninsula.

Simulium reptans (Linnaeus, 1758) and Simulium reptantoides Carlsson, 1962 are two species of the Simulium reptans group whose distribution is unclear because of their confusing taxonomy and systematics. Their genetic variability is well known for populations in northern and central Europe and shows that both species have two forms; however, the genetic variability of these species in southern and eastern Europe is unknown. To identify the status of these two species in southeast Europe, mtDNA was extracted from 19 individuals from 12 localities across the Balkan Peninsula. Phylogenetic analysis confirmed the existence of two species with 7.38-7.94% divergence. Each species was comprised of two clades, with 2.31% and 1.43% interclade divergence for S. reptans and S. reptantoides, respectively. This study revealed the presence of both species across the Balkans and that S. reptans occurs in this area in only one form (S. reptans B), while S. reptantoides is found in two genetic forms (A and B).

Phylogeography of Simulium Subgenus Wilhelmia (Diptera: Simuliidae)-Insights From Balkan Populations.

Many morphologically similar species of the simuliid (Diptera: Simuliidae) subgenus Wilhelmia, Enderlein are difficult to distinguish. Thus, the revision of the subgenus using various morphological, cytogenetic, and genetic analyses has been attempted. Neglected until now, the Balkan Peninsula, a crossroad between Europe and Anatolia, provides insight which could resolve problematic interrelationships of the taxa within this subgenus. To uncover the status and relations within the subgenus Wilhelmia, mtDNA was extracted from 47 individuals of six morphospecies: Simulium balcanicum (Enderlein, 1924), Simulium turgaicum Rubtsov, 1940, Simulium lineatum (Meigen, 1804), Simulium pseudequinum Séguy, 1921, Simulium equinum (Linnaeus, 1758), and Simulium paraequinum Puri, 1933 from 21 sites throughout the Balkan Peninsula. Phylogenetic analysis of the Wilhelmia species using mitochondrial DNA barcoding (COI) gene showed two major branches, the lineatum branch, which includes the lineages sergenti, paraequinum, and lineatum, and the equinum branch. In the equinum branch, the mtDNA sequences formed six clades, with high genetic distances, suggesting the existence of different species. Historically, the clades of the equinum branch appeared at numerous islands, perhaps as a result of allopatric speciation. The paraequinum lineage (lineatum branch) is composed of two species. However, six clades of the lineatum lineage overlapped with intra- and interspecific genetic distances. Our results revealed that the species S. balcanicum, S. pseudequinum B, and S. equinum were omnipresent in the Balkans. The results point to not only the fair diversity of Wilhelmia species in the Balkans, but also indicate that most Wilhelmia species live in sympatry.

Potentially toxic elements in muscle tissue of different fish species from the Sava River and risk assessment for consumers.

Fish from the Sava River are consumed daily by the local people: therefore, concern has been raised about the health implications of eating contaminated fish. In the present study, potentially toxic elements (PTE), such as Zn, Cu, Cr, Ni, Cd, Pb, As, Hg, and methylmercury (MeHg), were determined in ichthyofauna that are commonly consumed. PTE were determined in the fish muscle tissue. Fish were sampled at 12 locations from the source of the Sava River to its confluence with the Danube River during two sampling campaigns, namely; in 2014 under high water conditions and in 2015 under normal water conditions. Due to the different water regimes, different fish species were collected for chemical analysis. We observed that the concentrations of elements analysed in the fish muscle tissue were generally very low, except for those of Hg. Moreover, more than 90% of Hg present in the fish was in its most toxic form, namely MeHg. Especially in fish from the 2015 sampling campaign, Hg and MeHg concentrations increased with fish size, trophic level, and in the downstream direction. In addition, for Pb and As, and to some extent for Cd and Cr, spatial differences were detected in both years. The highest concentrations of PTE were detected in fish from sites with intensive industrial and agricultural activities. The consumption of fish in general does not pose a health risk for the PTE studied, except for Hg/MeHg at selected contaminated sites.

Long-term analysis of fish assemblage structure in the middle section of the Sava River - The impact of pollution, flood protection and dam construction.

At the beginning of the 20th century, the middle section of the Sava River in Croatia was unaffected by major human activities and rich in ichthyofauna. The Sava River was important for commercial and recreational fishing for the local population, which still remains today. However, the 1920s mining industry was established in Slovenia, which emitted carbon dust into the Sava River. At the same time, the construction of embankments to mitigate flooding started in the middle section. Furthermore, in the 1980s, the Krsko nuclear power plant (NPP), and in the 2010s, the Krsko hydropower plant (HPP) were built in Slovenia. These activities could have an impact on the composition of fish communities downstream from the major sources of disturbances. Therefore, the main aim of this paper were to analyze the changes in fish assemblages of the Sava River from 1978 to 2017, prior to and after the construction of Krsko NPP and HPP at the Medsave site on the Sava River, 20 km downstream from the major construction operations. Collected data were divided into four sampling periods (SP): SP1, from 1978 to 1980; SP2, from 1991 to 1994; SP3, from 2001 to 2006, and SP4 from 2011 to 2017. Besides alien fish species, water quality and hydromorphological modifications were identified as significant stressors. In SP1 and SP2 limnophilic and eurytopic fish groups were predominant, and 26 different fish species were identified, but in SP3 and SP4 rheophilic fish groups become dominant, and the diversity has declined to 21 species. Threatened species blageon, Telestes souffia seems to be missing from the main course of the Sava River in last 20 years. It can be concluded that disturbances in the fish assemblage pattern have coincided with the presence of multiple stressors of human origin.

Levels of regulated POPs in fish samples from the Sava River Basin. Comparison to legislated quality standard values.

Fish samples of different species (i.e. rainbow trout (Onchorhynchus mykiss), barbel (Barbus barbus) and European chub (Squalius cephalus)) were collected from the Sava River Basin for a preliminary investigation of the levels of PCDD/Fs, PCBs, PBDEs and PFAS as a whole. Concentrations of PCDD/Fs, in terms of pg WHO-TEQ/g ww, were below the maximum limit established at the Commission Regulation (EU) No 1259/2011. On the contrary, when DL-PCBs were also included, levels increase up to 11.7 pg WHO-TEQPCDD/Fs+DL-PCBs/g ww in a particular case, with two samples out of a total of ten exceeding the maximum set at this EU Regulation and the EQS established at the European Directive regarding priority substances in the field of water policy (0.0065 ng WHO-TEQPCDD/Fs+DL-PCBs/g ww). A similar trend was also observed for NDL-PCBs, whit the same two samples, from the lower stretch of the river basin, exceeding the maximum limit allowed at the EU Regulation (125 ng/g ww). For PBDEs, levels found in all the samples exceeded the EQS (0.0085 ng/g ww) up to more than a thousand times and 40% of the samples presented PFOS values above the EQS. Data from this study were compared to values reported at the literature for fish from other geographical areas.

Evaluation of genotoxic potential in the Velika Morava River Basin in vitro and in situ.

The Velika Morava River is the greatest national Serbian river and the significant tributary of the Danube River. The major problems in the Velika Morava River Basin (VMRB) represent untreated industrial and municipal wastewaters. In this study, the level of genotoxic potential at the sites along the VMRB was evaluated by parallel in vitro and in situ approach. Within in vitro testing, genotoxicity of native water samples collected from the sites in VMRB was evaluated by SOS/umuC test on Salmonella typhimurium TA1535/pSK1002 and by the comet assay on HepG2 cells. DNA damage in situ was assessed in bleak (Alburnus alburnus) erythrocytes by the comet (alkaline and Fpg-modified comet) and micronucleus assays. Additionally, the concentration of heavy metals in fish tissue was measured and this data, compiled with the data of the physico-chemical parameters measured in water, was used as a measure of the pollution pressure at the sites. Results showed that applied in vitro tests with native water samples are less sensitive in comparison with in situ tests and should be taken with precaution when making predictions on the status of the ecosystem. Within applied battery of in situ assays differential sensitivity of assays was observed where alkaline comet assay showed the highest potential in differentiation of the sites based on genotoxic potential. Integrated biomarker response showed that usage of the battery of bioassays provides better insight in a genotoxic effects in animals, and consequently, that the holistic approach is more suitable for this type of study.

The potential of chironomid larvae-based metrics in the bioassessment of non-wadeable rivers.

The chironomid community in non-wadeable lotic systems was tested as a source of information in the construction of biological metrics which could be used into the bioassessment protocols of large rivers. In order to achieve this, we simultaneously patterned the chironomid community structure and environmental factors along the catchment of the Danube and Sava River. The Self organizing map (SOM) recognized and visualized three different structural types of chironomid community for different environmental properties, described by means of 7 significant abiotic factors (a multi-stressor gradient). Indicator species analysis revealed that the chironomid taxa most responsible for structural changes significantly varied in their abundance and frequency along the established environmental gradients. Out of 40 biological metrics based on the chironomid community, the multilayer perceptron (MLP), an supervised type of artificial neural network, derived 5 models in which the abundance of Paratrichocladius rufiventis, Orthocladiinae, Cricotopus spp., Cricotopus triannulatus agg. and Cricotopus/Orthocladius ratio achieved a significant relationship (the r Pearson's linear correlation coefficient>0.7) with the multi stressor environment. The sensitivity analysis "partial derivatives" (PaD) method showed that all 5 biological metrics within the multi-stressor gradient were mostly influenced by dissolved organic carbon (DOC). Despite short and monotonous environmental gradients and the absence of reference conditions, the chironomid community structure and biological metrics predictably changed along the multistress range, showing a great potential for the bioassessment of large rivers.

Contribution to the improvement of diatom-based assessments of the ecological status of large rivers - The Sava River Case Study.

The Sava River Basin is a major drainage basin of southeastern Europe, significantly influenced by anthropogenic activities. Our study was focused on diatom communities as an indicator of the ecological status of running waters. We investigated over 937km of the Sava River at 19 sampling sites. Benthic diatom communities and 17 diatom indices were analyzed along with a large set of environmental parameters. CCA revealed that the most important elements along the spatial gradient were As and Si. Our results show that the species Navicula recens (Lange-Bert.) Lange-Bertalot and Eolimna minima (Grunow) Lange-Bertalot are very abundant at downstream localities where the highest concentrations of As were measured. The number of motile diatoms increased along the nutrient gradient, i.e. with Si availability. Correlations between diatom indices and selected environmental factors showed that temperature, As, Si and Fe are in significant negative correlation with most diatom indices. Analysis revealed the influence of As and metals in water on diatoms, although their concentrations did not exceed environmental standards. While our findings do not confirm that diatom indices reveal the intensity of pressures solely caused by nutrient and/or organic pollutants, they suggest that in moderately polluted large rivers benthic diatoms are good bioindicators of multiple pressures, and that diatom indices could serve as indicators of the level of overall degradation of an ecosystem.

The impact of multiple stressors on the biomarkers response in gills and liver of freshwater breams during different seasons.

Biomarkers attract increasing attention in environmental studies, as a tool for detection of exposure and effects of pollution, from both natural and anthropogenic sources. This study aims to assess the impact of multiple stressors during distinctive seasons, covering also extreme hydrological events (extensive flooding in the mid May 2014), on different levels of biological organization in the liver and gills of three closely related freshwater breams. Our previous study on DNA damage in blood cells of these specimens showed increased DNA damage in June 2014, one month after the flooding event. As a continuation of that research, the present study was conducted. As a biomarker of exposure DNA damage was measured by applying the alkaline comet assay, while histopathological alterations were monitored as a biomarker of effect. Additionally, concentrations of metals and metalloids in gills, liver and muscle were assessed. Sampling of fish tissues was performed in 2014, during winter (January and February), spring (March and early June) and summer (late June, July and August). Significant seasonal difference in DNA damage was observed for both tissues. During spring and summer the level of DNA damage in gills was significantly higher when compared to the liver. Histopathological analyses showed higher frequency of alterations in gills during spring, and in liver during summer, but without a significant seasonal difference. Gills had the highest concentration of metals and metalloids during the spring and summer, and liver during winter. Muscle was the least affected tissue during all three seasons. This study highlighted the importance of the multiple biomarker approach and the use of different fish tissues in assessment of surface water pollution.

River ecosystem processes: A synthesis of approaches, criteria of use and sensitivity to environmental stressors.

River ecosystems are subject to multiple stressors that affect their structure and functioning. Ecosystem structure refers to characteristics such as channel form, water quality or the composition of biological communities, whereas ecosystem functioning refers to processes such as metabolism, organic matter decomposition or secondary production. Structure and functioning respond in contrasting and complementary ways to environmental stressors. Moreover, assessing the response of ecosystem functioning to stressors is critical to understand the effects on the ecosystem services that produce direct benefits to humans. Yet, there is more information on structural than on functional parameters, and despite the many approaches available to measure river ecosystem processes, structural approaches are more widely used, especially in management. One reason for this discrepancy is the lack of synthetic studies analyzing river ecosystem functioning in a way that is useful for both scientists and managers. Here, we present a synthesis of key river ecosystem processes, which provides a description of the main characteristics of each process, including criteria guiding their measurement as well as their respective sensitivity to stressors. We also discuss the current limitations, potential improvements and future steps that the use of functional measures in rivers needs to face.

Characterization of sections of the Sava River based on fish community structure.

Sampling was undertaken, with the same fishing gear and along the Sava River, from its source to its confluence, in September 2014 and September 2015. In total, 44 fish species were identified, of which 37 were native species and 7 were alien. Fish samples revealed independence in terms of both species composition and their abundance under different hydrological conditions. During flooding and high water levels in 2014, pelagic fish species were sampled in greater proportion than at lower water levels in 2015 when benthic fish species were more abundant. The flood wave in 2014 was accompanied by catch of common carp, Cyprinus carpio, a typical lower rhithron fish species in the upper course, and of tench, Tinca tinca, a typical potamon fish species of backwaters, in the main channel of the lower Sava River. One specimen of bighead goby, Ponticola kessleri, which is common in the potamon fish community, was caught during the 2015 sampling close to the boundary between the upper and middle sections of the Sava. This is the first record of Ponto-Caspian gobies in the inland waters of Slovenia. Its finding far upstream indicates a strong effect of an as yet unidentified stress along the Sava River up to the spot where the bighead goby was sampled. Finally, these results indicate that pelagic fish species are more resistant to the stressful effect of flooding than benthic species, and that the structure of fish communities is influenced/affected by flooding as a short-term stressor. The progressively increasing number of alien fish species downstream in the Sava River point to the effects of long-term human-induced stressors in the area.

Conservation status of freshwater mussels in Europe: state of the art and future challenges.

Freshwater mussels of the Order Unionida provide important ecosystem functions and services, yet many of their populations are in decline. We comprehensively review the status of the 16 currently recognized species in Europe, collating for the first time their life-history traits, distribution, conservation status, habitat preferences, and main threats in order to suggest future management actions. In northern, central, and eastern Europe, a relatively homogeneous species composition is found in most basins. In southern Europe, despite the lower species richness, spatially restricted species make these basins a high conservation priority. Information on freshwater mussels in Europe is unevenly distributed with considerable differences in data quality and quantity among countries and species. To make conservation more effective in the future, we suggest greater international cooperation using standardized protocols and methods to monitor and manage European freshwater mussel diversity. Such an approach will not only help conserve this vulnerable group but also, through the protection of these important organisms, will offer wider benefits to freshwater ecosystems.

Evaluation of Genotoxic Pressure along the Sava River.

In this study we have performed a comprehensive genotoxicological survey along the 900 rkm of the Sava River. In total, 12 sites were chosen in compliance with the goals of GLOBAQUA project dealing with the effects of multiple stressors on biodiversity and functioning of aquatic ecosystems. The genotoxic potential was assessed using a complex battery of bioassays performed in prokaryotes and aquatic eukaryotes (freshwater fish). Battery comprised evaluation of mutagenicity by SOS/umuC test in Salmonella typhimurium TA1535/pSK1002. The level of DNA damage as a biomarker of exposure (comet assay) and biomarker of effect (micronucleus assay) and the level of oxidative stress as well (Fpg-modified comet assay) was studied in blood cells of bleak and spirlin (Alburnus alburnus/Alburnoides bipunctatus respectively). Result indicated differential sensitivity of applied bioassays in detection of genotoxic pressure. The standard and Fpg-modified comet assay showed higher potential in differentiation of the sites based on genotoxic potential in comparison with micronucleus assay and SOS/umuC test. Our data represent snapshot of the current status of the river which indicates the presence of genotoxic potential along the river which can be traced to the deterioration of quality of the Sava River by communal and industrial wastewaters. The major highlight of the study is that we have provided complex set of data obtained from a single source (homogeneity of analyses for all samples).