Issue |
Knowl. Manag. Aquat. Ecosyst.
Number 421, 2020
Topical Issue on Fish Ecology
|
|
---|---|---|
Article Number | 25 | |
Number of page(s) | 13 | |
DOI | https://doi.org/10.1051/kmae/2020016 | |
Published online | 05 June 2020 |
Research Paper
Genetic diversity of endangered Chondrostoma nasus in the River Rhine system: Conservation genetics considerations on stocking and reintroduction
Diversité génétique de Chondrostoma nasus, une espèce menacée, dans le bassin du Rhin : considérations de génétique de conservation sur l'empoissonnement et la réintroduction
1
iES Landau, Institute for Environmental Sciences, University of Koblenz-Landau, Fortstrasse 7, 76829 Landau, Germany
2
Bürogemeinschaft für fisch- & gewässerökologische Studien, Über dem Grund 1, 35041 Marburg-Michelbach, Germany
3
Aquatic Ecology and Evolution, Institute of Ecology and Evolution, University Bern, Baltzerstrasse 6, 3012 Bern, Switzerland
4
Department of Fish Ecology and Evolution, EAWAG, Seestrasse 79, 6047 Kastanienbaum, Switzerland
5
Eusserthal Ecosystem Research Station, University of Koblenz-Landau, Birkenthalstrasse 13, 76857 Eusserthal, Germany
* Corresponding author: wetjen@uni-landau.de
Received:
24
December
2019
Accepted:
30
April
2020
Reintroduction, stocking and translocation of freshwater fish are of growing concern given their importance for biodiversity conservation and ecosystem functioning. For successful management and stocking programmes, it is essential to incorporate genetics-based approaches. The nase (Chondrostoma nasus) constituted one of the most common fish species in European rivers. Its highly specialised and migratory nature exposed the species to human pressures, and thus, promoted its decline. Current knowledge of the genetic structure of C. nasus is considerably limited for Europe as a whole and for Germany specifically. To overcome this lack of information we present original data on C. nasus from different tributaries of the River Rhine. We analysed nine microsatellite markers and mtDNA Cytochrome b sequences to assess the distribution of genetic diversity and structure of this species across the study area. With the exception of the Lake Constance/Alpine Rhine population, C. nasus exhibited high gene flow within the Rhine system, and therefore, limited geographical genetic differences between populations where migration is not prevented by human intervention. The present study provides new insights into the levels of genetic variability of C. nasus in the Rhine system, providing useful information for guiding reintroduction and stocking programmes. Population genetic information will improve future preservation and management of this valuable freshwater fish species in Germany and beyond.
Résumé
La réintroduction, l'empoissonnement et le transfert de poissons d'eau douce sont de plus en plus préoccupants compte tenu de leur importance pour la conservation de la biodiversité et le fonctionnement des écosystèmes. Pour que les programmes de gestion et de repeuplement soient efficaces, il est essentiel d'intégrer des approches fondées sur la génétique. Le hotu (Chondrostoma nasus) constitue l'une des espèces de poissons les plus communes dans les rivières européennes. Sa biologie hautement spécialisée et migratoire a exposé l'espèce aux pressions humaines, et a donc favorisé son déclin. Les connaissances actuelles sur la structure génétique de C. nasus sont considérablement limitées pour l'Europe dans son ensemble et pour l'Allemagne en particulier. Pour pallier ce manque d'informations, nous présentons des données originales sur C. nasus provenant de différents affluents du Rhin. Nous avons analysé neuf marqueurs microsatellites et des séquences d'ADNmt du cytochrome b pour évaluer la répartition de la diversité génétique et la structure de cette espèce dans la zone d'étude. À l'exception de la population du lac de Constance et du Rhin alpin, C. nasus présentait un flux génétique élevé dans le système rhénan et, par conséquent, des différences génétiques géographiques limitées entre les populations où la migration n'est pas empêchée par l'intervention humaine. La présente étude apporte de nouvelles connaissances sur les niveaux de variabilité génétique de C. nasus dans le système rhénan, fournissant des informations utiles pour orienter les programmes de réintroduction et d'empoissonnement. L'information génétique des populations permettra d'améliorer la préservation et la gestion futures de cette précieuse espèce de poisson d'eau douce en Allemagne et au-delà.
Key words: genetic variability / cytochrome b / gene flow / source stock prioritisation / Cyprinidae
Mots clés : variabilité génétique / cytochrome b / flux génétique / priorisation du stock source / Cyprinidae
© M. Wetjen et al., Published by EDP Sciences 2020
This is an Open Access article distributed under the terms of the Creative Commons Attribution License CC-BY-ND (https://creativecommons.org/licenses/by-nd/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. If you remix, transform, or build upon the material, you may not distribute the modified material.
1 Introduction
Human interventions in freshwater ecosystems worldwide have affected the abundance and distribution of freshwater organisms with the consequences of causing taxonomic homogenisation by non-native species introduction (Villéger et al., 2011), population declines or even species losses (Vaughn, 2010). With at least 39% of freshwater fish being threatened, they are one of the most severely threatened taxonomic groups being assessed to date (Freyhof and Brooks, 2011). Even though fish restoration practices were realised throughout Germany and other countries, efforts often failed to establish self-sustaining and reproductive populations. Reasons of failure may be that management purposes often rest on private activities, are conducted at a local scale and still show deficits in the compilation, coordination and exchange of knowledge between scientists and managers (Lundmark et al., 2019). For many native freshwater fish species, neither national nor international, genetically informed management strategies exist. However, in Germany there are attempts to collate genetic information of freshwater fish e.g. in the AGRDEU database (https://agrdeu.genres.de/nationales-inventar-aqgr/) to provide information on the genetic diversity of hatchery strains and wild fish populations for political consulting, management purposes and information exchange. Genetic information is necessary for adopting appropriate conservation strategies, will be of help to maintain biodiversity and to prevent the loss of genetic resources and the eventual erosion of ecosystem function (Brodersen and Seehausen, 2014). Intraspecific diversity within a species may reflect either historical isolation (genetic drift, reinforced by limited gene flow between disconnected populations), patterns of local adaption (natural selection), or both. Either way, recognizing and preserving diversity is substantial for successful conservation on the long term. It is expected to stabilize ecosystem services and serve as a portfolio effect buffering metapopulations against pervasive effects of environmental change and perturbation (Schindler et al., 2010; Piccolo et al., 2018).
Even though the nase (Chondrostoma nasus L.) has a wide distribution area and was a common fish species in Europe (Flore and Keckeis, 1998), many populations have declined especially in Germany, Switzerland, Austria, Poland and the Czech Republic (Peňáz, 1996). The River Rhine system is the westernmost range limit of the species' natural distribution (Hudson et al., 2014), whereas nase has been introduced to further rivers in France, including the Rhone and Loire rivers. In the Rhine system and its tributaries such as the rivers Lahn, Neckar, Main, Nidda and Kinzig, nase was historically present in high abundances (Dußling and Berg, 2001; Hmuklv and Hessen-Forst Fena, 2014). However, in several surveys of the River Rhine basin in the 1990s, nase occurred only sparsely at low numbers and was considered extinct in the River Lahn (Freyhof, 1997, and cited literature within). In consequence of the drastic decline in the River Lahn, first attempts to reintroduce the species were implemented. In the years 1995–1997, 154.000 fish were stocked between the upper Lahn, near Kernbach and the middle reaches of the Lahn, near Steeden (Fig. 1). Stocked fish originated from a hatchery stock that derives from the River Nister in the River Rhine system (Fig. 2). Intense monitoring in subsequent years revealed no success (Schwevers and Adam, 1997) and have led to the conclusion that nase was in fact gone extinct in the entire Hessian catchment of the river Lahn. Sporadic and widely scattered encounters (N = 5 in total) were recorded and interpreted as remnants of specific local stocking efforts. However, no reproduction was detected (Schwevers and Adam, 1997). Only since 2006, fish were recorded five kilometres upstream Steeden in the middle reaches of the Lahn (Fig. 1). Since 2009 nase was caught frequently in the middle reaches (Hübner and Fricke, 2011, 2012, 2014) and a reproductive, self-sustaining population was re-established (Hübner et al., 2016). It remains uncertain whether fish in the middle reaches of the Lahn originated from the first reintroduction activities or from an autochthonous remnant stock of the watershed itself that has survived the collapse years ago and recovered again. However, in the upper River Lahn above Marburg, still no nase population exists. Further reintroduction attempts in the upper River Lahn in the years 2014 and 2015 used the same hatchery stock as in the 1990s, but failed again. Several weirs between the middle and upper reaches of the Lahn exist and may represent insurmountable obstacles for the upstream migration of nase (Fig. 1). In addition, different habitat conditions in the upper River Lahn (e.g., higher discharge, narrower riverbed) and/or improper genetic-based characteristics of the stocked fish could have prevented a re-establishment of nase in the upper Lahn. Since 2017, stocking strategies for the reintroduction of nase in the upper Lahn have been changed as part of the EU-LIFE project “Living Lahn” (Hübner et al., 2017). Since then, juveniles are produced by stripping spawning animals from the middle Lahn itself and the fertilized eggs are then raised up to stock size at a hatchery farm until they are released into the upper Lahn.
Genetic variability patterns are essential for a sustainable management and provide relevant information for future stocking programmes. Genetic information is needed to avoid inbreeding and to preserve the adaptive potential of a species since genetic diversity is linked to biological productivity and resilience to stressors. So far, the population genetic structure of the nase in Germany has been investigated to a very limited extent. To our knowledge, only one population genetics study investigated German nase populations (Hudson et al., 2014). In this context, the aims of the present study were (i) to assess the genetic diversity and variability of the widespread but locally threatened species in its historically native distribution range within the River Rhine system, (ii) to assess gene flow and identify possible distinct genetic populations, (iii) to ascertain the origin of the population found in the middle reaches of the River Lahn (stocked hatchery fish or survived autochthonous stock), and (iv) to identify the most suitable stock for future reintroduction of the nase in the upper River Lahn.
Based on genetic information, this study will contribute to improve our knowledge for the development of management strategies for the freshwater specialist C. nasus and its declining autochthonous populations in Germany.
Fig. 1 The course of the River Lahn and its obstacles impeding fish migration. Stocking sites of Chondrostoma nasus in the years 1995–1997 and 2014–2015 were located between the upper course (near Kernbach) and the middle course (near Steeden) of the River Lahn. |
Fig. 2 Sampling locations and observed haplotypes per sampling site of Chondrostoma nasus in the River Rhine basin. The size of the circles is proportional to the number of analysed individuals per sampling site. |
2 Material and methods
2.1 Sampling
Chondrostoma nasus individuals were collected from 10 locations in the River Rhine system (Tab. 1, Fig. 2) and from two hatchery stocks. The hatchery_1 stock (H1) was originally established from fish of the river Main near Obernburg, whereas the hatchery_2 stock (H2) was established from fish of the River Nister in 1984 and was used for the reintroduction of C. nasus in the River Lahn until 2015. The population named LA (Lahn) in this study originates from the middle reaches of the River Lahn. Several weirs mark the watercourse between the middle and upper River Lahn (Fig. 1) and it is uncertain whether nase is able to pass these upstreamwards. Investigated specimens from the River Murg, Wiese and Dornbirner Ach originate from the study of Hudson et al. (2014) but were reanalysed with a different mtDNA marker in this study. Wild fish were caught by electrofishing. From all fish, fin clips were taken in situ and stored in 96% ethanol at −20 °C until further processing. Fish were released at their home locality immediately after sampling.
Sampling sites and number of analysed Chondrostoma nasus individuals per population (Pop) at the sequence (NmtDNA) and microsatellite (NMS) level in the River Rhine system.
2.2 Sequence analysis
DNA extraction was performed following a modified (Wetjen et al., 2017, 2020) DNA salt-extraction protocol (Aljanabi and Martinez, 1997). A 1094 bp fragment of the mtDNA cytochrome b gene (Cytb) was obtained for 154 specimens using the primers Glu-L and Thr-H (Mäkinen and Merilä, 2008). Polymerase chain reactions (PCR) were performed using a Primus 96 Cycler (Peqlab Biotechnologie GmbH, Erlangen, Germany) under the following conditions: Denaturation (2 min, 95 °C) was followed by 25 cycles of 30 s at 94 °C, 30 s at 54 °C and 45 s at 72 °C, and a final elongation step (10 min, 72 °C). The company SeqIT GmbH & Co.KG (Kaiserslautern, Germany) conducted bidirectional sequencing of the DNA strands. Sequence segments were aligned using the programme Geneious v. 6.1.8 (Biomatters Ltd.). Generated sequences were compared with those available in GenBank using the ‘Basic alignment search tool' (BLAST).
2.3 Microsatellite analysis
Nine microsatellite loci were co-amplified for 284 individuals: LC27, LC290 (Vyskocilová et al., 2007), Lsou05, Lsou08, Lsou21 (Muenzel et al., 2007), SarN7F11b, SarN7F8, SarN7G5, SarN7K4 (Mesquita et al., 2003). Multiplex-PCRs were performed according to Hudson et al. (2014) with slight modifications resulting in a reduction of the used primer concentrations. Fragment analysis was performed on a Beckman Coulter CEQ 8000 eight capillary sequencer. Loci were scored using the software CEQ SYSTEM v. 9 (Beckman Coulter).
2.4 Data analysis
A median-joining (MJ) network (Bandelt et al., 1999) of the generated mtDNA sequences was constructed using the programme PopART (Leigh and Bryant, 2015). To obtain meaningful, unimpaired results, fish of the populations BA, LW and SM were only included in the MJ Network but excluded from any further analyses as only 2–3 individuals were available from these locations (Tab. 1). Haplotype (h) and nucleotide diversity (π) indices were estimated in DNAsp v.5.10.1. True diversity (DmtDNA) was calculated as the effective number of haplotypes per population (Jost, 2010). To analyse the demographic history and evidence of population expansion at the sequence level, Tajima's D (Tajima, 1989), Fu's FS (Fu, 1996), and Ramos-Onsins and Rozas R 2 (Ramos-Onsins and Rozas, 2002) statistics of neutrality were estimated in DNAsp v. 5.10.1. Significance was assessed from 1000 coalescent simulations.
MICRO-CHECKER v. 2.2.3 (Van Oosterhoudt et al., 2004) was used to test for the presence of null alleles, allele dropouts and scoring errors in the microsatellite dataset. Observed (H O) and expected (H E) heterozygosity, deviations from Hardy–Weinberg equilibrium (HWE), linkage equilibrium (LD), and the fixation index (F) were assessed using GenAlEx v.6.5 (Peakall and Smouse, 2012). Values of a fixation index close to zero are expected under random mating, whereas substantial positive values are indicative of inbreeding. True diversity (D MS) was estimated as the effective number of alleles per population (Jost, 2010). Allelic richness (AR, i.e. Petit et al., 1998) and private allelic richness (PAR, i.e. Kalinowski, 2004) were calculated using rarefaction in HP-RARE (Kalinowski, 2005). Effective population size (Ne) was estimated using NeEstimator v.2.1 (Waples et al., 2014) and the number of migrants (Nm) using GENEPOP v. 1.2 (Raymond and Rousset, 1995). Additionally, first generation migrants (F0) were identified in GeneClass2 (Piry et al., 2004) to evaluate migrants influence on genetic differentiation. To account for unsampled populations and maximising the analyses power, we used two test statistics (Lhome/Lmax and Lhome) to compute the likelihood of migrant detection (L) at the microsatellite level (Paetkau et al., 2004). Furthermore, a self-assignment computation was conducted for the individual genotypes using a frequency-based method with an assignment threshold of 0.05 in GeneClass2.
A Bayesian, Markov Chain Monte Carlo (MCMC) analysis was conducted to assess the extent to which allelic structures of the sampled populations could be clustered into the same group (K) utilizing the programme STRUCTURE v. 2.3.3 (Pritchard et al., 2000). An admixture model with correlated allele frequencies and prior information about individual populations was used to further shape the analysis. The analysis was run ten times for each K (1–16) with a burn-in period of 50,000 followed by 500,000 MCMC. The most probable K was ascertained following the delta K method (Evanno et al., 2005) using Structure Harvester (Earl and Von Holdt, 2012). The detection of one cluster would indicate high admixture between populations, whereas different geographic clusters would indicate that populations might have experienced genetic drift leading to genetic differentiation.
In addition, a distance-based principle component analyses (PCoA) was performed with data standardisation for nase populations to visualise the variance in microsatellite allelic structure within the River Rhine basin. PCoA was conducted using GenAlEx v.6.5 (Peakall and Smouse, 2012).
We computed measures of pairwise differentiation represented by the population level fixation indices F ST (microsatellites) and Φ ST (mtDNA sequences). Significance was assessed using 10,000 permutations. Hierarchical analysis of molecular variance (AMOVA) was additionally performed for both marker types to detect differentiation in genetic variability among populations. Both analyses were conducted in Arlequin v. 3.11 (Excoffier et al., 1992).
To assess the effective number of distinct populations, we calculated the between-population component ΔST for both marker types as described in Jost (2008). To determine the contribution of the individual populations to the total distinctiveness, ΔST was additionally calculated removing each individual population separately from the entire dataset (including all but one population). This allows for comparing the change in ΔST when one population is excluded, and hence, estimating its contribution to the overall effective number of distinct populations.
Demographic processes such as expansions and reductions in the effective population size affect genetic diversity. During a bottleneck, characteristic mode-shift distortion in the distribution of allele frequencies can be detected since allelic diversity is supposed to be reduced faster than heterozygosity (Luikart et al., 1998). While low frequency alleles will get lost, the remaining alleles maintain present in intermediate frequencies. Hence, it is expected that the observed heterozygosity is larger than the heterozygosity expected from the observed allele number at mutation-drift equilibrium (Cornuet and Luikart, 1996). We used the programme BOTTLENECK (Piry et al., 1999) to test for recent bottlenecks in individual populations applying the stepwise mutation model (SMM, Kimura and Ohta, 1978) and the two-phase mutation model (70% stepwise mutations, TPM). We used the sign test and the Wilcoxon signed-rank test to assess significant heterozygosity excess (Piry et al., 1999). Significance of multiple comparisons was determined by adjusting p-values for false discovery rates (FDR) following the method of Benjamini and Hochberg (1995).
3 Results
3.1 Haplotype diversity, distribution and demography
In total, 11 different mitochondrial Cytb haplotypes were observed in the River Rhine basin (Fig. 2). Haplotypes were differentiated by a maximum of 26 mutational steps (Fig. 3). The haplotype Hap_1 was most common (N = 100; 64.94%) and present in all populations except the Dornbirner Ach, a tributary of the Alpine Rhine upstream of Lake Constance. The haplotype Hap_1 was already identified in a previous study (Accession number: JQ652366.1) by Dubut et al. (2012), whereas the other ten haplotypes were newly identified for nase in this study (deposited in GenBank; Accession numbers: MN431659–MN431668). The highest number of private haplotypes was present in the River Lahn with PHap = 3. Diversity ranged from DmtDNA = 1.00 (DO) to DmtDNA = 2.88 (LA) at the sequence level (Tab. 2). Fish of the Dornbirner Ach exclusively carried a single haplotype (Hap_8). The haplotype Hap_8 was private for this site and is probably restricted to nase from the Alpine Rhine/Lake Constance basin. Noticeable is the haplotype Hap_11, which was only present in the Minthesee population with one individual. This haplotype was most distantly related to all other haplotypes by at least 15 mutational steps (Fig. 3). It was not found in GenBank but was closest related (98.99% identical; differentiated by 12 mutational steps) to a Chondrostoma angorense haplotype found in the Kizilirmark river, Turkey (Perea et al., 2010; Accession number: HM560078.1).
Neutrality tests revealed no recent population expansion within the studied populations (Tab. 3). Tajima's D were significantly negative only in the hatchery_2 and the Minthesee populations. However, Fu's FS and R 2 were positive and not significant for both these populations.
Fig. 3 Median-joining network of 154 individuals of Chondrostoma nasus of the River Rhine basin including reference haplotypes of Chondrostoma angorense (Hap_12) and Parachondrostoma toxostoma (Hap_13). The size of the circles is proportional to the frequency of fish carrying the respective haplotype. Colours refer to the origin-sampling site. Mutational steps are represented by intermediated dashes. Black dots represent missing haplotypes. |
Summary of genetic diversities within the investigated Chondrostoma nasus populations. For each population (Pop) true diversity (DMS), allelic richness (AR), private allelic richness (PAR), observed (HO) and expected heterozygosity (HE), fixation index (F), true diversity (DmtDNA), number of haplotypes (NHap), haplotype diversity index (h), number of private haplotypes (PHap), segregating sites (S) and the nucleotide diversity index (π) is reported.
Neutrality tests and evidence for recent population bottlenecks using the sign test and the Wilcoxon test of investigated populations (Pop) of Chondrostoma nasus. DT = Tajima's D; FS = Fu's Fs ; R 2 = Ramos-Onsins and Rozas R 2; TPM = two-phase model; SMM = stepwise mutation model. Significant values are in bold.
3.2 Microsatellite diversity
Data showed no evidence for null alleles, large allelic dropout or scoring errors. Genotype proportions for all populations conformed to HWE expectations (p > 0.05). Only few loci showed a deviation from HWE in some populations (Lsou05, SarN7G5 in DO; Lsou21, SarN7F8 in MI; LC27, LC290, Lsou08 in NK). No significant linkage disequilibrium between markers was detected. Allelic richness was highest in the Lahn (AR = 4.16) and the Minthesee (AR = 3.99) populations, whereas the hatchery_2 population (AR = 3.88) showed lowest mean allelic richness (Tab. 2). Private allelic richness was highest in the Lahn (PAR = 0.42) and Dornbirner Ach (PAR = 0.48) populations. True diversity was highest in the Murg (D MS = 3.21) and Lahn (D MS = 3.14). However, the hatchery_2 stock showed a comparatively lower overall diversity, indicating a slightly reduced genetic diversity in this stock (Tab. 2).
3.3 Gene flow and population structuring
Clustering of samples into population groups via the Bayesian clustering analysis gave the highest probability for K = 5 (Fig. 4). One cluster (K1) was formed by the hatchery_2 stock and a second cluster (K2) by the Lahn population. The third cluster (K3) was formed by the hatchery_1 stock and the fourth (K4) by fish from the Dornbirner Ach of the Alpine Rhine. All remaining populations (Nister, Wiese, Murg, Minthesee, and northern Ketscher) grouped predominantly into the fifth cluster (K5) but displayed high admixture, especially with individuals belonging to the clusters of both hatchery stocks and the Lahn population. We did not find any evidence for gene flow between the Dornbirner Ach and any of the Rhine populations from downstream of the Rhine Falls.
PCoA confirmed the unique genetic signature of the Dornbirner Ach population within the Rhine system (Fig. 5). In addition, the Lahn population and both hatchery stocks (H1 and H2) grouped further apart from the rest of the investigated Rhine populations.
The AMOVA (Tab. 4) revealed that the genetic variance was explained by the variance within populations by 70.39% and between populations by 29.61% at the sequence level. Microsatellites showed also that most variation was explained by the variation within the individual populations (95.45%), whereas only 4.55% of the variation was explained by the variation among the different populations.
Effective number of distinct populations resulted in Δ ST = 1.379 at the mtDNA level and Δ ST = 1.226 at the microsatellite level for the nine populations (Tab. 5). This indicates that we cannot find definite distinct structuring within the Rhine basin below the Rhine Falls. However, the contribution of the individual populations to the effective number of distinct populations was highest for Dornbirner Ach (mtDNAΔ ST = 1.104; MSΔ ST = 1.182) and Lahn (mtDNAΔ ST = 1.349; MSΔ ST = 1.212) and lowest for the hatchery_2 stock.
Estimated genetic differentiation in terms of F ST and ΦST were non-significant between most populations. Significant differentiations ranged from F ST = 0.015 (NI vs. MU) to F ST = 0.147 (DO vs. NI) at the microsatellite level and from ΦST = 0.118 (NI vs. MU) to ΦST = 0.942 (DO vs. H2) at the sequence level (Tab. 6). The only population from upstream the Rhine Falls, the Dornbirner Ach, was considerably more strongly differentiated from any of the other populations than the others were amongst each other. Assignment test indicated high self-assignment at the microsatellite level only for the populations Lahn (45.24%), Nister (47.06%), hatchery_1 (65.00%), hatchery_2 (78.00%) and Dornbirner Ach (100%). Correct self-assignment of individuals from the remaining populations ranged between 10-23% (Tab. 7).
Overall number of migrants was Nm = 6.47 (data not shown), indicating recent gene flow between many populations of the different tributaries of the Rhine. Highest Nm were observed between the northern Ketscher and the Lahn (Nm = 5.85), and the northern Ketscher and the Nister (Nm = 6.57) populations. However, Nm was explicitly lower between the Dornbirner Ach and all other investigated populations below the Rhine Falls and thus, gene flow from below the Rhine falls to the Constance/Alpine Rhine system is very unlikely. GeneClass2 identified seven and two first-generation migrants using the Lhome and Lhome/Lmax statistics, respectively (Tab. 7). These migrants where detected in the Lahn (N = 2), Murg (N = 3), northern Ketscher (N = 3), and Minthesee (N = 1).
Fixation indices were slightly negative in all populations, except for the Lahn where the highest value (F = 0.006) was obtained (Tab. 2). Hence, no evidence for significant inbreeding was observed at population level.
Estimated effective population size of Ne = 7.5 (CI: 3.5–13.3) was found for hatchery_2 stock. For all other populations Ne-values ranged from 78.6 (H1) to 364.3 (NK), whereas the populations Nister, Dornbirner Ach, Murg and Minthesee showed infinite Ne-values (data not shown). Infinite Ne-values in some populations are probably due to limited sample size. Infinite estimates are usually interpreted as no evidence for any disequilibrium caused by genetic drift (see Waples and Do, 2010; Marandel et al., 2019 for further interpretation of Ne estimates). However, obtained Ne-values are conforming to the tests for recent bottlenecks. Based on the microsatellite data, evidence for a recent bottleneck was only found in the hatchery_2 stock under both models and test statistics (Tab. 3). Instead, the Lahn and Dornbirner Ach displayed evidence of significant heterozygous deficiency under both models under the sign test. The Nister showed significant heterozygous deficiency only under the SMM in the sign test.
Fig. 4 Bayesian clustering analysis of the investigated Chondrostoma nasus populations: H2 = Hatchery_2, LA = Lahn, NI = Nister, DO = Dornbirner Ach, WI = Wiese, MU = Murg, MI = Minthesee, NK = northern Ketscher, H1 = Hatchery_1. Colours refer to the predicted clusters (1–5). Each individual is displayed by a vertical column with the length proportional to the individual's estimated affiliation to the respective cluster. |
Fig. 5 Plots of the first two axes of a principal component analysis (PCoA) of Chondrostoma nasus populations using GenAlEx v. 6.5. Each coloured dot represents one population: H1 = Hatchery_1, H2 = Hatchery_2, NI = Nister, LA = Lahn, NK = northern Ketscher, MI = Minthesee, MU = Murg, WI = Wiese, DO =Dornbirner Ach. |
Analysis of molecular variance (AMOVA) for the sequence (mtDNA Cytb) and microsatellite (MS) data of Chondrostoma nasus. Shown is the degree of freedom (df), the sum of squares (Sum of sq.) and the variation in percent for each source of variation.
Change in the effective number of distinct populations (ΔST ) excluding individual populations from the entire dataset (all Pops). ΔST and the change (difference) are shown for both, the sequence (mtDNA) and microsatellite (MS) data.
Pairwise ΦST- (above the diagonal) and F ST-values (below the diagonal) of the investigated Chondrostoma nasus populations (Pop). Significant values are in bold.
Results of the self-assignment test and first generation migrant (F0) analysis. Each row contains the individuals from one sampling location and the columns indicate the populations to which these individuals were assigned (highest likelihood). Total number of investigated individuals per population (Total N), percentage of correctly assigned individuals (Corr. Assign. %) and the number of first generation migrants with their most likely source population in brackets for the two tests statistics Lhome and Lhome/Lmax are given. For the first generation migrant analysis, both hatchery stocks (H1 and H2) were excluded from the analysis.
4 Discussion
We used two different molecular markers to investigate the genetic diversity and variability within the endangered fish species, Chondrostoma nasus, in the River Rhine district, Germany. Overall diversity was quite diverging between the individual populations. Ten out of 11 haplotypes recorded in total, of which seven were private haplotypes, were undescribed so far. This reflects the limited number of genetic studies of nase in Germany and entire Europe. These previously undetected haplotypes will be relevant for the further assessment of phylogenetic aspects of the species and are of concern for proper conservation strategies, not only within the Rhine basin but also at a larger scale. Furthermore, it will be relevant for further investigation of hybridisation of Parachondrostoma toxostoma, a threatened species endemic to France and Switzerland (Paz-Vinas et al., 2013), and C. nasus being invasive in France and Italy (e.g. Costedoat et al., 2007).
Overall, low to moderate levels of population differentiation were observed. Results confirm a genetically distinct population in the Dornbirner Ach only. Hence, the Rhine Falls act as a physical barrier and prevent genetic exchange even for migrating fish (Hudson et al., 2014). Furthermore, the Dornbirner Ach was represented by a haplotype (Hap_8; Fig. 2) that was not found in any other population downstream the Rhine Falls. Fish have colonized Lake Constance from two different glacial refugia (Vonlanthen et al., 2011), the Rhine and the Danube. Such a contact zone/co-occurrence of different phylogenetic lineages/clades was already detected e.g. for burbot, Lota lota (Barluenga et al., 2006; Wetjen et al., 2020) in Lake Constance. This may similarly reflect the distinctiveness of the Dornbirner Ach population in this study. Equivalent results were also apparent when the Dornbirner Ach was compared to other populations from Switzerland (Hudson et al., 2014). Specimens of the Dornbirner Ach were genetically even more similar to Danubian populations than to Rhine populations below the Rhine Falls (Vonlanthen et al., 2011). Such populations deserve high conservation priority by themselves.
The nase seems to disperse across large geographic distances (e.g. LA and MI; F0 migrants = 1; F ST = 0.028; ΦST = 0.094). This was already investigated for nase (> 400 km) and other cyprinids, e.g. Leuciscus idus (> 100 km) (Winter and Fredrich, 2003 and cited literature within) and is supported by the low, non-significant levels of genetic differentiation among sampling locations downstream the Rhine Falls. Additionally, it is confirmed by the self-assignment analysis. The majority of individuals of MI, MU, WI and NK were not assigned to their sampling population. A pronounced sub-structuring (despite DO and LA) of the populations was not present, and most genetic variation was found within rather than between populations. Gene flow exceeded by far the benchmark of one migrant per generation. Hence, long-distance dispersal prevents genetic drift and isolation as long as migration is not obstructed. Previous studies likewise detected distinct genetic groups of nase only between different major river basins where admixture is naturally impossible e.g. between the Rhine and the current Danube basin (Hudson et al., 2014) that were connected in former times. Consequently, the dramatic population declines of nase in the River Lahn and elsewhere in the recent past are likely an ecological consequence of habitat loss/changes and fragmentation leading to changes in demography, rather than a consequence of inbreeding and the loss of genetic variation alone.
The middle reaches of the River Lahn are about twice as wide as the upper River Lahn where runoff fluctuations are much greater. Hence, it is possible that the donor population was not locally sufficiently adapted for reintroduction in the upper Lahn. Therefore, artificial restoring of spawning grounds and their morphodynamics is important (e.g. Hauer et al., 2007; Kitada et al., 2019) in the upper River Lahn. Stable riverbed structures are needed for successful spawning. Furthermore, the permeability and oxygen supply of the interstitial are essential factors and higher fine sediment accumulation reduce successful hatching (Duerregger et al., 2018; Nagel et al., 2019). Discharge and water quantities differ within the River Lahn. Reduced input of gravel during floods, a deepened watercourse and a deficiency of lateral erosion result in a lack of suitable spawning grounds and growth habitats for nase in the upper River Lahn. Hence, in the context of the current EU-LIFE project further renaturation and habitat enhancement measures are implemented.
Moreover, a major point affecting the establishment of a self-reproducing population in the upper Lahn is the presence of weirs between the middle and upper River Lahn. Even though there have been several key improvements ensuring the passage for migrating fish in the tributaries of the River Rhine and the River Lahn itself, further removals or effective fish passages are needed. Migrating nase may not be able to pass several weirs upstream e.g. the “Nehmühle”, “Steinmühle”, “Am Grün”, and “Wehrda” weirs near Marburg (Fig. 1). When barriers remain impassable, nase may show atypical behaviour. Instead of migrating upstream, stocked fish may migrate further downstream to search for suitable spawning areas. Such a downstream migration was investigated in translocated nase in the River Amblève, Belgium (Ovidio et al., 2016). This might be an additional reason why nase was not noticed in the upper Lahn in recent years. Furthermore, the significant heterozygous deficits of the Lahn population indicate a recent population size expansion or an influx of alleles from different gene pools (Luikart and Cornuet, 1998), probably because of asymmetric (downstream-biased) gene flow (Paz-Vinas et al., 2013).
Reintroduction and stock enhancement are often seen as a means of recovering threatened or lost species populations. However, stocks that have been bred in captivity can have deleterious effects on wild populations (e.g. Araki et al., 2007; Christie et al., 2012). In general, donor populations of high genetic variation are suitable for reintroduction of a species in lost parts of a species' range (Drauch et al., 2007; Forsman, 2014) since heterozygosity is often related to fitness (García-Navas et al., 2014; Feiner et al., 2017). Interbreeding of wild populations with hatchery fish that were hatchery-bred over multiple generations or are of a different gene pool may cause changes in the genetic structure of wild populations due to genetic drift, loss of genetic variability, reduced effective population size, outbreeding depression and general reduced fitness (e.g. Araki et al., 2007; Christie et al., 2012; Naish et al., 2013). Our results show that the diversity was comparatively lower in the H2 stock. Furthermore, this stock passed through a recent bottleneck and had a comparatively low effective population size, which in in turn may result in inbreeding (Naish et al., 2013) and reduced genetic variability. Since the H2 stock was created in 1984, adaptation to captivity might have led to a reduced fitness of the individuals when they were released to the wild (Frankham, 2005) into the upper River Lahn. All these named aspects make the H2 stock rather unsuitable as a donor population for reintroduction purposes.
Even though, the H2 stock originated initially from the River Nister itself, the Lahn was stocked with H2 individuals, and the geographical proximity of LA to NI is more pronounced compared to LA and any other investigated population, these three populations showed significant genetic differentiation in terms of FST and Φ ST . If the existing Lahn population really consisted solely of the stocked H2 individuals, we would not have expected significant genetic differences to H2 and NI or would at least expect them to be more similar to each other than compared to the other populations further upstream the River Rhine. Increased genetic differentiation is predicted for declining, isolated populations. Therefore, genetic distance between LA and NI suggests that populations drifted apart and gene flow did not fully counterbalance this effect, which might be reinforced by migration barriers between the rivers. It is likely that the re-established Lahn population does not exclusively originate from former stocking activities but established from autochthonous remainders or a combination of both. Another option would be that the H2 stock was actually not solely established by individuals from the River Nister but was also influenced by another different genetic pool. Otherwise, effects of captive breeding (reduced effective population size, selection) and separate development (genetic drift) of the stocks since creation of the H2 stock in 1984 must have put forth these differences.
The LA population was one of the most diverse ones investigated in this study and exhibited a unique and healthy genetic signature. It is obvious that fish from above the Rhine Falls from Lake Constance should definitively not be used for stocking the River Lahn or other water bodies below the Rhine Falls. For future reintroduction, a high number of spawners, ideally from the water body itself should be used. Using specimens of the self-reproducing population of the middle reaches of the River Lahn ensures to maintain the genetic characteristics (e.g. private haplotypes and alleles) and allows for a similar spawning time of the individuals. Spawning time is an adaptation to the temperature regime and has a strong genetic component (Schneider, 2011). Individuals could then directly be released into the upper Lahn before spawning actually takes place. This was already tested for nase, where downstream populations are still abundant (Ovidio et al., 2016), and would be an alternative to stocking with hatchery-bred fish.
To further increase the number of specimens for reintroduction, juveniles may be artificially cultured by stripping the spawning animals directly at the middle Lahn. Since 2017, this method is exclusively used for reintroduction in the upper Lahn in the context of an EU-LIFE project (Hübner et al., 2017). Fertilized eggs are then raised up to stock size at a hatchery farm before they are released into the upper Lahn.
5 Conclusion
Chondrostoma nasus is a highly specialized freshwater species with a pronounced migration behaviour and high gene flow within the River Rhine system. Overall, low to moderate levels of population differentiation were estimated and highest genetic variance was found within rather than between the investigated populations. Distinct populations were solely found where migration is interrupted e.g., in the Dornbirner Ach upstream the Rhine Falls. However, individual populations (e.g. the River Lahn) downstream the Rhine Falls also account for specific genetic characteristics that should be regarded in future conservation strategies. Despite the limitation of an unambiguous determination of the origin of the existing nase population in the middle Lahn, several key recommendations can be derived from the results of this study to guide the restoration of the nase in the upper Lahn and/or other waters of the River Rhine basin. For a successful re-establishment of a self-reproducing nase population in the upper River Lahn on the long term, ideally spawners from the middle reaches of the water body itself should be used. However, it is not only important to carefully choose genetically sufficient stocking material, but also to restore habitats for spawning and growth, and to enhance river connectivity. Considering these aspects will allow for an establishment of a population robust enough in size and sufficiently diverse.
Acknowledgements
We are grateful to all collaborators, namely Theresa Graf, Roman Fricke, Thomas Schmidt, and especially Dr. Götz Kuhn, Pascal Vonlanthen, and Franck Bonell for collecting samples. Many thanks to Britta Wahl-Ermel for her assistance in the laboratory. The fishery department of the Regierungspräsidium Gießen (Hesse) sponsored genetic analyses of the populations Lahn, Nister and H2.
References
- Aljanabi SM, Martinez I. 1997. Universal and rapid salt-extraction of high quality genomic DNA for PCR-based techniques. Nucleic Acid Res 25: 4692–4693. [CrossRef] [PubMed] [Google Scholar]
- Araki H, Cooper B, Blouin MS. 2007. Genetic effects of captive breeding cause a rapid, cumulative fitness decline in the wild. Science 318: 100–103. [Google Scholar]
- Bandelt H, Forster P, Röhl A. 1999. Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol 16: 37–48. [CrossRef] [PubMed] [Google Scholar]
- Barluenga M, Sanetra M, Meyer A. 2006. Genetic admixture of burbot (Teleostei: Lota lota) in Lake Constance from two European glacial refugia. Mol Ecol 15: 3583–3600. [CrossRef] [PubMed] [Google Scholar]
- Benjamini Y, Hochberg Y. 1995. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Series B Stat Methodol 57: 289–300. [Google Scholar]
- Brodersen J, Seehausen O. 2014. Why evolutionary biologists should get seriously involved in ecological monitoring and applied biodiversity assessment programs. Evol Appl 7: 968–983. [CrossRef] [PubMed] [Google Scholar]
- Cornuet JM, Luikart G. 1996. Description and power analysis of two tests for detecting recent population bottlenecks from allele frequency data. Genetics 144: 2001–2014. [PubMed] [Google Scholar]
- Costedoat C, Pech N, Chappaz R, Gilles A. 2007. Novelties in hybrid zones: Crossroads between population genomic and ecological approaches. PLoS ONE 4: e357. [Google Scholar]
- Christie MR, Marine ML, French RA, Blouin MS. 2012. Genetic adaptation to captivity can occur in a single generation. PNAS 109: 238–242. [CrossRef] [Google Scholar]
- Dubut V, Fouquet A, Voisin A, Costedoat C, Chappaz R, Gilles E. 2012. From late miocene to holocene: Processes of differentiation within the telestes genus (Actinopterygii: Cyprinidae). PLoS ONE 7: e34423–03. [Google Scholar]
- Drauch AM, Fisher BE, Latch EK, Fike JA, Rhodes OE. 2007. Evaluation of a remnant lake sturgeon population's utility as a source for reintroductions in the Ohio River system. Conserv Genet 9: 1195–1209. [Google Scholar]
- Duerregger A, Pander J, Palt M, Mueller M, Nagel C, Geist J. 2018. The importance of stream interstitial conditions for the early-life-stage development of the European nase (Chrondrostoma nasus L.). Ecol Freshw Fish 00: 1–13. [Google Scholar]
- Dußling U, Berg R. 2001. Fische in Baden-Württemberg. Ministerium für Ernährung und Ländlichen Raum Baden-Württemberg. Stuttgart, Germany. 176 p. [Google Scholar]
- Earl DA, Von Holdt BM. 2012. STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv Genet Res 4: 359–361. [CrossRef] [Google Scholar]
- Evanno G, Regnaut S, Goudet J. 2005. Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol 14: 2611–1620. [CrossRef] [PubMed] [Google Scholar]
- Excoffier L, Smouse PE, Quattro JM. 1992. Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data. Genetics 131: 479–491. [PubMed] [Google Scholar]
- Feiner ZS, DeWoody JA, Breck JE, Höök TO. 2017. Influences of multilocus heterozygosity on size during early life. Ecol Evol 7: 2142–2154. [Google Scholar]
- Flore L, Keckeis H. 1998. The effect of water current on foraging behaviour of the rheophilic cyprinid Chondrostoma nasus (L.) during ontogeny: evidence of a trade-off between energetic gain and swimming costs. Regul Rivers: Res Manage 14: 141–154. [CrossRef] [Google Scholar]
- Forsman A. 2014. Effects of genotypic and phenotypic variation on establishment are important for conservation, invasion, and infection biology. Proc Natl Acad Sci 111: 302–307. [CrossRef] [Google Scholar]
- Frankham R. 2005. Stress and adaptation in conservation genetics. J Evol Biol 18: 750–755. [Google Scholar]
- Freyhof J. 1997. Remarks on the status of Chondrostoma nasus in the River Rhine. Folia Zool 46 (Suppl. 1): 61–66. [Google Scholar]
- Freyhof J, Brooks E. 2011. European Red List of Freshwater Fishes. Luxembourg: Publications Office of the European Union. [Google Scholar]
- Fu YX. 1996. New statistical tests of neutrality for DNA samples from a population. Genetics 143: 557–570. [CrossRef] [PubMed] [Google Scholar]
- García-Navas V, Cáliz-Campal C, Ferrer ES, Sanz JJ, Ortego J. 2014. Heterozygosity at a single locus explains a large proportion of variation in two fitness-related traits in great tits: a general or a local effect? J Evol Biol 27: 2807–2819. [Google Scholar]
- Hauer C, Unfer G, Schmutz S, Habersack H. 2007. The importance of morphodynamic processes at riffles used as spawning grounds during the incubation time of nase (Chondrostoma nasus). Hydrobiologia 579: 15–27. [Google Scholar]
- Hmuklv, Hessen-Forst Fena (Eds.) 2014. Atlas der Fische Hessens − Verbreitung der Rundmäuler, Fische, Krebse und Muscheln. In: FENA Wissen Ed. 2, Gießen, Wiesbaden, Germany. [Google Scholar]
- Hübner D, Fricke R. 2011. Pilotprojekt Aalmonitoring im Mittellauf der Lahn − Aufnahme von Aalbestand und Habitatbedingungen. Studie im Auftrag des Regierungspräsidiums Gießen − obere Fischereibehörde. 78 p. [Google Scholar]
- Hübner D, Fricke R. 2012. Pilotprojekt Aalmonitoring im Mittel- und Unterlauf der Lahn − Aufnahme von Aalbestand und Habitatbedingungen. Studie im Auftrag der Struktur und Genehmigungsdirektion Nord − obere Fischereibehörde. 73 p. [Google Scholar]
- Hübner H, Fricke R. 2014. Pilotprojekt Aalmonitoring. Optimiertes Aal-Besatzmanagement in der Lahn. Maßnahmen zur Steigerung des Besatzerfolges. Studie im Auftrag des Landes Hessen vertreten durch das Regierungspräsidium Gießen − obere Fischreibehörde. 57 p. [Google Scholar]
- Hübner D, Cramer C, Schmidt T. 2016. Wiederansiedlung der Nase (Chondrostoma nasus) im Oberlauf der Lahn. Studie im Auftrag des Regierungspräsidiums Gießen − obere Fischereibehörde. 46 p. [Google Scholar]
- Hübner D, Fricke R, Graf T. 2017. Maßnahmen zur Stützung der Bestände der kieslaichenden Fischarten Äsche und Nasen in der Oberen Lahn. EU-LIFE14 IPE/DE/022_C7B/D4.1. Zwischenbericht. Studie im Auftrag des Landes Hessen, Regierungspräsidium Gießen. 59 p. [Google Scholar]
- Hudson AG, Vonlanthen P, Seehausen O. 2014. Population structure, inbreeding and local adaptation within an endangered riverine specialist: the nase (Chondrostoma nasus). Conserv Genet 15: 933–951. [Google Scholar]
- Jost L. 2008. GST and its relatives do not measure differentiation. Mol Ecol 17: 4015–4026. [CrossRef] [PubMed] [Google Scholar]
- Jost L. 2010. The relation between evenness and diversity. Diversity 2: 207–232. [Google Scholar]
- Kalinowski ST. 2004. Counting alleles with rarefaction: Private alleles and hierarchical sampling designs. Conserv Genet 5: 539–543. [Google Scholar]
- Kalinowski ST. 2005. HP-RARE 1.0: a computer program for performing rarefaction on measures of allelic richness. Mol Ecol Notes 5: 187–189. [Google Scholar]
- Kimura M, Ohta T. 1978. Stepwise mutation model and distribution of allelic frequencies in a finite population. Proc Natl Acad Sci USA 75: 2868–2872. [CrossRef] [Google Scholar]
- Kitada S, Nakajima K, Hamasaki K, Shishidou H, Waples RS, Kishino H. 2019. Rigorous monitoring of a large-scale marine stock enhancement program demonstrates the need for comprehensive management of fisheries and nursery habitat. Sci Rep 9: 5290. [CrossRef] [PubMed] [Google Scholar]
- Leigh JW, Bryant D. 2015. PopART: Full-feature software for haplotype network construction. Methods Ecol Evol 6: 1110–1116. [Google Scholar]
- Luikart G, Allendorf FW, Cornuet JM, Sherwin WB. 1998. Distortion of allele frequency distributions provides a test for recent population bottlenecks. J Hered 89: 238–247. [Google Scholar]
- Luikart G, Cornuet JM. 1998. Empirical evaluation of a test for identifying recently bottlenecked populations from allele frequency data. Conserv Biol 12: 228–237. [Google Scholar]
- Lundmark C, Sandström A, Andersson C, Laikre L. 2019. Monitoring the effects of knowledge communication on conservation managers' perception of genetic biodiversity − A case study from the Baltic Sea. Mar Policy 99: 223–229. [Google Scholar]
- Mäkinen HS, Merilä J. 2008. Mitochondrial DNA phylogeography of the three-spined stickleback (Gasterosteus aculeatus) in Europe − evidence for multiple glacial refugia. Mol Phylogenet Evol 46: 167–182. [Google Scholar]
- Marandel F, Lorance P, Berthelé O, Trenkel VM, Waples RS, Lamy JB. 2019. Estimating effective population size of large marine populations, is it feasible? Fish Fisher 20: 189–198. [CrossRef] [Google Scholar]
- Mesquita N, Cunha C, Hänfling B, Carvalho GR, Zé-Zé L, Tenreiro R, Coelho MM. 2003. Isolation and characterization of polymorphic microsatellite loci in the endangered Portuguese freshwater fish Squalius aradensis (Cyprinidae). Mol Ecol Notes 3: 572–574. [Google Scholar]
- Muenzel FM, Santera M, Salzburger W, Meyer A. 2007. Microsatellites from the vairone Leuciscus souffia (pisces: Cyprinidae) and their application to closely related species. Mol Ecol Notes 7: 1048–1050. [Google Scholar]
- Nagel C, Pander J, Mueller M, Geist J. 2019. Substrate composition determines emergence success and development of European nase larvae (Chondrostoma nasus L.). Ecol Freshw Fish 00: 1–11. [Google Scholar]
- Naish KA, Seamons TR, Dauer MB, Hauser L, Quinn TP. 2013. Relationship between effective population size, inbreeding and adult fitness-related traits in a steelhead (Oncorhynchus mykiss) population released in the wild. Mol Ecol 22: 1295–1309. [CrossRef] [PubMed] [Google Scholar]
- Ovidio M, Hanzen C, Gennotte V, Michaux J, Benitez J-P, Dierckx A. 2016. Is adult translocation a credible way to accelerate the recolonization process of Chondrostoma nasus in a rehabilitated river? Cybium 40: 43–49. [Google Scholar]
- Paetkau D, Slade R, Burden M, Estoup A. 2004. Genetic assignment methods for the direct, real-time estimation of migration rate: a simulation-based exploration of accuracy and power. Mol Ecol 13: 55–65. [CrossRef] [PubMed] [Google Scholar]
- Paz-Vinas I, Comte L, Chevalier M, Dubut V, Veyssiere C, Grenouillet G, Loot G, Blanchet S. 2013. Combining genetic and demographic data for prioritizing conservation actions: insights from a threatened fish species. Ecol Evol 3: 2696–2710. [CrossRef] [PubMed] [Google Scholar]
- Paz-Vinas I, Quéméré E, Chikhi L, Loot G, Blanchet S. 2013. The demographic history of populations experiencing asymmetric gene flow: combining simulated and empirical data. Mol Ecol 22: 3279–3291. [CrossRef] [PubMed] [Google Scholar]
- Peakall R, Smouse PE. 2012. GENALEX 6.5: Genetic analysis in excel. Population genetic software for teaching and research − An update. Bioinformatics 28: 2537–2539. [CrossRef] [PubMed] [Google Scholar]
- Peňáz M. 1996. Chondrostoma nasus − its reproduction strategy and possible reasons for a widely observed population decline − a review. In: Kirchhofer A, Hefti D. (eds.). Conservation of Endangered Freshwater Fish in Europe. ALS Advances in Life Sciences. Birkhäuser Basel. [CrossRef] [Google Scholar]
- Perea S, Böhme M, Zupančič P, Freyhof J, Šanda R, Özuluğ M, Abdoli A, Doadrio I. 2010. Phylogenetic relationships and biogeographical patterns in circum-mediterranean subfamily Leuciscinae (Teleostei, Cyprinidae) inferred from both mitochondrial and nuclear data. BMC Evol Biol 10: 265. [CrossRef] [PubMed] [Google Scholar]
- Petit RJ, El Mousadik A, Pons O. 1998. Identifying populations for conservation on the bases of genetic markers. Conserv Biol 12: 844–855. [Google Scholar]
- Piccolo JJ, Unfer G, Lobón-Cerviá J. 2018. Why conserve native brown trout? In: Lobón-Cerviá J, Sanz N. (eds.) Brown Trout: Biology, Ecology and Management, First Edition, John Wiley & Sons Ltd. [Google Scholar]
- Piry S, Luikart G, Cornuet JM. 1999. BOTTLENECK: A computer program for detecting recent reductions in the effective population size using allele frequency data. J Hered 90: 502–503. [Google Scholar]
- Piry S, Alapetite A, Cornuet JM, Paetkau D, Baudouin L, Estoup A. 2004. GENECLASS2: A Software for genetic assignment and first-generation migrant detection. J Hered 95: 536–539. [Google Scholar]
- Pritchard JK, Stephens M, Donnelly P. 2000. Inference of population structure using multilocus genotype data. Genetics 155: 945–959. [Google Scholar]
- Ramos-Onsins SE, Rozas J. 2002. Statistical properties of new neutrality tests against population growth. Mol Biol Evol 19: 2092–2100. [CrossRef] [PubMed] [Google Scholar]
- Raymond M, Rousset F. 1995. GENEPOP (version 1.2): population genetics software for exact tests and ecumenicism. J Heredity 86: 248–249. [CrossRef] [Google Scholar]
- Schindler DE, Hilborn R, Chasco B, Boatright CP, Quinn TP, Rogers LA, Webster MS. 2010. Population diversity and the portfolio effect in an exploited species. Nature 465: 609–612. [CrossRef] [PubMed] [Google Scholar]
- Schneider J. 2011. Review of reintroduction of Atlantic salmon (Salmo salar) in tributaries of the Rhine River in the German Federal States of Rhineland-Palatinate and Hesse. J Appl Ichthyol 27: 24–32. [Google Scholar]
- Schwevers U, Adam A. 1997. Wiederansiedlung der Nase in der Lahn. Im Auftrag der IG-Lahn. 30 p. [Google Scholar]
- Tajima F. 1989. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics 123: 585–595. [PubMed] [Google Scholar]
- Van Oosterhoudt C, Hutchinson WF, Wills DPM, Shipley P. 2004. Micro-Checker: Software for identifying and correcting genotyping errors in microsatellite data. Mol Ecol Notes 4: 535–538. [Google Scholar]
- Vaughn CC. 2010. Biodiversity losses and ecosystem function in freshwaters: Emerging conclusions and research directions. BioScience 60: 25–35. [Google Scholar]
- Villéger S, Blanchet S, Beauchard O, Oberdorff T, Brosse S. 2011. Homogenization patterns of the world's freshwater fish faunas. PNAS 108: 18003–18008. [CrossRef] [Google Scholar]
- Vonlanthen P, Hudson A, Seehausen O. 2011. Genetische Differenzierung und lokale Anpassung der Nasenpopulationen in der Schweiz. Im Auftrag des Bundesamt für Umwelt (BAFU), Kastanienbaum, Switzerland. 42 p. [Google Scholar]
- Vyskocilová M, Simková A, Martin J-F. 2007. Isolation and characterization of microsatellites in Leuciscus cephalus (Cypriniformes, Cyprinidae) and cross-species amplification within the family Cyprinidae. Mol Ecol Notes 7: 1150–1154. [Google Scholar]
- Waples RS, Do C. 2010. Linkage disequilibrium estimates of contemporary Ne using highly variable genetic markers: a largely untapped resource for applied conservation and evolution. Evol Appl 3: 244–262. [CrossRef] [PubMed] [Google Scholar]
- Waples RS, Peel D, Mecbeth GM, Tillett BJ, Ovenden JR. 2014. NeEstimator v2: re-implementation of software for the estimation of contemporary effective population size (Ne) from genetic data. Mol Ecol Res 14: 209–214. [CrossRef] [Google Scholar]
- Wetjen M, Cortey M, Vera M, Schmidt T, Schulz R, García-Marín J-L. 2017. Occurrence of length polymorphism and heteroplasmy in brown trout. Gene Reports 6: 1–7. [Google Scholar]
- Wetjen M, Schmidt T, Schrimpf A, Schulz R. 2020. Genetic diversity and population structure of burbot Lota lota in Germany: Implications for conservation and management. Fish Manag Ecol 27: 170–184. [CrossRef] [Google Scholar]
- Winter HV, Fredrich F. 2003. Migratory behaviour of ide: a comparison between the lowland rivers Elbe, Germany, and Vecht, the Netherlands. J Fish Biol 63: 871–880. [Google Scholar]
Cite this article as: Wetjen M, Hübner D, Seehausen O, Schulz R. 2020. Genetic diversity of endangered Chondrostoma nasus in the River Rhine system: Conservation genetics considerations on stocking and reintroduction. Knowl. Manag. Aquat. Ecosyst., 421, 25.
All Tables
Sampling sites and number of analysed Chondrostoma nasus individuals per population (Pop) at the sequence (NmtDNA) and microsatellite (NMS) level in the River Rhine system.
Summary of genetic diversities within the investigated Chondrostoma nasus populations. For each population (Pop) true diversity (DMS), allelic richness (AR), private allelic richness (PAR), observed (HO) and expected heterozygosity (HE), fixation index (F), true diversity (DmtDNA), number of haplotypes (NHap), haplotype diversity index (h), number of private haplotypes (PHap), segregating sites (S) and the nucleotide diversity index (π) is reported.
Neutrality tests and evidence for recent population bottlenecks using the sign test and the Wilcoxon test of investigated populations (Pop) of Chondrostoma nasus. DT = Tajima's D; FS = Fu's Fs ; R 2 = Ramos-Onsins and Rozas R 2; TPM = two-phase model; SMM = stepwise mutation model. Significant values are in bold.
Analysis of molecular variance (AMOVA) for the sequence (mtDNA Cytb) and microsatellite (MS) data of Chondrostoma nasus. Shown is the degree of freedom (df), the sum of squares (Sum of sq.) and the variation in percent for each source of variation.
Change in the effective number of distinct populations (ΔST ) excluding individual populations from the entire dataset (all Pops). ΔST and the change (difference) are shown for both, the sequence (mtDNA) and microsatellite (MS) data.
Pairwise ΦST- (above the diagonal) and F ST-values (below the diagonal) of the investigated Chondrostoma nasus populations (Pop). Significant values are in bold.
Results of the self-assignment test and first generation migrant (F0) analysis. Each row contains the individuals from one sampling location and the columns indicate the populations to which these individuals were assigned (highest likelihood). Total number of investigated individuals per population (Total N), percentage of correctly assigned individuals (Corr. Assign. %) and the number of first generation migrants with their most likely source population in brackets for the two tests statistics Lhome and Lhome/Lmax are given. For the first generation migrant analysis, both hatchery stocks (H1 and H2) were excluded from the analysis.
All Figures
Fig. 1 The course of the River Lahn and its obstacles impeding fish migration. Stocking sites of Chondrostoma nasus in the years 1995–1997 and 2014–2015 were located between the upper course (near Kernbach) and the middle course (near Steeden) of the River Lahn. |
|
In the text |
Fig. 2 Sampling locations and observed haplotypes per sampling site of Chondrostoma nasus in the River Rhine basin. The size of the circles is proportional to the number of analysed individuals per sampling site. |
|
In the text |
Fig. 3 Median-joining network of 154 individuals of Chondrostoma nasus of the River Rhine basin including reference haplotypes of Chondrostoma angorense (Hap_12) and Parachondrostoma toxostoma (Hap_13). The size of the circles is proportional to the frequency of fish carrying the respective haplotype. Colours refer to the origin-sampling site. Mutational steps are represented by intermediated dashes. Black dots represent missing haplotypes. |
|
In the text |
Fig. 4 Bayesian clustering analysis of the investigated Chondrostoma nasus populations: H2 = Hatchery_2, LA = Lahn, NI = Nister, DO = Dornbirner Ach, WI = Wiese, MU = Murg, MI = Minthesee, NK = northern Ketscher, H1 = Hatchery_1. Colours refer to the predicted clusters (1–5). Each individual is displayed by a vertical column with the length proportional to the individual's estimated affiliation to the respective cluster. |
|
In the text |
Fig. 5 Plots of the first two axes of a principal component analysis (PCoA) of Chondrostoma nasus populations using GenAlEx v. 6.5. Each coloured dot represents one population: H1 = Hatchery_1, H2 = Hatchery_2, NI = Nister, LA = Lahn, NK = northern Ketscher, MI = Minthesee, MU = Murg, WI = Wiese, DO =Dornbirner Ach. |
|
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.