- Research article
- Open Access
Population connectivity: dam migration mitigations and contemporary site fidelity in arctic char
BMC Evolutionary Biologyvolume 11, Article number: 207 (2011)
Animal feeding and spawning migrations may be limited by physical barriers and behavioral interactions. Dam constructions (e.g. hydropower) commonly include gateways for fish migrations to sustain ecological connectivity. Relative genetic impacts of fish passage devices versus natural processes (e.g. hybrid inferiority) are, however, rarely studied. We examined genetic (i.e. microsatellite) population connectivity of highly migrating lake-dwelling Arctic char (Salvelinus alpinus), introduced 20 generations ago, across and within two subalpine lakes separated by a dam with a subterranean tunnel and spill gates after 7 generations. Due to water flow regime, the time window for fish migration is highly restricted.
Char populations, with similar genetic structuring and diversity observed across and within lakes, were admixed across the dam with fishways during feeding. For spawning, however, statistically significant, but very low population differentiation (θ; 0.002 - 0.013) was found in nine out of ten reproductive site comparisons, reflecting interactions between extensive migration (mean first generation (F0) = 10.8%) and initial site fidelity. Simulations indicated that genetic drift among relatively small effective populations (mean N e = 62) may have caused the observed contemporary differentiation. Novel Bayesian analyses indicated mean contributions of 71% F0 population hybrids in spawning populations, of which 76% had maternal or paternal native origin.
Ecological connectivity between lakes separated by a dam has been retained through construction of fishways for feeding migration. Considerable survival and homing to ancestral spawning sites in hybrid progeny was documented. Population differentiation despite preceding admixture is likely caused by contemporary reduced reproductive fitness of population hybrids. The study documents the beginning stages of population divergence among spatial aggregations with recent common ancestry.
Genetic differentiation, fundamental to population genetics, is initiated by restricted gene flow and reproductive isolation mechanisms within gene pools [1, 2]. Habitat and population fragmentation with concomitant disruption of ecological connectivity is threatening biota worldwide, on large scales determined by time since isolation and physical reproductive barriers among demes . On finer scales, genetic differentiation may additionally involve reproductive behaviors and social interactions [4, 5]. Addressing genetic effects of physical barriers and reproductive behaviors in novel populations may forecast evolutionary consequences of human intervention.
Natural ecological connectivity may fluctuate in time and space, but anthropogenic alterations (e.g. water regulation) typically exaggerate this variation. Human impacted water flows and water habitat fragmentation expanded greatly worldwide in the twentieth century . In Europe (excluding Russia) close to all large river systems are fragmented by dams , inevitably presenting impediments to migration. Construction of habitat corridors to sustain ecological connectivity may in theory counteract negative fragmentation effects , although empirical evidence is limited. A number of often retrospective fish passage improvements have been designed to reduce ecosystem wide effects of (hydropower) dams . Dams commonly allow passive or active downstream drift of fish, but obstruct upstream migration. World Commissions of Dams with contributing papers [10, 11] and recent reviews [12–14], highlight the need for empirical studies addressing long term effects of fish passage success. Migration barriers may obstruct passage directly or negate otherwise adaptive benefits of homing and assortative mating indirectly [15–17].
Genetic drift and behavior may induce reproductive isolation mechanisms operating on fine scales. Genetic microstructure even within continuous habitat is reported [18, 19], although the biological relevance of weak differentiation is debated [20–22]. Within-lake genetic structure for highly mobile fish without obvious migration barriers have similarly been documented [23–28] for a traditional evolutionary time perspective. Genetic structuring and adaptive microevolution may, however, also occur on contemporary time scales [29–33]. Recently founded or invasive species provide opportunities to study initial genetic divergence and evolution in the wild , but few studies have evaluated contemporary genetic structure in the invasive species [35–37]. Moreover, we are not aware of assessment of initial genetic structure in recently founded and highly migrating species that are temporally (i.e. seasonally) admixed. Recently founded salmonid species in freshwater lake habitats represent natural small scale experiments that provide opportunities to study the details of early stages of population divergence.
Micro-scale studies under initial admixture introduce challenges concerning sampling design and analytical approaches. For obvious reasons, sample sizes are often restricted both in numbers and temporal replicates. In continuous habitats sampling units may be difficult to delineate. Widely used Bayesian methods to infer cryptic population structure [38–40] may not, however, detect population substructure of recently established populations with low differentiation. Inclusion of population migrants and hybrids (i.e. admixed ancestry) in presumably distinct populations may lead to underestimates of differentiation and erroneous conclusions regarding reproductive isolation.
Reproductive isolation may be reinforced by homing, as documented within a wide range of species [41, 42]. Successful reproduction of migrants may be common, caused by relaxed selection against non-natives resulting in non-native offspring, or hybrid offspring with admixed ancestry which break down incipient population differentiation [5, 43]. Further differentiation requires reduced reproductive success of non-native and hybrid progeny. Therefore, studies of population admixture and hybrid contribution to population segregation under incomplete divergence are important and may possibly unveil biological relevant differentiation. Traditional measures of population differentiation combined with recently developed Bayesian statistics  now permit studies of how inter-population hybrids may affect population differentiation.
Here we study genetic differentiation of Arctic char (Salvelinus alpinus), introduced around 1920 (i.e. 20 generations ago), across and within two lakes: Pålsbufjorden (PAL) and the downstream Tunhovdfjorden (TUN) (Figure 1). The two lakes, naturally connected by a short river, were separated by a dam in 1946, creating two hydropower reservoirs. Both lakes have geographically distinct char spawning sites (PAL; P1 - 3, TUN; T1 - 2) for which we assess population connectivity implementing 10 microsatellites. The first objective was to examine genetic structure across lakes to quantify present effects of two migration enhancements at the dam; a subterranean tunnel and spill gates. An additional summer sample (X) of lake feeding and presumably admixed char from PAL was included to assist interpretation of differentiation. We expected to find greater genetic differentiation across lakes than across populations within lakes, and assignment of sample × to PAL populations. We also expected that passive migration downstream would exceed upward migration and result in observed asymmetric migration and associated genetic variation. The second was to examine migration patterns and population hybrid contributions potentially responsible for contemporary genetic substructure. Due to expectations of low population differentiation in the progress towards reproductive isolation, we intend to deliberate over the paradox of population divergence in the face of migration, rather than to state absolute assertions.
Amplification and allele calling of 190 individuals over 6 sites was obtained in 99.8% of the cases, and secondary amplifying and allele calling of a subsample gave consistent results. Number of alleles per locus averaged 15.4 (SD ± 7.2, range 2 - 26, Additional file 1). Quality control screening did not reveal indications of scoring error due to stuttering or evidence for large allele dropout. The test for null alleles assuming Hardy-Weinberg equilibrium (HW) suggested that null alleles may be present for locus Sco204 in sample site T2, by general excess of homozygotes. However, the combined probability for all homozygote class frequencies was not significant (P > 0.05). No deviations from HW were found, but nonrandom association of alleles from different loci was indicated in 12 and 9 (Sco 213/Sfo8, Smm17, Smm24, Str73; Sfo8/Smm17, Smm24, Str73; Smm24/Smm17, Str73) of 45 tests uncorrected and after false discovery rate (FDR) correction  at 5% level, respectively. Genetic variation revealed a mean value of local FST = 0.024 (range 0.013 - 0.036), and inbreeding F = 0.028 (range 0.003 - 0.053) over sites.
Dam barrier to migration and gene flow
Measured lake water levels and estimated velocities in the fishways indicated possible upstream migration only within very limited windows of time. During 2000 - 2010, migration through the tunnel was feasible in less than 3.6% of the year (mean 13 days per year, SD ± 20, min 0, max 55). Migration through the spill gates was possible in less than 1.7% of the year (mean 6 days per year, SD ± 4, min 1, max 16). Upstream migration was possible in all seasons through gates, but only through the tunnel in spring.
Observed heterozygosity and allelic richness was similar between lakes (Additional file 1). Number of private alleles found totaled 22, 6 in TUN and 16 in PAL (data not shown).
The total allele frequency variation consisted of 99.3% (P < 0.001) variation within populations, 0.7% (P = 0.001) variation among spawning sites within lakes and no variation between lakes (P = 0.499). Genetic differentiation of char from the two connected lakes PAL and TUN was, however, highly significant (P = 0.001) but very low (θ = 0.003, standardized (θ') = 0.014). Removal of first generation migrants (F0), amounting to 7 and 9% in PAL and TUN (not significant after correction, 170 tests, threshold P = 0.0003), respectively, did not affect this differentiation. STRUCTURE did not detect genetic structure among lakes or sites (K = 1).
Estimated FST and exact tests revealed that one spawning population (T2) was significantly different from a summer feeding sample (X) in PAL (Table 1). Sample × assigned to both lakes (TUN, 25%), and 30% was most similar to spawning site T1 on a population level (Table 2), indicating upstream migration.
Migration beyond F0 was assessed with a maximum likelihood approach for assignment and with Bayesian estimates for present and contemporary estimates. Self-assignment of individual char in GENECLASS2 revealed 93 and 46% correct assignments to Lake PAL and TUN, respectively, suggesting a downstream source-sink structure. Probability estimates of ancestral origin over the last generation estimated in the software program BIMR, also indicated high, but more symmetric migration rates. With the two lakes as populations, mode allele frequencies showed almost equal ancestral native origin for the two lakes (PAL; 56%, highest posterior density intervals (HPDI): 34 - 75, TUN; 59%, HPDI: 25 - 90), with estimated HPDI for non-native origin of 24 - 65% and 9 - 74%, in PAL and TUN respectively.
Mode probability of hybrid ancestral origin for the pooled samples within lakes estimated in BIMR were 50% in PAL (HPDI = 36 - 49) and also 50% (HPDI = 13 - 48) in TUN. Nevertheless, pure origin from the home lake was more common than pure origin from the other lake.
Lentic population migration and admixture
All except one test of genetic differentiation among spawning sites were significant (mean θ = 0.007, SD ± 0.003). Mean standardized differentiation (θ') was 0.030 (SD ± 0.016). Exact tests corroborated the pattern of significantly differentiated spawning sites (Table 1). Using spawning sites as prior populations, 37% of the char assigned to its sampled site, varying from 12 - 62% among sites (Table 2). Mean probability of best assignment, however, was only 51% (SD ± 28), indicating recent admixture, migration, or un-sampled populations. Exclusion of 18 char (11%) found to be F0 migrants across spawning sites (range 1 - 5, Table 2, not significant after correction, 170 tests, threshold P = 0.0003) increased population segregation and resulted in all spawning populations being significantly different (θ > 0.002, θ' > 0.008, P < 0.021). This finding corroborates the estimated migration rates. Analyses in IMMANC verified migrant detection (8%) with a power of 0.968. The Bayesian probability that the sampled alleles within spawning sites originated from the same site last generation (ancestral rate) were 0.47 across sites (SD ± 0.09, min 0.34, max 0.60), and higher than the probability of origin from any other population (mean 0.13, SD ± 0.10, min 0.02 max 0.29).
Mixed ancestry analyses indicated that population hybrids constituted a substantial proportion within spawning sites. The majority of hybrids within all sites were progeny of native and non-native char (mean 54%), whereas hybrids of two non-natives were less common (mean 17%, Figure 2). Similar to pooled lake estimates, the mode probability of home origin within sites was consistently higher than pure non-native origin within sites, and consequently higher for hybrids with partial native ancestry than for any other mixed ancestry across sites. Standard deviations of mixed ancestry estimates were generally similar to mode estimates.
Potential interpretation bias
Potential bias and biological relevance of observed data were addressed to assure sound interpretations. FST tests excluding loci with possible LD (Sco213 and Sfo8) revealed similar differentiations as for all loci combined (data not shown). The FST estimate per locus averaged 0.0064 (SD ± 0.006, min 0.001, max 0.023). Exclusion of single loci generally weakened spawning site differentiation, but did not change the above pattern, thus indicating neutral behavior of markers.
AMOVA analyses segregated based on year-classes revealed 100% variation within sites and a lack of variation among year-classes (FST = -0.001, P = 0.624), whereas hierarchical analyses based on spawning sites revealed that 99.45% of genetic variation stemmed from within-population variation and 0.55% (FST = 0.006, P = 0.002) from among-population variation. Thus, genetic variation was congruent with the high variability of microsatellites used, indicating temporal stability among consecutive year-classes. Pairwise year-classes were not significantly different within any spawning site (θ < 0.060, P > 0.013), or pooled across sites (θ < 0.004, P > 0.184). Removal of year-classes one-by-one changed the significance in 15.5% of the pairwise tests, but did not affect the pattern of differentiation among spawning sites.
We did not find evidence for family structures affecting population differentiations. Simulations indicated a power of 1.000 and 0.918 to discriminate unrelated from full sib and half sib, respectively. Only 0.1% of pairwise tests had relatedness coefficients ≥ 0.25, indicating half sib relationship, and half of these were across sites. Mean distribution of relatedness (LRM) for the pooled dataset was negative (μ = -0.003, SD ± 0.045, Figure 3), although positively skewed (1.640) and not normally distributed (P = 0.010). Sites individually revealed the same pattern (not shown), with lack of bimodal patterns that could have indicated kin-groups.
Effective population sizes (N e ) were estimated to be 52 - 71 individuals within spawning sites, 133 - 142 within lakes, and 233 in total (Table 3), supporting both across and within lake genetic structure. Estimates were not affected by prior values of N e . Estimates of effective parental population size (N b ) were generally higher but similar to N e , except negative and infinite for population P1. This may indicate sampling error, or no evidence for LD caused by genetic drift due to a finite number of parents. Tests for LD did not, however, give significant deviations for any loci combination at this site, after FDR correction. Recent bottlenecks over all loci were not confirmed for any of the tests (P > 0.080), assuming a mixed mutation model (TPM). Neither did tests for di-nucleotide and tetra-nucleotide loci separately, under TPM and a single-step mutation model (SMM) respectively, reveal signs of severe loss of allelic diversity (normal L-shaped distribution, P > 0.060), except for population P2 and tetra-nucleotide loci, where the Wilcoxon test for heterozygosity excess was significant (P = 0.030). Test for isolation by distance (IBD) was positive but not significant (r = 0.29, Pupper = 0.250, Plower = 0.760).
Forward simulated FST estimates (Figure 4), indicated that genetic drift alone may have caused the observed population differentiation. For instance, five populations with N e = 60 were significantly differentiated (mean θ = 0.005, SD ± 0.002, P < 0.047) after a mean migration rate of 0.3 in 20 generations, assuming initial admixture after introduction. Ten independent replicates of this scenario indicated similar differentiation (mean θ = 0.006, SD ± 0.001). IMMANC indicated detection power for migrants (8%) and hybrids (7%) in our sampled dataset to be 0.968 and 0.638, respectively.
A forward simulated alternative hypotheses of m = 0.01 revealed high differentiation (mean θ = 0.122, SD ± 0.019, P < 0.001) between reproductive sites, whereas m = 0.99 revealed lack of differentiation (mean θ = 0.001, SD ± 0.001, P > 0.176). A scenario without migration only indicated 5.0% (1.3% after correction) F0 migrants in GENECLASS, similar to expectations from type 1 errors. Similar analyses in IMMANC indicated 10.7% (1.0% after correction) F0 migrants and 10.0% (1.3% after correction) F1 hybrids, with mean power > 0.999. Coalescent simulations supported that drift alone may have caused contemporary site fidelity (Additional file 2).
Our results suggest that constructed gateways for migration ensured alternate habitat utilization and gene flow across the dam separating the two studied lakes. Significant differentiation was found between spawning char from the two interconnected lakes, but differentiation among populations within lakes equaled or exceeded that for lakes. Without assuming any direct relationship between FST and gene flow (below), these relative measures are low and indicate high connectivity across the dam. Char with genetic assignment to distinct populations in both lakes were caught during feeding migration in the upper lake in summer, indicating upstream migration across the dam. Corroborating this, gene diversity and genetic structure were similar across and within lakes.
Extensive migration reflecting seasonal foraging admixture across lakes and populations is independently corroborated by tagging studies [46, 47]. Barrier type, hydrology and life stages of focal species all interact in determining the impact of barriers . Tunnels associated with turbine outflow at dams are known to attract ascending salmonids [49–51]. The lack of turbines in the PAL dam tunnel makes migration through the tunnel possible. It is an open question, however, to what extent char in our study lakes utilize the spill gates or the tunnel for (upstream) migration. Ongoing studies in brown trout also reveal high dispersal across the focal dam, with unknown proportional use of the fishways (unpublished). Regardless, results indicate that char, in both cases, take advantage of the highly restricted time window suitable for upstream migration. Disentangling present from previous genetic structure and biological interpretation of low differentiation is, however, challenging (below), and results should therefore be interpreted with caution.
Contemporary genetic substructure
The present study documents that genetically differentiated sympatric populations can be established within decades. Pairwise genetic differentiation was found to be significant, but very low, both among lakes and among most spawning sites within lakes, indicating spatial reproductive isolation even within continuous habitat. Estimates of N e supported this. However, STRUCTURE did not reveal any population structuring. Performance of Bayesian clustering methods decreased for FST < 0.02 , and we propose that even the novel algorithm used in STRUCTURE  may not detect recently established populations with very low differentiation (FST < 0.01).
Low differentiation is expected due to the short time since species introduction, and the migratory behavior of char. Following the simplistic [53, 54] equation FST ≈ 1/4Nm + 1), our mean local FST estimates indicate long term Nm = 10.2, and migration rates of 0.17 assuming N e = 60. Fixed m = 0.5 would reveal the same differentiation if N e is as low as 20. As an example, 60 individuals and m = 0.5 result in the same differentiation as e.g. 1000 individuals with 3.0% gene flow. In contrast, if we assume F0 analyses (mean m = 0.11) are representative for present gene flow, and implement N e 's of 60, the present Nm is 6.6, likely overestimated due to common reproductive inferiority in migrants. Consequently, mean gene flow between populations may have decreased since introduction. Thus, present genetic structure does not reflect panmixia . Both estimates, however, indicate that gene flow is too small to bring about drift connectivity , but too large for Bayesian detection of subdivision . A denotation of the sampled reproductive entities as populations, according to the evolutionary and especially the ecological paradigm  is, however, debatable. Maximum FST estimates will be limited by heterozygosity , and estimated differentiation may underestimate true population differentiation , especially when gene diversity is high . Similar standardized population differentiation is previously found in char . Direct comparisons of observed θ are complicated, but similarly low and biological relevant differentiations are found in both mammals and fish [59, 60]. Interestingly, similarly low, but significant differentiations are also found among populations of char founded thousands of years ago , and population structure of lake-dwelling brook char (Salvelinus fontinalis) has been found on a smaller geographical scale under migration-drift equilibrium  than in our study. Low but significant pairwise differentiations are also found in invading catadromous crustaceans separated after recent colonization . This study also documented strong genetic drift simultaneously with significant differentiation (FST > 0.007), but among year classes, as opposed to our study.
True substructures in our study were corroborated by N e estimates, genetic variability, year-class stability and low number of sib-groups within sites. Invasive species, such as char in our lakes, are expected to express low genetic variation and N e , caused by founder events and population bottlenecks [63–65]. The high level of genetic variation within the examined loci in our populations does not suggest depauperate populations or bottlenecks on the surveyed scale, which would otherwise certainly support genetic differentiation caused by strong genetic drift in very small founding populations. The concept of founder populations, however, suggests that current structure is novel. We are not aware of evidence from the literature that fish populations bred in admixture may resume historic population structures in novel habitat. Simulations, however, documented that drift alone, even in larger populations, could cause rapid differentiation. Estimation of population size based on single samples relies on a number of simplifying assumptions [66, 67]. Our estimated population sizes should therefore be interpreted with caution. Echosounder and gillnetting experiments  suggest a population size of char around 9000 in PAL, returning N e /N ≈ 0.015, which is common for marine fish  but low compared to most salmonids [69, 70]. Nevertheless, estimated N e and N b gave similar results, and suggest population sizes sufficient to confer contemporary genetic stability as well, as opposed to a transitory or strict metapopulation structure . It is possible, however, that estimates are inflated by immigrants and sampling of hybrid zones, even though no deviations from HW were evident. Regardless, the question remains whether one is actually able to sample biological populations correctly, and if hybrid contribution represents true population structure. The presence of immigrants and hybrids, not contributing to the reproductive population gene pool, can easily complicate interpretations of genetic diversity, N e and genetic differentiation.
Initial genetic structure
The low level of interspecies competition and vacant niches fits well with rapid establishment of char in our lakes, and zooplanktivores such as char show a high rate of invasion success . Genetic differentiation by means of the evolution of reproductive isolation can be viewed as non-ecological modes of differentiation such as founder effect or genetic drift in small populations, or by divergent natural selection [73, 74]. Lag time between initial colonization and the onset of rapid population growth is expected in invasive species , in particular if evolutionary adaptation is important during colonization. Thus, the rapid expansion in the lakes surveyed  corroborates the independent simulations, pointing to genetic drift as an initial diversifying factor for observed population differentiation. However, theory [1, 77, 78], experiments  and empirical studies [80–85] all suggest a potential for rapid evolutionary changes involving adaptation in newly founded populations. Moderate levels of plasticity are also typical for char [86, 87], and may also have been optimal for fast genetic evolution [81, 88–90].
The combined genetic and ecological data suggest that the initial genetic structure and diversity within our lakes is a result of site fidelity and initial drift. Several evolutionary forces may, however, interact at the same time. Realistic and rapid drivers of cryptic kin selection  or discontinuous adaptive variation [32, 92–94] may add to incipient reproductive isolation. Once sub-structure is established, site fidelity, drift, and low effective migration rates may strengthen differentiation in time and space. Such isolation by adaptation  is likely facilitating drift in neutral loci by reduced gene flow as a general barrier mechanism , and even neutral markers may detect adaptation in the face of intermediate migration .
Differentiation despite extensive migration
Substantial migration among reproductive sites was found, without significant IBD, suggesting little present drift. Thus, although common in salmonids, equilibrium between gene flow and genetic drift is not present [98, 99]. Migration analyses may be biased by method assumptions, un-sampled populations, low differentiation and convergence problems [44, 100, 101]. However, all methods used to interpret migration in our study, indicated considerable migration among lakes and populations. Extensive migration among populations is often seen in salmonids . Homing and kin discrimination will, however, contribute to genetic structuring and is widely known in fish , particularly char [103–106]. Pairwise genetic differentiation of other neighbouring lacustrine char populations has been found to be highly significant, despite long term migration rates of 1.853 - 2.755 individuals per generation . Low levels of effective migration, i.e. gene flow among sympatric populations, contrary indicate breeding site fidelity in established populations [107, 108]. Reproductive divergence despite extensive migration in the study at hand indicates restricted effective gene flow between lakes and populations, even after 20 generations. Thus, at first glance, one may be tempted to conclude that results indicate reproductive selection against immigrants in newly established populations, along the lines of previous studies [77, 83, 109]. Alternatively, one may speculate that breakdown of non-native gene flow is caused by reduced reproductive fitness of population hybrids.
Hybrid contribution during naturalization
The hierarchical Bayesian analyses of mixed ancestry confirmed non-native spawning success. The large proportion of hybrids documents reproduction by immigrants, i.e. relaxed selection against non-natives early in differentiation and survival of hybrids. Selection against population hybrids and immigrants strengthen divergence [110, 111]. Hybrids of sympatric and closely related (FST = 0.070) perch (Perca fluviatilis) were found to have reduced fitness in early life stages in laboratory experiments, but the authors did not test for hybrid inferiority in vivo . Novel habitats without competition as in our study may have promoted survival of hybrids, despite potentially reduced fitness. However, hybrids may be less competitive in secondary stages of naturalization, as increased population divergence reduces hybrid fitness .
In our study, hybrids without parental origin from their spawning site (Figure 2; non-native hybrids) constituted a minor part of mature char caught at their reproductive site. The considerably larger proportion of individuals with partial native ancestral origin indicates hybrid homing. Reproductive units sustained by kin discrimination are found in a range of species [114–117], and hybrid juveniles of char will be in close proximity to kin when hatching. Population differentiation, however, could only be established if reproductive success of returning hybrids is low. Admixed individuals would otherwise have caused migration loads obscuring effects of selection and differentiation . Postmating reproductive isolation and hybrid inferiority have been addressed for a century [2, 111]. Both theoretical and empirical studies frequently reveal selection against hybrids during speciation, although its origin seems unclear [74, 96]. Sexual selection (mate choice) against hybrids is found across taxa [e.g. insects; , fish; ], and may also characterize char. Such mechanisms may evolve faster in small populations , as within our sites. Social recognition of relatives in fish using odour cues can induce assortative behavior in fish [121–123]. Finally, high dispersal rates may increase the absolute number of hybrids, but the increased competition reduces hybrid reproductive fitness, and consequently increases reproductive isolation of the resident population .
Migration corridors between populations separated by dams may be valuable for sustaining evolutionary potential. This study demonstrates that even temporally very limited connectivity through a subterranean tunnel and spill gates between two regulated lakes, likely counteract genetic isolation in char.
Initial stages of genetic divergence of subpopulations, despite high migration rates, are documented in this study. We are not aware of previous studies addressing hybrid migration and reproductive success in recently founded and admixed populations. Our study indicates that population differentiation may be detectable, despite inclusion of hybrid sub-populations, when assessing genetic structure among populations. Combined also with independent mark-recapture data by Aass [46, 47], our genetic approach indicates that hybrids contribute extensively to migration rates in recently founded populations. Population hybrid events may have been important in establishing the diverse population structure in the novel habitat . However, present differentiation implicates nascent non-native and hybrid inferiority, whether biologic or ethologic [32, 91]. Estimates generally gave wide HPDI, an inherent problem expected from low differentiation affecting accuracy , and results must be interpreted with caution.
Few population genetic software programs have been evaluated with respect to their performance in detecting low genetic differentiation, making interpretation of differentiation close to admixture difficult [53, 124]. However, restricted power in admixture analyses should not restrict studies of hybridization mechanisms, as they raise consequential questions in micro-evolution and behavioral ecology. Both small scale empirical and larger simulation studies of incomplete reproductive isolation may guide management of invasive and naturalized species, potentially unveiling initial population differentiation mechanisms.
The study species
The Arctic char is a salmonid fish with Holarctic distribution  showing extreme phenotypic and life history variation, exemplified by weight at maturation from 3 g to 12 kg . This reflects the species capability to evolve trophic polymorphism and possible genetic differentiations within drainages  and lakes [107, 127, 128]. While most populations of char result from natural postglacial colonization, high altitude populations are introduced by man. When in sympatry, Arctic char typically have a benthopelagic distribution in landlocked habitat. Young char typically feed in the littoral, shifting towards a highly migratory pelagic feeding pattern as adults. Lacustrine char are commonly 20 - 40 cm in length and a few years old when mature, and spawn in the littoral within a few weeks in autumn. Site fidelity is common, but reproductive ecology is poorly understood (see Jonsson et al. 2001 , Klemetsen et al. 2003  and Johnson 1980  for excellent reviews).
The study site
Lake Tunhovdfjord (TUN) and Lake Pålsbufjord (PAL) are part of a 35 km long sub alpine hydroelectric reservoir located in south-central Norway (48°E, 67°N, Figure 1), regulated first in 1919. A hydropower dam separating the two lakes was erected in 1946, restricting the previously free migration among the lakes to migration through spill gates in the dam (lower level 734.9 m a.s.l., c 10 m, c 3.1 m2, neutral gradient) and a subterranean anthropogenic branch (lower level 722.4 m a.s.l., 1300 m, 7.1 m2, neutral gradient) with outflow 600 m below the dam. Lake PAL now has a surface area of 5.25 - 19.5 km2 with maximum depth 25 m, 725.5 - 749.1 m a.s.l.. Lake Tunhovdfjord, located immediately downstream, has a surface area of 14 - 25 km2 and a maximum depth 70 m, 716.4 - 734.4 m a.s.l. Arctic char are sympatric with brown trout (Salmo trutta) and invasive European minnow (Phoxinus phoxinus) in both lakes. Brown trout have probably been native for > 6000 years [130, 131]. Minnow were introduced coincidentally around 1920.
In 1910, 10.000 fry from a natural char population in Lake Tinnsjø (187.2 - 191.2 m.a.s.l., 51 km2, mean depth 190 m), situated in a different watershed 80 km south-west of PAL and TUN were stocked in two alpine lakes above PAL [132, 133]. Char were first observed in PAL in 1919. The population expanded rapidly and was present in large numbers in both PAL and TUN after a decade . Densities have thereafter been highly variable, probably because of large water level fluctuations and associated egg mortality [68, 134]. Char in these lakes seldom become more than 8 years old and have a generation time of about 4 years (20 generations since introduction). They are usually < 200 g , although cannibalistic individuals of 3 kg have been reported . Extensive mark-recapture studies have revealed substantial char movements within and between the lakes , even after the construction of the dam in 1946. Char feeding migration intensity increases in late summer and autumn. The majority of adult char congregate in outlets and narrow parts of the lakes after spawning in autumn when water level drops, followed by c. 10 - 50% passive displacement downstream in the spring . Char emigrating PAL often attempt to re-enter the lake, but successful upstream migration is unknown .
Arctic char are low performance swimmers as opposed to brown trout . Maximum relative swimming speed (body length s-1) is about 2.8 , thus a water velocity of 150 cm s-1 is likely an absolute limit for upstream migration for studied char across the dam. The Mannings formula; v = M * Rh2/3 * √I , assuming M = 35, was used to estimate elevation heads associated with 150 cm s-1 water velocity in the fishways. The water level below the dam in TUN commonly increases (up-arches) to the lower level of the gates. There is also free entrance from TUN to the subterranean tunnel independent of water level. Thus, migration was estimated to be solely restricted by water levels in PAL at 734.9 - 735.2 m a.s.l. in the gates, and 722.4 - 726.4 m a.s.l. in the tunnel. The water level in PAL was registered from year 2000 - 2010.
A total of 190 mature char (mean total length 244 mm, SD ± 31, range 170 - 334) were sampled by gillnetting at four sites in PAL and two sites in TUN (Figure 1). Samples represented 54% males (i.e., presence of gametes) and the year classes 1999 - 2003 (N = 31, 37, 48, 37 and 2, respectively, as determined by scales and otoliths of 155 char). All char except those from T2 were aged, but these fish had similar size distribution. Samples from site P1 - 3 and T1 - 2 each included 34 char ready to spawn this season (gonadal stage 4 - 6) caught at traditional and geographically separated spawning sites below minimum water level (LRV) during spawning time, and are therefore treated as discrete populations, although true population boundaries are unknown. Additional spawning sites within the two lakes are not known. Site × (N = 20) is a randomly chosen mid-lake sample from late summer two months prior to spawning, representing admixtured non-spawning char.
Microsatellite genotyping and variability
Tissue samples of approximately 2 mm2 from the adipose fin were preserved in 96% ethanol in the field, and DNA was isolated using DNeasy kit (QIAGEN), following the manufacturer's guidelines. Microsatellite polymorphism was analyzed by means of 10 di- and tetra-nucleotide loci known to be polymorphic in S. alpinus; Mst-85 , Sco202, Sco204, Sco213 , Sfo-8, Sfo-23 , Smm-17, Smm-24 , Ssa-85  and Str73 . One primer for each locus was end-labeled with fluorescence (HEX, FAM and NED), and run partly as multiplex PCR reactions; multiplex A; Sfo-8, Smm-17, Smm-24, Str73 and Sco21, multiplex B; Sco202 and Sco204, and multiplex C; Sfo-23 and Ssa-85. Each PCR contained 2 μl genomic template DNA and 8 μl reaction mixtures containing 1 - 2 pmol primer, 50 mM KCl, 1.5 mM MgCl, 10 mM Tris-HCl, 0.2 mM dNTP and 0.25 U Taq polymerase enzymes (Ampliqon). Thermocycling parameters after denaturation at 95°C for 2 min were 24 - 34 cycles of 95°C for 30 sec, annealing temperature of 55°C for 30 sec, followed by an extension at 72°C for 45 sec. The last polymerization step was extended to 10 min. PCR products were added to buffer containing formamide and labeled standard (ROX Std 400, Applied Biosystems), and electrophoresed using an ABI Prism 3100 Genetic Analyzer (Applied Biosystems). The software GENEMAPPER v.3.7 (http://www.applied-biosystems.com) was used to score alleles, and all automated allele calling were controlled for by manual reading. Two positive controls were done in each run, and scoring was repeated twice for several loci, and for all private alleles, to check for consistency.
Quality control screenings were performed by testing for null alleles, large allele drop-outs and scoring errors in MICRO-CHECKER v.2.2.3 . The program TFPGA v.1.3  was used for descriptive statistics (number of alleles, observed and expected heterozygosity). Allelic richness was compared among lakes and spawning sites based on minimum sample size for each comparison separately in FSTAT v.184.108.40.206 . Departure from Hardy-Weinberg (HW) via separate one-tailed tests for heterozygote excess and deficiency for each locus in each site was tested by 60.000 randomizations, and a linkage disequilibrium (LD) test between pairs of loci across populations was performed by 450.000 permutations, both in FSTAT. Local FST (diversity standardized) and inbreeding coefficient F were calculated in BIMR v.1.1 , including 95% C.I., the F-model  and default values, except burn-in and sample size of 105 iterations, to ensure convergence.
Genetic structure and temporal stability
Possible genetic segregation among lakes and sites and their significance were tested with the FST analogue θ  and pairwise test of differentiation in FSTAT, and by exact test (Raymond and Rousset 1995) in TFPGA, all with 15.000 permutations. The estimated θ was chosen since it out-performs other FST analogues in detecting recently established reproductive isolation , and marker neutrality was addressed by jackknifing FST estimates over loci and populations. Standardized measures of genetic differentiation (θ' ) were calculated using the software RECODEDATA v.0.1 (http://www.bentleydrummer.nl/software), to address FST estimates dependence on the level of genetic variation. Finally, STRUCTURE v.2.3.2  was used to test whether population differentiation was detected in a Bayesian model based algorithm. The admixture locprior model with correlated allele frequencies was run for K = 1 - 10 with a burn in of 200.000 MCMC steps, followed by 500.000 steps in 10 replicates.
Temporal substructure stability was tested for by possible year-class variation, sib-ship analyses, effective size estimates and for isolation by distance (IBD). Allele frequency variation was tested in a hierarchical fashion (AMOVA; among versus within populations), and compared with an AMOVA with segregated year-classes (among versus within year-classes), to test for temporal stability within spawning sites where age was available (P1 - 3 and T1). Analyses were performed in ARLEQUIN v.3.1 , under the infinite alleles model (IAM), standard model, unknown gametic phase and 104 permutations. Year-class variation was further controlled for by testing pairwise differentiation among year-classes, and excluding one-by-one year-class from site wise θ tests (above). KININFOR v.1.0  simulated the power of applied markers and observed allele frequencies to discriminate between unrelated and full-/half-sib individuals. We ran the simulation with a prior Dirichlet distribution of 1 for Δ0, Δ1 and Δ2, with a 0.05 confidence level, a precision level of 0.01 and 106 simulated pairs of genotypes. Family structures within sites, possibly affecting population differentiation estimates, were evaluated in GENALEX v.6.2  by the algorithm of Lynch and Ritland [LRM; 154]. Normality of pairwise relatedness was evaluated in JMP v.8.0 .
Contemporary methods were used for estimating effective population size (N e ) and effective parental population size (N b ) within lakes and spawning sites from genetic data. Estimates of N e using summary statistics and approximate Bayesian statistics were computed in ONESAMP , after 50.000 iterations, using prior N e of 4 - 600 for individual sites and 10 - 2000 for the lake samples. Computations were repeated three times, reporting the median result among tests. Runs with prior values ranging from 2 - 1000 were compared to check for consistency and convergence. Parental N b was similarly estimated based on linkage disequilibrium using LDNE v.1.31 . The model is assuming closed populations, and is eliminating the possible bias on small sample sizes. We used 0.02 as the lowest allele frequency used in the computation to balance bias and precision , including a parametric 95% C.I. in the random mating model. The program BOTTLENECK, v.1.2.02  was used to test for recent bottlenecks using the Wilcoxon test and mode-shift tests with 104 replications, with the two phase mutation model (TPM, including 10% infinite alleles model) and the single step mutation model (SMM).
Presence of IBD, indicating equilibrium conditions, was addressed with a Mantel test of correlation between pairwise θ and linear geographic distance among sites in TFPGA, after 1000 iterations.
Migration and admixture
A combination of individual based maximum likelihood methods and population based Bayesian methods were employed to overcome the potential challenge of detecting gene flow under low differentiation. First generation migrants (F0) were in GENECLASS v.2.0  identified to evaluate migrants influence on genetic differentiation, using the test statistics L_home/L_max and L_home between lakes and spawning sites, respectively, to account for un-sampled populations and maximizing analyses power. Migrants between spawning sites were assigned to the population with highest (> 5%) self-assignment probability. Estimates were based on the assignment criteria of Rannala and Mountain  not assuming genetic equilibrium, and the re-sampling algorithm of Paetkau et al.  after 10.000 simulations and a threshold score of 5%. We applied IMMANC v.5.0 , with 10.000 iterations to verify F0 estimates and to assess the power of our dataset.
Population patterns of migration over the last generation within lakes and spawning sites were assessed in BIMR. The method implements estimation of inbreeding coefficients (F) to allow for departure from HW, and assumes sampling after reproduction, but before migration. Information used is gametic disequilibrium, and estimates are calculated using a Bayesian approach and MCMC technique, including 95% highest posterior density intervals (HPDI). Burn-in and sample size of 105, thinning 50, F-model and default values of pilot runs, priors and incremental values were used, reporting mode estimates from the run with lowest total deviance and acceptance rate of 25 - 45% after 10 replicates to ensure convergence [44, 100, 147]. Bayesian assignment tests were then performed on individual level among lakes and among all sampled spawning sites in GENECLASS (with above settings). The random mid-lake sample was similarly assigned to populations, to test for temporal continuous distribution.
Patterns of recent admixed ancestry were finally evaluated in BIMR to assess hybrid presence within spawning sites with above settings on a population level. Hybrid ancestry was only considered on population level, including HPDI estimates, as hybrid detection under admixture may be hampered by restricted power. IMMANC was used to assess F1 ancestral origin on an individual level. Multiplicity correction procedure of Benjamini and Hochberg was used, balancing the risk of Type 1 and Type 2 errors at α = 5% [FDR correction; 45].
Forward modeling of population differentiation was performed with EASYPOP v.2.0.1.  to evaluate if genetic drift could cause observed differentiation patterns alone, and for power evaluation. Simulations were performed with variation in gene flow (0.01 - 0.99), and N e (30 - 90), for 5 populations in 20 generations following panmixia (m = 1.0) with maximum variation and 15 allelic states, assuming random mating, equal sex ratio and an island migration model. Observed pairwise population differentiation from our sample was compared with matrix scenarios of migration and N e indicating similar FST estimates (above) from simulations. Ten independent replicates were obtained for the most likely model to address simulation variation. Similar backward modeling was run in SIMCOAL2  to test the likelihood of the drift model.
Fisher RA: The genetical theory of natural selection. 1930, Oxford: Oxford University Press
Dobzhansky T: Genetics and the origin of species. 1937, New York: Colombia University Press
Saunders DA, Hobbs RJ, Margules CR: Biological consequences of ecosystem fragmentation: a review. Conserv Biol. 1991, 5: 18-32. 10.1111/j.1523-1739.1991.tb00384.x.
Blais J, Plenderleith M, Rico C, Taylor MI, Seehausen O, Van Oosterhout C, Turner GF: Assortative mating among Lake Malawi cichlid fish populations is not simply predictable from male nuptial colour. BMC Evol Biol. 2009, 9:
Johnson JR, Fitzpatrick BM, Shaffer HB: Retention of low-fitness genotypes over six decades of admixture between native and introduced tiger salamanders. BMC Evol Biol. 2010, 10:
Dynesius M, Nilsson C: Fragmentation and flow regulation of river systems in the northern 3rd of the world. Science. 1994, 266: 753-762. 10.1126/science.266.5186.753.
Nilsson C, Reidy CA, Dynesius M, Revenga C: Fragmentation and flow regulation of the world's large river systems. Science. 2005, 308: 405-408. 10.1126/science.1107887.
Beier P, Noss RF: Do habitat corridors provide connectivity?. Conserv Biol. 1998, 12: 1241-1252. 10.1046/j.1523-1739.1998.98036.x.
Schilt CR: Developing fish passage and protection at hydropower dams. Appl Anim Behav Sci. 2007, 104: 295-325. 10.1016/j.applanim.2006.09.004.
World Commissions of Dams: Dams and development A new framework for decisions making. 2000, London: Earthscan Publications Ltd
Marmulla G: Dams, fish and fisheries Opportunities, challenges and conflict resolution. 2001, Rome: FAO Fisheries Technical Paper
Fullerton AH, Burnett KM, Steel EA, Flitcroft RL, Pess GR, Feist BE, Torgersen CE, Miller DJ, Sanderson BL: Hydrological connectivity for riverine fish: measurement challenges and research opportunities. Freshw Biol. 2010, 55: 2215-2237.
Roscoe DW, Hinch SG: Effectiveness monitoring of fish passage facilities: historical trends, geographic patterns and future directions. Fish Fish. 2010, 11: 12-33.
Kemp PS, O'Hanley JR: Procedures for evaluating and prioritising the removal of fish passage barriers: a synthesis. Fish Manag Ecol. 2010, 17: 297-322.
Taggart JB, McLaren IS, Hay DW, Webb JH, Youngson AF: Spawning success in Atlantic salmon (Salmo salar L.): a long-term DNA profiling-based study conducted in a natural stream. Mol Ecol. 2001, 10: 1047-1060. 10.1046/j.1365-294X.2001.01254.x.
Gerlier M, Roche P: A radio telemetry study of the migration of Atlantic salmon (Salmo salar L.) and sea trout (Salmo trutta trutta L.) in the upper Rhine. Hydrobiologia. 1998, 372: 283-293.
Thorstad EB, Okland F, Aarestrup K, Heggberget TG: Factors affecting the within-river spawning migration of Atlantic salmon, with emphasis on human impacts. Rev Fish Biol Fish. 2008, 18: 345-371. 10.1007/s11160-007-9076-4.
Knutsen H, Olsen EM, Jorde PE, Espeland SH, Andre C, Stenseth NC: Are low but statistically significant levels of genetic differentiation in marine fishes 'biologically meaningful'? A case study of coastal Atlantic cod. Mol Ecol. 2011, 20: 768-783. 10.1111/j.1365-294X.2010.04979.x.
Nielsen EE, Hemmer-Hansen J, Poulsen NA, Loeschcke V, Moen T, Johansen T, Mittelholzer C, Taranger GL, Ogden R, Carvalho GR: Genomic signatures of local directional selection in a high gene flow marine organism; the Atlantic cod (Gadus morhua). BMC Evol Biol. 2009, 9:
Waples RS: Separating the wheat from the chaff: Patterns of genetic differentiation in high gene flow species. J Hered. 1998, 89: 438-450. 10.1093/jhered/89.5.438.
Hedrick PW: Perspective: Highly variable loci and their interpretation in evolution and conservation. Evolution. 1999, 53: 313-318. 10.2307/2640768.
Waples RS, Gaggiotti O: What is a population? An empirical evaluation of some genetic methods for identifying the number of gene pools and their degree of connectivity. Mol Ecol. 2006, 15: 1419-1439. 10.1111/j.1365-294X.2006.02890.x.
Strange RM, Stepien CA: Genetic divergence and connectivity among river and reef spawning groups of walleye (Sander vitreus vitreus) in Lake Erie. Can J Fish Aquat Sci. 2007, 64: 437-448. 10.1139/f07-022.
Massicotte R, Magnan P, Angers B: Intralacustrine site fidelity and nonrandom mating in the littoral-spawning northern redbelly dace (Phoxinus eos). Can J Fish Aquat Sci. 2008, 65: 2016-2025. 10.1139/F08-114.
Seehausen O, Terai Y, Magalhaes IS, Carleton KL, Mrosso HDJ, Miyagi R, van der Sluijs I, Schneider MV, Maan ME, Tachida H, Imai H, Okada N: Speciation through sensory drive in cichlid fish. Nature. 2008, 455: 620-626. 10.1038/nature07285.
Bergek S, Bjorklund M: Cryptic barriers to dispersal within a lake allow genetic differentiation of Eurasian perch. Evolution. 2007, 61: 2035-2041. 10.1111/j.1558-5646.2007.00163.x.
Heggenes J, Roed KH, Jorde PE, Brabrand A: Dynamic micro-geographic and temporal genetic diversity in vertebrates: the case of lake-spawning populations of brown trout (Salmo trutta). Mol Ecol. 2009, 18: 1100-1111. 10.1111/j.1365-294X.2009.04101.x.
Pavey SA, Nielsen JL, Hamon TR: Recent ecological divergence despite migration in sockeye salmon (Oncorhynchus nerka). Evolution. 2010, 64: 1773-1783. 10.1111/j.1558-5646.2009.00927.x.
Hendry AP, Kinnison MT: Perspective: The pace of modern life: Measuring rates of contemporary microevolution. Evolution. 1999, 53: 1637-1653. 10.2307/2640428.
Reznick DN, Ghalambor CK: The population ecology of contemporary adaptations: what empirical studies reveal about the conditions that promote adaptive evolution. Genetica. 2001, 112: 183-198.
Hendry AP, Nosil P, Rieseberg LH: The speed of ecological speciation. Funct Ecol. 2007, 21: 455-464. 10.1111/j.1365-2435.2007.01240.x.
Hendry AP: Ecological speciation! Or the lack thereof?. Can J Fish Aquat Sci. 2009, 66: 1383-1398. 10.1139/F09-074.
Rolshausen G, Segelbacher G, Hobson KA, Schaefer HM: Contemporary evolution of reproductive isolation and phenotypic divergence in sympatry along a migratory divide. Current Biology. 2009, 19: 2097-2101. 10.1016/j.cub.2009.10.061.
Carroll SP: Facing change: forms and foundations of contemporary adaptation to biotic invasions. Mol Ecol. 2008, 17: 361-372. 10.1111/j.1365-294X.2007.03484.x.
Herborg LM, Weetman D, Van Oosterhout C, Hanfling B: Genetic population structure and contemporary dispersal patterns of a recent European invader, the Chinese mitten crab, Eriocheir sinensis. Mol Ecol. 2007, 16: 231-242.
Dupont L, Jollivet D, Viard F: High genetic diversity and ephemeral drift effects in a successful introduced mollusc (Crepidula fornicata : Gastropoda). Mar Ecol Prog Ser. 2003, 253: 183-195.
Dupont L, Viard F, Dowell MJ, Wood C, Bishop JDD: Fine- and regional-scale genetic structure of the exotic ascidian Styela clava (Tunicata) in southwest England, 50 years after its introduction. Mol Ecol. 2009, 18: 442-453. 10.1111/j.1365-294X.2008.04045.x.
Corander J, Siren J, Arjas E: Bayesian spatial modeling of genetic population structure. Comput Stat. 2008, 23: 111-129. 10.1007/s00180-007-0072-x.
Hubisz MJ, Falush D, Stephens M, Pritchard JK: Inferring weak population structure with the assistance of sample group information. Mol Ecol Resour. 2009, 9: 1322-1332. 10.1111/j.1755-0998.2009.02591.x.
Chen C, Durand E, Forbes F, Francois O: Bayesian clustering algorithms ascertaining spatial population structure: a new computer program and a comparison study. Mol Ecol Notes. 2007, 7: 747-756. 10.1111/j.1471-8286.2007.01769.x.
Stabell OB: Homing and olfaction in salmonids: A critical review with special reference to the Atlantic salmon. Biol Rev Camb Philos Soc. 1984, 59: 333-388. 10.1111/j.1469-185X.1984.tb00709.x.
Wallraff HG: Navigation by homing pigeons: updated perspective. Ethol Ecol Evol. 2001, 13: 1-48. 10.1080/08927014.2001.9522786.
de Leaniz CG, Fleming IA, Einum S, Verspoor E, Jordan WC, Consuegra S, ubin-Horth N, Lajus D, Letcher BH, Youngson AF, Webb JH, Vøllestad LA, Villanueva B, Ferguson A, Quinn TP: A critical review of adaptive genetic variation in Atlantic salmon: implications for conservation. Biol Rev. 2007, 82: 173-211.
Faubet P, Gaggiotti OE: A new Bayesian method to identify the environmental factors that influence recent migration. Genetics. 2008, 178: 1491-1504. 10.1534/genetics.107.082560.
Benjamini Y, Hochberg Y: Controlling the False Discovery Rate: A practical and powerful approach to multiple testing. J R Stat Soc Series B Methodol. 1995, 57: 289-300.
Aass P: Rødingvandringar i Pålsbufjord och Tunhovdfjord. 1965, Drottningholm
Aass P: Management and utilization of Arctic charr in Norwegian hydroelectric reservoirs. Biology of the Arctic charr Proceedings of the International Symposium on Arctic charr. Edited by: Johnson L, Burns B. 1984, Winnipeg: Univ. Manitoba Press, 277-291.
Northcote TG: Migratory strategies and production in freshwater fishes. Ecology of Freshwater Fish Production. Edited by: Gerking SD. 1978, Oxford: Blackwell Scientific, 326-359.
Arnekleiv JV, Ronning L: Migratory patterns and return to the catch site of adult brown trout (Salmo trutta L.) in a regulated river. River Research and Applications. 2004, 20: 929-942. 10.1002/rra.799.
Rivinoja P, McKinnell S, Lundqvist H: Hindrances to upstream migration of Atlantic salmon (Salmo salar) in a northern Swedish river caused by a hydroelectric power-station. Regulated Rivers-Research & Management. 2001, 17: 101-115. 10.1002/rrr.607.
Thorstad EB, Okland F, Kroglund F, Jepsen N: Upstream migration of Atlantic salmon at a power station on the River Nidelva, Southern Norway. Fish Manag Ecol. 2003, 10: 139-146. 10.1046/j.1365-2400.2003.00335.x.
Latch EK, Dharmarajan G, Glaubitz JC, Rhodes OE: Relative performance of Bayesian clustering software for inferring population substructure and individual assignment at low levels of population differentiation. Conserv Genet. 2006, 7: 295-302. 10.1007/s10592-005-9098-1.
Lowe WH, Allendorf FW: What can genetics tell us about population connectivity?. Mol Ecol. 2010, 19: 3038-3051. 10.1111/j.1365-294X.2010.04688.x.
Whitlock MC, McCauley DE: Indirect measures of gene flow and migration: F-ST not equal 1/(4 Nm+1). Heredity. 1999, 82: 117-125. 10.1038/sj.hdy.6884960.
Mallet J, Meyer A, Nosil P, Feder JL: Space, sympatry and speciation. J Evol Biol. 2009, 22: 2332-2341. 10.1111/j.1420-9101.2009.01816.x.
Hendry AP, Castric V, Kinnison MT, Quinn TP: The evolution of philopatry and dispersal: homing versus straying in salmonids. Evolution illuminated: Salmon and their relatives. Edited by: Hendry AP, Stearns SC. 2004, Oxford: Oxford University Press, 52-91.
Jost L: G(ST) and its relatives do not measure differentiation. Mol Ecol. 2008, 17: 4015-4026. 10.1111/j.1365-294X.2008.03887.x.
Gomez-Uchida D, Knight TW, Ruzzante DE: Interaction of landscape and life history attributes on genetic diversity, neutral divergence and gene flow in a pristine community of salmonids. Mol Ecol. 2009, 18: 4854-4869. 10.1111/j.1365-294X.2009.04409.x.
Petit E, Balloux F, Goudet J: Sex-biased dispersal in a migratory bat: A characterization using sex-specific demographic parameters. Evolution. 2001, 55: 635-640. 10.1554/0014-3820(2001)055[0635:SBDIAM]2.0.CO;2.
Hauser L, Carvalho GR: Paradigm shifts in marine fisheries genetics: ugly hypotheses slain by beautiful facts. Fish Fish. 2008, 9: 333-362.
Bernatchez L, Rhydderch JG, Kircheis FW: Microsatellite gene diversity analysis in landlocked Arctic char from Maine. Trans Am Fish Soc. 2002, 131: 1106-1118. 10.1577/1548-8659(2002)131<1106:MGDAIL>2.0.CO;2.
Dynes J, Magnan P, Bernatchez L, Rodriguez MA: Genetic and morphological variation between two forms of lacustrine brook charr. J Fish Biol. 1999, 54: 955-972. 10.1111/j.1095-8649.1999.tb00850.x.
Nei M, Maruyama T, Chakraborty R: Bottleneck effect and genetic variability in populations. Evolution. 1975, 29: 1-10. 10.2307/2407137.
Dlugosch KM, Parker IM: Founding events in species invasions: genetic variation, adaptive evolution, and the role of multiple introductions. Mol Ecol. 2008, 17: 431-449. 10.1111/j.1365-294X.2007.03538.x.
Roman J, Darling JA: Paradox lost: genetic diversity and the success of aquatic invasions. Trends Ecol Evol. 2007, 22: 454-464. 10.1016/j.tree.2007.07.002.
Tallmon DA, Koyuk A, Luikart G, Beaumont MA: ONeSAMP: a program to estimate effective population size using approximate Bayesian computation. Mol Ecol Resour. 2008, 8: 299-301. 10.1111/j.1471-8286.2007.01997.x.
Waples RS, Do C: LDNE: a program for estimating effective population size from data on linkage disequilibrium. Mol Ecol Resour. 2008, 8: 753-756. 10.1111/j.1755-0998.2007.02061.x.
Brabrand Å, Bremnes T, Saltveit SJ, Aass P, Wollebæk J, Heggenes J, Røed K: Fiskeribiologiske undersøkelser i Pålsbufjorden Hovedrapport. 2008, Oslo
Shrimpton JM, Heath DD: Census vs. effective population size in chinook salmon: large- and small-scale environmental perturbation effects. Mol Ecol. 2003, 12: 2571-2583. 10.1046/j.1365-294X.2003.01932.x.
Frankham R: Effective population size adult population size ratios in wildlife - a review. Genet Res. 1995, 66: 95-107. 10.1017/S0016672300034455.
Levins R: Some demographic and genetic consequences of environmental heterogeneity for biological control. Bull Entomol Soc. 1969, 15: 237-240.
Moyle PB, Light T: Biological invasions of fresh water: Empirical rules and assembly theory. Biol Conserv. 1996, 78: 149-161. 10.1016/0006-3207(96)00024-9.
Schluter D: Ecological speciation in postglacial fishes. Philos Trans R Soc Lond, B, Biol Sci. 1996, 351: 807-814. 10.1098/rstb.1996.0075.
Rundle HD, Nosil P: Ecological speciation. Ecol Lett. 2005, 8: 336-352. 10.1111/j.1461-0248.2004.00715.x.
Sakai AK, Allendorf FW, Holt JS, Lodge DM, Molofsky J, With KA, Baughman S, Cabin RJ, Cohen JE, Ellstrand NC, McCauley DE, O'Neil P, Parker IM, Thompson JN, Weller SG: The population biology of invasive species. Annu Rev Ecol Syst. 2001, 32: 305-332. 10.1146/annurev.ecolsys.32.081501.114037.
Huitfeldt-Kaas H: Der Einfluss der Gewasserregelungen auf den Fischbestand in Binnenseen. 1935, Oslo
Hendry AP: Selection against migrants contributes to the rapid evolution of ecologically dependent reproductive isolation. Evol Ecol Res. 2004, 6: 1219-1236.
Templeton AR: The reality and importance of founder speciation in evolution. Bioessays. 2008, 30: 470-479. 10.1002/bies.20745.
Conover DO, Munch SB: Sustaining fisheries yields over evolutionary time scales. Science. 2002, 297: 94-96. 10.1126/science.1074085.
Stockwell CA, Hendry AP, Kinnison MT: Contemporary evolution meets conservation biology. Trends Ecol Evol. 2003, 18: 94-101. 10.1016/S0169-5347(02)00044-7.
Hendry AP: Adaptive divergence and the evolution of reproductive isolation in the wild: an empirical demonstration using introduced sockeye salmon. Genetica. 2001, 112: 515-534.
Kinnison MT, Bentzen P, Unwin MJ, Quinn TP: Reconstructing recent divergence: evaluating nonequilibrium population structure in New Zealand chinook salmon. Mol Ecol. 2002, 11: 739-754. 10.1046/j.1365-294X.2002.01477.x.
Taylor EB: A review of local adaptation in Salmonidae, with particular reference to Pacific and Atlantic salmon. Aquaculture. 1991, 98: 185-207. 10.1016/0044-8486(91)90383-I.
Koskinen MT, Haugen TO, Primmer CR: Contemporary fisherian life-history evolution in small salmonid populations. Nature. 2002, 419: 826-830. 10.1038/nature01029.
Ayllon F, Davaine P, Beall E, Garcia-Vazquez E: Dispersal and rapid evolution in brown trout colonizing virgin Subantarctic ecosystems. J Evol Biol. 2006, 19: 1352-1358. 10.1111/j.1420-9101.2005.01075.x.
Jonsson B, Jonsson N: Polymorphism and speciation in Arctic charr. J Fish Biol. 2001, 58: 605-638. 10.1111/j.1095-8649.2001.tb00518.x.
Klemetsen A, Amundsen PA, Dempson JB, Jonsson B, Jonsson N, O'Connell MF, 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-59. 10.1034/j.1600-0633.2003.00010.x.
Price TD, Qvarnstrom A, Irwin DE: The role of phenotypic plasticity in driving genetic evolution. Proc R Soc Lond, B, Biol Sci. 2003, 270: 1433-1440. 10.1098/rspb.2003.2372.
Gilchrist GW, Lee CE: All stressed out and nowhere to go: does evolvability limit adaptation in invasive species?. Genetica. 2007, 129: 127-132. 10.1007/s10709-006-9009-5.
Ghalambor CK, Mckay JK, Carroll SP, Reznick DN: Adaptive versus non-adaptive phenotypic plasticity and the potential for contemporary adaptation in new environments. Funct Ecol. 2007, 21: 394-407. 10.1111/j.1365-2435.2007.01283.x.
Hatchwell BJ: Cryptic kin selection: Kin structure in vertebrate populations and opportunities for kin-directed cooperation. Ethology. 2010, 116: 203-216. 10.1111/j.1439-0310.2009.01732.x.
Puebla O: Ecological speciation in marine v. freshwater fishes. J Fish Biol. 2009, 75: 960-996. 10.1111/j.1095-8649.2009.02358.x.
Berner D, Grandchamp AC, Hendry AP: Variable progress toward ecological speciation in parapatry: Stickleback across eight lake-stream transitions. Evolution. 2009, 63: 1740-1753. 10.1111/j.1558-5646.2009.00665.x.
Labonne J, Hendry AP: Natural and sexual selection giveth and taketh away reproductive barriers: Models of population divergence in Guppies. Am Nat. 2010, 176: 26-39. 10.1086/652992.
Nosil P, Funk DJ, Ortiz-Barrientos D: Divergent selection and heterogeneous genomic divergence. Mol Ecol. 2009, 18: 375-402. 10.1111/j.1365-294X.2008.03946.x.
Thibert-Plante X, Hendry AP: Five questions on ecological speciation addressed with individual-based simulations. J Evol Biol. 2009, 22: 109-123. 10.1111/j.1420-9101.2008.01627.x.
Thibert-Plante X, Hendry AP: When can ecological speciation be detected with neutral loci?. Mol Ecol. 2010, 19: 2301-2314. 10.1111/j.1365-294X.2010.04641.x.
Hutchison DW, Templeton AR: Correlation of pairwise genetic and geographic distance measures: Inferring the relative influences of gene flow and drift on the distribution of genetic variability. Evolution. 1999, 53: 1898-1914. 10.2307/2640449.
Whiteley AR, Spruell P, Rieman BE, Allendorf FW: Fine-scale genetic structure of bull trout at the southern limit of their distribution. Trans Am Fish Soc. 2006, 135: 1238-1253. 10.1577/T05-166.1.
Faubet P, Waples RS, Gaggiotti OE: Evaluating the performance of a multilocus Bayesian method for the estimation of migration rates. Mol Ecol. 2007, 16: 1149-1166. 10.1111/j.1365-294X.2007.03218.x.
Wilson GA, Rannala B: Bayesian inference of recent migration rates using multilocus genotypes. Genetics. 2003, 163: 1177-1191.
Griffiths SW: Learned recognition of conspecifics by fishes. Fish Fish. 2003, 4: 256-268.
Olsen KH, Grahn M, Lohm J, Langefors A: MHC and kin discrimination in juvenile Arctic charr, Salvelinus alpinus L. Anim Behav. 1998, 56: 319-327.
Svenning MA, Grotnes P: Stationarity and homing ability of landlocked Arctic charr. Nord J Freshw Resour. 1991, 66: 36-43.
Fraser DJ, Duchesne P, Bernatchez L: Migratory charr schools exhibit population and kin associations beyond juvenile stages. Mol Ecol. 2005, 14: 3133-3146. 10.1111/j.1365-294X.2005.02657.x.
Nordeng H: Char ecology. Natal homing in sympatric populations of anadromous Arctic char Salvelinus alpinus L.: Roles of pheromone recognition. Ecol Freshw Fish. 2009, 18: 41-51. 10.1111/j.1600-0633.2008.00320.x.
Adams CE, Fraser D, Wilson AJ, Alexander G, Ferguson MM, Skulason S: Patterns of phenotypic and genetic variability show hidden diversity in Scottish Arctic charr. Ecol Freshw Fish. 2007, 16: 78-86. 10.1111/j.1600-0633.2006.00182.x.
Adams CE, Hamilton DJ, Mccarthy I, Wilson AJ, Grant A, Alexander G, Waldron S, Snorasson SS, Ferguson MM, Skulason S: Does breeding site fidelity drive phenotypic and genetic sub-structuring of a population of Arctic charr?. Evol Ecol. 2006, 20: 11-26. 10.1007/s10682-005-2489-4.
Bolnick DI, Nosil P: Natural selection in populations subject to a migration load. Evolution. 2007, 61: 2229-2243. 10.1111/j.1558-5646.2007.00179.x.
Barton NH, Hewitt GM: Analysis of hybrid zones. Annu Rev Ecol Syst. 1985, 16: 113-148. 10.1146/annurev.es.16.110185.000553.
Nosil P, Vines TH, Funk DJ: Perspective: Reproductive isolation caused by natural selection against immigrants from divergent habitats. Evolution. 2005, 59: 705-719.
Behrmann-Godel J, Gerlach G: First evidence for postzygotic reproductive isolation between two populations of Eurasian perch (Perca fluviatilis L.) within Lake Constance. Front Zool. 2008, 5:
Edmands S: Between a rock and a hard place: evaluating the relative risks of inbreeding and outbreeding for conservation and management. Mol Ecol. 2007, 16: 463-475.
Brown GE, Brown JA: Kin discrimination in salmonids. Rev Fish Biol Fish. 1996, 6: 201-219.
Ward AJW, Webster MM, Magurran AE, Currie S, Krause J: Species and population differences in social recognition between fishes: a role for ecology?. Behav Ecol. 2009, 20: 511-516. 10.1093/beheco/arp025.
Gerlach G, Hodgins-Davis A, Avolio C, Schunter C: Kin recognition in zebrafish: a 24-hour window for olfactory imprinting. Proc R Soc, B, Biol Sci. 2008, 275: 2165-2170. 10.1098/rspb.2008.0647.
Mehlis M, Bakker TCM, Frommen JG: Smells like sib spirit: kin recognition in three-spined sticklebacks (Gasterosteus aculeatus) is mediated by olfactory cues. Animal Cognition. 2008, 11: 643-650. 10.1007/s10071-008-0154-3.
Naisbit RE, Jiggins CD, Mallet J: Disruptive sexual selection against hybrids contributes to speciation between Heliconius cydno and Heliconius melpomene. Proc R Soc Lond, B, Biol Sci. 2001, 268: 1849-1854. 10.1098/rspb.2001.1753.
van der Sluijs I, Van Dooren TJM, Hofker KD, van Alphen JJM, Stelkens RB, Seehausen O: Female mating preference functions predict sexual selection against hybrids between sibling species of cichlid fish. Philos Trans R Soc, B, Biol Sci. 2008, 363: 2871-2877. 10.1098/rstb.2008.0045.
Nei M, Maruyama T, Wu CI: Models of evolution of reproductive isolation. Genetics. 1983, 103: 557-579.
Ward AJW, Webster MM, Hart PJB: Social recognition in wild fish populations. Proc R Soc, B, Biol Sci. 2007, 274: 1071-1077. 10.1098/rspb.2006.0231.
Olsen KH, Grahn M, Lohm J: Influence of MHC on sibling discrimination in Arctic char, Salvelinus alpinus (L.). J Chem Ecol. 2002, 28: 783-795. 10.1023/A:1015240810676.
Courtenay SC, Quinn TP, Dupuis HMC, Groot C, Larkin PA: Factors affecting the recognition of population-specific odours by juvenile coho salmon. J Fish Biol. 1997, 50: 1042-1060.
Kalinowski ST: The computer program STRUCTURE does not reliably identify the main genetic clusters within species: simulations and implications for human population structure. Heredity. 2011, 106: 625-632. 10.1038/hdy.2010.95.
Behnke RJ: Systematics of salmonid fishes of recently glaciated lakes. J Fish Res Board Can. 1972, 29: 639-671. 10.1139/f72-112.
Primmer CR, Aho T, Piironen J, Estoup A, Cornuet JM, Ranta E: Microsatellite analysis of hatchery stocks and natural populations of Arctic charr, Salvelinus alpinus, from the Nordic region: implications for conservation. Hereditas. 1999, 130: 277-289.
Brunner PC, Douglas MR, Osinov A, Wilson CC, Bernatchez L: Holarctic phylogeography of Arctic charr (Salvelinus alpinus L.) inferred from mitochondrial DNA sequences. Evolution. 2001, 55: 573-586. 10.1554/0014-3820(2001)055[0573:HPOACS]2.0.CO;2.
Wilson AJ, Gislason D, Skulason S, Snorrason SS, Adams CE, Alexander G, Danzmann RG, Ferguson MM: Population genetic structure of Arctic Charr, Salvelinus alpinus from northwest Europe on large and small spatial scales. Mol Ecol. 2004, 13: 1129-1142. 10.1111/j.1365-294X.2004.02149.x.
Johnson L: The Arctic charr, Salvelinus alpinus. Charrs, Salmonid fishes of the genus Salvelinus. Edited by: Balon EK. 1980, The Hague: Dr. W. Junk Publishers, 15-98.
Huitfeldt-Kaas H: Ferskvandsfiskens utbredelse og indvandring i Norge med et tillæg om krebsen. 1918, Kristiania: Centraltrykkeriet
Indrelid S: De første bosetterne. Hardangervidda. Edited by: Barth EK. 1985, Oslo: Luther forlag, 97-111.
Aagaard B: Utsettelse og gjenfangst av røie i Bredevandet, Hol, Buskeruds amt. Norsk Fiskeritidende. 1915, 34: 225-226.
Aass P: The winter migration of charr, Salvelinus alpinus L, in the hydroelectric reservoirs Tunhovdfjord and Pålsbufjord, Norway. 1970, Drottningholm
Aass P: Some effects of lake impoundments on salmonids in Norwegian hydroelectric reservoirs. Dissertation. PhD Thesis. 1973, University of Uppsala, Faculty of Science
Behlke CE, Kane DL, McLean RF, Travis MD: Fundamentals of culvert design for passage of weak-swimming fish. 1991, Alaska
Jones DR, Kiceniuk JW, Bamford OS: Evaluation of swimming performance of several fish species from Mackenzie river. J Fish Res Board Can. 1974, 31: 1641-1647. 10.1139/f74-206.
Tritton DJ: Physical fluid dynamics. 1988, Oxford: Clarendon Press
Presa P, Guyomard R: Conservation of microsatellites in three species of salmonids. J Fish Biol. 1996, 49: 1326-1329.
Dehaan PW, Ardren WR: Characterization of 20 highly variable tetranucleotide microsatellite loci for bull trout (Salvelinus confluentus) and cross-amplification in other Salvelinus species. Mol Ecol Notes. 2005, 5: 582-585. 10.1111/j.1471-8286.2005.00997.x.
Angers B, Bernatchez L, Angers A, Desgroseillers L: Specific microsatellite loci for brook charr reveal strong population subdivision on a microgeographic scale. J Fish Biol. 1995, 47: 177-185.
Crane PA, Lewis CJ, Kretschmer EJ, Miller SJ, Spearman WJ, DeCicco AL, Lisac MJ, Wenburg JK: Characterization and inheritance of seven microsatellite loci from Dolly Varden, Salvelinus malma, and cross-species amplification in Arctic char, S alpinus. Conserv Genet. 2004, 5: 737-741. 10.1007/s10592-004-1853-1.
OReilly PT, Hamilton LC, McConnell SK, Wright JM: Rapid analysis of genetic variation in Atlantic salmon (Salmo salar) by PCR multiplexing of dinucleotide and tetranucleotide microsatellites. Can J Fish Aquat Sci. 1996, 53: 2292-2298.
Estoup A, Presa P, Krieg F, Vaiman D, Guyomard R: (Ct)(N) and (Gt)(N) microsatellites: a new class of genetic markers for Salmo trutta L. (brown trout). Heredity. 1993, 71: 488-496. 10.1038/hdy.1993.167.
Van Oosterhout C, Hutchinson WF, Wills DPM, Shipley P: MICRO-CHECKER: software for identifying and correcting genotyping errors in microsatellite data. Mol Ecol Notes. 2004, 4: 535-538. 10.1111/j.1471-8286.2004.00684.x.
Miller MP: Tools for population genetic analyses (TFPGA) ver. 1.3. 1997, Department of Biological Sciences, Nortern Arizona University, Flaggstaff, Arizona
Goudet J: FSTAT (Version 1.2): A computer program to calculate F-statistics. J Hered. 1995, 86: 485-486.
Gaggiotti OE, Foll M: Quantifying population structure using the F-model. Mol Ecol Resour. 2010, 10: 821-830. 10.1111/j.1755-0998.2010.02873.x.
Weir BS, Cockerham CC: Estimating F-statistics for the analysis of population structure. Evolution. 1984, 38: 1358-1370. 10.2307/2408641.
Sefc KM, Payne RB, Sorenson MD: Genetic differentiation after founder events: an evaluation of F-ST estimators with empirical and simulated data. Evol Ecol Res. 2007, 9: 21-39.
Meirmans PG: Using the AMOVA framework to estimate a standardized genetic differentiation measure. Evolution. 2006, 60: 2399-2402.
Excoffier L, Laval G, Schneider S: Arlequin ver. 3.0: An integrated software package for population genetics data analyses. Evol Bioinform Online. 2005, 1: 47-50.
Wang JL: Informativeness of genetic markers for pairwise relationship and relatedness inference. Theor Popul Biol. 2006, 70: 300-321. 10.1016/j.tpb.2005.11.003.
Peakall R, Smouse PE: GENALEX 6: genetic analysis in Excel. Population genetic software for teaching and research. Mol Ecol Notes. 2006, 6: 288-295. 10.1111/j.1471-8286.2005.01155.x.
Ritland K: Marker-inferred relatedness as a tool for detecting heritability in nature. Mol Ecol. 2000, 9: 1195-1204. 10.1046/j.1365-294x.2000.00971.x.
JMP: JMP v.8.0. (8.0). 2008, SAS Institute Inc
Waples RS, Do C: Linkage disequilibrium estimates of contemporary N-e using highly variable genetic markers: a largely untapped resource for applied conservation and evolution. Evol Appl. 2010, 3: 244-262. 10.1111/j.1752-4571.2009.00104.x.
Cornuet JM, Luikart G: Description and power analysis of two tests for detecting recent population bottlenecks from allele frequency data. Genetics. 1996, 144: 2001-2014.
Piry S, Alapetite A, Cornuet JM, Paetkau D, Baudouin L, Estoup A: GENECLASS2: A software for genetic assignment and first-generation migrant detection. J Hered. 2004, 95: 536-539. 10.1093/jhered/esh074.
Rannala B, Mountain JL: Detecting immigration by using multilocus genotypes. Proc Natl Acad Sci USA. 1997, 94: 9197-9201. 10.1073/pnas.94.17.9197.
Paetkau D, Slade R, Burden M, Estoup A: Genetic assignment methods for the direct, real-time estimation of migration rate: a simulation-based exploration of accuracy and power. Mol Ecol. 2004, 13: 55-65. 10.1046/j.1365-294X.2004.02008.x.
Balloux F: EASYPOP (version 1.7): A computer program for population genetics simulations. J Hered. 2001, 92: 301-302. 10.1093/jhered/92.3.301.
Laval G, Excoffier L: SIMCOAL 2.0: a program to simulate genomic diversity over large recombining regions in a subdivided population with a complex history. Bioinformatics. 2004, 20: 2485-2487. 10.1093/bioinformatics/bth264.
The present study is financed by the Norwegian Ministry of Education and Research. Freshwater Ecology and Inland Fisheries Laboratory, University of Oslo, is acknowledged for assistance in providing samples and age analyses. For help in the laboratory we are in debt to Liv Midthjell.
All authors; JW, JH and KH conceived the study and participated in the project design. JW conducted the field and laboratory work, and did the genetic analyses. All authors wrote and approved the final manuscript.