Research Article |
|
Corresponding author: Erin S. McCallum ( erin.mccallum@slu.se ) Academic editor: Markéta Ondračková
© 2025 Erin S. McCallum, Kristina M. Sefc, Tomas Brodin, Patricia Burkhardt-Holm, Karen Bussmann-Charran, Ann-Britt Florin, Juergen Geist, Michal Janáč, Pavel Jurajda, Jake M. Martin, Joachim Pander, Aneesh P. H. Bose.
This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Citation:
McCallum ES, Sefc KM, Brodin T, Burkhardt-Holm P, Bussmann-Charran K, Florin A-B, Geist J, Janáč M, Jurajda P, Martin JM, Pander J, Bose APH (2025) Round goby population differentiation across river barriers in Central Europe. Aquatic Invasions 20(3): 355-370. https://doi.org/10.3391/ai.2025.20.3.152950
|
River barriers such as hydropower dams and weirs can negatively affect river ecosystems by disrupting connectivity and reducing biodiversity. However, such barriers could also limit the spread of invasive species. Here, we used a spatial population genetics approach to test whether river barriers act as a hindrance to gene flow in the invasive round goby (Neogobius melanostomus Pallas, 1814). We sampled gobies from four different rivers across their invasive range in Central Europe (the Danube, Dyje, Morava, and Rhine rivers), with locations on either side of eight major river barriers. Using microsatellite genotyping, we found that round goby populations were differentiated with increasing number of river barriers and with increasing distance between sampling sites, depending on the river system in focus. We found significant population differentiation across three individual barriers, but no clear indication that this was related to barrier type as barriers were highly diverse. We also found reduced genetic diversity in populations that were more recently established. Our findings suggest that successive river barriers can sometimes slow the spread of round goby. Further research on the features of barriers that hinder round goby movement will help to design barrier passage solutions that will both limit spread of this invasive species and maintain connectivity for the native fauna.
Aquatic habitat invasion, dispersal, genetic differentiation, hydropower dam, invasive species, Neogobius melanostomus, river connectivity
Aquatic invasive species are a significant driver of global biodiversity loss, especially in freshwaters (
The round goby is considered a prime model of a highly invasive species (
To assess the extent to which river barriers restrict round goby dispersal, we have used a spatial population genetics approach and estimated pair-wise population differentiation as a proxy of gene flow. This approach captures the effects of migration (or lack thereof) over long time periods, given that cross-barrier movements are probably too infrequent to be tracked in real-time for individual fish. We sampled round goby populations at sites throughout their invasive range in Central Europe where round goby already had established populations on either side of river barriers. We then related genetic population differentiation to geographic distance and the number of river barriers (e.g., weirs, hydropower dams). We predicted that population structure would follow a pattern of isolation by resistance, which would be evidenced by higher levels of population differentiation across multiple river barriers relative to barrier-less river stretches. This would indicate either a limited capacity for the fish to bypass barriers of their own accord and/or limited human-facilitated transfer (e.g., bait fish dumping, boat traffic). The barriers studied here from river systems in the wild are diverse, and they vary in type and potential impact on fish dispersal. We can therefore only preliminarily explore whether certain river barrier types appear to be more impactful. Finally, this data allowed us to investigate if more upstream populations showed reduced genetic variation in relation to downstream populations, as expected when gene flow towards the invasion front is restricted.
We sampled round goby from four European rivers (Rhine, Danube, Dyje and Morava rivers) in three countries, and all field collections were carried out under the appropriate licenses issued within each country (see section "Ethics and permits" at the end of the publication for details). Sampling sites in the Danube, Dyje, and Morava rivers are thought to have a single invasion front for round goby that is progressing upstream, to the best of our knowledge. In the Rhine river, sampling sites were located between two commercial harbours that receive ships from different shipping corridors, thereby increasing the chances of two potential invasion fronts (implications of which are detailed in the discussion). Sampling sites were located upstream and downstream of pre-existing and diverse river barriers types such as weirs and hydropower dams, some of which were outfitted with fish passages (Table
Description of river barriers and the sampling environments, including the estimated year of the arrival of round goby at that site.
| River, country | Barrier site | Barrier | Fish pass description | Lock system | River width | Cargo boats | Sampling year, Estimated year of goby arrival |
|---|---|---|---|---|---|---|---|
| Danube, DE | Poikam | Dam | No fish pass | Yes | 120–180 m | Yes | 2022, 2009 |
| Danube, DE | Bergheim | Dam | Nature-like fish pass | No | 75–240 m | No | 2022, 2020 |
| Dyje, CZ | Bulhary | Weir | Boulder ramp | No | 50–60 m | No | 2022, 2012 |
| Dyje, CZ | Lednice | Weir | Boulder ramp | No | 50–60 m | No | 2022, 2009 |
| Dyje, CZ | Břeclav | Weir | Boulder ramp | No | 50–60 m | No | 2022, 2008 |
| Morava, CZ | Morava | Weir | No fish pass | No | 30–35 m | No | 2022, 2008 |
| Rhine, CH | Birsfelden | Dam | Vertical slot | Yes | 200–250 m | Yes | 2022, 2013 |
| Rhine, CH | Augst-Wyhlen | Dam | Vertical slot & fish lifta | Yes | 200–250 m | Yes | 2022, 2015 |
Pairwise population differentiation within each river system. Each river barrier has both an upstream (US) and downstream (DS) site and a sample size of fish collected from each site. FST values are reported along with corresponding P-values (applying a Benjamini-Hochberg correction for multiple comparisons within each river), and cells are bolded when they indicate significant differentiation between populations on either side of a river barrier.
| Danube river | ||||||
| Bergheim US | Bergheim DS | Poikam US | Poikam DS | |||
| N = 50 | N = 50 | N = 50 | N = 50 | |||
| Bergheim US | NA | |||||
| Bergheim DS | -0.0043 | NA | ||||
| P = 0.91 | ||||||
| Poikam US | 0.069 | 0.058 | NA | |||
| P < 0.0001 | P < 0.0001 | |||||
| Poikam DS | 0.11 | 0.093 | 0.013 | NA | ||
| P < 0.0001 | P < 0.0001 | P = 0.0094 | ||||
| Dyje river | ||||||
| Bulhary US | Bulhary DS | Lednice US | Lednice DS | Břeclav US | Břeclav DS | |
| N = 45 | N = 32 | N = 41 | N = 50 | N = 48 | N = 50 | |
| Bulhary US | NA | |||||
| Bulhary DS | 0.0010 | NA | ||||
| P = 0.41 | ||||||
| Lednice US | 0.025 | 0.0078 | NA | |||
| P < 0.0001 | P = 0.11 | |||||
| Lednice DS | 0.041 | 0.020 | 0.0014 | NA | ||
| P < 0.0001 | P = 0.0092 | P = 0.40 | ||||
| Břeclav US | 0.034 | 0.022 | 0.0012 | -0.0026 | NA | |
| P < 0.0001 | P < 0.0001 | P = 0.40 | P = 0.78 | |||
| Břeclav DS | 0.049 | 0.021 | 0.0063 | 0.0084 | 0.012 | NA |
| P < 0.0001 | P = 0.0021 | P = 0.11 | P = 0.039 | P = 0.0021 | ||
| Morava river | ||||||
| Morava US | Morava DS | |||||
| N = 44 | N = 50 | |||||
| Morava US | NA | |||||
| Morava DS | 0.00009 | NA | ||||
| P = 0.45 | ||||||
| Rhine river | ||||||
| Augst US | Augst DS | Birsfelden US | Birsfelden DS | |||
| N = 18 | N = 28 | N = 49 | N = 33 | |||
| Augst US | NA | |||||
| Augst DS | 0.0054 | NA | ||||
| P = 0.35 | ||||||
| Birsfeld US | 0.0030 | 0.0092 | NA | |||
| P = 0.40 | P = 0.090 | |||||
| Birsfeld DS | 0.0090 | 0.013 | 0.0092 | NA | ||
| P = 0.16 | P = 0.047 | P = 0.047 | ||||
In May 2022, we sampled round goby from sites along the Dyje and Morava rivers in the Czech Republic via electrofishing (Fig.
Maps of sections of the Dyje, Rhine, Danube, and Morava rivers where we sampled round goby populations (green circles) on either side of river barriers (yellow rectangles). Estimates of population differentiation, FST, and the respective P-values (after applying a Benjamini-Hochberg correction for multiple comparisons within each river) are given for consecutive sampling sites along each river (white boxes), some of which are separated by a barrier (described in beige boxes with site name). Note that in the Danube, between Bergheim and Poikam, there are two other hydropower dams, each at Vohburg and Ingolstadt (with boat locks).
Considering previous findings on limited movement of gobies across the banks of larger rivers (
DNA was extracted from fin clips using Qiagen DNeasy blood and tissue kits and diluted 1:30 with 10 mM Tris-HCl pH 8, following manufacturer instructions. Each fish was then genotyped at ten microsatellite loci using two multiplex reactions. PCR volumes consisted of 5 μl of Qiagen HotStarTaq Master Mix, 2 μl of template DNA, 1.5 μl of primer mix (0.3 μM concentration per forward and reverse primer), and 1.5 μl of ddH2O. Forward primers were labelled with one of the following fluorescent dyes: FAM, ATTO550, ATTO565, ATTO532. We used the following PCR programme settings: denaturation at 95°C for 15 min, followed by 35 cycles at 94°C for 30 s, annealing at 56°C (for both multiplexes) for 90 s, extension at 72°C for 60 s and a final extension at 72°C for 30 min. We scored allele sizes against an internal standard (GeneScan 500 LIZ; Applied Biosystems) in an automatic sequencer (ABI3730 Genetic Analyzer; Applied Biosystems).
Genotyping revealed several round goby individuals that had three allele peaks at certain loci. This occurred in seven fish sampled by the Poikham hydropower dam in the Danube river, Germany (two fish had three alleles at locus Nme3, two fish had three alleles at locus Ame17, and three fish had three alleles at locus Ame133). We omitted the allelic information at these loci for these fish.
However, our genotyping also revealed three alleles at locus Ame133 for 17 fish sampled by the Augst and Birsfelden hydropower dams in the Rhine river in Switzerland. In these cases, the third allele was always associated with the genotype 195/219/X (fragment sizes in base pairs [bp]), with X being the third allele. Furthermore, the allele size 219 bp always co-occurred with allele size 195 bp, but was rare amongst the individuals that did not possess a third allele (found in four such individuals). This pattern suggests that a chromosomal region bearing locus Ame133 may have been duplicated in the Swiss population, giving rise to a haplotype carrying alleles 195 bp and 219 bp. We represented this haplotype by the 195 bp allele; that is, we omitted the 219 bp allele from the dataset, which resulted in genotypes 195/219/X being transformed into 195/X, and 195/219 being transformed into 195/195 for locus Ame133 (Note that analyses that exclude this locus entirely from the Rhine sites produce qualitatively similar results).
We used Arlequin to test for Hardy-Weinberg equilibrium (HWE) for each microsatellite marker in each population sample (ten markers each tested in 16 population samples = 160 HWE tests). Several deviations were detected after Benjamini-Hochberg corrections, however, given the recent colonization of these populations, some deviations from HWE are to be expected due to the relatively stronger effects of genetic drift at small (presumed) population sizes (Suppl. material
We first used Arlequin to calculate pairwise differentiation, FST, between populations within each river. P-values were calculated based on 1,023 permutations of the data, and adjusted for multiple comparisons within each river using the Benjamini-Hochberg method. We also used Arlequin to calculate expected heterozygosity (He) as a measure of genetic diversity within populations at each of our sampling sites, as well as to calculate the modified Garza-Williamson index (M) as a measure of recent population bottlenecks (
Within each river system, some pairs of sampling sites straddled river barriers (in relatively close proximity to the barriers), while others were separated by longer stretches of river that may have been free of barriers or may have had several barriers (Fig.
We next tested whether He declined as we sampled further upstream from our most downstream sites, and as we sampled across consecutive river barriers. Expected heterozygosity is predicted to decline further upstream as sampling sites approach the invasion front where the round goby population is more recently founded. Expected heterozygosity is also predicted to decline as additional river barriers are surpassed by small subsets of the downstream populations. We therefore fit a linear mixed effects model (LMM) to the He data (means across loci for each population) using the ‘glmmTMB’ R package, v 1.1.8 (
Our MRM models indicated that population pair-wise FST in the Danube was marginally associated with geographic distance between sampling sites (P = 0.05; Fig.
Pairwise FST values plotted against A distance between sampling sites, or B number of river barriers between sampling sites, paneled by river system. * shows a significant (P < 0.05) result from the MRM models. Slopes derived from a linear regression to aid in visualizing the MRM model results, coloured data points show raw values. Note, the Morava river was not included in the MRM analyses because only two sites were sampled.
Significant population differentiation was detected across three out of the eight barriers: Poikam hydropower dam in the Danube, Birsfelden hydropower dam in the Rhine river, and Břeclav weir in the Dyje river (Fig.
Within-population genetic diversity represented by He declined as we sampled further upstream (LMM; est. ± std. error = –0.0014 ± 0.0005, z = –2.84, P = 0.0046), but was not clearly associated with the number of barriers between the sampling site and the most downstream site (est. ± std. error = –0.010 ± 0.0086, z = –1.18, P = 0.24, Fig.
Model prediction (fitted line with 95% confidence interval) showing how the expected heterozygosity (He) of round goby populations decreases as sampling progressed upriver and approached the various rivers’ invasion fronts. * shows a significant (P < 0.05) result for geographic separation from the LMM. Coloured data points show raw He values.
We investigated whether river barriers slowed the dispersal of invasive round goby using a spatial population genetics approach, sampling round goby populations on either side of barriers in four different rivers in their invasive range in Europe. We found that populations were significantly differentiated over increasing distances and increasing numbers of barriers, depending on the river system. Starting broadly, round goby population differentiation was marginally associated with distance and significantly associated with barrier number in the Danube River (Germany); with increasing number of barriers, but not with distance in the Dyje River (Czech Republic); and with neither distance nor number of barriers in the Rhine River (Switzerland). Our findings that populations differentiate with distance (albeit marginally after accounting for barrier number) are in line with previous work at similar spatial scales. For instance
River barriers are highly diverse, ranging in size and structure and differing in their implementation of fish passage solutions (Table
We observed significant upstream-downstream differentiation across three individual barriers, one within each of the three river systems. The rivers and barriers studied here were highly diverse, reflecting the diversity of barrier types that have been installed in European river systems. Therefore, we explored our data to investigate whether there were any consistent patterns of differentiation associated with one barrier type. However, given the low level of replication within barrier types in our dataset, this investigation should be treated as preliminary. Overall, however, we saw no clear indiciation that any particular river barrier type (e.g., weir, dam – with or without fish passage solutions) was related to round goby population differentiation. In two of the cases where we observed significant cross-barrier differentiation, it occurred at small spatial scales with less than 2 km between the upstream and downstream sites. The Poikam dam in the Danube has no fish passage solution (e.g., fish ladders), which may hinder round goby dispersal across the barrier. However, at the same time there is a large lock system allowing ships to pass by the dam, which fish can presumably use as well. Meanwhile, the Birsfelden dam in the Rhine also showed genetic differentiation and has a vertical slot fish passage. Vertical slot passages may be more challenging for round goby to pass if the flow is high, the passage is steep, or the passage contains no bottom structure on which the goby can anchor itself (
It should be noted that several mechanisms can contribute to our observed population structures, including stochasticity in the genetic makeup of the founder population, gene flow, genetic drift (particularly in small populations), mutations, and different selective regimes across populations. Our measures of population differentiation, however, ought to be most strongly influenced by founder effects, gene flow (or hinderances to gene flow) and somewhat less by genetic drift, as our study populations are extremely young (founded within the past 8–15 years) (
We also found that round goby showed less genetic variation in the furthest upstream—and therefore presumably most recently invaded—sampling sites relative to the downstream sites within the same river system. Lower genetic diversity at the invasion front can occur in newly invaded areas and diversity can increase later as alleles accumulate over time (
In conclusion, we have used a spatial genetics sampling approach to test the hypothesis that river barriers would slow round goby dispersal as indicated by population differentiation. This work was motivated by an increasing recognition of the benefits of river connectivity, which may also come with the cost of increased spread of invasive species, like the round goby (
ESM: Conceptualization, Funding acquisition, Investigation, Writing original draft, Writing – review & editing Visualization; KMS: Formal analysis, Writing - review & editing; PB: Investigation, Resources, Writing – review & editing; TB: Resources; KB: Investigation, Resources, Writing – review & editing; AF: Funding acquisition, Writing – review & editing; JG: Investigation, Resources, Validation, Writing – review & editing; MJ: Investigation, Resources, Validation, Writing – review & editing; PJ: Investigation, Resources, Validation, Writing – review & editing; JP: Investigation, Resources, Validation, Writing – review & editing; JMM: Investigation, Validation, Writing – review & editing; APHB: Conceptualization, Funding acquisition, Formal Analysis, Investigation, Writing original draft, Writing – review & editing, Visualization.
This research was supported by funding from the Swedish Environmental Protection Agency (Naturvårdsverket, 2020-00051 to A.F. and E.S.M.) the Kempe Foundation (JCK22-0043 to E.S.M. and A.B.), and the Swedish Research Council Formas (2022-02796/2023-01253 to J.M.M). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
We sampled round goby from four European rivers in three countries, and all field collections were carried out under the appropriate licenses issued within each country (Czech Republic: JMK 153264/2020; Switzerland: authorizations for special fishing catches 01/02/2022 by Authority for Environment and Energy; Germany: electrofishing activities carried out under the license number 31-7562 issued by the district office of Freising, Bavaria, Germany).
Species georeferenced records are available at the European Alien Species Information Network: https://easin.jrc.ec.europa.eu/easin/RJD/Download/5c243fea-47d6-489d-8d37-129fb35ff91b.
Besides the first, second, and last authors, the remaining authors are ordered alphabetically and contributions should be considered relatively equal. We would like to thank the anonymous reviewers for the feedback on this publication.
Supplementary data
Data type: xlsx
Explanation note: Supplementary data file containing the data and the Arelquin analysis output that were used in the analyses of the data presented here.
Supplementary figure and table
Data type: docx
Explanation note: Word document containing a supplementary figure and table. figure S1 shows the model prediction of how the expected heterozygosity (He) of round goby populations decreases as sampling progresses upstream and approaches various invasion fronts. table S1 shows the marker polymorphism of 10 microsatellites used in this study for each of the 16 populations samples.