Skip to main content

Long-term monitoring of a brown trout (Salmo trutta) population reveals kin-associated migration patterns and contributions by resident trout to the anadromous run

Abstract

Background

In species showing partial migration, as is the case for many salmonid fishes, it is important to assess how anthropogenic pressure experienced by migrating individuals affects the total population. We focused on brown trout (Salmo trutta) from the Guddal River in the Norwegian Hardanger Fjord system, which encompasses both resident and anadromous individuals. Aquaculture has led to increased anthropogenic pressure on brown trout during the marine phase in this region. Fish traps in the Guddal River allow for sampling all ascending anadromous spawners and descending smolts. We analyzed microsatellite DNA markers from all individuals ascending in 2006–2016, along with all emigrating smolts in 2017. We investigated (1) if there was evidence for declines in census numbers and effective population size during that period, (2) if there was association between kinship and migration timing in smolts and anadromous adults, and (3) to what extent resident trout were parents of outmigrating smolts.

Results

Census counts of anadromous spawners showed no evidence for a decline from 2006 to 2016, but were lower than in 2000–2005. Estimates of effective population size also showed no trends of declines during the study period. Sibship reconstruction of the 2017 smolt run showed significant association between kinship and migration timing, and a similar association was indicated in anadromous spawners. Parentage assignment of 2017 smolts with ascending anadromous trout as candidate parents, and assuming that unknown parents represented resident trout, showed that 70% of smolts had at least one resident parent and 24% had two resident parents.

Conclusions

The results bear evidence of a population that after an initial decline has stabilized at a lower number of anadromous spawners. The significant association between kinship and migration timing in smolts suggests that specific episodes of elevated mortality in the sea could disproportionally affect some families and reduce overall effective population size. Finally, the results based on parentage assignment demonstrate a strong buffering effect of resident trout in case of elevated marine mortality affecting anadromous trout, but also highlight that increased mortality of anadromous trout, most of which are females, may lower overall production in the system.

Peer Review reports

Introduction

Individuals within species can exhibit different life history strategies which are often associated with important phenotypic variation, can differ between sexes and overall have pervasive ecological implications [1]. Hence, individuals representing different life history types may differentially allocate their amount of available energy between maintenance and reproduction functions to maximise their fitness. This polymorphism in life history strategies is maintained within species because their costs and benefits vary according to the environmental contexts [1]. As an example, within many species known to undertake migrations, some individuals migrate while others from the same population remain on the same site across their lifespan, referred to as partial or facultative migration [2, 3].

Among fishes, many salmonid species show anadromous life history forms, which means that juveniles hatch in freshwater and undertake feeding migrations at sea before returning to freshwater for spawning [4]. Their populations often include both sea-migratory and resident individuals that remain in freshwater, therefore showing facultative anadromy [5]. Coexistence between resident and migratory life-history strategies involves a fine balance between their respective costs and benefits [6]. Increased food availability in marine environments may lead to better growth and higher fecundity of anadromous individuals [7, 8]. On the other hand, residency can be advantageous when costs of migration become higher than benefits, due to factors such as predation risk, additional exposure to pathogens and parasites, or energetic costs for the migration process itself [6, 9].

Facultative anadromy is usually considered a quantitative trait, controlled by the action of multiple genes and their interaction with environmental factors [10, 11]. However, recent studies have shown that traits related to migration and life history in some salmonid species can be under control of single genes [12,13,14], whereas other studies point towards important elements of epigenetic regulation [15]. It is furthermore noteworthy that proportions of migrants and residents within a population may vary across years according to environmental factors or anthropogenically induced disturbances [16,17,18].

Migratory species, and not least salmonids, may be particularly susceptible to anthropogenic impact due to their dependence on several different habitats and connectivity between them [19]. For instance, fishing pressure in the sea and decreased access to marine environments due to dams represent important issues [20, 21]. Moreover, emerging threats related to climate change altering marine temperature regimes and adverse effects of marine aquaculture, such as accummulation of parasites that subsequently infect wild populations have become increasingly important [22,23,24,25]. In addition to general population declines, increased mortality at sea could also disproportionally affect the resident or migratory components of populations showing facultative anadromy. Also, in the case of sex-ratio differences between life-history types [26], increased mortality of one sex could reduce the total effective population size, resulting in a lower ratio between effective and census population size and leading to increased inbreeding and loss of genetic variation [27]. Finally, some studies have suggested that related fish tend to group together during migration [28, 29]. In addition to active association of kin [29], this could also reflect the mere fact that closely related individuals, especially smolts, may be “physiologically timed” to migrate at the same time [30]. If related individuals migrate together and are subjected to specific incidences of e.g., exposure to parasites, this could lead to high variance of survival among families and ultimately increase the variance of their reproductive success, a factor also leading to decreased effective population size [27].

The brown trout (Salmo trutta) is a species often showing facultative anadromy within populations [26, 31]. Recent studies have shown that environmental factors such as water temperature and food availability for brown trout juveniles can have contrasting effects on their migration tendency, with food limitation generally favoring anadromy and increasing temperature favoring residency [32,33,34]. These environmental factors interact with inherited genetic factors to shape life history strategies [33, 35, 36]. Brown trout shows a sex ratio typically skewed towards females among anadromous individuals and towards males among resident individuals, the latter including precocious mature male parr that can successfully fertilize eggs by adopting a sneaking behaviour [5, 26, 37,38,39]. Fjords in Norway have experienced an increased establishment of Atlantic salmon (Salmo salar) farms since the 1990s [40]. The aquaculture industry poses major problems for wild salmon and brown trout during marine migration [41]. Hence, the high concentration of farmed fish attracts and accummulates parasites such as sea lice Lepeophtheirus salmonis that subsequently infect wild salmon and trout passing by, often in lethal doses [23, 42,43,44,45,46]. Anadromous brown trout in Norway typically do not migrate beyond a ca. 80 km range from their native river before returning to spawn [47]. Compared to Atlantic salmon, their home range is therefore more restricted to the fjords, where sea lice are concentrated. In some regions, low marine survival of brown trout has indeed been recorded during the last decades and ascribed to increasing exposure to salmon lice [43]. There is also evidence for genetic variation in susceptibility of brown trout populations to salmon lice infestation, further underpinning salmon lice as an important factor in marine mortality [48].

The present study focuses on the brown trout population of the River Guddal, located in the central region of the Norwegian Hardanger Fjord (Fig. 1). A trap facility encompassing two types of traps was installed in the early 2000’s after a reported decrease in the number of sea trout in the fjord [49]. It allows for a full monitoring of ascending adults and descending smolts, from which phenotypic and genetic data have been secured annually. In turn, this provides a unique long-term data set on the down- and up-stream migratory patterns of brown trout in this river, located in a major farming region where marine survival is known to be low [49].

Fig. 1
figure1

The trap facility is located in the river Guddal, in the central part of the Hardanger Fjord. The black dot indicates the location of the Guddal river mouth. Inset map shows location of the study site in Norway

We analyzed microsatellite DNA variation in all anadromous trout ascending the river from 2006 to 2016 and in all smolts emigrating from the river in 2017. We assessed temporal trends in the number of ascending spawners across years and used genetic data to estimate effective population sizes. We tested the hypothesis that census and effective population sizes had declined over the time span studied. Second, as migratory behaviour in general may involve genetic components [33, 35, 36] and specifically as some heritability has been shown in the timing of migration in other salmonid species [28, 50], we used genetically based parentage assignment and sibship reconstruction to test whether related individuals are more likely to migrate together than non-related individuals. Third, using parentage assignment of smolts from 2017 we assessed the degree to which anadromous and resident trout contributed to the smolt run, thus providing information on the possibility that resident trout can compensate for increased mortality in the marine environment.

Material and methods

Study site

The river Guddal (59° 58′ N, 6° 00′ E) is 13.5 km long and located in the central region of the Hardanger Fjord (Fig. 1), the second longest fjord in Norway and showing the highest concentration of Atlantic salmon farms [51]. Survival of wild anadromous trout in the system has been found to be low [43], presumably as a consequence of exposure to sea lice derived from salmon farms, as has recently been determined with genetic studies on lice [46, 52]. The river Guddal is partly fed by melting water running from the Folgefonna glacier, with relatively cold mean summer temperatures ranging from 8.5 to 11.9 °C between 2007 and 2016.

Sampling

Each year, all smolts descending the river for their first seaward migration are caught in a Wolf trap [53] that covers the whole river transect, located at about 100 m from the tidal zone. The trap is operated annually from March or early April, depending on the water discharge, until the end of the smolt run. In principle, the trap should sample all downstream migrating smolts, but in some years there are shorter periods ranging from a few hours to a couple of days, where water discharge prevents full sampling, so that the trap is estimated to catch 90% of the migrating smolts. Smolts caught in the trap are anaesthetised with benzocaine, measured, weighed, adipose fin-clipped, and tagged with a passive integrated transponder (PIT tags 12 × 2 mm) since 2007 [49, 54] before being released downstream. A specific test showed that among 600 smolts retained in a prolonged recovery, only 4 individuals (0.7%) died (Ø. Skaala, unpublished results). Hence, we do not expect that handling has severely affected the smolts, and any effects should be randomly distributed among individuals and families. In addition, an upstream trap is operated each year during the upstream migration to capture all anadromous trout ascending the river. Each fish is anaesthetised, fin-clipped, checked for PIT tag, measured and weighed before being released upstream to continue its migration. Tissue samples for both smolts and ascending fish were stored in 95% ethanol. For this study, we used data on 711 ascending fish collected between 2006 and 2016 (Table 1), representing all the adults entering the system in this period. Census numbers (but not genetic data) of ascending spawners between 2000 and 2005 were also available [43]. Furthermore, we analysed the entire smolt run of 2017, representing 965 fish (Table 1).

Table 1 Summary of brown trout trap sampling and genotyping, both for sea trout ascending the Guddal River in 2006–2016 and for the smolts emigrating to sea in 2017

Genotyping

DNA was isolated and extracted using Qiagen’s DNeasy® Blood & Tissue kit (Qiagen Inc.), following the manufacturer’s recommendations. DNA concentration was measured with a NanoDrop 1000 Spectrophotometer (Thermo Fisher Scientific Inc.). Dilutions were subsequently conducted using a Freedom EVOware® robot (Tecan Inc.), yielding ca. 16.6 ng/µL of DNA per sample.

All samples were genotyped at 21 microsatellite loci (Additional file 1: Appendix 1), divided into three multiplexes for the Polymerase Chain Reaction (PCR). The third multiplex included the sex-specific markers, Exon 2 and Exon 4 [55] from the sdY sequence, a male-specific-Y-chromosome gene that is highly conserved in salmonids [56]. Although not always 100% accurate, due to the occasional occurrence of autosomal pseudocopies of sdY, at least in brown trout´s closest relative, Atlantic salmon [57], this provides more accurate sexing than phenotypic sex as morphological dimorphism can be difficult to ascertain in the field [58]. PCR was conducted using a Verity 96 well thermal cycler and GeneAmp PCR system 9700 (Applied Biosystems). Conditions for the cycling reactions and primer mixes used for the 3 multiplexes are detailed in Additional file 1: Appendix 2. PCR products were then diluted at 1:15 and separated by capillary electrophoresis using a 3730 DNA Analyzer (Applied Biosystems). Alleles were scored using GeneMapper 5 with GeneScan 500 LIZ size standards (Applied Biosystems). A total of 30 adults and 94 smolts were identified as Atlantic salmon as they showed several allele sizes that do not correspond to alleles known from brown trout, which left a total of 1560 trout sampled (711 adult ascenders and 851 smolt, see Table 1).

Two microsatellite loci (Ssa85 and Ssa197) were represented twice in separate multiplexes in order to provide a control of the allele scoring (Additional file 1: Appendix 1). PCR was repeated for each sample until data were obtained for most loci. In the whole dataset, 3 individuals had missing data at 3 loci, 4 at 2 loci and 56 at 1 locus, out of the 21 loci genotyped. When several PCRs had been conducted for the same individuals, they were compared and the one with the most clear-cut signal was considered right in case of conflict. The locus BG935488 was excluded from the dataset due to difficulties with reliable scoring of alleles. As loci MHC-I and Sasa-TAP2A are closely linked to loci associated with immune responses and have previously been suggested to be under diversifying selection among trout populations from the Hardanger Fjord [59]; they were omitted from analyses of effective population size but were used for sibship reconstruction (see below).

Tissue samples were accidentally missing for 20 of the ascending trout (Table 1). However, the PIT tag number from a missing fish of 2016 had a match with one captured in 2015, allowing to find its genotype based on the previous sampling.

Individual identification

The “matches” algorithm implemented in GenAlEx 6.5 [60] was used to identify identical genotypes (allowing for one mismatching locus), corresponding to individuals that had ascended the river Guddal several times among or within years. This procedure was implemented as PIT tags can be lost between captures [54], particularly in females during spawning, thus compromising individual identification.

Genetic diversity and effective population size

Total number of alleles and allelic richness (Ar) per locus were computed with the R package diveRsity [61], using R 3.6.1 [62]. Observed (HO) and unbiased expected heterozygosity (HE), as well as inbreeding coefficient (FIS) were computed with GenAlEx 6.5. Tests for Hardy–Weinberg equilibrium were performed with Genepop 4.7 [63] using the following Markov chain parameters: 10,000 dememorisation steps, 1000 batches and 10,000 iterations per batch. False discovery rate (FDR) correction was applied to account for multiple testing [64], using the method by Y Benjamini and Y Hochberg [65] implemented in Myriads 1.1 [66].

Effective population sizes (Ne) for the ascending sea trout and the smolts were estimated with the linkage disequilibrium method implemented in NeEstimator 2.1 [67], using allele frequencies higher than 0.05. As several cohorts are represented among spawning individuals, we assume that the estimates measure Ne rather than Nb (the effective number of breeders in a single breeding event [68,69,70]). A smaller number of cohorts were expected to be represented in the sample of smolts and hence estimates could be shifted more towards estimating Nb.

Sibship and parentage assignment

COLONY 2.0.5.1 [71] was used to infer full and half sibships both in the 2006–2016 ascending trout and the smolt 2017 datasets. This software infers all possible relationships (siblings and parentage) of all individuals (all offspring and all candidate parents) simultaneously in a full likelihood framework. Analyses of the ascending trout dataset were conducted with no information on parental genotypes, assuming both male and female polygamy as well as possible inbreeding. The full-likelihood method was used at very high precision for the full likelihood calculation and medium run length together with the options sibship scaling, no updated allele frequencies and no sibship priors. In the analysis of the smolt dataset, the trout ascending in or after 2010 were used as candidate parents. We assumed that when fish were last caught before 2010, they were unlikely to be candidate parents, as smolts at these latitudes leave the river between age 2 and 7 years [72]. Similarly, it is unlikely that reproduction in 2016 would result in smolts sampled in 2017, but several trout ascending in 2016 had also ascended in previous years (see Additional file 1: Appendix 2) and could consequently be parents of smolts as a result of previous spawning events. We made use of this information for providing an empirical assessment of the quality of parentage assignment. We thus predicted that anadromous individuals that had ascended the river in 2016 and were identified as parents of smolts from 2017 should also have ascended the river in previous years.

Three COLONY runs with different random number seeds were used to check the reliability of the results. Only the relationships found in the 3 runs were kept in the inferred pedigree. Allelic dropout rate in the input file was estimated with the PopGenReport R package [73]. Genotyping error rate was set to 0.01 for all loci. Each time, full and half-sibships were inferred from the “best configuration” COLONY output files, which is more accurate than the results found by the pairwise analysis in the files full and halfsib dyads [74].

More than 50% of male parr in a population can be precocious [75,76,77]. As it would be infeasible to sample all resident trout and not least precocious male parr from the system, we inferred the parentage of resident trout indirectly. Hence, when none or only one of the two parents of a smolt was identified among the anadromous candidate parents, it was assumed that missing parents were part of the resident proportion of the population. This was a realistic assumption as all ascending anadromous spawners are assumed to have been caught in the trap and genotyped (only 8 adults caught in the trap were accidentally not sampled for genotyping between 2010 and 2016, Table 1). Moreover, given sex ratio differences between anadromous and resident trout, we further corroborated this assumption by expecting a higher number of supposedly resident males than resident females.

Sibship and timing of migration events

To test whether related individuals tend to migrate at the same time (both for ascending adults and descending smolts), Mantel tests were conducted between a distance matrix comprising the sibship previously found between pairs of individuals, coded as 1 for unrelated, 0.75 for halfsibs and 0.5 for fullsibs, and a matrix with the distance in days between their date of capture in the traps. The R package ecodist [78] was used for these analyses, and significance was assessed using 10,000 permutations together with 1000 bootstraps. The test was conducted separately for each year in the ascending spawners. Fish that did not have any siblings within the samples of the year or in the smolt sample were removed from the analyses, in order not to bias the results.

Results

Numbers of ascending anadromous spawners

The number of anadromous trout ascending the river Guddal showed pronounced variation among years between 2006 and 2016, with a mean of 64.9 ± 22.9 individuals (± SD), ranging from 25 in 2010 to 96 in 2015, but with no general tendency for a decline over this time period (linear regression: F(1,9) = 0.29, R2adj = -0.07, P = 0.61; see also Table 1). However, compared to the mean number of anadromous trout recorded in the trap in in 2000–2005 in the same river (100 ± 43 (mean ± SD), [43]) there were indications of a decline over a longer time span, although for the entire period of 2000–2016 the result remained non-significant (linear regression: F(1,15)=2.90, R2adj = 0.11, P = 0.11).

After genotyping, we identified 591 individuals among the 711 tissue samples from returning adult sea trout collected between 2006 and 2016. The matching genotypes showed that 82 individuals ascended the river more than once in the study period. Eight individuals were found returning to the river twice within the same year, and one returned three times in the same year (Table 1). Assignment tests using the simulation option in GeneClass2 [79] provided no evidence that the latter individuals were strayers from different populations, except for a single individual whose multilocus genotype was unlikely to occur based on the allele frequencies observed in the Guddal population (data not shown).

Individuals returned to the river on average 1.19 times, ranging from 1 to 6 times (recaptures are illustrated in Additional file 1: Appendix 3). Individuals ascending multiple times returned within a period of 2.7 ± 1.0 (± SD) years on average (Additional file 1: Appendix 3). During the study period, most of the sea trout ascended the river between mid-July and early September, but the spawning run had its earliest start in 2016 on the 14th of May and the latest end in 2010 on the 25th of November (Additional file 1: Appendix 4). The spawning run could be roughly divided into two peaks of migration for most of the years, one in summer and one in early autumn, but the pattern was not clear-cut (Additional file 1: Appendix 4).

A χ2 goodness of fit test revealed that the number of females exceeded number of males among anadromous trout across all years, with sex ratios (females:males) ranging from 1.04 in 2013 to 2.13 in 2009 (Table 2, χ2 = 28.2, df = 10, P = 0.002).

Table 2 Summary of the population genetic statistics, both for anadromous trout ascending the river Guddal in 2006–2016 and for smolts sampled in 2017

Genetic variation and effective population size estimated from ascending spawners

A total of three deviations from Hardy–Weinberg equilibrium were found in the microsatellite DNA dataset for annually ascending trout after Benjamini and Hochberg [80] False Discovery Rate correction (Table 2). As they concerned three different loci that were not deviating in other annual samples, they were all retained in the analyses. Mean unbiased expected heterozygosity (HE) was high and stable across years, ranging from 0.74 to 0.76, and FIS was not significantly different from 0, except in 2012 and 2013 where it nevertheless remained low (Table 2).

Estimates of mean effective population size (Ne) showed considerable variation, ranging between 83 in 2014 and 3646 in 2010 (Table 2). The highest Ne estimate coincided with the lowest sample size (N = 25) and showed a very wide 95% confidence interval (94–∞), indicating that the point estimate is not informative. Considering only the years with sample sizes > 50 (a total of 8 years), the Ne point estimates ranged from 83 to 548, with a mean of 223.1. Similar to the number of spawners per year, there was no evidence for a temporal decline of Ne.

Genetic variation and effective population size estimated from 2017 smolts

Mean unbiased expected heterozygosity (HE) was similar to that observed in anadromous spawners (Table 2). In contrast, however, 17 deviations from Hardy–Weinberg equilibrium were found among the 18 loci analysed (Table 2), likely reflecting the large number of full- and half-sibs (see below) and essentially violating the criterion for Hardy–Weinberg equilibrium of infinite population size, combined with very high statistical power due to the sample size of N = 851. Sex-ratio in the total sample of downstream migrating smolts was more biased towards females than for the anadromous spawners (2.15 against a mean of 1.48 in the spawners, Table 2). The point estimate of effective population size (Ne) was 51, much lower than estimates based on anadromous spawners (Table 2).

Relatedness and timing of migration events

Among the 851 smolts migrating to the sea in 2017, the consensus pedigree achieved after three different COLONY runs identified 3198 full and 11,065 halfsib dyads. This represents 317 different fullsib families, with a mean of 2.7 offspring per family (ranging between 1 and 42). The Mantel test between distance matrices composed of distance in Julian date of seaward migration and distance in inferred sibship relation yielded a significant positive correlation for the 815 smolts which had siblings in the run (r = 0.027, P = 0.0002, Table 3).

Table 3 Mantel test results between sibship and distance in week of ascendance for ascending trouts between 2006–2016 and between Julian date of seaward migration for 2017 smolt run

For the samples of anadromous spawners, the numbers of individuals with half- or full-sibs in the spawning run of the same year ranged from 0 in 2010 to 56 in 2015 (Table 3). Mantel tests for association of sibship and migration timing conducted for ascending anadromous trout year by year used the week of ascendance in the matrix, as the time window for the ascending migration is larger than for descending migration of smolts. The results were significant for the adults ascending in 4 years (2006: r = 0.112, P = 0.0040; 2007: r = 0.225, P = 0.0425; 2008: r = 0.092, P = 0.0400, and 2015: r = 0.055, P = 0.0405; Table 3), but not for individuals ascending in the other years.

Resident trout contribution to anadromous trout run

Parents of the 2017 smolt run were found among the anadromous trout ascending the river Guddal between 2011 and 2016. Among 140 candidate anadromous fathers and 175 candidate anadromous mothers, 33 males and 52 females contributed to the total smolt sample. Among these individuals, four males and 12 females ascended the river in 2016, but as stated previously, the 2017 smolt run would be unlikely to include offspring from the previous year. In accordance with this, all but one of these individuals had in fact been recorded as ascending the river and presumably spawning in previous years, thus lending indirect empirical support to the robustness of parentage assignment. The putative parent recorded in 2016 but not found to ascend the river in previous years was a male. In addition to the possibility of incorrect parentage assignment, it is possible that the male reproduced as mature male parr and smolted and undertook migration to the sea afterwards. It is also possible that this male represents some of the individuals that were accidentally not genotyped (see Table 1).

Based on the assumption that a parent not represented by anadromous spawners caught in the trap must have been a resident individual, the freshwater resident part of the population was inferred to contribute to 70% of the 2017 smolt run (Fig. 2A). Identified anadromous females contributed to 72% of the smolt sample. The opposite pattern was observed for males, with putatively resident males (including precocious male parr) contributing to 66% of the smolt sample, which was significantly more than anadromous males (binomial tests, N = 851, P < 0.001, Fig. 2B). These patterns are in accordance with the skewed sex ratio observed among anadromous spawners (Table 2). A total of 46% of the smolts had one putatively resident parent while the other was a migrant, 24% had two resident parents and 30% had two anadromous parents (Fig. 2A). Among these matings involving two different life history strategy types, 91% took place between a migrant female and a putatively resident male (Fig. 2A).

Fig. 2
figure2

The parental origin of smolts in the 2017 run based on parentage assignment using COLONY 2.0.5.1 [71]. A Percentage of smolts according to their parental origin. B Contribution of sexes of anadromous and resident trout to the smolt run. Error bars represent the 95% confidence interval of the percentage. Significance of the binomial tests to compare contribution of anadromous and resident individuals of the same sex, ***: P < 0.001

Discussion

Our study represents a unique long term monitoring effort of a salmonid population in an environment experiencing significant anthropogenic pressure. In addition to providing insights into demographic trends, the results raised the interesting possibility of concordance between migration timing and kinship in smolts and adult spawners. Finally, and most importantly, parentage assignment of an entire cohort of migrating smolts showed that resident parents contribute substantially to the anadromous run. We discuss these findings and their conservation implications in more detail below.

Temporal trends of census and effective population size

The number of anadromous trout ascending the Guddal River varied considerably among years, but did not show a general tendency for decline over the study period, although numbers of spawners were lower in 2006–2016 than in 2000–2006. Salmon farming in the Hardangerfjord accellerated during the 1990s, so it is expected that the most severe population declines occurred during that period and have subsequently stabilized, as indicated by our results.

Similar to census sizes of anadromous spawning run, effective population size estimates also showed no general trends towards declines during the study period. This comes with the caveat that estimates of Ne based on linkage disequilibrium are sensitive to low sample sizes [81], which for some years/samples resulted in low precision of estimates. It may seem surprising that the Ne estimate based on smolts (51) was much lower than any of the estimates based on anadromous spawners. However, as noted previously this estimate is not completely comparable to estimates from adult anadromous spawners that represent different numbers of sea-winters before they return to the river; the smolt run in a given year is expected to represent fewer cohorts compared to adult spawners and would therefore represent something in between Ne and Nb (the effective number of breeders in a single breeding event [68, 69]).

The Ne estimates of a few hundreds are comparable to most other estimates found in anadromous brown trout populations using temporal or LD-based methods [59, 82,83,84], but higher than most estimates from strictly resident populations [85,86,87,88]. Whereas Ne in the Guddal population is lower than the 500 or even 1000 assumed to be required for maintaining evolutionary potential [89,90,91,92], it should be noted that it is part of a larger system in the Hardanger Fjord where gene flow occurs among populations [59]. In general, anadromous brown trout populations have been found to exhibit a hierarchical genetic structure shaped by both geographical distance between populations and environmental parameters, with low genetic differentiation among local populations resulting from gene flow [93]. Evolutionary potential should thus be considered across several neighbouring populations, where total Ne is expected to be higher [94].

Concordance of migration timing of related individuals

Association of kin along with the possibility of kin selection has been studied intensively in salmonids [29, 95,96,97,98,99]. As sibs during the earliest life stages are situated in the same spawning redds, spatial association of kin would be expected to occur immediately after hatching, whereas subsequent drift and dispersal would lead to decreased kin association over time, unless active kin recognition and association takes place [99]. The significant association betwen migration timing and relatedness as observed in smolts in the present study can hardly be ascribed to reminiscent patterns of association of kin several years back in time in their spawning redds, but could reflect: (1) active aggregation of kin, (2) genetic components in the timing of smolt migration, and/or (3) similarity in size of sibs and thereby propensity for migrating at the same time, simply because sibs hatched and emerged from the same redds at the same time. Whereas the study does not allow for distinguishing between these possibilities, we note that (2) and (3) are the most parsimonious explanations and also indirectly supported by empirical results [47, 100, 101], including data from Atlantic salmon demonstrating clear genetic components in migration timing [96, 102].

Interestingly, our results are at odds with those from a different study on migration timing and kinship in Atlantic salmon, which found no significant association between kinship and schooling and migration timing in smolts [103]. Part of the reason for the discrepancy of results could lie in different experimental set-ups. The study by Fernandes et al. [103] was based on experimental full-sib families stocked into a natural environment at the same point in time, whereas our study encompassed the total smolt run composed of families naturally spawned and hatched over an extended period of time. This would leave more statistical power in our study for detecting association between kinship and migration timing resulting from similar hatching time and size of sibs, without necessarily involving active kin aggregation or genetic components in migration timing.

We also found some support for association between kin and timing of upstream migration among anadromous trout, although significant associations were observed in only 4 out of 11 years. Few studies of possible kin-biased aggregation of adult salmonid fishes have been conducted, undoubtedly due to challenges with sampling. However, one study found kin associations to occur at the adult stage outside the spawning period in brook char (Salvelinus fontinalis) inhabiting a large freshwater lake [29]. Nevertheless, kinship analysis of spawners in a tributary to the same lake provided no evidence for association of kin [104]. As noted by the authors, this could be an effect of accummulated mortality over time, leaving few surviving related individuals at the time of spawning. This could also be the case in our present study, where numbers of anadromous spawners per year were overall low.

In total, there was evidence for association between kinship and migration time in smolts, and also evidence, albeit less consistent, for such an association in spawners returning to the river. The association found in smolts raises the possibility that episodes of increased marine mortality, e.g., due to salmon lice exposure or fluctuating aggregations of predators, could potentially increase variance in mortality among families, which could again lead to higher variance in reproductive success among families and lower effective population size. Sibs were found among anadromous spawners in all years except for 2010 (with only 25 ascending anadromous spawners), but whether this reflects a disproportionally high variance in mortality among families compared to undisturbed conditions cannot be assessed. This would require comparable data from the system before major environmental disturbance of the Hardangerfjord system took place.

To what extent does the resident stock of Guddal brown trout population contribute to the sea run?

In systems like the Hardanger Fjord with increased marine mortality [43,44,45, 49], it is important to assess to which extent the resident part of the population can compensate for recruitment in the case of a reduced number of anadromous spawners. Moreover, it is important to consider to what extent this will drive changes in anadromy.

We did not genotype candidate parents among the resident trout and made the assumption that non-genotyped parents corresponded to resident individuals. This is an important limitation of the study and for instance precludes distinguishing between parents that are mature male parr and adult resident trout. However, the trap in which ascending spawners were sampled is of a construction that makes it unlikely that individuals can escape further upstream without being registered. In very few instances tissue samples were by mistake not taken, but this is unlikely to account for all the parental genotypes not represented among the ascending anadromous spawners. Moreover, our findings of parentage are in accordance with expectations given the skewed sex-ratio observed among anadromous and resident spawners [5, 26, 37, 39], providing further confidence in our results. Hence, only 4% of all smolts in 2017 had an anadromous father and putatively resident mother, whereas 42% had a putatively resident father and anadromous mother. In total, 70% of the smolts had either one or two inferred resident parents, with the latter category accounting for 24% of all smolts.

In the case of resident males, it is likely that many of them are in fact precocious male parr, as studies of both brown trout and other salmonids have shown that they can be both numerous and have significant reproductive success [6, 38, 39, 105]. On the other hand, the contribution to the 2017 smolt run was higher for anadromous than putatively resident females (72% versus 28%), which could reflect a higher number of anadromous relative to resident females and/or the fact that it is more advantageous as a female to migrate to sea in order to maximize body weight and egg production and thereby reproductive success [31, 106]. In the context of elevated mortality rates at sea, these results demonstrate that resident trout may indeed have some buffering effects towards a decline of anadromous spawners and that a sizeable proportion of smolts in fact have two resident parents. Given the high contribution of anadromous females to the smolt run it is nevertheless also evident that strong declines of the anadromous portion of the population would likely have significant negative demographic consequences for the total population. It should also be stressed that high marine mortality would lead to reduced gene flow among populations, hence reducing overall effective metapopulation size (as in for instance the entire Hardanger Fjord system) and potentially leading to inbreeding and loss of variation in individual populations [94].

Could long-term elevated marine mortality select against anadromy and ultimately remove it from the population? If we assume that anadromy is a quantitative trait with an environmentally-cued threshold [10], then in this case high levels of genetic variation could be preserved even under directional selection acting against it [107]. This way, even if migration costs are increasing, the propensity to migrate within a population may persist. This finds empirical support in studies of brown trout [84] and other salmonid species that have been landlocked for centuries (e.g., by dams) such as brook charr [108], rainbow trout Oncorhynchus mykiss [109] and bull trout Salvelinus confluentus [110], but where migratory behaviour is retained. Nevertheless, recent studies have demonstrated major quantitative trait loci for important life history and migratory traits in e.g., Chinook Oncorhynchus tshawytscha [12] and Atlantic salmon [13]. If such genetically based variation also exists in brown trout, then increased marine mortality could exert strong selection and lead to genetic and phenotypic changes in populations, even if anadromy per se is retained.

Conclusions

The unique long-term monitoring of ascending anadromous spawners in the Guddal River allowed us to track both census and effective population size over an extended time period coinciding with adverse anthropogenic conditions in the sea. We did not observe general trends of declines during this period, suggesting that the population had stabilized after initial declines prior to the study. We found a significant association between kinship of smolts and their timing of emigration from the river, which raises the possibility that periodically adverse conditions in the sea could disproportionally affect some families and potentially lead to decreased effective population size. A similar association between kinship and migration timing was also indicated in ascending anadromous spawners but was not consistent across years. It is possible that accummulated mortality until this life stage would decrease the number of surviving sibs and thereby weaken signals of an otherwise genuine association. Finally, using parentage assignment of the total smolt run within a year, we estimated that 70% of smolts had at least one resident parent, and in 24% of cases both parents were inferred to be resident. Hence, the resident proportion of the population played a major role in recruitment of anadromous trout and would be expected to provide some buffer against elevated marine mortality. Nevertheless, as the majority of smolts had an anadromous mother, it is also envisaged that elevated marine mortality would have important negative consequences for production of the population as a whole and could also lead to altered selection pressure for important life history traits. In total, the study thus provides important new information about recruitment and dynamics of populations showing partial migration, and how this may interact with anthropogenic environmental disturbance.

Availability of data and materials

The data that support the findings of this study have been deposited at the Institute of Marine Research electronic archive at: https://hdl.handle.net/11250/2740709. It consists of three text files: a genepop-format file of all the individuals genotyped at 20 microsatellite loci, and two files containing the phenotypic information for adults and smolts, respectively.

Abbreviations

Ar :

Allelic richness

HO :

Observed heterozygosity

HE :

Expected heterozygosity

FIS :

Fixation index in subpopulations (“inbreeding coefficient”)

FDR:

False discovery rate

Ne :

Effective population size

References

  1. 1.

    Gross MR. Alternative reproductive strategies and tactics: diversity within sexes. Trends Ecol Evol. 1996;11(2):92–8.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  2. 2.

    Dingle H. Migration: the biology of life on the move. Oxford: Oxford University Press; 1996.

    Google Scholar 

  3. 3.

    Chapman BB, Brönmark C, Nilsson JÅ, Hansson LA. The ecology and evolution of partial migration. Oikos. 2011;120(12):1764–75.

    Article  Google Scholar 

  4. 4.

    McDowall R. Diadromy: origins and definitions of terminology. Copeia. 1992;1992(1):248–51.

    Article  Google Scholar 

  5. 5.

    Jonsson B, Jonsson N. Partial migration: niche shift versus sexual maturation in fishes. Rev Fish Biol Fisher. 1993;3(4):348–65.

    Article  Google Scholar 

  6. 6.

    Hendry AP, Bohlin T, Jonsson B, Berg OK. To sea or not to sea? Anadromy versus non-anadromy in salmonids. In: Hendry AP, Stearns SC, editors. Evolution illuminated: salmon and their relatives. New York: Oxford University Press; 2004. p. 92–125.

    Google Scholar 

  7. 7.

    Gross MR, Coleman RM, McDowall RM. Aquatic productivity and the evolution of diadromous fish migration. Science. 1988;239(4845):1291–3.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  8. 8.

    Fleming IA, Gross MR. Breeding competition in a Pacific salmon (coho: Oncorhynchus kisutch): measures of natural and sexual selection. Evolution. 1994;48(3):637–57.

    PubMed  PubMed Central  Google Scholar 

  9. 9.

    Halttunen E, Gjelland KØ, Hamel S, Serra-Llinares RM, Nilsen R, Arechavala-Lopez P, Skarðhamar J, Johnsen IA, Asplin L, Karlsen Ø. Sea trout adapt their migratory behaviour in response to high salmon lice concentrations. J Fish Dis. 2017;41:953–67.

    PubMed  Article  PubMed Central  Google Scholar 

  10. 10.

    Ferguson A. Genetics of sea trout, with particular reference to Britain and Ireland. In: Harris G, Milner N, editors. Sea trout: biology, conservation and management. Oxford: Blackwell publishing; 2006. p. 157–82.

    Google Scholar 

  11. 11.

    Ferguson A, Reed TE, McGinnity P, Prodöhl PA. Anadromy in brown trout (Salmo trutta): a review of the relative roles of genes and environmental factors and the implications for management and conservation. In: Harris GS (Ed.). Sea Trout: from Science to Management. Self-published; 2016.

  12. 12.

    Thompson TQ, Bellinger MR, O’Rourke SM, Prince DJ, Stevenson AE, Rodrigues AT, Sloat MR, Speller CF, Yang DY, Butler VL, et al. Anthropogenic habitat alteration leads to rapid loss of adaptive variation and restoration potential in wild salmon populations. P Natl Acad Sci USA. 2019;116(1):177–86.

    CAS  Article  Google Scholar 

  13. 13.

    Barson NJ, Aykanat T, Hindar K, Baranski M, Bolstad GH, Fiske P, Jacq C, Jensen AJ, Johnston SE, Karlsson S, et al. Sex-dependent dominance at a single locus maintains variation in age at maturity in salmon. Nature. 2015;528:405–8.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  14. 14.

    Czorlich Y, Aykanat T, Erkinaro J, Orell P, Primmer CR. Rapid sex-specific evolution of age at maturity is shaped by genetic architecture in Atlantic salmon. Nat Ecol Evol. 2018;2(11):1800–7.

    PubMed  PubMed Central  Article  Google Scholar 

  15. 15.

    Baerwald MR, Meek MH, Stephens MR, Nagarajan RP, Goodbla AM, Tomalty KMH, Thorgaard GH, May B, Nichols KM. Migration-related phenotypic divergence is associated with epigenetic modifications in rainbow trout. Mol Ecol. 2016;25(8):1785–800.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  16. 16.

    Wilcove DS, Wikelski M. Going, going, gone: is animal migration disappearing? PLoS Biol. 2008;6(7):1361–4.

    CAS  Article  Google Scholar 

  17. 17.

    Wysujack K, Greenberg L, Bergman E, Olsson I. The role of the environment in partial migration: food availability affects the adoption of a migratory tactic in brown trout Salmo trutta. Ecol Freshw Fish. 2009;18(1):52–9.

    Article  Google Scholar 

  18. 18.

    Shaw AK. Drivers of animal migration and implications in changing environments. Evol Ecol. 2016;30(6):991–1007.

    Article  Google Scholar 

  19. 19.

    Webster MS, Marra PP, Haig SM, Bensch S, Holmes RT. Links between worlds: unraveling migratory connectivity. Trends Ecol Evol. 2002;17(2):76–83.

    Article  Google Scholar 

  20. 20.

    Rand PS, Goslin M, Gross MR, Irvine JR, Augerot X, McHugh PA, Bugaev VF. Global assessment of extinction risk to populations of sockeye salmon Oncorhynchus nerka. PLoS ONE. 2012;7(4):e34065.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  21. 21.

    Limburg KE, Waldman JR. Dramatic declines in North Atlantic diadromous fishes. Bioscience. 2009;59(11):955–65.

    Article  Google Scholar 

  22. 22.

    Shelton AO, Sullaway GH, Ward EJ, Feist BE, Somers KA, Tuttle VJ, Watson JT, Satterthwaite WH. Redistribution of salmon populations in the northeast Pacific ocean in response to climate. Fish and Fisheries 2021, in press.

  23. 23.

    Krkošek M, Lewis MA, Volpe JP. Transmission dynamics of parasitic sea lice from farm to wild salmon. Proc Royal Soc London B: Biol Sci. 2005;272(1564):689–96.

    Google Scholar 

  24. 24.

    Torrissen O, Jones S, Asche F, Guttormsen A, Skilbrei OT, Nilsen F, Horsberg TE, Jackson D. Salmon lice—impact on wild salmonids and salmon aquaculture. J Fish Dis. 2013;36(3):171–94.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  25. 25.

    Taranger GL, Karlsen O, Bannister RJ, Glover KA, Husa V, Karlsbakk E, Kvamme BO, Boxaspen KK, Bjorn PA, Finstad B, et al. Risk assessment of the environmental impact of Norwegian Atlantic salmon farming. ICES J Mar Sci. 2015;72(3):997–1021.

    Article  Google Scholar 

  26. 26.

    Jonsson B. Life history patterns of freshwater resident and sea-run migrant brown trout in Norway. T Am Fish Soc. 1985;114(2):182–94.

    Article  Google Scholar 

  27. 27.

    Palstra FP, Ruzzante DE. Genetic estimates of contemporary effective population size: what can they tell us about the importance of genetic stochasticity for wild population persistence? Mol Ecol. 2008;17(15):3428–47.

    PubMed  Article  PubMed Central  Google Scholar 

  28. 28.

    Bentzen P, Olsen J, McLean J, Seamons T, Quinn T. Kinship analysis of Pacific salmon: insights into mating, homing, and timing of reproduction. J Hered. 2001;92(2):127–36.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  29. 29.

    Fraser DJ, Duchesne P, Bernatchez L. Migratory charr schools exhibit population and kin associations beyond juvenile stages. Mol Ecol. 2005;14(10):3133–46.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  30. 30.

    Killen SS, Marras S, Nadler L, Domenici P. The role of physiological traits in assortment among and within fish shoals. Philos Trans R Soc B. 2017;372(1727):20160233.

    Article  Google Scholar 

  31. 31.

    Elliott JM. Quantitative ecology and the brown trout. Oxford: Oxford University Press; 1994.

    Google Scholar 

  32. 32.

    Archer LC, Hutton SA, Harman L, McCormick SD, O’Grady MN, Kerry JP, Poole WR, Gargan P, McGinnity P, Reed TE. Food and temperature stressors have opposing effects in determining flexible migration decisions in brown trout (Salmo trutta). Global Change Biol. 2020;26(5):2878–96.

    Article  Google Scholar 

  33. 33.

    Archer LC, Hutton SA, Harman L, O'Grady MN, Kerry JP, Poole WR, Gargans P, McGinnity P, Reed TE. The interplay between extrinsic and intrinsic factors in determining migration decisions in brown trout (Salmo trutta): an experimental study. Front Ecol Evol 2019; 7.

  34. 34.

    Nevoux M, Finstad B, Davidsen JG, Finlay R, Josset Q, Poole R, Hojesjo J, Aarestrup K, Persson L, Tolyanen O, et al. Environmental influences of life history strategies in partial anadromous brown trout (Salmo trutta, Salmonidae). Fish Fish. 2019;20(6):1051–82.

    Article  Google Scholar 

  35. 35.

    Lemopoulos A, Uusi-Heikkila S, Huusko A, Vasemagi A, Vainikka A. Comparison of Migratory and resident populations of brown trout reveals candidate genes for migration tendency. Genome Biol Evol. 2018;10(6):1493–503.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  36. 36.

    Lemopoulos A, Uusi-Heikkila S, Hyvarinen P, Alioravainen N, Prokkola JM, Elvidge CK, Vasemagi A, Vainikka A. Association mapping based on a common-garden migration experiment reveals candidate genes for migration tendency in brown trout. G3-Genes Genom Genet. 2019;9(9):2887–96.

    CAS  Google Scholar 

  37. 37.

    Bekkevold D, Hansen MM, Mensberg KLD. Genetic detection of sex-specific dispersal in historical and contemporary populations of anadromous brown trout Salmo trutta. Mol Ecol. 2004;13(6):1707–12.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  38. 38.

    Klemetsen A, Amundsen PA, Dempson J, Jonsson B, Jonsson N, O’connell M, Mortensen E. Atlantic salmon Salmo salar L., brown trout Salmo trutta L. and Arctic charr Salvelinus alpinus (L.): a review of aspects of their life histories. Ecol Freshw Fish. 2003;12(1):1–59.

    Article  Google Scholar 

  39. 39.

    Cucherousset J, Ombredane D, Charles K, Marchand F, Baglinière J-L. A continuum of life history tactics in a brown trout (Salmo trutta) population. Can J Fish Aquat Sci. 2005;62(7):1600–10.

    Article  Google Scholar 

  40. 40.

    Glover KA, Solberg MF, McGinnity P, Hindar K, Verspoor E, Coulson MW, Hansen MM, Araki H, Skaala O, Svasand T. Half a century of genetic interaction between farmed and wild Atlantic salmon: status of knowledge and unanswered questions. Fish Fish. 2017;18(5):890–927.

    Article  Google Scholar 

  41. 41.

    Forseth T, Barlaup BT, Finstad B, Fiske P, Gjoaester H, Falkegard M, Hindar A, Mo TA, Rikardsen AH, Thorstad EB, et al. The major threats to Atlantic salmon in Norway. ICES J Mar Sci. 2017;74(6):1496–513.

    Article  Google Scholar 

  42. 42.

    Thorstad EB, Finstad B. Impacts of salmon lice emanating from salmon farms on wild Atlantic salmon and sea trout. NINA Report. 2018;1449:1–22.

    Google Scholar 

  43. 43.

    Skaala Ø, Kålås S, Borgstrøm R. Evidence of salmon lice-induced mortality of anadromous brown trout (Salmo trutta) in the Hardangerfjord, Norway. Marine Biol Res. 2014;10(3):279–88.

    Article  Google Scholar 

  44. 44.

    Johnsen IA, Stien LH, Sandvik AD, Asplin L, Oppedal F. Optimal estimation of lice release from aquaculture based on ambient temperatures. Aquacult Env Interac. 2020;12:179–91.

    Article  Google Scholar 

  45. 45.

    Harvey AC, Quintela M, Glover KA, Karlsen O, Nilsen R, Skaala O, Saegrov H, Kalas S, Knutar S, Wennevik V. Inferring Atlantic salmon post-smolt migration patterns using genetic assignment. Roy Soc Open Sci. 2019;6(10):190426.

    CAS  Article  Google Scholar 

  46. 46.

    Fjørtoft HB, Nilsen F, Besnier F, Stene A, Bjorn PA, Tveten AK, Aspehaug VT, Finstad B, Glover KA. Salmon lice sampled from wild Atlantic salmon and sea trout throughout Norway display high frequencies of the genotype associated with pyrethroid resistance. Aquacult Env Interac. 2019;11:459–68.

    Article  Google Scholar 

  47. 47.

    Thorstad EB, Todd CD, Uglem I, Bjørn PA, Gargan PG, Vollset KW, Halttunen E, Kålås S, Berg M, Finstad B. Marine life of the sea trout. Mar Biol. 2016;163(3):47.

    Article  Google Scholar 

  48. 48.

    Glover KA, Skaala O, Nilsen F, Olsen R, Teale AJ, Taggart JB. Differing susceptibility of anadromous brown trout (Salmo trutta L.) populations to salmon louse (Lepeophtheirus salmonis (Kroyer, 1837)) infection. ICES J Mar Sci. 2003;60(5):1139–48.

    Article  Google Scholar 

  49. 49.

    Skaala Ø. A summary of 20 years (1998–2017) of scientific work on genetics and survival in anadromous brown trout (Salmo trutta L) and Atlantic salmon (S. salar L) in the river Guddalselva western Norway. Fisken Og Havet. 2017;5:46.

    Google Scholar 

  50. 50.

    Dodson JJ, Aubin-Horth N, Thériault V, Páez DJ. The evolutionary ecology of alternative migratory tactics in salmonid fishes. Biol Rev. 2013;88(3):602–25.

    PubMed  Article  PubMed Central  Google Scholar 

  51. 51.

    Halttunen E, Gjelland KØ, Hamel S, Serra-Llinares RM, Nilsen R, Arechavala-Lopez P, Skarðhamar J, Johnsen IA, Asplin L, Karlsen Ø. Sea trout adapt their migratory behaviour in response to high salmon lice concentrations. J Fish Dis. 2017;00:1–15.

    Google Scholar 

  52. 52.

    Fjørtoft HB, Besnier F, Stene A, Nilsen F, Bjørn PA, Tveten A-K, Finstad B, Aspehaug V, Glover KA. The Phe362Tyr mutation conveying resistance to organophosphates occurs in high frequencies in salmon lice collected from wild salmon and trout. Sci Rep-Uk. 2017;7(1):14258.

    Article  CAS  Google Scholar 

  53. 53.

    Wolf P. A trap for the capture of fish and other organisms moving downstream. T Am Fish Soc. 1951;80(1):41–5.

    Article  Google Scholar 

  54. 54.

    Gibbons JW, Andrews KM. PIT tagging: simple technology at its best. Bioscience. 2004;54(5):447–54.

    Article  Google Scholar 

  55. 55.

    Eisbrenner WD, Botwright N, Cook M, Davidson EA, Dominik S, Elliott NG, Henshall J, Jones SL, Kube PD, Lubieniecki KP. Evidence for multiple sex-determining loci in Tasmanian Atlantic salmon (Salmo salar). Heredity. 2014;113(1):86.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  56. 56.

    Yano A, Nicol B, Jouanno E, Quillet E, Fostier A, Guyomard R, Guiguen Y. The sexually dimorphic on the Y-chromosome gene (sdY) is a conserved male-specific Y-chromosome sequence in many salmonids. Evol Appl. 2013;6(3):486–96.

    CAS  PubMed  Article  Google Scholar 

  57. 57.

    Ayllon F, Solberg MF, Besnier F, Fjelldal PG, Hansen TJ, Wargelius A, Edvardsen RB, Glover KA. Autosomal sdY pseudogenes explain discordances between phenotypic sex and DNA marker for sex identification in Atlantic salmon. Front Genet. 2020;11:544207.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  58. 58.

    Quéméré E, Perrier C, Besnard A-L, Evanno G, Baglinière J-L, Guiguen Y, Launey S. An improved PCR-based method for faster sex determination in brown trout (Salmo trutta) and Atlantic salmon (Salmo salar). Conserv Genet Resour. 2014;6(4):825–7.

    Article  Google Scholar 

  59. 59.

    Hansen MM, Skaala O, Jensen LF, Bekkevold D, Mensberg KLD. Gene flow, effective population size and selection at major histocompatibility complex genes: brown trout in the Hardanger Fjord, Norway. Mol Ecol. 2007;16(7):1413–25.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  60. 60.

    Peakall R, Smouse P. GenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and researchd-an update. Bioinformatics. 2012;28:2537–9.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  61. 61.

    Keenan K, McGinnity P, Cross TF, Crozier WW, Prodöhl PA. diveRsity: an R package for the estimation and exploration of population genetics parameters and their associated errors. Methods Ecol Evol. 2013;4(8):782–8.

    Article  Google Scholar 

  62. 62.

    RCoreTeam: R: A language and environment for statistical computing. R foundation for Statistical Computing, Vienna. 2018.

  63. 63.

    Rousset F. GENEPOP’007: a complete re-implementation of the genepop software for Windows and Linux. Mol Ecol Resour. 2008;8(1):103–6.

    PubMed  Article  PubMed Central  Google Scholar 

  64. 64.

    Waples RS. Testing for Hardy-Weinberg proportions: have we lost the plot? J Hered. 2014;106(1):1–19.

    PubMed  Article  PubMed Central  Google Scholar 

  65. 65.

    Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Roy Stat Soc: Ser B (Methodol). 1995;57:289–300.

    Google Scholar 

  66. 66.

    Carvajal-Rodriguez A. Myriads: P-value-based multiple testing correction. Bioinformatics. 2017;1:3.

    Google Scholar 

  67. 67.

    Do C, Waples RS, Peel D, Macbeth G, Tillett BJ, Ovenden JR. NeEstimator v2: re-implementation of software for the estimation of contemporary effective population size (Ne) from genetic data. Mol Ecol Resour. 2014;14(1):209–14.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  68. 68.

    Ferchaud AL, Perrier C, April J, Hernandez C, Dionne M, Bernatchez L. Making sense of the relationships between Ne, Nb and Nc towards defining conservation thresholds in Atlantic salmon (Salmo salar). Heredity (Edinb). 2016;117(4):268–78.

    CAS  Article  Google Scholar 

  69. 69.

    Waples RS. A generalized approach for estimating effective population size from temporal changes in allele frequency. Genetics. 1989;121(2):379–91.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  70. 70.

    Waples RS. Genetic estimates of contemporary effective population size: to what time periods do the estimates apply? Mol Ecol. 2005;14(11):3335–52.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  71. 71.

    Jones OR, Wang J. COLONY: a program for parentage and sibship inference from multilocus genotype data. Mol Ecol Resour. 2010;10(3):551–5.

    PubMed  Article  PubMed Central  Google Scholar 

  72. 72.

    Økland F, Jonsson B, Jensen A, Hansen L. Is there a threshold size regulating seaward migration of brown trout and Atlantic salmon? J Fish Biol. 1993;42(4):541–50.

    Article  Google Scholar 

  73. 73.

    Adamack AT, Gruber B. PopGenReport: simplifying basic population genetic analyses in R. Methods Ecol Evol. 2014;5(4):384–7.

    Article  Google Scholar 

  74. 74.

    Wang JL. Estimating genotyping errors from genotype and reconstructed pedigree data. Methods Ecol Evol. 2018;9(1):109–20.

    Article  Google Scholar 

  75. 75.

    L’Abee-Lund JH, Jensen AJ, Johnsen BO. Interpopulation variation in male parr maturation of anadromous brown trout (Salmo trutta) in Norway. Can J Zool. 1990;68(9):1983–7.

    Article  Google Scholar 

  76. 76.

    Dellefors C, Faremo U. Early sexual maturation in males of wild sea trout, Salmo trutta L, inhibits smoltification. J Fish Biol. 1988;33(5):741–9.

    Article  Google Scholar 

  77. 77.

    Birnie-Gauvin K, Thorstad EB, Aarestrup K. Overlooked aspects of the Salmo salar and Salmo trutta lifecycles. Rev Fish Biol Fisher. 2019;29(4):749–66.

    Article  Google Scholar 

  78. 78.

    Dray S, Dufour A-B. The ade4 package: implementing the duality diagram for ecologists. J Stat Softw. 2007;22(4):1–20.

    Article  Google Scholar 

  79. 79.

    Piry S, Alapetite A, Cornuet J-M, Paetkau D, Baudouin L, Estoup A. GENECLASS2: a software for genetic assignment and first-generation migrant detection. J Hered. 2004;95(6):536–9.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  80. 80.

    Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Royal Stat Soc Series B (Methodological). 1995;57:289–300.

    Article  Google Scholar 

  81. 81.

    Wang J. Estimation of effective population sizes from data on genetic markers. Phil Trans Royal Soc B: Biol Sci. 2005;360(1459):1395–409.

    CAS  Article  Google Scholar 

  82. 82.

    Hansen MM, Ruzzante DE, Nielsen EE, Bekkevold D, Mensberg KLD. Long-term effective population sizes, temporal stability of genetic composition and potential for local adaptation in anadromous brown trout (Salmo trutta) populations. Mol Ecol. 2002;11(12):2523–35.

    PubMed  Article  PubMed Central  Google Scholar 

  83. 83.

    Hansen MM, Fraser DJ, Meier K, Mensberg KLD. Sixty years of anthropogenic pressure: a spatio-temporal genetic analysis of brown trout populations subject to stocking and population declines. Mol Ecol. 2009;18(12):2549–62.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  84. 84.

    Hansen MM, Limborg MT, Ferchaud AL, Pujolar JM. The effects of Medieval dams on genetic divergence and demographic history in brown trout populations. BMC Evol Biol. 2014;14:122.

    PubMed  PubMed Central  Article  Google Scholar 

  85. 85.

    Charlier J, Palme A, Laikre L, Andersson J, Ryman N. Census (NC) and genetically effective (Ne) population size in a lake-resident population of brown trout Salmo trutta. J Fish Biol. 2011;79(7):2074–82.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  86. 86.

    Serbezov D, Jorde PE, Bernatchez L, Olsen EM, Vollestad LA. Short-term genetic changes: evaluating effective population size estimates in a comprehensively described brown trout (Salmo trutta) population. Genetics. 2012;191(2):579–92.

    PubMed  PubMed Central  Article  Google Scholar 

  87. 87.

    Jorde PE, Ryman N. Demographic genetics of brown trout (Salmo trutta) and estimation of effective population size from temporal change of allele frequencies. Genetics. 1996;143(3):1369–81.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  88. 88.

    Charlier J, Laikre L, Ryman N. Genetic monitoring reveals temporal stability over 30 years in a small, lake-resident brown trout population. Heredity. 2012;109(4):246–53.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  89. 89.

    Franklin IR. Evolutionary change in small populations. In: Soulé ME, Wilcox BA, editors. Conservation biology: an evolutionary-ecological perspective. Sunderland: Sinauer; 1980. p. 135–49.

    Google Scholar 

  90. 90.

    Frankham R, Bradshaw CJ, Brook BW. Genetics in conservation management: revised recommendations for the 50/500 rules, Red List criteria and population viability analyses. Biol Conserv. 2014;170:56–63.

    Article  Google Scholar 

  91. 91.

    Franklin IR, Frankham R. How large must populations be to retain evolutionary potential? Anim Conserv. 1998;1(1):69–70.

    Article  Google Scholar 

  92. 92.

    Reed DH, Briscoe DA, Frankham R. Inbreeding and extinction: the effect of environmental stress and lineage. Conserv Genet. 2002;3(3):301–7.

    CAS  Article  Google Scholar 

  93. 93.

    Bekkevold D, Hojesjo J, Nielsen EE, Aldven D, Als TD, Sodeland M, Kent MP, Lien S, Hansen MM. Northern European Salmo trutta (L.) populations are genetically divergent across geographical regions and environmental gradients. Evol Appl. 2020;13(2):400–16.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  94. 94.

    Fraser DJ, Hansen MM, Østergaard S, Tessier N, Legault M, Bernatchez L. Comparative estimation of effective population sizes and temporal gene flow in two contrasting population systems. Mol Ecol. 2007;16(18):3866–89.

    PubMed  Article  PubMed Central  Google Scholar 

  95. 95.

    Griffiths SW, Armstrong JD. Kin-biased territory overlap and food sharing among Atlantic salmon juveniles. J Anim Ecol. 2002;71(3):480–6.

    Article  Google Scholar 

  96. 96.

    Olsén KH, Petersson E, Ragnarsson B, Lundqvist H, Järvi T. Downstream migration in Atlantic salmon (Salmo salar) smolt sibling groups. Can J Fish Aquat Sci. 2004;61(3):328–31.

    Article  Google Scholar 

  97. 97.

    Hoare D, Krause J, Peuhkuri N, Godin JG. Body size and shoaling in fish. J Fish Biol. 2000;57(6):1351–66.

    Article  Google Scholar 

  98. 98.

    Killen SS, Marras S, Nadler L, Domenici P. The role of physiological traits in assortment among and within fish shoals. Phil Trans R Soc B. 2017;372(1727):20160233.

    PubMed  PubMed Central  Article  Google Scholar 

  99. 99.

    Carlsson J, Carlsson JEL, Olsen KH, Hansen MM, Eriksson T, Nilsson J. Kin-biased distribution in brown trout: an effect of redd location or kin recognition? Heredity. 2004;92(2):53–60.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  100. 100.

    Bohlin T, Dellefors C, Faremo U. Date of smolt migration depends on body-size but not age in wild sea-run brown trout. J Fish Biol. 1996;49(1):157–64.

    Article  Google Scholar 

  101. 101.

    Griffiths SW, Brockmark S, Höjesjö J, Johnsson J. Coping with divided attention: the advantage of familiarity. Proc Royal Soc B: Biol Sci. 2004;271(1540):695.

    CAS  Article  Google Scholar 

  102. 102.

    Skaala O, Besnier F, Borgstrom R, Barlaup B, Sorvik AG, Normann E, Ostebo BI, Hansen MM, Glover KA. An extensive common-garden study with domesticated and wild Atlantic salmon in the wild reveals impact on smolt production and shifts in fitness traits. Evol Appl. 2019;12(5):1001–16.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  103. 103.

    Fernandes WPA, Ibbotson AT, Griffiths SW, Maxwell DL, Davison PI, Riley WD. Does relatedness influence migratory timing and behaviour in Atlantic salmon smolts? Anim Behav. 2015;106:191–9.

    Article  Google Scholar 

  104. 104.

    Meli A, Fraser DJ. Kinship analysis of brook trout Salvelinus fontinalis during their breeding migration. J Fish Biol. 2013;82(5):1514–22.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  105. 105.

    Martinez J, Moran P, Perez J, De Gaudemar B, Beall E, Garcia-Vazquez E. Multiple paternity increases effective size of southern Atlantic salmon populations. Mol Ecol. 2000;9(3):293–8.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  106. 106.

    Goodwin JC, Andrew King R, Iwan Jones J, Ibbotson A, Stevens JR. A small number of anadromous females drive reproduction in a brown trout (Salmo trutta) population in an English chalk stream. Freshwater Biol. 2016;61(7):1075–89.

    Article  Google Scholar 

  107. 107.

    Pulido F. Evolutionary genetics of partial migration–the threshold model of migration revis (it) ed. Oikos. 2011;120(12):1776–83.

    Article  Google Scholar 

  108. 108.

    Curry RA, Bernatchez L, Whoriskey F, Audet C. The origins and persistence of anadromy in brook charr. Rev Fish Biol Fisher. 2010;20(4):557–70.

    Article  Google Scholar 

  109. 109.

    Holecek DE, Scarnecchia DL, Miller SE. Smoltification in an impounded, adfluvial redband trout population upstream from an impassable dam: does it persist? T Am Fish Soc. 2012;141(1):68–75.

    Article  Google Scholar 

  110. 110.

    Quinn TP, Bond MH, Brenkman SJ, Paradis R, Peters RJ. Re-awakening dormant life history variation: stable isotopes indicate anadromy in bull trout following dam removal on the Elwha River, Washington. Environ Biol Fish. 2017;100(12):1659–71.

    Article  Google Scholar 

Download references

Acknowledgements

We acknowledge Britt Iren Østebø for extensive field work and data punching, Bjørghild Breistein Seliussen for her help with the genotyping process, and Dr. Baocheng Guo and two anonymous reviewers for many insightful comments and suggestions that significantly improved the paper.

Funding

The study was funded by the Norwegian Ministry of Trade, Industry and Fisheries.

Author information

Affiliations

Authors

Contributions

ØS, ED, KAG and MMH conceived and planned the study. ED, MQ, GD, AD and VW genotyped the samples and conducted the statistical analyses. ED, MMH, ØS, and KAG wrote the manuscript with input from the other authors. All authors read and approved the final version of the manuscript.

Corresponding authors

Correspondence to Eloïse Duval or Øystein Skaala or Michael M. Hansen.

Ethics declarations

Ethics approval and consent to participate

All experimental protocols were in accordance with permits obtained by the Hordaland County Governor, and protocols approved by the Norwegian Environment Agency and Norwegian Food Safety Authority. Specifically, this involves permits to capture, tag, and release smolt and recapture adult fish. Permit numbers are 07/13020-1, 2011/88155-1 and 17/15858-1. All applications and protocols for animal handling have been approved by the Institute of Marine Research's local animal welfare unit before further processing by the Norwegian Food Safety Authority.

Consent for publication

Not applicable.

Competing interests

The authors declare no conflict of interest.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

Supplementary material on microsatellite loci used (Appendix 1), PCR conditions (Appendix 2), overview of anadromous trout returning to Guddal River several times during the study period (Appendix 3), and proportion of anadromous trout ascending the river in different weeks in specific years (Appendix 4).

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Duval, E., Skaala, Ø., Quintela, M. et al. Long-term monitoring of a brown trout (Salmo trutta) population reveals kin-associated migration patterns and contributions by resident trout to the anadromous run. BMC Ecol Evo 21, 143 (2021). https://doi.org/10.1186/s12862-021-01876-9

Download citation

Keywords

  • Partial migration
  • Salmo trutta
  • Life-history types
  • Parentage assignment
  • Sibship reconstruction
  • Migration timing
  • Effective population size