Temporal analysis shows relaxed genetic erosion following improved stocking practices in a subarctic transnational brown trout population

Maintaining standing genetic variation is a challenge in human-dominated landscapes. We used genetic (i.e., 16 short tandem repeats) and morphological (i.e., length and weight) measurements of 593 contemporary and historical brown trout (Salmo trutta) samples to study fine-scale and short-term impacts of different management practices. These had changed from traditional breeding practices, using the same broodstock for several years, to modern breeding practices, including annual broodstock replacement, in the transnational subarctic Pasvik River. Using population genetic structure analyses (i.e., Bayesian assignment tests, DAPCs, and PCAs), four historical genetic clusters (E2001A-D), likely representing family lineages resulting from different crosses, were found in zone E. These groups were characterized by consistently lower genetic diversity, higher within-group relatedness, lower effective population size, and significantly smaller body size than contemporary stocked (E2001E) and wild fish (E2001F). However, even current breeding practices are insufficient to prevent genetic diversity loss and morphological changes as demonstrated by on average smaller body sizes and recent genetic bottleneck signatures in the modern breeding stock compared to wild fish. Conservation management must evaluate breeding protocols for stocking programs and assess if these can preserve remaining natural genetic diversity and morphology in brown trout for long-term preservation of freshwater fauna.

www.nature.com/scientificreports/ locally adapted populations 5,8 . Further, changes in population genetic structure caused by release of translocated broodstocks, stocking, and accidental escapes from fish farms has been in the focus of fishery research 5,8 . Other anthropogenic impacts, like overharvesting and the destruction of natural spawning and nursery areas, can lower reproduction resulting in decreasing population sizes and thus, decreased genetic diversity over time 4 . In addition, hydroelectric dams increase isolation and genetic divergence between river zones 4 . Negative impacts of anthropogenic infrastructure may have been unforeseen in some cases, but often, such as with construction of hydroelectric dams, they were to some extent foreseen. Moreover, at least in Norway, dam construction projects were often coupled with positivistic stocking programs as mitigation and replacement for lost breeding potentials as part of the concession 9 . Different hatchery breeding strategies likely have different impacts on the genetic composition of released fish. For example, rearing repeatedly from a single captive broodstock without replacing individuals, referred to in this paper as traditional multi-generational supportive breeding, can result in inbreeding patterns and genetic drift because the same individuals are crossed and their offspring released to the wild where they can breed with each other 10 . By contrast, modern single-generational supportive breeding, with annual replacement of the parental fish with wild-caught individuals that are used only once for breeding purposes (i.e., one generation), is expected to have fewer negative impacts because the number of breeders should theoretically be larger and therefore, genetic diversity should be better preserved. Furthermore, the fact that the fish are not as long in captivity should reduce domestication effects and selection for hatchery conditions, which in turn may lead to lower survival and reduced fitness of released fish in the wild 11,12 . Brown trout in the subarctic Pasvik River in northern Norway. The Pasvik River is a transnational subarctic river system ( Fig. 1) known for its large-growing brown trout (Salmo trutta L. 1758). As shown by 4 , the natural genetic integrity of the population is threatened by habitat destruction, fragmentation caused by hydropower developments, and the long-term practice of supportive stocking. Different management practices in Norway and Russia have led to distinct genetic diversity patterns in this population. Mainly, river zones situated in Russia are not stocked, whereas Norwegian-Russian zones have been stocked with ~ 5000 offspring of local specimens annually since the 1970s to supplement natural recruitment 13,14 . Stocking was introduced because seven hydroelectric dams without fish passes were built between 1932 and 1978. The breeding program was initiated as a requirement for building the last dam (i.e., Melkefoss, Fig. 1). The dams fragmented the river habitat and destroyed natural spawning and nursery grounds, reducing brown trout recruitment potential, and creating barriers to migration and gene flow among river zones. The Norwegian-Russian river zones nowadays contain an estimated 70-90% stocked fish [15][16][17] . Furthermore, some of these stocked river zones show signs of recent bottlenecks, suggesting that they are losing genetic variation 4 .

Figure 1.
Overview of study area. The different river parts that are separated by hydroelectric dams (orange) are denoted by different letters (A-J). This is a modified map from 4  www.nature.com/scientificreports/ The wild large-growing brown trout population in this watercourse may be facing two major challenges concerning the preservation of genetic integrity: reduced effective population size and potential mating with stocked fish 18,19 . Both factors can lead to allele frequency shifts and losses along with loss of local population structure. Different management and breeding practices can also manifest themselves in morphological changes with different outcomes. Generally, the rate of morphological change introduced by human activity outpaces natural causes 20,21 . In fish species, multigenerational breeding of the same broodstock may lead to stronger deviations from the wild phenotype than single-year breeding programs because of domestication effects 22 , potentially including rapid evolutionary effects 23 . Further, breeding of closely related specimens, either due to a small start population or interbreeding of generations, may lead to inbreeding depression, for instance, for adult size and weight 24 . Inbreeding depression can negatively affect long-term survival of populations and thereby compromises supportive breeding programs.
To test this, we applied a temporal study design to estimate genetic and demographic changes during 17 years of stocking (2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015)(2016)(2017)(2018). Hypotheses included declining genetic diversity and/or allele frequency changes as well as bottlenecks and reductions in effective population sizes through time due to increasing stocking impacts (i.e., homogenization effect and/or replacement). Furthermore, we investigated the differential effects of traditional and modern supportive breeding practices, which shifted in the mid-1990s from a multi-to a single-generational (annually renewed) broodstock of local trout specimens. Overall, the aim was to study how different breeding practices have influenced genetic variation and demographics through time in brown trout.

Methods
Sample collection. In total, the data set consisted of 593 tissue samples of brown trout from the Pasvik River. From two zones we had time series data including 303 historical samples collected in 2001 and 2007-08 from river zones E and I (Fig. 1 Fig. 1 (e.g., A, B, C, etc.) while for historical samples, we additionally indicate the collection year as mentioned above. The study was performed in strict accordance with Norwegian legislation. Fish were euthanized by means of cerebral concussion prior to sample collection. A fishing license is required from the fishing right owner. Accordingly, we obtained fishing permissions for the Pasvik watercourse from the County Governor of Finnmark with legal authority through LOV 1992-05-15 nr 47, §13. No ethical permission is required from the Norwegian Animal Research Authority for collection with gill nets and the associated sacrifice of fish (FOR 1996-01-15 nr 23, Norwegian Ministry of Agriculture and Food). Additional samples were obtained by opportunistic collection of samples from local anglers with permission to fish in the river for personal sustenance. Molecular methods. All samples were genotyped at 16 short-tandem repeat (STR) loci and PCR conditions, electrophoresis with an Applied Biosystems 3730xl Genetic Analyzer (Applied Biosystems, UK), and allele scoring are detailed in Supplementary Material S1 and Klütsch et al. 4 .  26 using the Markov chain method with 10,000 dememorization steps, 5000 batches, and 10,000 iterations to estimate exact P-values. The same settings were employed to test for linkage disequilibrium between pairs of loci. Standard summary statistics, including allelic proportions per STR locus for the time series data, were calculated with GenAlEx 6.51b2 27,28 . ADZE 1.0 29 was used to estimate allelic richness (A R ) and private allelic richness (A RP ) corrected for differences in sample sizes when considering the different zones separately (standardized sample size of 10) and at the group level based on genetic clusters, stocking history, and different time periods (standardized sample size of 50). For the second analysis, six major groups were therefore identified: (i) the Russian zones B-C as no supportive breeding occurs there, (ii) E2008/2007 because these samples are from an intermediate time period, (iii) E2001A-E2001D as representative of the traditional breeding practice (i.e., multigenerational), (iv) E2001E and E2001F as these represent the modern stock and the stock that most closely resembles wild trout in the river system (we did not separate these two as some admixture occurs between them), (v) zones E-I as these zones are regularly stocked, and (vi) I2007 because these samples represent an intermediate time period.
Relatedness. Since stocking is typically performed with offspring from a small subset of individuals, it is expected that stocked fish show closer family relationships. In addition, in salmonid species, the juveniles are often found in shoals/schools of related individuals 30 and it has been shown that detection of genetic clusters in STRU CTU RE 31 and DAPC 32 can be caused by groups of closely related individuals in the data set 30 . Although we only sampled adult fish 27-70 cm in length, encompassing several age classes, we estimated relatedness based on allele frequencies with the R package RELATED v.1.0 33 to test whether a higher-than-expected relatedness level is present in any of the groups studied. First, we used the compareestimator function to compare the performance of four different relatedness estimators. The most fitting relatedness estimator based on the highest correlation coefficient (R) was used to estimate relatedness within spatial and temporal groups and to estimate pairwise relatedness within those groups that showed higher-than-expected relatedness.
Demography. We used the program BOTTLENECK 1.2.02 34 to assess whether recent genetic bottlenecks, linked to breeding practices over the last 20-30 years, could be found. The underlying algorithm assumes that allelic diversity is lost faster than heterozygosity. Hence, it tests for heterozygosity excess compared to expectations at mutation-drift equilibrium 35 . We tested two mutation models, the infinite-alleles-model (IAM) and the two-phase-model (TPM) following the recommendations by 36 . The latter encompasses two mutation models, the IAM and the stepwise mutation model (SMM), and different proportions of STRs can display either one of these mutation models. The TPM model was run three times for each population with the percentage of stepwise mutations being 20%, 50%, and 70%, respectively. We used the 1-way Wilcoxon sign-rank test 37 and 10,000 iterations to determine significance.
Effective population size. The effective population size may decrease because many individuals from the same parents are released relative to the natural recruitment in the river. To test for differences of effective population size in a spatio-temporal context as well as investigate the effects of different breeding practices (i.e., multiversus single-generational), we used the program NeEstimator 2.1 38 . We calculated the effective population size (Ne) with three point estimators (i.e., Linkage Disequilibrium method (LD; 39 ), Heterozygote Excess method (HE; 40 , and the Molecular Coancestry method (MCA; 41 ). In all instances, we used allele frequencies ≥ 0.02 as recommended by 38 . We obtained confidence intervals with the jackknife method of 39 . We added the sibship frequency (SF) estimator as a fourth effective population size estimator and used COLONY 2.0.6.5 42 with default priors, because 43 showed through simulations that the most widely used LD method overestimates Ne when sample sizes are considerably smaller than the actual Ne, and that it can be inaccurate when the assumption of random mating is violated.
Population-genetic differentiation. Pairwise population genetic differentiation was estimated with G ST 44 and Jost's D 45 using GENALEX 6.51b2 with 9,999 random permutations 27,28 . To correct for multiple testing, the false discovery rate as described by Benjamini and Hochberg ( 46 -BH-FDR) was calculated with the p.adjust function in R version 3.6.0 47 to minimize Type I and II errors, while controlling FDR by using an ɑ-level of 0.01 48 .
STRU CTU RE 2.3.4 31 was run with the admixture model with correlated allele frequencies 49 for several datasets. An initial run included all samples from the entire river and all time periods. In addition, the two river sections with temporal samples (i.e., zones E and I) were separately analysed to check whether additional substructuring could be found and to confirm the initially identified genetic clusters. Using two sets of runs, one with the LocPrior option and one without, K = 1-10 was tested with 40 replicates each. Run parameters further included 1,000,000 MCMC steps and a burn-in period of 100,000. The LocPrior option allows for detection of weak population-genetic structure 50 . All STRU CTU RE analyses were carried out on the CIPRES Portal v3.3 (https:// www. phylo. org/; 51 that allows for parallelised calculation with the R package PARALLELSTRU CTU RE 52 . To summarize replicates and determine the most likely number of genetic clusters in the data set, the program STRU CTU RESELECTOR 53 was utilized to account for uneven sample sizes in the data set. Four different estimators introduced by 54 were calculated: the median of means (MedMeaK), maximum of means (MaxMeaK), median of medians (MedMedK), and maximum of medians (MaxMedK). Summary bar plots were prepared with the program CLUMPAK 55  www.nature.com/scientificreports/ As initial tests indicated that highly related individuals were present in some of the identified genetic clusters, we ran a second set of Bayesian assignment analyses with a reduced data set to account for potential false identification of genetic clusters by the presence of closely related specimens in the data set 56 . This second set of Bayesian assignment analyses had closely related individuals randomly reduced by ≥ 50% following recommendations by 30 , whose simulations showed that total purging of related individuals also leads to biases 30 .
In addition, a Discriminant Analysis of Principal Components (DAPC) was conducted to test whether temporal populations within zones E and I are genetically different using the ADEGENET package 32 in R version 3.6.0 47 . To avoid over-fitting, a cross-validation function with 100 replicates was used to determine the optimal number of principal components to be retained. Training versus test set proportions were 0.1 and 0.9, respectively.
Finally, a Principal Component Analysis (PCA) was carried out in the ADEGENET package 32 in R version 3.6.0 47 , which is interfaced with the package ADE 57 . This was done as an additional verification that identified groups in DAPC can also be found without pre-defined groupings.

Morphological variation.
For the historical samples in E2001, measures of body size, weight, and a morphology-based classification into stocked and wild fish were available (Table 2), which were used to test whether genetic clusters showed differences in these morphological traits and whether there was a correlation between genetic diversity and the origin of the fish. The classification into stocked versus wild fish was done by the fishermen, based on characteristic fin damage in stocked fish 14 . The accuracy of this classification was evaluated by by otolith analysis 17 of a large subsample (N = 202), resulting in the same classification as the fin damage method in 97.5% of the cases.
Initial statistical tests included tests for multicollinearity assessed by a Pearson correlation, homogeneity of variances assessed by Levene's test 58 , and tests for outliers with the identify_outliers() function in the rstatix package 59 . Normality was evaluated by plotting the correlation between the data and a normal distribution per genetic cluster in a QQ plot, because with high sample sizes, the typically applied Shapiro-Wilk test becomes sensitive to even minor deviations from normality. Subsequently, we used a Welch one-way ANOVA plus Games-Howell post hoc tests for all pairwise comparisons. P-values were adjusted for multiple testing using the Tuk-eyHSD() function in the rstatix package 59 .

Results
Genetic variation and summary statistics. Initial tests showed strong substructure within the historical population E2001, which was ultimately sorted into six subpopulations (referred to as E2001A to E2001F), based on high assignment scores, to break down linkage disequilibrium. In the final population set, linkage disequilibrium was significant (p < 0.05) in only 7/2051 pairwise comparisons after Bonferroni correction. Deviations from Hardy-Weinberg equilibrium were observed in one locus (SsaD157) in two populations (zone I and I2007) after Bonferroni correction.
Heterozygosity values (Table 1) suggested that historical samples from zone E showed negative F IS estimates. Allelic richness and private allelic richness, as well as allelic proportions per locus were particularly low in the historical genetic clusters, except for E2001E and E2001F, which corresponded to the contemporary populations (Table 1, S15-S30). Four alleles were only found in E2001E and F and not in any other historical or contemporary groups.
MICRO-CHECKER identified only one locus (SsaD157) in one group (I2007), showing a significant presence of null alleles, and found no signs of large allele dropout at any locus. Furthermore, only one population (zone E) showed potential scoring errors due to stutter at one locus (OMM1152). Hence, all loci were retained for further analysis.
Tests for genetic clustering and differentiation. For all samples combined, STRU CTU RE and STRU CTU RESELECTOR determined K = 8 as the likely number of genetic clusters (Figs. 2 and 3). Spatially, genetic differentiation occurred between contemporary non-stocked and stocked river zones, consistent with 4 (Fig. 3). Temporally, elevated genetic structure occurred in historical compared to contemporary samples (Figs. 2 and  3), with three out of six genetic clusters found only among historical samples (i.e., E2001A, C, and D) (Fig. 3). Table 2. Overview of length and weight data for the historical samples in zone E, sorted according to genetic subgroups E2001A-F. In addition, the morphological classification into wild and stocked fish based on fin damage is given in percent and total numbers. For each measurement or classification, the number of individuals is given (i.e., N LENGTH , N WEIGHT , N CLASSIFICATION ).  www.nature.com/scientificreports/ Analysing only river zones with both historical and contemporary samples (i.e. zones E and I) confirmed these results, providing evidence that three genetic clusters went extinct during the 17-year study period (Fig. 3). For zone I, the historical sample was weakly genetically differentiated from the contemporary sample (Fig. 3). The DAPC analysis largely confirmed the STRU CTU RE results (Figs. 4a-f, 5), identifying three clearly differentiated historical populations (i.e., E2001A, E2001C, and E2001D; Fig. 4). Also, E2001B was differentiated but showed some overlap with both E2001E and E2001F (Fig. 4d-f). Analysing zone E groups alone, E2001A-D were genetically differentiated from both more recent historical samples E2007/E2008 and contemporary samples ( Fig. 4d-f), whereas E2001E-F clustered with contemporary populations. In zone I, the historical sample I2007 was genetically different from the contemporary sample (Fig. 5). Re-running the analyses with the reduced data set did not change the results, confirming that the observed patterns were not caused by high relatedness in some groups (Supplementary Material Fig. S1-S8).
The results of the PCA confirmed the DAPC results (Supplementary Material Fig. S9-S14), although distances between groups were smaller, as expected.
Estimates of population-genetic differentiation using G ST and D EST were largely significant in pairwise assessments (Supplementary Material Tables S1 and S2). Historical (stocked) genetic clusters E2001A-D were genetically differentiated from the contemporary stocked populations. In comparison, the historical groups E2001E-F showed relatively low genetic differentiation both from the more recent historical samples E2007/E2008/I2007 and most of the contemporary populations. Finally, the historical sample from I2007 was weakly but significantly differentiated from the contemporary sample in the same zone.
Relatedness, bottlenecks, and effective population size. For our data set, the relatedness estimator (rw) by 60 showed the highest correlation coefficient between observed and expected values (r = 0.839). Historical samples E2001A-D (Table 1) displayed considerably higher-than-expected relatedness estimates. Randomly removing 50% of the closely related individuals and re-analysis in STRU CTU RE and ADEGENET led to no loss of population structure (Supplementary Material Tables S1 and S2, Supplementary Material Figures S5-S8).
Both bottleneck tests ( Table 3) and estimates of effective population size (Table 4) pointed to reduced genetic variation and effective population sizes in the historical genetic clusters that went extinct (E2001A-D), compared to those that survived (E2001E-F). Of those that survived, E2007F showed no signs of bottlenecks and had the largest effective population size, whereas E2001E showed signs of a recent bottleneck and an intermediate effective population size. The results also confirmed previous findings 4 that non-stocked zones in Russia (i.e., B and C) do not show signs of demographic bottlenecks while stocked zones G, H, and I showed consistent signs Here, we show only runs without use of the LOCPRIOR option; runs for which the LOCPRIOR option has been used can be found in the supplementary material. On the left-hand side, STRU CTU RE bar plots for all investigated river zones (i.e., A-J) plus temporal groups (i.e., E2008, E2007, E2001A-E2001F, and I2007), are shown. On the right-hand side, additional STRU CTU RE bar plots for zone E and I are displayed that include temporal samples.  (Table 3). Regarding effective population sizes, the non-stocked Russian zones B and C showed either comparable (SF estimator) or higher estimates (LD estimator) than the stocked Norwegian-Russian zones E-I (Table 4).  www.nature.com/scientificreports/ Morphological data. Exploratory statistical tests showed strong multicollinearity between the variables weight and length (r = 0.9, p < 0.0001), but also differences in variances among the groups defined based on genetics (p < 0.0001). Therefore, we used a Welch one-way ANOVA test that does not assume homogeneity of variances in combination with a Games-Howell post hoc test for each variable separately. Prior to the analyses, outliers were removed from the data sets; however, for weight, after the first round of outlier removal, more outliers were found that remained in the analysis. For both length and weight, the ANOVA was significant (length: p < 0.0001; weight: p < 0.0001). For length, the pairwise post hoc tests showed that E2001E and E2001F were significantly different from all other genetic clusters, but not from each other (Fig. 6a). For weight, a similar result was found, but E2001E was not significantly different from E2001B and E2001C (Fig. 6b). Genetic clusters E2001A, E2001B, E2001C, and E2001D did not significantly differ in length and weight from one another.
Regarding the morphological classification of caught fish into wild and stocked, E2001F consisted of 63% wild fish, followed by E2001E with ~ 16% wild fish while genetic clusters E2001A-D showed very low proportions of wild fish (Table 2).

Discussion
Our results revealed that significant temporal changes in genetic structure, effective population size, and allelic composition occurred in the trout population of the Pasvik River during 17 years (i.e., 2001-2018). Three unique genetic clusters could only be found in historical samples within zone E, likely being the result of traditional breeding practices leading to considerable allele frequency shifts. A fourth genetic cluster became rare with time, although still present in low numbers and as admixed individuals at the end of the 17-year study period. In all these four genetic clusters that either disappeared or declined during the study period, genetic diversity measures and effective population sizes were distinctly lower than in the modern breeding stock and the wild population, which was consistent with indications of recent bottlenecks and considerably lower effective population sizes. Overall, these results suggest a loss of genetic diversity and a homogenization effect as well as isolation by time.
Interbreeding between the historical groups was marginal, except for 2001E and F. Previous studies in salmonids found similarly low introgression levels from domestic fish into the native population 61,62 . Potential causes include both lower rates of survival and reproduction and changes in the timing of maturity and spawning introduced by captive-breeding programs 11,[63][64][65][66][67] . Shifted spawning timing may also contribute to non-random mating or reproductive isolation by time (temporal assortative mating; 68 ), resulting in low admixture with wild populations 64 . If new parents are caught every year for breeding, these effects may be less severe, but they are still present 11,65 . For example, in a meta-analysis, relative reproductive success of early-generation hatchery-reared salmonids in the wild was only 50% on average compared to wild fish 65 . Importantly, genetic effects seem to be present even in F2 generation specimens that have lived their entire lives in natural environments 11 . Further, mortality rates may be higher in captivity-bred than wild fish because of less efficient foraging behaviour, leading to smaller body sizes and poor physical condition, and more risk-prone or exploratory behaviour that may result in higher predation and angling pressure 67 . The progeny of hatchery-raised trout may potentially also suffer from fewer resources at egg hatching in the wild, as early juvenile survival is positively correlated to egg size. Egg size typically increases with mother size and may be smaller than in wild fish for two reasons; (1) there may be a hatchery selection towards smaller egg size 23 , and (2) hatchery fish grow less than wild fish. The latter will also reduce the mothers' ability to make deep and well-protected nests 69 . These factors may select particularly against domestic trout that grew up in artificial environments and have reduced fitness in comparison to wild trout 11 , but also be detrimental to wild fish if the characters are introgressed back to the wild populations.
Introgression between domestic and native trout can be density-dependent, increasing with the relative size of the domestic population 8 . Although this may partly explain the resilience of the native group with the highest proportion of wild fish (E2001F), it does not explain why domestic fish did not interbreed. Christie et al. 10 pointed out that there is a trade-off between the number of stocked fish that should be allowed on the spawning Table 4. Effective population-size estimates for all zones and temporal groups identified. We used four different estimators: linkage disequilibrium method (LD 38,39 ), heterozygote excess method (HE 40 , molecular coancestry method (MCA 41 ), and the sibship frequency estimator (SF 43 ). For all estimators, the point estimate as well as 95% confidence intervals (CI) are given. In addition, we provide the effective population size estimates for the groups that had 50% of closely related individuals removed to check whether this has any considerable effect on these estimates; these estimates are labelled with the suffix _R. www.nature.com/scientificreports/ grounds because releasing very few might not help the goal of supportive breeding whereas allowing too many will reduce the effective population size considerably. There were no obvious geographical barriers that can explain the strong separation of genetic clusters within zone E. The relaxation of linkage disequilibrium observed when sorting into different groups and the elevated relatedness of individuals in four out of six groups (i.e., E2001A-D), suggest non-random mating. Ecological, behavioural, historical, and anthropogenic factors may explain genetic sub-structuring of fish populations 70 . We can exclude the possibility of food partitioning between juveniles and adults as only adults were collected. Furthermore, as coregonids (mainly vendace Coregonus albula) are the main prey, constituting > 90% of the diet of adult brown trout in the study system 14,15 , we can also exclude the possibility of trophic polymorphism as a mechanism for population sub-structuring.
Our results are consistent with these studies, given that the four genetic clusters E2001A-E2001D represent domestic multi-generational broodstocks. The clear genetic separation and low admixture rates of the genetic clusters in combination with the phenotypic classification of being stocked fish support this. Further, the high relatedness and low allelic richness give further credence to this interpretation. Differences in morphological traits, like body size and weight, between domestic and wild fish have also been found to be linked to hatcherybreeding 63 . However, size variation within those four groups suggests the presence of multiple generations. Thus, either a low reproduction rate was present, or several age classes were released; however, the latter seemed to have been stopped sometime in the 1980s (Pasvik Kraft, personal communication). In any case, reproductive success, www.nature.com/scientificreports/ particularly interbreeding with wild trout, was low. Also, stocked fish from groups E2001A-D did not interbreed with each other, suggesting that all released stocks had very low rates of survival and reproductive success. Our results suggest that particularly traditional stocking practices can lead to severe and rapid reductions in effective population sizes and allelic diversity, strong linkage disequilibrium, and increased relatedness in multiple-year breeding programs. Small numbers of breeders plus variation in reproductive success contribute to this pattern. This is consistent with findings by others (e.g., 10 ). The most likely explanation for the four observed genetic clusters is that they correspond to four brood years (i.e., four different instances of releases of hatchery-bred fish) with very low reproduction levels in the wild. These traditional breeding practices have also led to significantly smaller body sizes/weight than found in wild fish. Together with the common signatures of recent bottlenecks and low genetic diversity, inbreeding depression was a likely factor in the traditional breeding program.
The results indicated that the change to supportive breeding practices, using new parental fish every year to reduce the number of crossings per parent, has led to a 'modern' stocked group (i.e., E2001E), which showed higher allelic richness and effective population sizes than the previously stocked groups and which co-exists and interbreeds with wild brown trout in the river. Nevertheless, E2001F, the group that most closely resembled wild brown trout at the time, harboured private alleles that were not found in any modern samples, suggesting that genetic diversity was lost through time. Similarly, bottlenecked modern river zones indicate an ongoing loss of genetic diversity 4 . In addition, also population E2001E, suggestively representing the modern breeding practices, showed signs of bottlenecks. Thus, the effective size of the modern breeding stock is larger than in historical practices, but it is not large enough to prevent further loss of genetic diversity.
One remaining potential problem in this system is that annual catches for the broodstock come from the same zone (i.e., H) and that due to the dams, natural migration is limited or non-existent. This probably results in a situation in which the local fish population encompasses a high proportion of stocked fish, as shown by the current study and others 15,16 . Using the Bayesian assignment scores and under the assumption that the Russian non-stocked zones are genetically representative of wild fish, the estimated proportions of wild fish in the contemporary stocked Norwegian-Russian zones equalled 47.8% (zone E), 32.0% (zone F), 13.0% (zone G), 3.9% (zone H), and 8.9% (zone I). This gives an average of 16.2% of wild fish in the stocked zones, but with considerable variation between zones. For example, zone H showed the lowest proportion of wild fish; however, this is the zone in which breeding specimens are caught annually. To preserve natural genetic diversity, breeding practices will need to be refined to increase the number of unrelated breeders and the proportion of wild specimens in the breeding program. Four out of the six zone E populations (E2001B-E) showed largely consistent significant bottlenecks. The main currently surviving wild population in zone E (i.e., E2001F) showed neither elevated relatedness estimates and nor genetic bottlenecks. The other population (E2001E) corresponded to the main contemporary population of the Norwegian-Russian part of the river where stocking and genetic bottlenecks are common 4 . Several factors may explain these results. The bottlenecked historical populations from E2001B-D could already back then have suffered from considerable population decline and ultimately gone partially extinct, whereas the still occurring contemporary population E2001E may be impacted by stocking practices. Stocking and destruction of spawning habitats occurred also before 2001 and may have led to population declines. The only other surviving genetic cluster E2001F shows partially admixed ancestry with E2001E and this may explain the non-significant bottleneck tests. Effective population sizes in non-stocked Russian zones seem to be comparable to or higher than in stocked Norwegian-Russian zones. Thus, stocking does not appear to increase effective population sizes in any significant way, and in worst case scenarios, decreases them considerably 10 .
In conclusion, our study documents large impacts on the genetic structure and genetic diversity of brown trout populations from annual compensatory stockings following construction of hydropower dams. Four genetic clusters were morphologically identified as stocked fish and showed consistently lower genetic diversity, higher within-group relatedness, lower effective population sizes, and significantly smaller body sizes than the other two genetic clusters, which were genetically more similar to contemporary stocked and wild fish. Over a period of 17 years, three of these genetic clusters within one of the studied river sections went extinct and a fourth cluster distinctly declined, indicating that reproduction success of these groups was low. A change in breeding practices from multi-to single-generational breeding of the locally caught broodstock resulted in the maintenance of higher genetic diversity, but even these modernized breeding practices are not enough to prevent genetic diversity loss when compared to wild fish (see also 4 ). Hence, conservation management should carefully assess future breeding protocols to preserve remaining natural genetic diversity in brown trout in the Pasvik River and other rivers where compensatory stockings are implemented. Morphological and molecular screening including the determination of relatedness levels before breeding may be a viable option for quickly implementable conservation measures in this system (e.g., 71 ). However, this cannot replace other efforts like habitat restoration to increase natural reproduction. Alternatives include the abandonment of the breeding program altogether. In many cases this would be very drastic, and potentially detrimental to the fish population and human interest as well if not followed up by other measures. We strongly advise starting by improving conditions for natural reproduction, including increased connectivity by opening fish passages around hydropower installations to facilitate fish migration. Our findings have important implications for stocking program management and indicate that a detailed evaluation will be needed to find the best way forward for stocking practices where such programs exist.