- Research article
- Open Access
The effects of Medieval dams on genetic divergence and demographic history in brown trout populations
BMC Evolutionary Biology volume 14, Article number: 122 (2014)
Habitat fragmentation has accelerated within the last century, but may have been ongoing over longer time scales. We analyzed the timing and genetic consequences of fragmentation in two isolated lake-dwelling brown trout populations. They are from the same river system (the Gudenå River, Denmark) and have been isolated from downstream anadromous trout by dams established ca. 600–800 years ago. For reference, we included ten other anadromous populations and two hatchery strains. Based on analysis of 44 microsatellite loci we investigated if the lake populations have been naturally genetically differentiated from anadromous trout for thousands of years, or have diverged recently due to the establishment of dams.
Divergence time estimates were based on 1) Approximate Bayesian Computation and 2) a coalescent-based isolation-with-gene-flow model. Both methods suggested divergence times ca. 600–800 years bp, providing strong evidence for establishment of dams in the Medieval as the factor causing divergence. Bayesian cluster analysis showed influence of stocked trout in several reference populations, but not in the focal lake and anadromous populations. Estimates of effective population size using a linkage disequilibrium method ranged from 244 to > 1,000 in all but one anadromous population, but were lower (153 and 252) in the lake populations.
We show that genetic divergence of lake-dwelling trout in two Danish lakes reflects establishment of water mills and impassable dams ca. 600–800 years ago rather than a natural genetic population structure. Although effective population sizes of the two lake populations are not critically low they may ultimately limit response to selection and thereby future adaptation. Our results demonstrate that populations may have been affected by anthropogenic disturbance over longer time scales than normally assumed.
Nearly all natural ecosystems, and the species that they encompass, experience anthropogenic pressure resulting particularly from habitat destruction, climate change, overharvesting and introduction of exogenous species [1–4]. Land use leading to habitat fragmentation is considered one of the biggest threats to biodiversity , also encompassing short- and long-term negative genetic consequences [5–7]. Fragmentation problems are particularly severe in freshwater systems , where the one-dimensional structure of river systems causes dams, weirs and other human-made obstructions to represent impassable barriers to many fishes and invertebrates. Accordingly, several studies have documented negative effects of dams on freshwater fish populations that depend on migration between different spawning, nursery and foraging habitats in different parts of their life cycle. In some cases this had led to extirpation of entire populations [9–12]. In other cases fragmentation due to dams has been shown to exert negative genetic effects caused by restricted gene flow and declining population sizes [13–17].
Populations within a species are often distributed across habitats that exhibit different environmental conditions. These differences may lead to variation in local selection regimes and potentially local adaptation , the degree of which is determined by the opposing forces of local selection on the one side and gene flow and random genetic drift on the other [19, 20]. Whereas the limiting factors of gene flow and drift under pristine conditions are influenced by natural barriers in the landscape and habitat size and quality [21, 22], anthropogenic habitat fragmentation leading to decreased gene flow and effective population sizes may substantially shift migration-selection-drift equilibria . The outcome for local adaptation is not trivial [24, 25]. Reduced effective population sizes may cause drift to overwhelm directional selection  and limit evolutionary potential . Conversely, reduced gene flow may at the same time increase possibilities for local adaption, as directional selection becomes a stronger evolutionary force than immigration [19, 28], though in the longer term gene flow into small populations may benefit adaptive responses by introducing new variation. In total, when investigating adaptive divergence between populations in human-altered environments the question arises if local adaption reflects historical or more recent anthropogenically modified selection regimes and demographic parameters.
Most salmonid fish species form local, genetically differentiated populations, many of which are locally adapted [29, 30]. This is also the case for the brown trout (Salmo trutta), where several recent studies based on genetic markers and outlier scans [31–33], transcriptomics  and quantitative genetics experiments [35–38] have provided evidence for local adaptation. Specifically, in a common garden experiment a lake-dwelling brown trout population from Lake Hald, Central Jutland, Denmark (Figure 1) showed evidence for being adapted to higher water temperatures during incubation of eggs and larvae as compared to other populations . In a second study, both the Lake Hald and another lake-dwelling trout population from Lake Mossø (Figure 1) showed several outlier loci in microsatellite DNA outlier scans encompassing anadromous, landlocked and hatchery trout populations .
Both Lake Hald and Lake Mossø are part of the same river system, the Gudenå River (Figure 1). There are no natural barriers to migration within the system and historically anadromous trout would have access to all parts of the river. However, over the past ca. 800 years impassable dams have increasingly been established. During the Medieval water mills, often owned and managed by monasteries, were established at many rivers for manufacturing of textiles and grinding of flour. Whereas the first types of water mills built from ca. 1200 may not have completely and permanently blocked the rivers, from ca. 1400 new types of mills became established that involved permanent damming and created impassable barriers . Lake Hald and Lake Mossø became isolated from the headwaters of the Gudenå River system during that time and must be assumed to have been inaccessible to anadromous brown trout from the downstream part of the river system ever since. Hence, the question arises as to how historical anthropogenic fragmentation has affected the genetic structure, demographic parameters and potentially adaptive divergence of these populations. Did the lake trout populations diverge naturally thousands of years ago, potentially as far back in time as the end of the last Glaciation ca. 12,000 years bp? Or is habitat fragmentation due to establishment of dams the major factor shaping the current genetic population structure? What are the genetic consequences of several centuries of reproductive isolation of the populations?
In the present study we analyze variation at 44 microsatellite loci in brown trout populations from Lake Hald, Lake Mossø and anadromous trout from the lower part of the Gudenå River (the Lilleå River tributary). For reference, we additionally include 10 other anadromous populations and two hatchery strains that have been used for stocking in the region. Using individual Bayesian clustering  we first assess if stocking has significantly affected the current Gudenå River populations and genetic structure. Next, using analyses based on Approximate Bayesian Computation (ABC)  and an Isolation-with-migration model  we estimate divergence time between anadromous trout from the lower Gudenå River and the two lake populations. If the populations already were significantly reproductively isolated through natural processes thousands of years ago, we expect this to be reflected in the divergence time estimates. Conversely, if anthropogenic fragmentation is the primary determinant shaping the genetic structure of the populations, we expect estimates of divergence time to coincide with the establishment of dams ca. 500–800 years ago. Finally, based on results from the same analyses along with estimates of effective population size based on a linkage disequilibrium method  we assess the impact of dams on demographic parameters of the populations and discuss the consequences for extant and future adaptive responses.
Lake Hald (in the following abbreviated HAL); see Figure 1) covers an area of 3.3 km2 and is inhabited by brown trout that spawn in the tributaries. The outlet is Non Mølleå River, the Danish name of which indicates that nuns from the nearby Asmild Monastery maintained a water mill at the river. Damming of the outlet and isolation of HAL from the Gudenå River system is assumed to have taken place sometime in the period from 1400 to 1500, although older mills have been constructed even earlier. Stocking of adult hatchery strain trout directly in the lake took place in the 1970s-1980s, but a study based on allozymes suggested little if any genetic impact .
Lake Mossø (in the following abbreviated MOS; see Figure 1) covers an area of 16.9 km2 and is also inhabited by brown trout that spawn in the tributaries. The nearby Øm Monastery established a water mill and permanent dam at Rye Mølle downstream of the lake around 1500. However, remnants have been found of older water mills constructed 100–200 years earlier. There has been limited stocking of hatchery strain trout in some of the tributaries, and transplants from the lower part of the Gudenå River could potentially also have affected the genetic composition.The Lilleå River (see Figure 1; abbreviated LIL) is a major tributary of the lower Gudenå River. Access to the Gudenå River for spawning anadromous trout was significantly reduced after construction of the Tange hydropower plant (see Figure 1) in 1920. The majority of the present anadromous trout population of the Gudenå River now spawn in LIL downstream of the hydropower plant, and we considered trout from this river as representative of the trout population downstream HAL and MOS. The population could potentially be affected by stocked hatchery strain trout, in particular indirectly due to stocking elsewhere in the system.
To investigate if stocked hatchery strain trout genetically affected the populations, which could influence estimation of divergence time and demographic parameters, we included two hatchery strains in the analyses, Vork Hatchery (VOR) and Hårkær Hatchery (HAR). Stocking with hatchery strain trout has not been permitted in Denmark since 2003, but was wide-spread previously. More than 80% of all stocked fish were derived from four strains that share a common history and have been shown to exhibit close genetic relationships . VOR and HAR belong to this group of strains, which was originally founded by wild spawners from the Vejle and Kolding River on the East coast of Jutland.
Finally, to assess possible sources of gene flow into the LIL population and to provide a comparison to estimates of demographic parameters and to populations known to be introgressed by stocked hatchery strain trout, we included samples from 6 anadromous trout populations on the Jutland West coast (Storå River, STO; Skjern River, SKJ; Varde River, VAR; Sneum River, SNE; Kongeå River, KON; Ribe River, RIB), two populations from the Limfjord (Karup River, KAR; Skals River, SKA) and two populations from the Jutland East coast (Villestrup River, VIL; Kolding River, KOL) (see Figure 1). The West coast rivers, particularly SKJ and VAR, have previously been found to be strongly admixed with hatchery strain trout . In contrast, KAR has been heavily stocked with hatchery trout, but nevertheless shows very low admixture . SKA has not been stocked and to our knowledge the same applies to the VIL and KOL populations. However, KOL is known to be a main source for the founding of the quantitatively most important hatchery strains, including VOR and HAR, and could for that reason be expected to show close genetic relationships to these strains.
All samples were collected by electrofishing conducted from 1999–2009 and encompassed multiple cohorts. Electrofishing was conducted by technical staff at the Danish Institute for Fisheries Research (now National Institute of Aquatic Resources, technical University of Denmark), who had all necessary permits from the Danish Ministry of Food, Agriculture and Fisheries for electrofishing and sampling of tissue and followed all required regulations. The fish were anaesthetized and a small piece of adipose fin was removed and stored in 96% ethanol, following which the fish were released. Sampled localities, sample abbreviations, ecotypes (lake-dwelling, anadromous, hatchery trout) and sample sizes are summarized in Table 1.
The study was based on a subset of the data set from a previous publication . This paper analyzed 74 microsatellite markers and used outlier scans and landscape genomics approaches for detecting loci under possible selection. For the present study we omitted loci that were previously found to be outliers in terms of genetic differentiation, showed significant association with environmental parameters, were known to be linked to functional loci, and/or showed distributions of allele sizes that would indicate strong deviation from a stepwise mutation pattern. In total, we analyzed 44 loci, which are listed in Additional file 1: Table S1. For technical details on DNA extraction and molecular analyses we refer to a previous publication .
Analysis of genetic variation and differentiation
Exact tests for deviations from Hardy-Weinberg equilibrium, observed (HO) and expected (He) heterozygosity and allelic richness was quantified with FSTAT 2.9.1 . Genetic differentiation between populations was quantified with an estimator of F ST  and significance tested by permuting alleles (104 times) among populations using MSA 4.05 .
Bayesian cluster analysis
Bayesian clustering implemented in STRUCTURE 2.3.4 [51, 52] was used for estimating the number of populations/groups represented by the sampled individuals (k) and for estimating individual and population-level admixture proportions. The analysis assumed correlated allele frequencies and was furthermore based on the LOCPRIOR model , where information on the sample origin of individuals was used as prior information. For estimating the most likely k we conducted runs assuming k of 1 through 16. Each run consisted of a burn-in of 105 MCMC steps, followed by 2×105 steps. Ten replicates were conducted for each k. We plotted the probability of the data [(P(D)] and the ad hoc statistic ΔK (Evanno et al. 2005) which measures the steepest increase of the probability of k, using STRUCTURE HARVESTER . CLUMPP  was used to find the consensus configuration based on the ten replicates for a given k, implementing the LargeKGreedy algorithm. Finally, membership proportions of individuals and samples to the identified clusters were visualized using DISTRUCT .
Estimation of contemporary effective population size
We estimated the contemporary effective population size of all populations using the linkage disequilibrium (or more precisely gametic phase disequilibrium) method [43, 56]. For this purpose we used the software LDNE  and excluded alleles occurring at frequencies below 0.05. Given that the samples analyzed included multiple age classes, we assume that the estimates are closer to representing Ne (the effective population size per generation) than Nb (the effective population size per cohort).
Estimation of divergence time and demographic parameters
We used Approximate Bayesian Computation (ABC) as implemented in DIYABC 2.0  to assess posterior likelihoods of divergence time between LIL and HAL and between LIL and MOS, along with effective population sizes. It should be noted that DIYABC assumes that no gene flow occurs after populations have split. We assumed a generation time of 3.5 years . We used wide and flat priors of N e for all populations [50; 10,000] and for divergence time, t (30; 3,000 generations, corresponding to 105 and 10,500 years, respectively). We assumed a generalized stepwise mutation model  with a uniform prior distribution of mean mutation rate from 10-4 to 10-3, a prior distribution of individual locus mutation rates from 10-5 to 10-2 following a Gamma distribution with mean determined by the mean mutation rate across loci. The number of repeats per mutational event across loci was assumed to follow a geometric distribution with a uniform prior for the parameter P ranging from 0.1 to 0.6, whereas P for individual loci followed a Gamma distribution with mean determined by mean across loci and prior ranging from 10-2 to 9 × 10-1. A number of summary statistics can be chosen for estimating posterior distributions of parameters, some of which are, however, partly redundant. We followed the approach of a previous study  and conducted three different sets of analyses for each scenario based on combinations of summary statistics found to be useful in previous studies; 1) mean number of alleles across loci within populations, mean expected heterozygosity within populations, mean value of M within populations, F ST between populations and across loci, and index of individual assignment  across populations ; 2) mean number of alleles across loci, mean expected heterozygosity, mean allele size variance within populations within populations, and the same statistics across populations ; 3) F ST between pairs of populations, mean individual assignment likelihoods of individuals collected in one of a pair of populations and assigned to the other , the mean number of alleles per locus between pairs of populations, mean expected heterozygosity between pairs of populations and mean variance of allele size  between pairs of populations .
We pre-evaluated each scenario and set of priors by performing 104 simulations, conducting a Principal Component Analysis (PCA) based on the summary statistics and checking if the observed data set and the cloud of simulated data were congruent. The analyses were based on simulating 106 data sets, and the posterior distribution of parameters was estimated using the logit approach based on the 104 (1%) data sets closest to the observed data. Following estimation of parameters we checked the model by simulating 103 data sets based on the predictive posterior distribution of parameters but using a new set of summary statistics, and conducting a PCA based on the summary statistics of these simulated data sets along with the 104 data sets simulated based on the prior distribution (see above) and the summary statistics from the observed data. These model checks were conducted by switching summary statistics among analyses. Hence, for analyses based on summary statistics set a) we checked the model using summary statistics set b), for the analyses using set b) we checked the model using set c), and for the analyses using set c) we checked the model using set a).
We further estimated divergence time, effective population sizes and gene flow using the Bayesian, coalescence and Markov Chain Monte Carlo (MCMC) based method IMa . The method is based on an isolation-with-gene-flow model, where a single population back in time splits into two, which are subsequently connected by some gene flow. The following parameters are estimated: t, the point in time when the ancestral population split into two, q A , the effective population size in the ancestral population prior to splitting; q 1 and q2, the effective population size of population 1 and 2 after divergence; m 1 and m 2 , the migration rate from population 2 into population 1 and from population 1 into population 2. The parameters are scaled by mutation rate. To provide unscaled parameters we assumed the microsatellite mutation rate estimated using DIYABC (see above). Based on initial trial runs we found that a heated chain approach encompassing 30 chains (parameters g1 = 0.7 and g2 = 0.4), an initial burn-in of 2 × 106 MCMC steps followed by 2 × 106 steps yielded convergence. We assumed a generation time of 3.5 years (see above) and the following upper bounds on scaled priors: q1 = q2 = 10, qA = 100, m1 = m2 = 175 and t = 2. We conducted three replicate runs for each population pair.
Genetic variation and differentiation
Summary statistics for all loci in all populations (expected and observed heterozygosity, allelic richness, tests for Hardy-Weinberg equilibrium) are listed in Additional file 2: Table S2. Pairwise F ST ranged from 0.002 (between the neighbouring populations RIB and KON) to 0.078 between HAL and the hatchery strain HAR. All pairwise F ST estimates were statistically significant except for that between RIB and KON. The two lake populations HAL and MOS were significantly differentiated from all other populations with F ST ranging from 0.025 to 0.062. Differentiation between HAL and the downstream anadromous LIL population within the Gudenå River system was 0.033, whereas F ST between MOS and LIL was 0.037. Finally, F ST between HAL and MOS was 0.053 (Additional file 3: Table S3).
Bayesian cluster analysis
Analysis of the STRUCTURE output showed that ΔK was highest for k = 2, but with additional peaks for k = 4, 7 and 12 (Additional file 4: Figure S1). At k = 2 one cluster was prevalent (>90%) in the two hatchery strains, VOR and HAR, whereas the other cluster was prevalent (>90%) in KAR, SKA, MOS and HAL. The other populations showed significant proportions of both clusters (data not shown). Assuming k = 7 provided the most detailed separation of populations (see below), whereas assuming k = 12 did not provide further biologically meaningful separation (data not shown).
At k = 7, the hatchery strains VOR and HAR were characterized primarily by the “red” cluster in Figure 2, which also occurred at high frequency in the West coast populations, consistent with previous findings of admixture with hatchery strain trout . The cluster was virtually absent from KAR, SKA, VIL, HAL and MOS, but was present at low frequency (0.075) in LIL and high frequency (0.423) in KOL. This could reflect a common population history, as the two hatchery strains were founded by trout from the two East coast populations KOL and the neighbouring Vejle River, or it could represent admixture due to stocking. Comparison of population and individual level cluster membership proportions (Figure 2) showed a distinct pattern for the stocked West coast populations, where some individuals showed strong admixture proportions of the “red” hatchery-specific cluster, whereas other were virtually non-admixed. In contrast, all individuals in LIL and KOL showed approximately equal proportions of the “red” cluster, indicating that a common population history rather than admixture underlies the results. We therefore assume that LIL is not admixed with hatchery trout and is representative of indigenous anadromous Gudenå River trout.HAL and MOS were characterized by two different distinct clusters (“orange” and “yellow”) whereas the two Limfjord populations KAR and SKA were characterized by a “green” cluster (Figure 2). The southern populations on the Jutland West coast (SNE, KON, RIB) showed high proportions of a “dark blue” cluster that was near absent in other populations. Finally, a “light blue” cluster was found in high proportion in the East coast populations VIL, LIL and KOL and to a smaller extent in some of the West coast populations. In total, the results suggest that the most likely sources of gene flow into LIL encompass other East coast populations and that the gene pool of this population, along with MOS and HAL are unlikely to be strongly affected by stocked hatchery strain trout.
Contemporary effective population size
Estimates of contemporary effective population size (N e ) using LDNE generally ranged from ca. 250 to > 1,000 in the anadromous trout populations, the exception being SKA, where the point estimate of Ne was 74 (Table 2). In RIB, N e was too high to be estimated (denoted by ∞, meaning that sampling variance exceeded the signal from drift), but the lower 95% confidence interval was 480. In four anadromous trout populations (including RIB) the upper 95% confidence interval could not be determined. Specifically for LIL, representing anadromous trout from the lower Gudenå River system, N e was 288. A comparable N e estimate was obtained for MOS (252), whereas the estimate for HAL was lower (153).
The BOTTLENECK tests  provided no evidence for recent bottlenecks in neither LIL, HAL nor MOS.
Divergence time and demographic parameters
The analyses using DIYABC yielded point estimates of mean mutation rate ranging from 2.90 × 10-4 to 3.85 × 10-4 (Table 3). Estimates of divergence time for LIL-HAL ranged from 602 to 742 years, assuming different sets of summary statistics, whereas for LIL-MOS it ranged from 599 to 641 years, thus providing strong support for divergence caused by the establishment of water mills and associated impassable dams during the Medieval. Estimates of N e generated by DIYABC were substantially higher than those obtained using LDNE. Hence, point estimates of N e in LIL ranged from 5790 to 7470 across different analyses, in HAL it ranged from 1980–3470 and in MOS from 3420–4070 (Table 3).
For the IMa analyses we assumed a mutation rate of 3.0 × 10-4 corresponding to the estimates provided by DIYABC. The results showed highly consistent outcomes across the three replicate runs for LIL-HAL and LIL-MOS, both concerning estimates of divergence time (Figure 3) and effective population size and migration rate (Table 4). The mean divergence time across the three runs was 789 years for LIL-HAL and 696 years for LIL-MOS (Table 4), corresponding well to the estimates obtained using DIYABC. There was virtually no statistical support for divergence having occurred naturally further back in time (>1,500 years; Figure 3). Estimates of N e of the two isolated lake populations MOS and HAL were lower than for the downstream LIL population (Table 4). Average modes across three runs yielded estimates of 304 in MOS and 369 in HAL, whereas average mode of N e in LIL across six runs was 1585. Hence, the estimates were higher than those obtained using LDNE, but considerably lower than suggested by DIYABC. Estimates of N e of the ancestral population prior to divergence were high, generally on the order of 40,000 (Table 4). Gene flow estimates from LIL and upstream above dams to HAL and MOS were, as expected, virtually zero. Also, there was limited support for gene flow downstream from HAL to LIL, whereas low but non-negligible gene flow from MOS to LIL was indicated (Table 4).
Divergence time estimates obtained using two different methods based on Approximate Bayesian Computation (DIYABC) and an Isolation-with-migration model (IMa) led us to reject the hypothesis that divergence of HAL and MOS from the anadromous LIL population is natural and occurred several thousand years ago. Instead, divergence time coincided remarkably well with establishment of water mills and impassable dams in the Medieval from ca. 1200 to 1500. There was less congruence regarding the demographic impact of the establishment of dams, as DIYABC consistently provided higher estimates of effective population sizes than did IMa and estimates based on the linkage disequilibrium method (LDNE).
Bayesian cluster analysis showed that the two lake populations HAL and MOS were genetically distinct, whereas the anadromous trout populations showed a regional genetic structure, reflecting ongoing gene flow and isolation-by-distance, as documented in other studies [32, 46]. Importantly, these results also showed a minimal genetic contribution by stocked trout to the focal populations LIL, MOS and HAL, which could otherwise complicate estimation of divergence time and demographic parameters.
In the following, we first discuss the reliability of the estimates provided by DIYABC and IMa, and subsequently the impact of the impassable dams on genetic population structure and adaptive divergence.
Reliability of estimates
The conclusion that the two lake populations diverged from anadromous trout due to habitat fragmentation during the Medieval relies on the results using the methods IMa  and DIYABC . In both cases we used wide and flat priors that should not strongly influence posterior distributions. The mutation rate assumed, however, is a factor of uncertainty. For the IMa analyses we assumed a rate of 3.0 × 10-4, estimated using DIYABC. This is quite close to estimates of mutation rate at dinucleotide loci in e.g. humans of 2.73 × 10-4 and 5.56 × 10-4 in common carp (Cyprinus carpio) . Hence, the mutation rate assumed must be considered realistic, and it would require a much lower mutation rate (2.01 × 10-5) to change the estimate of divergence time for LIL-HAL from 789 years to, say, 12,000 years, coinciding with the end of the last Glaciation.
Whereas the congruence of results from different replicated runs of IMa suggests good convergence, a potential problem consists in the simplified model assumed, where LIL exchanges migrants with either HAL or MOS. In reality, both populations could simultaneously contribute to gene flow, but even more importantly LIL would be subject to gene flow from other anadromous populations. We are unable to resolve the magnitude of this potential problem, but we note that the results obtained generally make sense. For both lake populations the estimate of gene flow from LIL was virtually zero, as would be expected given the problems of passing the dams upstream. In contrast, there was a signal of downstream gene flow into LIL from MOS. Gene flow in this direction is certainly a realistic possibility, as fish can pass the dams downstream by simply being flushed with the current.
N e estimates of HAL and MOS obtained using IMa were of the same order of magnitude, albeit higher than estimates obtained independently using LDNE. In contrast, estimates of N e in LIL obtained by IMa were more than 5 times times higher than those derived from LDNE. Although recent bottlenecks could explain this result, there was no evidence for this to have taken place. Alternatively, the result could reflect gene flow from other anadromous populations into LIL, thereby inflating historical N e estimates (as in IMa) compared to contemporary estimates (as in LDNE) that pertain to the current or only a few generations back in time. The complexities of estimating N e in populations showing some degree of geographical continuity are increasingly acknowledged , and we tentatively suggest that the N e estimate in LIL based on IMa should be interpreted more broadly as encompassing LIL and neighboring populations contributing to gene flow.
The DIYABC analyses appeared robust towards the choice of summary statistics, as evidenced by similar outcomes of the different analyses. However, unlike IMa it is assumed that there is no gene flow following splitting of populations. In the present case, this is probably not a major problem, as IMa suggested low gene flow between LIL and HAL/MOS, and divergence time estimates were quite similar between methods. Effective population size estimates, however, differed strongly between IMa and DIYABC, with the latter providing 5–10 times higher estimates than IMa. In both cases the estimates are “historical”, i.e. encompassing many generations, although IMa estimates N e from the ancestral population prior to divergence and N e of the separate populations after divergence, the latter thereby giving more weight to the recent past (in the present case the centuries since establishment of dams). Most summary statistics applied in DIYABC, such as expected heterozygosity, F ST and numbers of alleles should primarily reflect historical N e over a very long time scale , with the exception of M which detects recent bottlenecks. Differences of the time scale over which N e is estimated could therefore explain the discrepancies of results. We emphasize the N e estimates obtained by IMa as being most realistic and relevant in a conservation context; they are generally congruent with estimates of contemporary N e obtained using LDNE and with estimates from Danish brown trout populations obtained using temporal methods, that for anadromous populations in Danish rivers generally range from ca. 250 to > 1,000 for large river systems [46, 61].
In total, we find the estimates of divergence time obtained by IMa and DIYABC to be robust, whereas N e estimates differ, probably reflecting the different time scales that they apply to.
Impact on genetic population structure and adaptive divergence
Our results suggest that the two lake populations HAL and MOS were previously part of a coherent anadromous trout population inhabiting the Gudenå River system. The Gudenå River is the largest watershed in Denmark, and the effective population size of the ancestral trout population prior to fragmentation has presumably been high; > 1,000, given the estimates obtained from the KON and RIB populations (see Table 2 and ). Gene flow has occurred to and from other such populations enabling introduction of new adaptive variation and potentially both enhancing and limiting local adaptation, depending on local selection regimes and migration-drift-selection equilibrium [24, 25]. Following the establishment of water mills, HAL and MOS have been disconnected from these dynamics for several centuries.
What is the primary impact of this fragmentation? Contemporary N e was estimated to 153 and 252 for HAL and MOS, respectively. Although these values cannot be considered high, they should also not be cause for immediate conservation concern. Indeed, they are higher than several estimates of N e in undisturbed resident brown trout populations from Sweden [75, 76] that have been isolated naturally presumably for even longer time spans that the MOS and HAL populations, and are in fact comparable to N e in some of the anadromous populations included in the study (Table 2). It should be considered, however, that several of these latter populations included for reference have likely declined recently, as documented by comparison of historical and contemporary samples from the populations .
We suggest that the most important consequence of the dams concerns local adaptation and evolutionary potential. Even if trout ascend the dams downstream and migrate to sea, they will be unable to return to spawn. This could impose strong selection against anadromy, although the response to selection depends on heritability of the trait. There are few studies available that have estimated heritability of anadromy in salmonids. However, a study based on pedigreeing a wild brook trout (Salvelinus fontinalis) population found heritability of life-history tactics (anadromy versus residency) as high as 0.52-0.56 . If this is also the case in brown trout, then the 600–700 years since establishment of dams corresponding to 170–200 generations should have left ample opportunities for selective responses to occur. On the other side, a recent study of trout in HAL documented that 15% of individuals aged between 1 and 3 years and with a length exceeding 12 cm left the lake and would potentially undertake migration to the sea, whereas 40% migrated into the lake and 44% remained in the tributaries . Hence, a potential for long-distance migration involving anadromy seems still to be present in the population, although it is unknown if a larger proportion of individuals would have left the lake prior to the establishment of dams.
The lack of immigration from other populations into the lakes should in the short term reduce influx of locally maladaptive alleles thereby shifting the migration-selection balance in favor of local adaptation. A common garden experiment including both LIL and HAL trout demonstrated significantly different temperature-related reaction norms for early life history traits, with HAL showing adaptation to higher incubation temperatures during winter owing to the spawning tributaries being fed by ground water . Whereas this selection regime would have been similar prior to the establishment of dams, local adaptation would be expected to be reduced depending on the rate of immigration from other populations and their degree of maladaption .
In the long term, it is expected that reproductive isolation of the lakes combined with the relatively low effective population sizes should impose limits to the response to selection, due to the fact that loss of potentially adaptive variation is not counteracted by introduction of new variation through migration. This is all the more serious, as anthropogenic pressure increases the need for populations to adapt to environmental change [2, 4, 80], and specifically for brown trout it is predicted that its future distribution will become significantly reduced due to loss of suitable habitat . The requirements for maintaining evolutionary potential have been expressed in the classical “500 rule”  (but see also  and  for discussion), stating that N e should be 500 or more for new mutations to balance loss of quantitative variation by drift. This criterion is not fulfilled in HAL and MOS. Assuming no mutation, the limit to response to selection, R(∞) should be given by the expression R(∞) = 2N e R(1), where R(1) denotes the initial response to selection. Moreover, the time until 50% of the response, t50% is roughly equal to 1.4N e . If we assume that HAL became isolated 700 years ago, corresponding to 200 generations and that N e is 153, then t50% is 214 generations, approximately equal to the number of generations that have elapsed since isolation of the population. Although this is based on simplified assumptions, it nevertheless suggests that there should still be adaptive potential within HAL and the presumably larger MOS population. Of course, different conclusions could apply to other anthropogenically isolated populations with much lower N e .
Landscapes are presently being strongly modified by humans, but in some regions, such as large parts of Europe this is a process that has already been ongoing for centuries and even millennia. Hence, it is often difficult to ascertain that a given genetic population structure is natural rather than the result of anthropogenic modifications. In this study, we show that genetic divergence of lake-dwelling trout in two Danish lakes reflects establishment of water mills and impassable dams ca. 600 years ago rather than a natural genetic population structure. Also, the populations have historically been part of a larger system of anadromous brown trout. The reproductive isolation of the lakes is likely to have affected adaptive divergence among populations, in the short term by shifting the migration-selection balance. However, in the long term lack of gene flow combined with relatively low effective population size may compromise evolutionary potential.
The results demonstrate that caution is warranted when analyzing genetic population structure and interpreting it as “natural”. Even though anthropogenic modification of habitats has accelerated during the past century, it may nevertheless be a process that has been ongoing over a much longer time scale with significant consequences for the affected populations. Our results exemplify the use of molecular dating for assessing anthropogenic factors affecting wild populations. This is similar to other recent studies demonstrating that population decline in Borneo orang-utans (Pongo pygmaeus) is more recent and drastic than previously assumed , and conversely that Californian fishers (Martes pennanti) became fragmented and experienced population declines prior to European settlement .
Availability of supporting data
The data set supporting the results of this article is available from Dryad: http://datadryad.org/resource/doi:10.5061/dryad.d364t at http://datadryad.org .
Sala OE, Chapin FS, Armesto JJ, Berlow E, Bloomfield J, Dirzo R, Huber-Sanwald E, Huenneke LF, Jackson RB, Kinzig A, Leemans R, Lodge DM, Mooney HA, Oesterheld M, Poff NL, Sykes MT, Walker BH, Walker M, Wall DH: Biodiversity - Global biodiversity scenarios for the year 2100. Science. 2000, 287: 1770-1774.
Parmesan C: Ecological and evolutionary responses to recent climate change. Annu Rev Ecol Evol Syst. 2006, 37: 637-669.
Myers RA, Worm B: Rapid worldwide depletion of predatory fish communities. Nature. 2003, 423: 280-283.
Hansen MM, Waller DM, Olivieri I, Nielsen EE, The Genetic Monitoring Group: Monitoring adaptive genetic responses to environmental change. Mol Ecol. 2012, 21: 1311-1329.
Athrey G, Barr KR, Lance RF, Leberg PL: Birds in space and time: genetic changes accompanying anthropogenic habitat fragmentation in the endangered black-capped vireo (Vireo atricapilla). Evol Appl. 2012, 5: 540-552.
Bijlsma R, Loeschcke V: Genetic erosion impedes adaptive responses to stressful environments. Evol Appl. 2012, 5: 117-129.
Rosvold J, Roed KH, Hufthammer AK, Andersen R, Stenoien HK: Reconstructing the history of a fragmented and heavily exploited red deer population using ancient and contemporary DNA. BMC Evol Biol. 2012, 12: 191-
Zwick P: Stream habitat fragmentation - a threat to biodiversity. Biodivers Conserv. 1992, 1: 80-97.
Morita K, Yamamoto S: Effects of habitat fragmentation by damming on the persistence of stream-dwelling charr populations. Conserv Biol. 2002, 16: 1318-1323.
Gustafson RG, Waples RS, Myers JM, Weitkamp LA, Bryant GJ, Johnson OW, Hard JJ: Pacific salmon extinctions: Quantifying lost and remaining diversity. Conserv Biol. 2007, 21: 1009-1020.
Iwamoto EM, Myers JM, Gustafson RG: Resurrecting an extinct salmon evolutionarily significant unit: archived scales, historical DNA and implications for restoration. Mol Ecol. 2012, 21: 1567-1582.
Billard R, Lecointre G: Biology and conservation of sturgeon and paddlefish. Rev Fish Biol Fish. 2001, 10: 355-392.
Alo D, Turner TF: Effects of habitat fragmentation on effective population size in the endangered Rio Grande silvery minnow. Conserv Biol. 2005, 19: 1138-1148.
Horreo JL, Martinez JL, Ayllon F, Pola IG, Monteoliva JA, Heland M, Garcia-Vazquez E: Impact of habitat fragmentation on the genetics of populations in dendritic landscapes. Freshw Biol. 2011, 56: 2567-2579.
Letcher BH, Nislow KH, Coombs JA, O’Donnell MJ, Dubreuil TL: Population response to habitat fragmentation in a stream-dwelling brook trout population. PLoS ONE. 2007, 2: e1139-
Meldgaard T, Nielsen EE, Loeschcke V: Fragmentation by weirs in a riverine system: A study of genetic variation in time and space among populations of European grayling (Thymallus thymallus) in a Danish river system. Conserv Genet. 2003, 4: 735-747.
Faulks LK, Gilligan DM, Beheregaray LB: The role of anthropogenic vs. natural in-stream structures in determining connectivity and genetic diversity in an endangered freshwater fish, Macquarie perch (Macquaria australasica). Evol Appl. 2011, 4: 589-601.
Kawecki TJ, Ebert D: Conceptual issues in local adaptation. Ecol Lett. 2004, 7: 1225-1241.
Endler JA: Natural Selection in the Wild. 1986, Princeton: Princeton University Press
Lenormand T: Gene flow and the limits to natural selection. Trends Ecol Evol. 2002, 17: 183-189.
Vaha JP, Erkinaro J, Niemela E, Primmer CR: Life-history and habitat features influence the within-river genetic structure of Atlantic salmon. Mol Ecol. 2007, 16: 2638-2654.
Orsini L, Vanoverbeke J, Swillen I, Mergeay J, De Meester L: Drivers of population genetic differentiation in the wild: isolation by dispersal limitation, isolation by adaptation and isolation by colonization. Mol Ecol. 2013, 22: 5983-5999.
Frankham R, Ballou JD, Briscoe DA: Introduction to Conservation Genetics. 2002, Cambridge: Cambridge University Press
Garant D, Forde SE, Hendry AP: The multifarious effects of dispersal and gene flow on contemporary adaptation. Funct Ecol. 2007, 21: 434-443.
Yeaman S, Whitlock MC: The genetic architecture of adaptation under migration-selection balance. Evolution. 2011, 65: 1897-1911.
Maruyama T: Distribution of gene frequencies in a geographically structured finite population. 1. Distribution of neutral genes and of genes with small effect. Ann Hum Genet. 1972, 35: 411-423.
Franklin IR: Evolutionary change in small populations. Conservation biology: an evolutionary-ecological perspective. Edited by: Soulé ME, Wilcox BA. 1980, Sunderland: Sinnauer, 135-150.
Haldane JBS: A mathematical theory of natural and artificial selection. Proc Camb Philos Soc. 1930, 26: 220-230.
Fraser DJ, Weir LK, Bernatchez L, Hansen MM, Taylor EB: Extent and scale of local adaptation in salmonid fishes: review and meta-analysis. Heredity. 2011, 106: 404-420.
Taylor EB: A review of local adaptation in Salmonidae, with particular reference to Pacific and Atlantic salmon. Aquaculture. 1991, 98: 185-207.
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: 1413-1425.
Meier K, Hansen MM, Bekkevold D, Skaala O, Mensberg KLD: An assessment of the spatial scale of local adaptation in brown trout (Salmo trutta L.): footprints of selection at microsatellite DNA loci. Heredity. 2011, 106: 488-499.
Keller I, Taverna A, Seehausen O: Evidence of neutral and adaptive genetic divergence between European trout populations sampled along altitudinal gradients. Mol Ecol. 2011, 20: 1888-1904.
Meier K, Hansen MM, Normandeau E, Mensberg KLD, Frydenberg J, Larsen PF, Bekkevold D, Bernatchez L: Local adaptation at the transcriptome level in brown trout: evidence from early life history temperature genomic reaction norms. PLoS ONE. 2014, 9: e85171-
Rogell B, Dannewitz J, Palm S, Dahl J, Petersson E, Laurila A: Adaptive divergence in body size overrides the effects of plasticity across natural habitats in the brown trout. Ecol Evol. 2013, 3: 1931-1941.
Rogell B, Dannewitz J, Palm S, Petersson E, Dahl J, Prestegaard T, Jarvi T, Laurila A: Strong divergence in trait means but not in plasticity across hatchery and wild populations of sea-run brown trout Salmo trutta. Mol Ecol. 2012, 21: 2963-2976.
Jensen LF, Hansen MM, Pertoldi C, Holdensgaard G, Mensberg KLD, Loeschcke V: Local adaptation in brown trout early life-history traits: implications for climate change adaptability. Proc Royal Soc London Series B: Biol Sci. 2008, 275: 2859-2868.
Westley PAH, Ward EJ, Fleming IA: Fine-scale local adaptation in an invasive freshwater fish has evolved in contemporary time. Proc Royal Soc London Series B: Biol Sci. 2013, 280: 1-9.
Sand-Jensen K: Mennesket og de ferske vande. Naturen i Danmark De Ferske Vande. Edited by: Sand-Jensen K. 2013, Copenhagen, Denmark: Gyldendal, 33-46.
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.
Cornuet JM, Ravigne V, Estoup A: Inference on population history and model checking using DNA sequence and microsatellite data with the software DIYABC (v1.0). BMC Bioinformatics. 2010, 11: 401-
Hey J, Nielsen R: Integration within the Felsenstein equation for improved Markov chain Monte Carlo methods in population genetics. Proc Natl Acad Sci U S A. 2007, 104: 2785-2790.
Waples RS, Do C: Linkage disequilibrium estimates of contemporary Ne using highly variable genetic markers: a largely untapped resource for applied conservation and evolution. Evol Appl. 2010, 3: 244-262.
Hansen MM, Loeschcke V, Rasmussen G, Simonsen V: Genetic differentiation among Danish brown trout Salmo trutta populations. Hereditas. 1993, 118: 177-185.
Hansen MM, Nielsen EE, Bekkevold D, Mensberg KLD: Admixture analysis and stocking impact assessment in brown trout (Salmo trutta), estimated with incomplete baseline data. Can J Fish Aquat Sci. 2001, 58: 1853-1860.
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: 2549-2562.
Hansen MM: Estimating the long-term effects of stocking domesticated trout into wild brown trout (Salmo trutta) populations: an approach using microsatellite DNA analysis of historical and contemporary samples. Mol Ecol. 2002, 11: 1003-1015.
Goudet J: FSTAT (Version 1.2): A computer program to calculate F-statistics. J Hered. 1995, 86: 485-486.
Weir BS, Cockerham CC: Estimating F-statistics for the analysis of population structure. Evolution. 1984, 38: 1358-1370.
Dieringer D, Schlotterer C: Microsatellite Analyser (MSA): A platform independent analysis tool for large microsatellite data sets. Mol Ecol Notes. 2003, 3: 167-169.
Pritchard JK, Stephens M, Donnelly P: Inference of population structure using multilocus genotype data. Genetics. 2000, 155: 945-959.
Falush D, Stephens M, Pritchard JK: Inference of population structure using multilocus genotype data: linked loci and correlated allele frequencies. Genetics. 2003, 164: 1567-1587.
Earl DA, Vonholdt BM: STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv Genet Resour. 2012, 4: 359-361.
Jakobsson M, Rosenberg NA: CLUMPP: a cluster matching and permutation program for dealing with label switching and multimodality in analysis of population structure. Bioinformatics. 2007, 23: 1801-1806.
Rosenberg NA: DISTRUCT: a program for the graphical display of population structure. Mol Ecol Notes. 2004, 4: 137-138.
Waples RS: A bias correction for estimates of effective population size based on linkage disequilibrium at unlinked gene loci. Conserv Genet. 2006, 7: 167-184.
Waples RS, Do C: LDNE: a program for estimating effective population size from data on linkage disequilibrium. Mol Ecol Resour. 2008, 8: 753-756.
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, Luikart G, Cornuet JM: BOTTLENECK: A computer program for detecting recent reductions in the effective population size using allele frequency data. J Hered. 1999, 90: 502-503.
Dirienzo A, Peterson AC, Garza JC, Valdes AM, Slatkin M, Freimer NB: Mutational processes of simple-sequence repeat loci in human populations. Proc Natl Acad Sci U S A. 1994, 91: 3166-3170.
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: 2523-2535.
Estoup A, Jarne P, Cornuet JM: Homoplasy and mutation model at microsatellite loci and their consequences for population genetics analysis. Mol Ecol. 2002, 11: 1591-1604.
Stockwell CA, Heilveil JS, Purcell K: Estimating divergence time for two evolutionarily significant units of a protected fish species. Conserv Genet. 2013, 14: 215-222.
Garza JC, Williamson EG: Detection of reduction in population size using data from microsatellite loci. Mol Ecol. 2001, 10: 305-318.
Rannala B, Mountain JL: Detecting immigration by using multilocus genotypes. Proc Natl Acad Sci U S A. 1997, 94: 9197-9201.
Miller N, Estoup A, Toepfer S, Bourguet D, Lapchin L, Derridj S, Kim KS, Reynaud P, Furlan L, Guillemaud T: Multiple transatlantic introductions of the western corn rootworm. Science. 2005, 310: 992-992.
Beaumont MA: Joint determination of topology, divergence time and immigration in population trees. Simulations, genetics and human prehistory. Edited by: Matsumura S, Forster P, Renfrew C. 2008, Cambridge: McDonald Institute for Archeological Research, 135-154.
Pascual M, Chapuis MP, Mestres F, Balanya J, Huey RB, Gilchrist GW, Serra L, Estoup A: Introduction history of Drosophila subobscura in the New World: a microsatellite-based survey using ABC methods. Mol Ecol. 2007, 16: 3069-3083.
Estoup A, Beaumont M, Sennedot F, Moritz C, Cornuet JM: Genetic analysis of complex demographic scenarios: Spatially expanding populations of the cane toad, Bufo marinus. Evolution. 2004, 58: 2021-2036.
Guillemaud T, Beaumont MA, Ciosi M, Cornuet JM, Estoup A: Inferring introduction routes of invasive species using approximate Bayesian computation on microsatellite data. Heredity. 2010, 104: 88-99.
Sun JX, Helgason A, Masson G, Ebenesersdottir SS, Li H, Mallick S, Gnerre S, Patterson N, Kong A, Reich D, Stefansson K: A direct characterization of human mutation based on microsatellites. Nat Genet. 2012, 44: 1161-
Yue GH, David L, Orban L: Mutation rate and pattern of microsatellites in common carp ( Cyprinus carpio L.). Genetica. 2007, 129: 329-331.
Neel MC, McKelvey K, Ryman N, Lloyd MW, Bull RS, Allendorf FW, Schwartz MK, Waples RS: Estimation of effective population size in continuously distributed populations: there goes the neighborhood. Heredity. 2013, 111: 189-199.
Nei M: Molecular Evolutionary Genetics. 1987, New York: Columbia University Press
Palme A, Laikre L, Ryman N: Monitoring reveals two genetically distinct brown trout populations remaining in stable sympatry over 20 years in tiny mountain lakes. Conserv Genet. 2013, 14: 795-808.
Palm S, Laikre L, Jorde PE, Ryman N: Effective population size and temporal genetic change in stream resident brown trout (Salmo trutta, L.). Conserv Genet. 2003, 4: 249-264.
Theriault V, Garant D, Bernatchez L, Dodson JJ: Heritability of life-history tactics and genetic correlation with body size in a natural population of brook charr (Salvelinus fontinalis). J Evol Biol. 2007, 20: 2266-2277.
Boel M, Aarestrup K, Baktoft H, Larsen T, Sondergaard Madsen S, Malte H, Skov C, Svendsen JC, Koed A: The physiological basis of the migration continuum in brown trout (Salmo trutta). Physiol Biochem Zool. 2014, 87: 334-345.
Moore JS, Gow JL, Taylor EB, Hendry AP: Quantifying the constraining influence of gene flow on adaptive divergence in the lake-stream threespine stickleback system. Evolution. 2007, 61: 2015-2026.
Hoffmann AA, Sgro CM: Climate change and evolutionary adaptation. Nature. 2011, 470: 479-485.
Filipe AF, Markovic D, Pletterbauer F, Tisseuil C, De Wever A, Schmutz S, Bonada N, Freyhof O: Forecasting fish distribution along stream networks: brown trout (Salmo trutta) in Europe. Divers Distrib. 2013, 19: 1059-1071.
Lynch M, Lande R: The critical effective size for a genetically secure population. Anim Conserv. 1998, 1: 70-72.
Franklin IR, Frankham R: How large must populations be to retain evolutionary potential?. Anim Conserv. 1998, 1: 69-70.
Robertson A: A theory of limits in artificial selection. Proc Royal Soc London Series B: Biol Sci. 1960, 153: 235-249.
Goossens B, Chikhi L, Ancrenaz M, Lackman-Ancrenaz I, Andau P, Bruford MW: Genetic signature of anthropogenic population collapse in orang-utans. PLoS Biol. 2006, 4: 285-291.
Tucker JM, Schwartz MK, Truex RL, Pilgrim KL, Allendorf FW: Historical and contemporary DNA indicate fisher decline and isolation occurred prior to the European settlement of California. PLoS ONE. 2012, 7: e52803-
Hansen MM, Limborg MT, Ferchaud A, Pujolar J: Data from: The effects of Medieval dams on genetic divergence and demographic history in brown trout populations. Dryad Digital Repository. 2014, doi:10.5061/dryad.d364t at http://datadryad.org
We thank the Danish Natural Science Research Council for funding (Grant no. 272-05-0202 to MMH and a Sapere Aude Career Grant no. 12–126687 to MTL), Karen-Lise D. Mensberg and Kristian Meier for technical assistance and two anonymous reviewers for constructive comments.
The authors declare that they have no competing interests.
MMH and MTL conceived the study. MMH and MTL analyzed the data, with inputs and suggestions from ALF and JMP. MMH drafted the paper with contributions from MTL, ALF and JMP. All authors contributed to, read and approved the final manuscript.
Electronic supplementary material
Additional file 4: Figure S1: Assessment of the number of groups (k) represented by the sampled individuals, based on STRUCTURE [51, 52] analysis of multilocus genotypes. (DOCX 96 KB)
About this article
Cite this article
Hansen, M.M., Limborg, M.T., Ferchaud, A. et al. The effects of Medieval dams on genetic divergence and demographic history in brown trout populations. BMC Evol Biol 14, 122 (2014) doi:10.1186/1471-2148-14-122
- Approximate Bayesian Computation
- Divergence time
- Effective population size
- Habitat fragmentation
- Isolation-with-gene-flow model
- Microsatellite DNA