- Research article
- Open Access
Genetic variation in brown trout Salmo truttaacross the Danube, Rhine, and Elbe headwaters: a failure of the phylogeographic paradigm?
BMC Evolutionary Biology volume 13, Article number: 176 (2013)
Brown trout Salmo trutta have been described in terms of five major mtDNA lineages, four of which correspond to major ocean basins, and one, according to some authors, to a distinct taxon, marbled trout Salmo marmoratus. The Atlantic and Danubian lineages of brown trout meet in a poorly documented contact zone in Central Europe. The natural versus human mediated origin of the Atlantic lineage in the upper Danube is a question of both theoretical and practical importance with respect to conservation management. We provide a comprehensive population genetic analysis of brown trout in the region with the aim of evaluating the geographic distribution and genetic integrity of these two lineages in and around their contact zone.
Genetic screening of 114 populations of brown trout across the Danube/Rhine/Elbe catchments revealed a counter-intuitive phylogeographic structure with near fixation of the Atlantic lineage in the sampled portions of the Bavarian Danube. Along the Austrian Danube, phylogeographic informative markers revealed increasing percentages of Danube-specific alleles with downstream distance. Pure Danube lineage populations were restricted to peri-alpine isolates within previously glaciated regions. Both empirical data and simulated hybrid comparisons support that trout in non-glaciated regions north and northeast of the Alps have an admixed origin largely based on natural colonization. In contrast, the presence of Atlantic basin alleles south and southeast of the Alps stems from hatchery introductions and subsequent introgression. Despite extensive stocking of the Atlantic lineage, little evidence of first generation stocked fish or F1 hybrids were found implying that admixture has been established over time.
A purely phylogeographic paradigm fails to describe the distribution of genetic lineages of Salmo in Central Europe. The distribution pattern of the Atlantic and Danube lineages is extremely difficult to explain without invoking very strong biological mechanisms.
The peri-alpine distribution of relict populations of pure Danubian lineage brown trout implies that they colonized headwater river courses post-glacially ahead of the expansion of the Atlantic lineage. The recognition of natural as opposed to anthropogenic introgression of the Atlantic lineage into Danubian gene pools is of fundamental importance to management strategies.
An overwhelming majority of species in nature exhibit some level of phylogeographic structure and this structure is very often shown to correspond to major landscape features born out of paleo-environmental processes. The phylogeographic revolution  has promoted a paradigm, especially in Europe, whereby mountain ranges, peninsulas, and for freshwater organisms, river catchments are a priori landscape units expected to correspond to some level of genetic subdivision within a species. To the extent that studies support these expectations, the evolutionary significance of this non-biological source of genetic structure remains controversial . Moreover, broad-scale phylogeographic studies of European fishes almost invariable invoke complex within-drainage or between drainage lineage structure, which is presumed to be the result of paleo-hydrological re-arrangements of river networks or for some managed species, the result of human-mediated transport and release. Several studies support a clear distinction between genetic lineages of cold-tolerant fishes across adjacent headwater regions of the Rhine and Danube catchments such as for Telestes souffia and European grayling Thymallus thymallus, but a comparably low level of distinction for Cottus gobio. In contrast, lineage sharing between the Danube and Rhine catchments presumably based on natural colonization patterns has been demonstrated in European perch Perca fluviatilis and European chub Leuciscus cephalus and between the Danube and Vistula drainages for Barbus carpathicus. The pre-molecular view of the Ponto-Caspian (including the Danube) basin in aquatic zoogeography was that it served as the major refuge for the post-glacial re-colonization of Central and Northern Europe . While molecular studies highlight the evolutionary importance of the Ponto-Caspian basin in terms of genetic diversity, there is often little evidence or a lack of clarity concerning its role as a source of post-glacial expansion for specific species. None of the above studies were explicitly aimed at investigating the Danube/Rhine/Elbe headwater region with respect to these considerations. Moreover, the population-level knowledge base of the investigated species does not approach that of the brown trout.
The brown trout was one of the first non-model vertebrates to attract the attention of population geneticists in Europe [10–12] and subsequently one of the first organisms to be analyzed with molecular markers for a pan-European phylogeographic structure . Five major mtDNA lineages of the species were identified corresponding to four major catchment areas (Danube, Atlantic, Mediterranean & Adriatic) with the fifth lineage associated with the marbled phenotype found in the Adriatic basin. This structure has proliferated in the literature as the “five major evolutionary lineages” of the species . Subsequent studies based on intensive sampling deliver a more complicated perspective than that of major mtDNA lineages associated with major river basins (e.g. [15–17]). An ancient split within the Atlantic lineage was reported  and described as a sixth (Duero) major mtDNA lineage . The allopatric origin of the Mediterranean lineage was questioned , and the “Adriatic” lineage was considered a misnomer  as it is extensively distributed from Iberia to Turkey, with no evidence that its origin is in the Adriatic. Although the marbled trout phenotype is limited to portions of Adriatic drainage (Po drainage in Italy, Slovenia, Bosnia-Herzegovina & Montenegro), the so-called marble trout mtDNA lineage has been found in Albania , central Italy , Greece  and Corsica (Snoj A. personal communication, unpublished data). Thus, the mtDNA lineage appears to have an independent origin with respect to the taxon’s phenotypic divergence.
In Central Europe, the notion of a Danube lineage of brown trout has been well-established and unquestioned within management and conservation circles. The lineage, however, exhibits extensive life-history diversity and stretches across the entire Ponto-Caspian basin, exhibiting its deepest split across the Tigris-Euphrates catchment , leaving some uncertainty about where the lineage actually arose as well as its meaning related to evolutionary or taxonomic debate. Recently, populations of the upper Danube basin have been assigned a new name Salmo labrax, although it is unclear what the distribution of this taxon should be, or on what basis if any it should be genetically delineated. Weiss et al.  reported on the extensive sympatric occurrence of Atlantic and Danube lineage mtDNA in Austrian and Bavarian streams of the upper Danube, and argued that although introduced hatchery strains may account for much of this pattern, the statistical distribution of Atlantic lineage mtDNA suggests a degree of natural occurrence on the north slopes of the Alps. Watershed patterns of alleles at the LDH-C1 locus reflected some known river capture events across adjoining headwater tributaries of the Rhine and Danube in Baden-Württemberg . A review of molecular genetic data including extensive allozyme studies  concluded that there is no modern genetic signal of brown trout gene flow from the Ponto-Caspian Basin into north draining river systems of Central Europe (e.g. Rhine & Elbe), despite the generally accepted zoogeographic model of this basin serving northern post-glacial expansion and the fact that the Danube basin has repeatedly lost (through river capture) area to the north-flowing Rhine and Elbe systems over the last few ice ages . Nonetheless, no study to date has systematically evaluated the genetic structure of brown trout across the Danube/Rhine/Elbe interface.
In this paper we characterize the population genetic structure in the upper Danube in the Alpine and sub-Alpine regions of Austria and Bavaria including comparison to samples from adjacent areas of the Rhine and Elbe drainages (Figure 1). This study was broadly designed to: 1) evaluate the hypothesis of  and  suggesting that the upper Danube in Austria and Germany may have been naturally colonized by the Atlantic lineage; 2) identify pure Danubian lineage populations for conservation purposes, and 3) to question the utility of river basin-based phylogeographic expectations versus biological hypotheses for explaining the extant distribution of genetic variation across large-scale aquatic landscapes. More specific hypotheses concerning the potential hybrid origins of groups of populations are presented below.
MtDNA and LDH-C1 analysis
Atlantic lineage mtDNA was fixed in the Rhine and Elbe drainages and there was a general decreasing trend in its mean occurrence from west to east (i.e. downstream) within the Danube basin ranging from a high of 93.5% in the Bavarian Danube to a low of 32.1% in the Austrian Drau, south of the Alps (Table 1; Figure 2). Although the percentage of Atlantic mtDNA varied widely among sample sites within the Danube catchment, in our data only previously glaciated regions within the Austrian Mur and Drau drainages revealed some populations (N = 11) fixed for Danubian mtDNA (Table 1). An analogous west–east trend in decreasing mean frequency was seen for the Atlantic basin LDH-C1 90* allele (Table 1; Figure 1), with near fixation in the Elbe (97.8%) and Rhine (95.5%) basins and the lowest occurrence in the Austrian Drau (37.3%). The alternate LDH-100* allele was fixed in five (ETR, MOO, BIS, SAG, WOL) of the 11 populations (ETR, KOT, MOO, TIE, ANR, BIS, SAG, SCH, WOL, ANL and WIS) that were fixed for the Danubian mtDNA lineage, but nowhere else (Table 1). Fixation of the LDH-C1 90* allele was found in only two populations in the Danube catchment (in Bavaria), and both were fixed for Atlantic mtDNA (Ilz & Lech, Table 1). There was a highly significant correlation (Kendall’s tau 0.609; P < 0.001) between the frequencies of Atlantic mtDNA and the LDH-C1 90* allele across all populations. Correcting for multiple testing, there was no deviation from HWE at the population level for the LDH-C1 locus. For some comparisons in the text we will refer to populations fixed for Danube mtDNA as “pure Danubian”.
The mean allelic richness across all loci at the population level ranged from 2.5 (BIS, Drau drainage) to 7.9 (SUL, Mur drainage, Additional file 1). This richness was significantly different (Kruskal-Wallis test; P = 0.018) among drainages with the highest values found in the Austrian Danube (6.4) and the lowest in the Austrian Drau (5.1). There was a highly significant difference (Mann–Whitney U, P < 0.001) in mean allelic richness between populations fixed for Danubian mtDNA (3.8) and the remaining Danubian populations of the Mur and Drau drainages (6.2). HE likewise ranged widely across populations (0.374 to 0.836) and showed wholly analogous patterns of statistical significance within and among major basins (data not shown). There was a significant positive correlation between HE and the percent occurrence of Atlantic basin mtDNA across populations in the Austrian Mur (R2 = 0.488, P < 0.01) and the Austrian Drau (R2 = 0.391, P < 0.05) drainages supporting the influence of recent admixture; in contrast, no such correlation was found across populations in the Bavarian or the Austrian Danube.
Most populations (N = 74, 80%) showed no deviation from HWE, not considering table-wide corrections for multiple testing. The remaining samples that showed significantly negative FIS-values stem either from hatcheries or small headwater streams. Samples that showed significantly positive FIS-values are known to be either moderately or heavily managed. After correcting for multiple testing only a single heavily managed population (Schwarza, lower Austria) deviated significantly from HWE, with a positive FIS-value, presumably reflecting the presence of both wild and hatchery-origin fish.
The global FST across all 97 populations was 0.116 indicating a moderate level of differentiation among populations. Among drainage (Elbe, Rhine, Bavarian Danube, Austrian Danube, Mur and Drau) differentiation was considerably lower (FST = 0.031) while within drainage differentiation varied widely with the lowest values revealed among populations within the Bavarian Rhine (FST = 0.044) and Elbe (FST = 0.068), and the highest values found among populations within the Mur (FST =0.1465) and Drau (FST =0.188) where pure Danubian populations were found. A maximum FST value of 0.228 was observed among pure Danubian populations, presumably reflecting long periods of isolation relative to other among population comparisons within the global data set.
Using analog F-statistics (Φ), the AMOVA drainage model (Rhine, Elbe, Bavarian Danube, Austrian Danube, Mur and Drau) revealed 2.6% of the genetic variance distributed among drainages (ΦCT), 10% among populations within drainages (ΦSC), and 87.4% within populations (ΦST), all of which were highly significant (P < 0.0001). The portion of among drainage variance (ΦCT) rose to 4.8% when using R-statistics reflecting the influence of mutation on the among drainage variance. A higher (ΦCT) value (13.8%) was seen based on a reduced model for differentiation between native Danube and remaining Danubian populations reflecting the relatively high divergence of pure Danubian populations within the global data set.
The allele size permutation test further supported the influence of mutation on phylogeographic structure (and thus longer periods of isolation) for global comparisons (e.g. among all populations, or among all drainages) as well as those involving among population differentiation of the native Danube populations (Figure 3). In contrast, non-mutational mechanisms (drift and gene flow) were presumably dominant in comparisons involving Bavarian populations, supporting a surprising lack of geographic structure across the sampled sites of the three major drainages (Rhine, Elbe, Danube) in this region.
The principal component analysis of populations based on microsatellites revealed one axis (PC-1, 28.2%) generally reflecting population differentiation corresponding to the two mtDNA (Atlantic and Danubian) lineages (Figure 4). Populations fixed for, or revealing high percentages of Danubian mtDNA occur to the right along the x axis, whereas those with very low frequencies are to the left. The second axis (PC-2, 7.2%) tends to distinguish among populations within the Atlantic mtDNA lineage; hatchery populations are at the low end of the axis, whereby Bavarian populations (red text) are at the high end and Austrian Danube populations (black text) are generally intermediate (Figure 4). The location of the OBE (Oberallacher Bach) population from the Drau drainage at the low end of this axis is concordant with fixation of Atlantic lineage mtDNA despite its south Alpine location. Subsequent to our genetic screening of this population, local authorities told us that this river stretch was fishless (above a waterfall) and was stocked by an unknown source. Our data imply that the source of stocking was a pure Atlantic hatchery strain. All additional PCA factors reflected only minor portions of variance in the data set, and were neither biologically or geographically interpretable at the global scale (Additional file 2).
Two of the ten hatcheries (B &C) have an intermediate position. The first is known to be supplemented with locally caught fish, while the second is of Czech origin and presumed to be founded with fish from the Elbe drainage. A few Austrian Danube populations, all north of the Danube (e.g. KK, THA, GAL, VIE), are located intermediate along the X-axis, suggesting some admixture with pure Danube-lineage fish. A large number of Austrian Danube populations are found around the center of both axes potentially reflecting admixture between one or more lineages.
Structure analysis revealed three main clusters primarily corresponding to: (1) Bavarian populations irrespective of drainage; (2) the so-called pure Danubian populations from the Mur and Drau; and (3) remaining populations primarily from the Danube catchment (data not shown). However, only the pure Danubian populations showed a high mean proportion of self-assignment (based on mean q-values) (Table 2) (Additional file 3). Within other geographic groups, the population-level q-values (data not shown) varied widely, and some individual populations were assigned to groups outside of their region. For example, several Austrian populations north of the Danube could be assigned to the Bavarian group (e.g. the mean Q value for VIE = 0.707, HM = 0.699, KK = 0.779, and HOL = 0.702). Likewise, a number of Bavarian populations, including one from the Elbe (LAM) could be assigned to the third cluster representing hatcheries. While K = 3 showed the highest delta K, K = 4 also revealed a biologically meaningful group consisting of five Austrian populations (KK, HM, KM, GAL & VIE) all from the granitic regions (Waldviertal & Muhlviertal) north of the Danube, thus providing some level of distinction from Bavarian populations.
Statistically significant but weak isolation-by-distance signals were seen for Bavaria (R2 = 0.029, P = 0.004), the Austrian Danube (R2 = 0.075, P = 0.026) and pure Danubian populations (R2 = 0.438, P = 0.003) as well as the Drau drainage alone (Drau R2 = 0.172, P = 0.026), suggesting that natural as opposed to anthropogenic mechanisms are still operative, if not dominant, in determining among population differentiation. The population tree highly supported both the distinction of Bavaria relative to Danube populations (including hatcheries), and the two members of a south alpine group in the Danube drainages (Mur and Drau). In addition, it provided moderate support for a group containing the Rhine and Elbe catchments (Figure 5). A weak but statistically significant signal for demographic expansion was seen for the Rhine basin (P = 0.04), as well as the Rhine and Elbe basins combined (P = 0.04) using the within-locus k test. This test is noted for being more sensitive to recent as opposed to ancient expansion events compared to the inter-locus g test, which did not reveal any statistically significant signal.
Global patterns of admixture
To generate hypothetical sources of admixture, we simulated six different crosses between four population groups (identified with Structure): 1) Hatcheries, 2) pure Danubian, 3) Bavarian, 4) Northeast Austrian granitic group (KK, HM, KM, GAL & VIE). There were three presumably admixed population groups of interest: 1) the Austrian Danube, 2) the Mur, and 3) the Drau. For each of these three groups, we applied Bayesian clustering tests to each of the simulated hybrid crosses, resulting in 18 tests. Across all 18 tests, the posterior probably of individual assignment to an F1 hybrid was negligible (<1%) (Figures 6 & 7). However, meaningful levels of individual assignment to subsequent classes (i.e. F2 and backcrosses) were seen for the Austrian Danube group (Figure 6) as well as the Mur and Drau (Figure 7). The percentage assignment to a post-F1 class for Austrian Danube individuals reached 28 and 29% for the two simulated crosses involving hatcheries and other Austrian population groups (i.e. native Danube & granitic North). For the Mur and Drau individuals, assignment to a post-F1 class reached a high of 53% for the simulated cross between hatcheries and the pure Danubian populations (Figure 7).
Multiple marker systems (mtDNA, microsatellites, LDH-C1) strongly support the natural occurrence of the Atlantic lineage of brown trout throughout large areas of the upper Danube catchment, whereas there is no evidence of the Danube lineage in the sampled populations of the Bavarian Rhine and Elbe catchments and very little trace (< 7% mtDNA) of the Danube-lineage in the sampled portion of the Bavarian Danube. The cause of this large-scale pattern is natural, reflecting the admixed faunal character of the upper Danube established over both geological and glacial time scales . These results are largely concordant with the distribution of the Atlantic lineage of brown trout in Czech and Slovakian tributaries of the upper Danube , but these authors did not argue for natural colonization. While the timing and colonization corridors of especially north-to-south colonizing fish taxa (i.e. from Atlantic to the Danube basin) are poorly understood, it is clear that the genera Salvelinus and Coregonus have arrived in the upper Danube relatively recently during the last (Würm) or penultimate (Riss) ice age. As divergence for brown trout across the Rhine/Danube/Elbe contact zones are minimal to non-existent, we presume the species arrived in the upper Danube from the north during these ice age transitions. However, we emphasize that described river capture events are primarily in the opposite direction (i.e. from the Danube into the Atlantic basin), and this corresponds well with the paleo-hydrological perspective that the Danube basin has been continuously loosing area to the Rhine and Elbe basins across the ice ages .
While a relatively clear gradient of increasing Danube-specific alleles is found with increasing distance downstream (west-to-east) in the Danube catchment, the source of admixture displays a geographic as well as glacial pattern. Hatchery Atlantic strains appear to be the source of admixture in glaciated regions, especially south and southeast of the Alps, whereas natural sources of Atlantic lineage colonization from the upper Bavarian Danube appear to be the primary source of admixture north and northeast of the Alps, especially in unglaciated regions. Nonetheless, the influence of stocking in the region can be evaluated at a population level, as hatchery and wild Atlantic genotypes can be distinguished, but the overall presence and dominance of Atlantic genotypes is based on natural colonization.
Regardless of the orientation of the Alps, pure Danube lineage populations are exceedingly rare and appear limited to isolated high elevation populations in previously glaciated areas. This inference is further supported by an earlier study reporting 72% Danube lineage mtDNA and nine populations (three reported herein) fixed for Danube-lineage mtDNA in the upper Inn, Lech, Salzach and Drau drainages in Austria . Curiously, all occur within previously glaciated areas and all over 1000 meters above sea level (MASL) (Table 1; Figure 1). The distribution of populations with 100% Danube-linage mtDNA form a peri-alpine ring that is difficult to explain based on either purely anthropogenic influence or any simple natural mechanism of colonization. Since the analysis of this data set, three additional pure Danubian-lineage populations have been identified, all on the south slopes of the Alps in Carinthia, and all over 1000 MASL (Weiss, unpublished data).
Brown trout undoubtedly occupied large areas of Austria even during the height of the Würm glaciation, as glaciers did not reach the Danube from the south. No real glaciation occurred either north of the Danube or across large areas east and southeast of the Alps. In the mountainous region of Salzburg, for example, glaciers reached a maximum extension of about 5–600 MASL approximately 24–22,000 years ago . Based on accessibility, it is assumed that cold-tolerant fish re-colonized the Salzach River catchment post-glacially from the Danube and its tributaries to a maximum extent ca.11,000 year ago  during the early Holocene warm period. Between 11 and 4,000 years ago depositional and erosional processes created numerous impassible barriers isolating cold-tolerant fishes, namely brown trout and sculpin (Cottus gobio) in some river reaches found between 1 and 2,000 MASL today. We assume that this scenario was similar for many other glaciated river valleys of the Austrian Alps. Thus, pure Danubian lineage brown trout are primarily, if not exclusively found in headwater tributaries that were accessible in the early Holocene, but are now physically isolated from downstream colonization.
For regions south and southeast of the Alps, Atlantic lineage trout apparently exist exclusively through human-mediated mechanisms. That the admixture in these regions is more recent than to the north is strongly supported by the elevated levels of allelic diversity in the presence of both mtDNA strains in the Mur and Drau compared to the lack of such a correlation in the north. However, elsewhere the existence of pure Danubian lineage brown trout in isolated headwater systems means that they colonized these systems in the early stages of the current interglacial in the absence of the Atlantic lineage. This is difficult to explain as the Atlantic lineage not only dominates brown trout stocks in the upper Danube in Bavaria, but also reveals no signal of recent (i.e. post-glacial) expansion and thus must have been available for post-glacial colonization of the Alps. There is also no indication of a breeding barrier between Danubian and Atlantic lineages. Thus, the post-glacial expansion into alpine rivers of only one of these two co-existing lineages strongly implies some level of pre-glacial physical separation within the upper Danube. Interestingly, some authors have considered the Atlantic lineage brown trout stocks in the region as Salmo trutta, but the Danubian linage brown trout stocks in the upper Danube as Salmo labrax admitting a hybrid zone in the upper Danube . Regardless of one’s taxonomic viewpoint, we postulate that the Danube lineage of brown trout must have existed in numerous small-scale peri-Alpine refugia and were the initial colonizers of at least some Alpine streams following glacier retreat. Such a scenario begs the question of whether the Danube lineage is more adapted to steeper gradient or glacially influenced habitats. Increasing attention is being paid to finer than basin-scale local adaption of brown trout . Keller et al.  investigated the possibility of an elevation gradient being associated with outlier loci showing some signs of divergent selection within drainages but the association with elevation were weak, and the Danube lineage was essentially not involved in their work. The existing isolated populations of the pure Danubian lineage in high elevation habitats of Austria provide an intriguing source of material for future experiments on their potential for adaptive advantage compared to the Atlantic lineage. In contrast, one could pose the hypothesis of whether or not the Atlantic lineage has an adaptive advantage over the Danube lineage in the lower elevation river courses.
Regardless of such experimental avenues, it is apparent that a phylogeographic paradigm for European brown trout is failing to provide both researchers and managers with the proper framework for providing biological insights. From a purely management standpoint, the identification of a “Danubian” or an “Atlantic” genetic lineage is useless if not misleading in many (but not all) areas of the Danube basin, because the so-called Atlantic lineage is the dominant native trout lineage in much if not all of the Bavarian Danube and overlaps extensively in many drainages of Austria as well. Moreover, a homogeneous view of the Danubian lineage for management purposes might very well be counter-productive if strong selective mechanisms are operative in higher elevation populations. Indeed, elevation as a proxy for specific environmental conditions may turn out to be an important functional characteristic at the population or meta-population level for conservation purposes. Discarding or lumping admixed populations might also be counter-productive if some regions have arrived at this admixture through human-mediated events, while others have been admixed for thousands of years, representing unique and well-adapted gene pools. Our results and baseline data allow us to distinguish not only between native and anthropogenic introgression within the Atlantic lineage, but also the relict status of peri-alpine Danubian populations providing a valuable resource for future conservation planning.
The initial naming and description of five major mtDNA lineages in brown trout  was a breakthrough for European phylogeographic studies. Numerous subsequent studies have added detail and some controversy as to the origin of these lineages, but to date no study has provided a single biologically relevant character that is fixed or even predominant in one or more of these so-called lineages. Considering the high-level of physical fragmentation of all brown trout populations, it may be time to question the usefulness of the mtDNA phylogeographic lineages in management schemes, but also in the logical construction of basic science research programs. The universality of this concern related to other phylogeographically circumscribed lineages of plants and animals can be debated, but for highly fragmented species with large ranges it might indeed be broadly applicable, as recently suggested for European grayling , another pan-European salmonid species.
The Atlantic lineage of brown trout is native to large areas of the upper Danube in Bavaria where it predominates, but also in numerous drainages of the Austrian Danube. The lineage has been a very successful post-glacial invader of river courses of the Austrian Danube, but to date there is no evidence of natural occurrence south and southeast of the Alps in the provinces of Styria and Carinthia. Pure Danubian lineage populations are found primarily in previously glaciated regions at higher elevations, in a peri-alpine distribution implying that the lineage colonized headwater river courses post-glacially ahead of the expansion of the Atlantic lineage. The distribution pattern of these two lineages is difficult if not impossible to describe without invoking strong biological mechanisms, meaning that a purely phylogeographic paradigm for these lineages fails to explain their distribution in the upper Danube basin. More recently, anthropogenic activities have aided the spread of the Atlantic lineage and there is no clear evidence of niche segregation or a breeding barrier where the two lineages are found in sympatry today.
A total of 2568 brown trout were sampled from 97 populations (including 10 hatcheries) across Austrian and Bavarian reaches of the upper Danube basin, as well as Bavarian reaches of the upper Elbe and Rhine catchments (Table 1; Figure 1). Fish were sampled with certified back-pack electric fishing generators, licensed operators, and with written permission from the local (district) level authorities as required by law. Small (1.2 mm) samples were cut from regenerative caudal fin tissue and stored in 96% ethanol. All fish were released unharmed back into their respective habitats. A wide variety of natural and relatively intact habitats were sampled including small headwaters with little or no known history of stocking as well as moderate to heavily managed fisheries. For post-hoc descriptive purposes, each population was given a three scale rating in relation to known management history: 1) no known history of stocking and unlikely affected by dispersal from managed waters; 2) at least one known event of stocking and/or lack of isolation from managed waters: 3) known history of stocking and management. Genomic DNA was extracted from fin clips using a high-salt extraction technique . For some analyses and inferences concerning mtDNA only, published data  from an additional 17 populations primarily from the upper Inn drainage in western Austria were also used (Table 1; Figure 1). Thus the total number of populations integrated into the inferences of this study is 114, whereas results of 97 of these are reported for the first time.
Microsatellite analysis was based on a two reaction, 12 locus multi-plex assay specifically developed and optimized for the region as previously reported . SSR profiles were recorded using the GeneMapper Software v3.7 (Applied Biosystems). Following evaluation with the program Micro-Checker 2.2 , one locus (Ssa85) was removed from the analysis due to the unambiguous presence of null alleles. Variation at two additional markers (mtDNA and LDH-C) was evaluated due to their broad-scale phylogeographic information content. Raw allelic data for all loci and populations is provided [Additional file 4].
As both the Atlantic and Danubian mtDNA lineages of browns trout exist in the region [13, 14, 25], we developed an allele-specific PCR assay to assign an individual’s mtDNA to one of these two lineages. A multi-plex allele-specific assay was developed to screen two diagnostic positions in the control region (positions 26 & 389 in ). One primer (Da26 –GACTTTTCAGCTATGTACAATAACAAA) was combined with three published primers (LN19 and HN20 ; and 28ribaF ) to produce diagnostic bands for the Atlantic (452 bp) and Danubian (1035 bp) lineages. The PCR was performed in a total volume of 10 μl with 1X reaction buffer (PeqLab), 0.2 mM of each dNTP, 0.06 μM of primer At389Rd, 0.4 μM of primer Da26 + 2Fb, 0.15 μM of primer CytR, 0.03 μM of primer 28ribaF, 1.1 mM of MgCl2, 0.4 U Taq polymerase (PeqLab) and 20–50 ng genomic DNA. The reaction consisted of 3 min denaturation at 94°C, 32 cycles of 45 s at 94°C, 15 s at 61°C, 30 s at 72°C, and a final extension step of 7 min at 72°C. The reactions were loaded on a 2% agarose gel.
We developed an allele-specific duplex assay to replace the PCR-RFLP assay from  in order to efficiently screen the bi-allelic phylogeographically informative LDH-C1 locus . We selected two primers (Ldhxon4F-100 and Ldhxon4R-At) with the 3′end located on the allele-defining substitution, which in combination with Ldhxon3F and Ldhxon4R  allowed detection of all three relevant genotypes in one PCR. PCR was performed in a total volume of 10 μl with 1X reaction buffer (PeqLab), 0.2 mM of each dNTP, 0.8 μM of primer Ldxon4F-100 (ATTGTTCTCCCACGGTCAGA), 0.1 μM of primer Ldxon4R-At (GTTCGCCGTCACAGAGTAGC), 0.03 μM of primer Ldhxon3F, 0.03 μM of primer Ldhxon4R, 1.8 mM of MgCl2, 0.8 U Taq polymerase (PeqLab) and 20–50 ng genomic DNA. The reaction consisted of 3 min denaturation at 94°C, 32 cycles of 45 s at 94°C, 10 s at 70°C, 30 s at 72°C, and a final extension step of 7 min at 72°C. The reactions were loaded on a 3.8% NuSieve GTG agarose gel (Cambrex Bio Science). Amplified fragments of 340 bp and 100 bp characterize the *90 and *100 alleles, respectively, while an additional 440 bp fragment spanning regions of exons 3 and 4, as described in  is also observed.
The number of alleles per locus, allelic richness, observed (HO) and expected (HE) heterozygosity and tests for deviations from Hardy-Weinberg expectations (HWE) were performed with FSTAT 22.214.171.124 . The distribution of genetic variation across major drainage basins was evaluated with an hierarchical analysis of genetic variation (AMOVA)  using Arlequin v 3.11  whereby a number of geographic partitions were evaluated using both FST- and RST-statistics. The contribution of microsatellite allele sizes (i.e. a measure reflecting mutation) to the distribution of this variation was examined using the allele size permutation test  implemented in SPAGeDi 1.2 . The distribution of RST values from 10 000 permutations (pRST) was compared to the observed RST, which is analogous to FST when mutations are not contributing to genetic differentiation, and significantly greater than FST when they do. Such an evaluation should allow more temporal insight into the comparison of pairwise divergence estimates among groups of populations and whole basins. To gain perspective on the overall genetic relationships among the large number of populations screened we also carried out a Principal Component Analysis (PCA) on gene frequency data using the software PCAGEN 1.2 . Populations were also grouped into drainages and unrooted Neighbor-Joining trees were constructed using the program Populations 1.2.30 . Population trees were constructed based on the Dc chord distance . Node confidence was evaluated by bootstrapping over loci (500 replicates).
An a priori evaluation of genetic structure across the global data set was carried out with the program Structure v2.1  assuming admixture and independent allele frequencies. For each simulated K, the first 20,000 steps were discarded as burn-in, followed by 100,000 iterations to collect the data. We derived the most likely value of K using the second-order rate of change L”(K) following  using the on-line tool Structure Harvester (http://taylor0.biology.ucla.edu/struct_harvest). Structure analysis was repeated separately for a number of population groups defined by major drainages or sub-drainages (i.e. Rhine, Elbe, Bavarian Danube, Austrian Danube, Mur & Drau).
Correlations between pairwise geographic (GGD) and genetic (FST) distances were examined using the Mantel test implemented in the GeneAlEx 6 software . To assess whether significant drainage-wide demographic events may affect our global inferences we tested for both expansion and bottlenecks across groups of populations. Demographic expansion within drainages was assessed using the within-locus k test and the inter-locus g test  using the Excel (Microsoft) Macro developed by . The significance of k was evaluated with a one-tailed binomial distribution, whereas the significance of g was based on cut-off values as described in .
Finally, to evaluate the potential that genotypic profiles within the Austrian Danube have a hybrid origin, we combined a series of simulated hybrid data sets based on hypothetical hybrid scenarios with Bayesian assessment of membership of individual multi-locus genotypes using the software NewHybrids Version 1.1 . Genotypes from the Austrian Danube basin were tested for membership to various hybrid classes using six different hybrid scenarios, whereby at least one parental source of Atlantic origin was used in five of six scenarios and Atlantic lineage genotypes of both natural and hatchery origin were considered separately. Considering the historical complexity of these systems as well as the size of the data set, our goal with this analysis was to gain global insight on the plausibility of various hybrid scenarios without elaborating on the statistical confidence of individual genotype assignment or hybrid architecture of individual populations. As the precise procedure was based to some extent on other major results of our study, more detailed description of these analyses is integrated into the Results section.
Availability of supporting data
The data sets supporting the results of this article are included within the article and its additional files.
Avise JC: Twenty-five key evolutionary insights from the phylogeographic revolution in population genetics. Phylogeography of southern european refugia. Edited by: Weiss S, Ferrand N. 2006, Dordrecht, The Netherlands: Springer, 7-21.
Moritz C, McGuigan K, Bernatchez L: Conservation of freshwater fishes: integrating evolution and genetics with ecology. Freshwater Fish Conservation: options for the future. Edited by: Collares-Pereira MJ, Coelho MM, Cowx IG. 2002, Oxford: Blackwell Science, 293-310.
Salzburger W, Brandstatter A, Gilles A, Parson W, Hempel M, Sturmbauer C, Meyer A: Phylogeography of the vairone (Leuciscus souffia, Risso 1826) in Central Europe. Mol Ecol. 2003, 12 (9): 2371-2386. 10.1046/j.1365-294X.2003.01911.x.
Gross R, Kühn R, Baars M, Schröder W, Stein H, Rottmann O: Genetic differentiation of European grayling populations across the Main, Danube and Elbe drainages in Bavaria. J Fish Biol. 2001, 58 (1): 264-280. 10.1111/j.1095-8649.2001.tb00513.x.
Englbrecht CC, Freyhof J, Nolte A, Rassmann K, Schliewen U, Tautz D: Phylogeography of the bullhead Cottus gobio (Pisces: Teleostei: Cottidae) suggests a pre-pleistocene origin of the major central European populations. Mol Ecol. 2000, 9 (6): 709-722. 10.1046/j.1365-294x.2000.00912.x.
Nesbo CL, Fossheim T, Vollestad LA, Jakobsen KS: Genetic divergence and phylogeographic relationships among European perch (Perca fluviatilis) populations reflect glacial refugia and postglacial colonization. Mol Ecol. 1999, 8 (9): 1387-1404. 10.1046/j.1365-294x.1999.00699.x.
Durand JD, Persat H, Bouvet Y: Phylogeography and postglacial dispersion of the chub (Leuciscus cephalus) in Europe. Mol Ecol. 1999, 8 (6): 989-997. 10.1046/j.1365-294x.1999.00654.x.
Konopinski MK, Amirowicz A, Kukula K: Probable direction of the postglacial colonization of rivers on northern slopes of the Carpathian Ridge by Barbus carpathicus (Teleostei: Cyprinidae) evidenced by cline of genetic variation. J Fish Biol. 2007, 70: 406-415.
Bânârescu P: Zoogeography of fresh waters, Volume 2: distribution and dispersal of freshwater animals in North America and Eurasia. 1991, Wiesbaden: AULA-Verlag
Allendorf F, Ryman N, Stennek A, Stahl G: Genetic-variation in Scandinavian brown trout (Salmo trutta L) - evidence of distinct sympatric populations. Hereditas. 1976, 83 (1): 73-82.
Ryman N, Allendorf FW, Stahl G: Reproductive isolation with little genetic divergence in sympatric populations of brown trout (Salmo trutta). Genetics. 1979, 92 (1): 247-262.
Ferguson A, Mason FM: Allozyme evidence for reproductively isolated sympatric populations of brown trout Salmo trutta L in Lough Melvin, Ireland. J Fish Biol. 1981, 18 (6): 629-642. 10.1111/j.1095-8649.1981.tb03805.x.
Bernatchez L, Guyomard R, Bonhomme F: DNA sequence variation of the mitochondrial control region among geographically and morphologically remote European brown trout Salmo trutta populations. Mol Ecol. 1992, 1 (3): 161-173. 10.1111/j.1365-294X.1992.tb00172.x.
Bernatchez L: The evolutionary history of brown trout (Salmo trutta L.) inferred from phylogeographic, nested clade, and mismatch analyses of mitochondrial DNA variation. Evolution. 2001, 55 (2): 351-379.
Cortey M, Pla C, Garcia-Marin JL: Historical biogeography of Mediterranean trout. Mol Phylogenet Evol. 2004, 33 (3): 831-844. 10.1016/j.ympev.2004.08.012.
Sušnik S, Snoj A, Wilson IF, Mrdak D, Weiss S: Historical demography of brown trout (Salmo trutta) in the Adriatic drainage including the putative. Mol Phylogenet Evol. 2007, 44 (1): 63-76. 10.1016/j.ympev.2006.08.021.
Apostolidis AP, Madeira MJ, Hansen MM, Machordom A: Genetic structure and demographic history of brown trout (Salmo trutta) populations from the southern Balkans. Freshwater Biol. 2008, 53 (8): 1555-1566. 10.1111/j.1365-2427.2008.01988.x.
Suárez J, Bautista JM, Almodóvar A, Machordom A: Evolution of the mitochondrial control region in Palaearctic brown trout (Salmo trutta) populations: the biogeographical role of the Iberian Peninsula. Heredity. 2001, 87: 198-206. 10.1046/j.1365-2540.2001.00905.x.
Martinez P, Bouza C, Castro J, Hermida M, Pardo BG, Sanchez L: Analysis of a secondary contact between divergent lineages of brown trout Salmo trutta L. from Duero basin using microsatellites and mtDNA RFLPs. J Fish Bio. 2007, 71B: 195-213.
Snoj A, Maric S, Bajec SS, Berrebi P, Janjani S, Schoffmann J: Phylogeographic structure and demographic patterns of brown trout in North-West Africa. Mol Phylogenet Evol. 2011, 61 (1): 203-211. 10.1016/j.ympev.2011.05.011.
Splendiani A, Giovannotti M, Cerioni PN, Caniglia ML, Caputo V: Phylogeographic inferences on the native brown trout mtDNA variation in central Italy. Ital J Zool. 2006, 73 (2): 179-189. 10.1080/11250000600679751.
Apostolidis AP, Triantaphyllidis C, Kouvatsi A, Economidis PS: Mitochondrial DNA sequence variation and phylogeography among Salmo trutta L (Greek brown trout) populations. Mol Ecol. 1997, 6 (6): 531-542. 10.1046/j.1365-294X.1997.d01-176.x.
Sušnik S, Schoffmann J, Weiss S: Genetic verification of native brown trout from the Persian Gulf (Catak Cay River, Tigris basin). J Fish Biol. 2005, 67 (3): 879-884. 10.1111/j.0022-1112.2005.00780.x.
Kottelat M, Freyhof J: Handbook of European Freshwater Fishes. 2007, Cornol: Kottelat, Berlin: Freyhof
Weiss S, Schlotterer C, Waidbacher H, Jungwirth M: Haplotype (mtDNA) diversity of brown trout Salmo trutta in tributaries of the Austrian Danube: massive introgression of Atlantic basin fish - by man or nature?. Mol Ecol. 2001, 10 (5): 1241-1246. 10.1046/j.1365-294X.2001.01261.x.
Riffel M, Storch V, Schreiber A: Allozyme variability of brown trout (Salmo trutta L.) populations across the Rhenanian-Danubian watershed in southwest Germany. Heredity. 1995, 74 (3): 241-249. 10.1038/hdy.1995.37.
Weiss S, Antunes A, Schlotterer C, Alexandrino P: Mitochondrial haplotype diversity among Portuguese brown trout Salmo trutta L. populations: relevance to the post-Pleistocene recolonization of northern Europe. Mol Ecol. 2000, 9 (6): 691-698. 10.1046/j.1365-294x.2000.00910.x.
Hantke R: Flußgeschichte Mitteleuropas – Skizzen zu einer, Vegetations- und Klimageschichte der letzten 40 Millionen Jahre. 1993, Stuttgart: Ferdinand Enke Verlag
Schliewen U, Englbrecht C, Rassmann K, Miller M, Klein L, Tautz D: Veränderungen der genetischen Vielfalt: Molekulare und populations-ökologische Charkterisierung autochoner und durch Besatz beeinflusster Salmoniden-Populationen (Bachforelle, Alpen-Seesaibling) in Bayern. 2001, Berlin: Umweltbundesamt, 206-Available through http://www.umweltdaten.de/publikationen/fpdf-l/2014.pdf
Baric S, Riedl A, Meraner A, Medgyesy N, Lackner R, Pelster B, Via JD: Alpine headwater streams as reservoirs of remnant populations of the Danubian clade of brown trout. Freshwater Biol. 2010, 55 (4): 866-880. 10.1111/j.1365-2427.2009.02318.x.
van Husen D: Die Ostalpen in den Eiszeiten: palägeographischen Karte der letzen Eiszeit 1:5000.000. 1987, Wien: Geologische Bundesanstallt
Kohout J, Jaškova I, Papoušek I, Šedivá A, Šlechta V: Effects of stocking on the genetic structure of brown trout, Salmo trutta, in Central Europe inferred from mitochondrial and nuclear DNA markers. Fisheries Manag Ecol. 2012, 19: 252-263. 10.1111/j.1365-2400.2011.00828.x.
Winkler KA, Pamminger-Lahnsteiner B, Wanzenbock J, Weiss S: Hybridization and restricted gene flow between native and introduced stocks of Alpine whitefish (Coregonus sp.) across multiple environments. Mol Ecol. 2011, 20 (3): 456-472. 10.1111/j.1365-294X.2010.04961.x.
Fink J: Die Paläogeographie der Donau. Limnologie der Donau. vol. Teil III. 1966, Stuttgart: Schweizerbar’scher Verlagsbuchhandlung, 1-50.
Ibetsberger H, Jäger P, Häupl M: Der Zerfall des Salzachgletschers und die nacheiszeitliche Entwicklung des Salzburger Gewässersystems aus der Sicht der Wiederbesiedlung der Salzburger Gewässer mit Fischen. Nacheiszeitliche Entwicklung des Salzburger Gewässernetzes. 2010, Land Salzburg: Österreich, 7-54.
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 (3): 488-499. 10.1038/hdy.2010.164.
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 (9): 1888-1904. 10.1111/j.1365-294X.2011.05067.x.
Weiss S, Kopun T, Sušnik Bajec S: Assessing natural and disturbed population structure in European grayling, Thymallus thymallus (Salmonidae): melding phylogeographic, population genetic, and jurisdictional perspectives for conservation planning. J Fish Biol. 2013, 82 (2): 505-521. 10.1111/jfb.12007.
Miller SA, Dykes DD, Polesky HF: A simple salting out procedure for extracting DNA from human nucleated cells. Nucleic Acids Res. 1988, 16 (3): 1215-1215. 10.1093/nar/16.3.1215.
Lerceteau-Köhler E, Weiss S: Development of a multiplex PCR microsatellite assay in brown trout Salmo trutta, and its potential application for the genus. Aquaculture. 2006, 258 (1–4): 641-645.
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 (3): 535-538. 10.1111/j.1471-8286.2004.00684.x.
McMeel OM, Hoey EM, Ferguson A: Partial nucleotide sequences, and routine typing by polymerase chain reaction-restriction fragment length polymorphism, of the brown trout (Salmo trutta) lactate dehydrogenase, LDH-C1*90 and*100 alleles. Mol Ecol. 2001, 10 (1): 29-34. 10.1046/j.1365-294X.2001.01166.x.
Hamilton KE, Ferguson A, Taggart JB, Tomasson T, Walker A, Fahy E: Post-glacial colonization of brown trout, Salmo trutta L - LDH-5 as a phylogeographic marker locus. J Fish Biol. 1989, 35 (5): 651-664. 10.1111/j.1095-8649.1989.tb03017.x.
Goudet J: FSTAT, version 2.9.3: a program to estimate and test gene diversities and fixation indices. [Available via http://www2.unil.ch/popgen/softwares/fstat.htm]
Excoffier L, Smouse PE, Quattro JM: Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data. Genetics. 1992, 131 (2): 479-491.
Excoffier L, Laval G, Schneider S: Arlequin ver. 3.0: an integrated software package for population genetics data analysis. Evol Bioinform Online. 2005, 1: 47-50.
Hardy OJ, Charbonnel N, Freville H, Heuertz M: Microsatellite allele sizes: a simple test to assess their significance on genetic differentiation. Genetics. 2003, 163 (4): 1467-1482.
Hardy OJ, Vekemans X: SPAGeDi: a versatile computer program to analyse spatial genetic structure at the individual or population levels. Mol Ecol Notes. 2002, 2 (4): 618-620. 10.1046/j.1471-8286.2002.00305.x.
Goudet J: PCAGEN 1.2. [Available via http://www2.unil.ch/popgen/softwares/pcagen.htm]
Langella O: Populations 1.2.30. [Available http://bioinformatics.org/~tryphon/populations]
Cavalli-Sforza LL, Edwards AWF: Phylogenetic analysis: models and estimation procedures. Am J Hum Genet. 1967, 19: 233-257.
Pritchard JK, Stephens M, Donnelly P: Inference of population structure using multilocus genotype data. Genetics. 2000, 155 (2): 945-959.
Evanno G, Regnaut S, Goudet J: Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol. 2005, 14 (8): 2611-2620. 10.1111/j.1365-294X.2005.02553.x.
Peakall R, Smouse PE: GenAlEX 6: genetic analysis in excel. Population genetic software for teaching and research. Mol Ecol Notes. 2006, 6 (1): 288-295. 10.1111/j.1471-8286.2005.01155.x.
Reich DE, Goldstein DB: Genetic evidence for a paleolithic human population expansion in Africa (vol 95, pg 8119, 1998). Proc Natl Acad Sci U S A. 1998, 95 (18): 11026-11026.
Bilgin R: Kgtests: a simple excel macro program to detect signatures of population expansion using microsatellites. Mol Ecol Notes. 2007, 7 (3): 416-417. 10.1111/j.1471-8286.2006.01671.x.
Reich DE, Feldman MW, Goldstein DB: Statistical properties of two tests that use multilocus data sets to detect population expansions. Mol Biol Evol. 1999, 16 (4): 453-466. 10.1093/oxfordjournals.molbev.a026127.
Anderson EC, Thompson EA: A model-based method for identifying species hybrids using multilocus genetic data. Genetics. 2002, 160 (3): 1217-1229.
The work was primarily funded by the provincial governments of lower Austria and Styria as well as the lower Austrian fishery association under the so-called “TROUCTCHECK” project, done in close cooperation with Günter Unfer of the Institute of Hydrobiology and Aquatic Ecosystem Management of the University of Natural Resources and Life Sciences, Vienna. Additional funding was obtained through Klaus Kugi of the Nature Protection Agency of the Carinthia. Numerous people have assisted in the sample collection and several students including Ana Sprajc and Stephan Kurat contributed to the laboratory work. A special thanks goes to Dr. Otto Schwomma of the lower Austrian fishery association for his promotion and sincere interest of this project throughout all phases.
The authors declare that they have no competing interests.
EL-K carried out the majority of the laboratory work as well as nearly all statistical analyses and contributed to the first draft of the manuscript as a post-doc employed in the research group of SW. US contributed to the design and development of this research study as well as to the collection and maintenance of all Bavarian sample material. TK contributed substantially to the laboratory work and analysis. SW conceived the design of the project and supervised the research study at all stages, contributed some statistical analyses and designed and wrote the first draft of the manuscript. All authors read, contributed revisions and approved the final manuscript.
Electronic supplementary material
Additional file 1: Population Genetic Statistics of all Populations. Population genetic statistics for all sampled populations. Shown is the full name of the population, a three letter code, AR = Allelic Richness, HE = Expected Heterozygosity, HO = Observed Heterozygosity, F IS = Inbreeding Coefficient, Sig. = P-value for the inbreeding coefficient, whereby NS is not significant. (DOCX 25 KB)
Additional file 2: Additional PCA factors. Scatterplots of PCA factor 1 (X Axis) against factors 3 through 6 (Y Axis), based on a PCA of microsatellite allele frequencies in all 97 newly typed populations. The additional factors further support several outlier populations as well as the uniqueness of several Pure Danubian populations (e.g. ANR and BIS) and the differentiation between these and Bavarian populations. Several of these populations defining the Y – axis are labeled with the population code for reference. There is however no general pattern or inference that can be drawn from these additional vectors of variation. (DOCX 115 KB)
Additional file 3: Q-Values of Pure Danubian Populations. Percentage Self Assignment (Q-values from STRUCTURE analysis) of all pure Danubian populations found in this study. One known stock transport from Anrasersee to the Anlaufbach/Winbach drainage was made and is clearly evidenced here. Additionally, fish from these streams were released into Fuscher Ache by local authorities. (PDF 1 MB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Lerceteau-Köhler, E., Schliewen, U., Kopun, T. et al. Genetic variation in brown trout Salmo truttaacross the Danube, Rhine, and Elbe headwaters: a failure of the phylogeographic paradigm?. BMC Evol Biol 13, 176 (2013) doi:10.1186/1471-2148-13-176