Postglacial recolonizations, watershed crossings and human translocations shape the distribution of chub lineages around the Swiss Alps
© The Author(s). 2016
Received: 11 April 2016
Accepted: 19 August 2016
Published: 9 September 2016
Distributions of European fish species were shaped by glaciations and the geological history of river networks until human activities partially abrogated the restrictions of biogeographical regions. The nearby origins of the Rhine, Rhone, Danube and Po rivers in the Swiss Alps allow the examination of historical and human-influenced patterns in fish genetic structure over a small geographic scale. We investigated these patterns in the widespread European chub (Squalius cephalus) from the Rhone, Rhine and Danube catchments and its proposed southern sister species Italian chub (Squalius squalus) from the Po catchment.
A phylogenetic tree constructed from mitochondrial Cytochrome b and COI sequences was consistent with earlier work in that it showed a separation of European chub and Italian chub, which was also reflected in microsatellite allele frequencies, morphological traits and shape differences quantified by geometric morphometrics. A new finding was that the predominant mitochondrial haplotype of European chub from the Rhine and Rhone catchments was also discovered in some individuals from Swiss populations of the Italian chub, presumably as a result of human translocation. Consistent with postglacial recolonizations from multiple refugia along the major rivers, the nuclear genetic structure of the European chub largely reflected drainage structure, but it was modified by watershed crossings between Rhine and Rhone near Lake Geneva as well as between Danube and Rhine near Lake Constance.
Our study adds new insights into the cyprinid colonization history of central Europe by showing that multiple processes shaped the distribution of different chub lineages around the Swiss Alps. Interestingly, we find evidence that cross-catchment migration has been mediated by unusual geological events such as drainage captures or watershed crossings facilitated by retreating glaciers, as well as evidence that human transport has interfered with the historical distribution of these fish (European chub haplotypes present in the Italian chub). The desirable preservation of evolutionarily distinct lineages will thus require the prevention of further translocations.
Strong climatic fluctuations during the Pleistocene induced extensive population contractions and expansions in the European flora and fauna . The geographic situation of Europe, surrounded by sea and divided by east–west mountain barriers, repeatedly trapped species in glacial refuges until recolonization during the warm interglacial periods. The repeated separation into glacial refuges led to genetic subdivision in many species, with hybrid zones forming where they came into contact again following recolonization . In freshwater fishes, populations are additionally constrained by the geological history of a region. They remain restricted to their hydrographic basins unless new interconnections or chance dispersal over land allow further expansion. Accordingly, the highest species diversity is found in the historically ice-free but isolated river catchments of Peri-Mediterranean and Ponto-Caspian Europe, and the lowest diversity in northern and central Europe ( and references therein).
Human activities have partially abrogated the borders between biogeographical regions. Today, shipping waterways such as the Rhine-Main-Danube canal or the Rhone-Rhine canal connect most important watersheds throughout Europe, and commercial and recreational fisheries spread species of interest either actively by stocking or accidentally by live bait release, increasing the potential for hybridization between related but geographically separated taxa. Indeed, non-native fish have become one of the major threats for native species of the Mediterranean region by enabling hybridization and introgression, but also by the transmission of novel parasites and diseases or competition and predation .
Here we investigate the potentially human-influenced current biogeography of European chub from the Squalius cephalus complex in Central Europe, with a particular focus on Switzerland, where four major river systems originate in close proximity (Rhine, Rhone, Danube and Po). The European chub originated during the Pliocene (3–2.5 Myr ago) in the Tigris-Euphrates basin, where secondary contacts and hybridization with other Squalius species of previous waves of colonization played an important role in the evolutionary history of the group . The invasion was at the end of an invasion period of several fish genera with Asian origin, which lasted from the Oligocene to the Pliocene (44–2.5 Myr ago), composing the large majority of today’s European fish species . A rapid radiation of the chub from its Mesopotamian origin led to four contemporary phylogroups that are represented by four major clades in a mitochondrial cytochrome b (Cyt b) phylogeny, namely a Western, Adriatic, Aegean and Eastern lineage [7, 8]. Collectively they cover almost the entire European continent . Italian and Greek haplotypes found by Durand et al.  are frequently considered a separate species, the Italian chub S. squalus [9, 10], although this is not universally adopted . Western Europe appears to have been colonized out of the Danubian basin before the last Würm glaciation, probably during the Riss - Würm interglacial period (about 100,000 years ago), and holocene recolonization of the Western lineage presumably started from multiple Würm refuges (Rhone, Rhine, Danube, and tributaries of the Black Sea) [7, 8].
Mitochondrial genetic variation
The distinct status of Italian chub from the Po catchment was also supported by quantifying sequence divergence at the barcoding gene COI only. In total we could distinguish 8 COI haplotypes that fell into the same two major clades with a sequence divergence of 3.8 % between clades (Kimura’s two-parameter model ) and a mean sequence divergence of only 0.5 % within European chub sequences and 0.2 % within Italian chub sequences. Separate phylogenies for the two mitochondrial gene sequences are provided in Additional file 1: Figure S1.
Of the European chub haplotypes (A-K), haplotype E was by far the most common. It was the predominant haplotype in the Rhine catchment, virtually the only haplotype present in the Rhone catchment, where only one European chub out of 56 had a different haplotype (D), and it was also found in two chub from the Danube catchment (Fig. 1). Surprisingly, haplotype E was also discovered in some Italian chub from two sites south of the Alps in the Po catchment (Fig. 1). This is the only haplotype shared between European and Italian chub, which are otherwise clearly distinct mitochondrially (Fig. 2).
Nuclear genetic variation
Genetic diversity measures for European chub (Squalius cephalus) from the Danube, Rhine and Rhone catchments, and for the Italian chub (S. squalus) from the Po catchment
Danube (n = 16)
Rhine (n = 74)
Rhone (n = 57)
Po (n = 20)
Genetic differentiation estimated as F ST (above diagonal) and D est (below diagonal) among chub from the Rhine, Danube and Rhone catchments (Squalius cephalus) and the Po catchment (S. squalus)
Because the presence of the mitochondrial haplotype E in Italian chub from canton Ticino indicated the possible introduction of European chub in the Po catchment, we subjected all microsatellite genotypes to an analysis with the software NewHybrids 1.1 Beta3 , but none of the individuals were identified with any confidence as either F1or F2 hybrids or backcrosses.
Differences in body shape expressed as Procrustes distances (above diagonal) and associated p-values (below diagonal) among chub from the Danube, Rhine, Rhone (Squalius cephalus) and Po catchments (S. squalus)
Pairwise tests for differences in morphology between chub from the four drainages, based on MANOVAs on scores from the first 14 PCs (explaining 80 % of the variance) from a PCA on all traditional mophological traits (see Additional file 4: Table S2)
Pillai’s Trace = 0.33 p = 0.012
Pillai’s Trace = 0.46 p = 0.003
Pillai’s Trace = 0.91 p < 0.001
Pillai’s Trace = 0.39 p = 0.003
Pillai’s Trace = 0.87 p < 0.001
Pillai’s Trace = 0.90 p < 0.001
Clear mitochondrial separation of European and Italian chub with evidence for human translocation
We observed a deep split between two major mitochondrial clades, one associated with chub from the Danube, Rhine and Rhone catchments (assigned to the Western lineage by Durand et al. ), and one associated with chub from the Po catchment (assigned to the Adriatic lineage by Durand et al. ). Sequence divergence between these clades at the DNA barcoding gene COI was 3.8 %, which corresponds to species-level divergence by commonly applied standards , although we acknowledge that gene trees can and should not be equated with species trees [21, 22]. The distinct status of chub from south of the Alps was also evident from the nuclear microsatellite data and from morphological variation. Two microsatellite loci did not amplify in most individuals from the Po catchment, suggesting high-frequency null alleles (i.e. mutations in the priming sites) at these loci, and the other loci showed strong differences in allele frequencies, resulting in highly significant genetic differentiation. The seemingly low F ST values for such a deep split are a consequence of the dependency of F ST on within-population diversity, which biases differentiation towards zero in highly polymorphic markers such as microsatellites . Using Jost’s D , which does not suffer from this problem, estimates of genetic differentiation among catchments increased and were highest between the Italian chub from the Po catchment and European chub from the other catchments (Table 2), although the differentiation between European chub from the Rhone and Danube catchments was quite high as well. Italian chub also exhibited a distinct body shape (Fig. 4), and differed significantly from European chub in the ensemble of metric, meristic and qualitative morphological traits. Taking together the mitochondrial and nuclear genetic differentiation as well as the morphological differences we observed, our data do not contradict the treatment of S. cephalus and S. squalus as distinct species [9, 10], and for practical purposes, fin color remains the easiest trait to distinguish them. Squalius squalus reliably lacked the characteristic red hue present particularly on ventral and anal fins of S. cephalus. However, if S. squalus is given species status, the same might be warranted for chub belonging to the Aegean and Eastern lineages (sensu Durand et al. ), because at least mitochondrially, these are even more differentiated from the Western lineage [7, 8]. While distinct species names have been introduced for geographically restricted taxa that, according to , fall into the Aegaen lineage (e.g. S. laietanus for chub in northern Spain and southern France , or S. orpheus for chub in northern Greece and Bulgaria ), the mitochondrially most distinct chub from the Eastern lineage remain part of S. cephalus in the 'Handbook of European Freshwater Fishes' , which appears inconsistent. Thus we agree with Seifertova et al.  that a systematic revision of the S. cephalus complex is urgently needed.
Given the clear differences between European and Italian chub and the deep split between the two mitochondrial haplogroups, we were surprised to find mitochondrial haplotype E, the dominant haplotype in European chub from the Rhine and Rhone catchments, in some Italian chub from canton Ticino (in 6 individuals out of 21). It is theoretically possible that this reflects incomplete lineage sorting , i.e. that we recovered a mitochondrial polymorphism that pre-dates the evolutionary split between European and Italian chub and has persisted in the Italian chub. However, we do not consider this explanation likely because despite considerable haplotype diversity in both mitochondrial clades (Fig. 2), we did not discover some additional haplotype related to the European chub haplogroup in Italian chub from canton Ticino, we discovered a haplotype that is an exact match to the dominant haplotype in European chub of the geographically closest European chub populations from the Rhine and Rhone catchments (Fig. 1). Also, earlier studies did not observe any shared haplotypes between chub from the Po drainage and chub from north of the Alps [7, 8]. We consider it more likely, therefore, that this finding reflects a mitochondrial introgression after human-mediated transport across the Alps. Although the European chub is not comprised in the list of 38 introduced fish species in Italy , unrecorded releases cannot be excluded, and because our study region in the Swiss canton of Ticino is well-frequented by tourists, translocation as live fishing bait is a likely route of introduction. It is maybe suggestive that four Italian chub possessing the E haplotype showed somewhat more signs of nuclear admixture than other Italian chub in our Bayesian clustering analysis (Fig. 3), yet the analysis with the software NewHybrids detected neither first or second generation hybrids nor backcrosses. That does not exclude the possibility of nuclear admixture, however, because our power to detect hybrid individuals was very limited with only 8 loci. For efficient detection of F1 hybrids, 12–24 loci are recommended , even more are required for reliable detection of backcrosses, and the additional substructure among European chub from different catchments could have additionaly complicated hybrid detection. Furthermore, the suspected translocations may well be less recent. Introgression from the European chub has already been reported in the Catalan chub (S. laietanus) in southern France . Depending on how frequently it occurs, such introgression may also be of concern for the lineage integrity of the Italian chub in Switzerland, where it has recently been taken up in the fisheries legislation as a distinct species, in accordance with the modification of the IUCN species list .
Evidence of watershed crossings in European chub
For the chub's colonization of Europe north of the Alps, a two-step scenario was proposed by Durand et al.  and Seifertova et al. : An initial Pleistocene colonization via the Danube, and a rapid post-Pleistocene expansion into Rhine and Rhone from western refugia following the last glaciation. Our results are consistent with this scenario. The Rhine and Rhone populations share the same dominant mitochondrial haplotype and possess very low mitochondrial diversity compared to European chub from the Danube, the nuclear genetic structure largely reflects drainage structure, and especially the chub from the Rhine catchment also exhibit comparatively low allelic richness (AR) at nuclear microsatellite loci (Table 1). Because AR is more sensitive than heterozygosity to founder effects , this is suggestive of a recent population expansion in this drainage.
However, the Bayesian clustering analysis of microsatellite genotypes also provided interesting evidence of genetic exchange between catchments. Firstly, chub from the Lake Geneva area of the Rhone catchment were assigned to a separate cluster and seemed more closely related to chub from the nearby parts of the upper Rhine catchment than to other chub from the Rhone catchment. This pattern is remarkable because it is also observed in the brown trout, Salmo trutta  in the European grayling, Thymallus thymallus  and in the bullhead, Cottus gobio [32, 33]. It has been explained by a postglacial watershed crossing, facilitated by the retreating Rhone glacier after the Würm glaciation . It appears that the same process has also led to the colonization of Lake Geneva by chub from the Rhine catchment, which is the first example of this for a cyprinid fish for this region. These chub remained to some extent distinct because upstream colonization from the lower Rhone was impeded by the “Pertes du Rhône” that was only submerged since the construction of Lake Genissiat in 1948 (see Background).
Secondly, many chub from the Rhine catchment are assigned to the Danube cluster with high probability, indicating extensive admixture between these two catchments. This admixture may reflect the historical opportunity for genetic exchange between the two catchments during the last retreat of the Rhine glacier (see Background), it may be mediated by the extant drainage capture via the Danube sinkhole, or nowadays it may also be facilitated by the Rhine-Main-Danube canal. Considering that admixture is most evident in the sites Rhine 3–4 and Aach (Fig. 2), all located in the vicinity of Lake Constance (Fig. 1), it is more likely that it reflects a legacy of exchange during glacial retreat or that the ongoing drainage capture through the Danube Sinkhole plays a significant role, because water from the Danube enters the Rhine drainage at the western end of Lake Constance via the Radolfzeller Aach. However, this water has to travel for several days and up to 19 km through a complex underground karst system , suggesting it is unpassable for fish. Thus, exchange via meltwater streams from the retreating Rhine glacier to the Danube at the end of the Würm glacial period seems the more likely explanation. Of course, an additional role of human-mediated transport, e.g. via the release of bait fish, cannot be excluded.
Finally, many chub from the river Neckar, which is part of the Rhine catchment, were also assigned with high probability to the genetic cluster associated with the Danube catchment. A possible explanation for this pattern is yet another drainage capture. The same observation was made in the bullhead (Cottus gobio), in which populations from particular tributaries of the Neckar also show a genetic affinity with Danubian bullheads rather than bullheads from the Rhine drainage system . The authors of that study suggest that this may be related to the fact that the river Neckar gradually eroded backward during the Pleistocene, resulting in the capture of several tributaries that were previously part of the Danubian drainage system [35, 36]. A similar scenario for the chub could also help explain the relatively high mitochondrial haplotype diversity in the Neckar compared to other sites from the Rhine catchment (Fig. 1). However, this has to remain speculative for the moment and further work is required to establish whether that is indeed why chub from the river Neckar are genetically more similar to chub from the Danube than to chub from the nearest site in the Rhine.
Four major European river catchments originate in close proximity around the Swiss Alps: The Rhone, the Rhine, the Danube and the Po catchment. A combined genetic and morphological analysis of chub from these catchments identifies the chub from the Po drainage, which are part of the Adriatic lineage sensu Durand et al. , as clearly distinct. Our analyses are thus in agreement with its treatment as a distinct species, the Italian chub S. squalus. While this is commonly (not universally) adopted [8, 9, 37], it would presumably require elevating other phylogenetic clades to species level as well to restore consistency, highlighting the need for a systematic revision of the S. cephalus complex. The other drainages contain European chub, S. cephalus, belonging to the Western lineage as defined by Durand et al. . The higher mitochondrial diversity in the Danube catchment and a nuclear genetic structure coarsely reflecting catchment structure is generally consistent with the proposed scenario of a Pleistocene colonization of central Europe from east to west via the Danube [7, 8], followed by post-glacial recolonizations from multiple refugia such as the lower Danube, the lower Rhone and possibly also the Rhine. However, we found evidence that interesting additional processes have shaped the phylogeography of the European chub around the Swiss Alps. These include a watershed crossing between the upper parts of the Rhine catchment and the Lake Geneva area of the Rhone catchment, which has also been inferred from phylogeographic patterns observed in other fish species [32, 33], as well as between the Rhine and Danube catchments near Lake Constance, which has also been inferred previously for perch . An important new finding of the present study was the discovery of the predominant European chub haplotype from the Rhone and Rhine drainages in some Italian chub from southern Switzerland. This is difficult to explain by any process other than human-mediated transport. Although we do not currently have any evidence for that in the Italian chub, there are deterrent examples of endemic fish species losing their ‘identity’ by hybridization, e.g. the Italian barbel Barbus plebejus suffering from introgressive hybridization with introduced B. barbus , or Adriatic trout getting dissolved in large populations of stocked Atlantic trout . The desirable preservation of distinct evolutionary lineages will thus require the prevention of human translocations.
Chub were collected from a total of 23 sites, three from the Danube (sites nr. 1–3), 11 from the Rhine (sites nr. 4–14), six from the Rhone (sites nr. 15–20) and three from the Po catchment (sites nr. 21–23) (Fig. 1, Additional file 1: Figure S1). Fish were caught by electro fishing with a backpack generator (FEG 1700, EFKO comm., Leutkirch, Germany) or with rod and line. The adult fish designated for morphometric analysis were anesthetized with clove oil and killed according to the animal protection laws by gill cut before photographing and taking fin clips as tissue samples. Juvenile fish were fin clipped (approx. 1 mm2) for genotyping and released thereafter. Fin clips were stored in 99 % ethanol until DNA extraction.
We used the salting-out DNA extraction protocol developed by Sunnucks and Hales , adapted to a 96 deep well plate format. Fin clips were first air dried in 8-strip microtubes. Thereafter, 300 μl of TNES buffer (50 mM Tris, pH 7.5, 400 mM NaCl, 20 nM EDTA, 0.5 % SDS) and 5 μl of 10 mg/ml proteinase K (Roche Inc., Basel, Switzerland) was added, followed by incubation at 55 °C on a shaker (Thermomixer Comfort, Eppendorf Inc., Hamburg, Germany) for 60 min at 300 rpm. Protein precipitation was performed by adding 85 μl of 5 M NaCl and shaking for 10 s. After the proteins were pelleted in a centrifuge at 4700 rpm for 10 min (Heraeus Megafuge 40R, Thermo Fisher Scientific inc, Waltham, MA, USA), the clear supernatants were transferred into a 96 deep well block. The DNA was precipitated by adding 400 μl of ice cold 100 % ethanol and pelleted by centrifugation for 10 min at 4700 rpm. The DNA pellet was washed with 700 ml of 70 % ethanol and air dried. For storage at −20 °C the DNA was resuspended in 100 μl of 1× TE buffer (100 mM Tris–HCl, 10 mM EDTA).
Two multiplex PCR reactions were used to genotype the individuals at ten microsatellite loci: LceA149, LceC1, LceCb , N7G5, N7K4  and LC128, LC27, LC290, LC32, LC93 . Amplifications were performed in a total reaction volume of 10 μl, with primer concentrations and cycling conditions as described in .
The phylogenetic relationships among different chub populations were investigated by sequencing parts of two mitochondrial genes, the cytochrome oxidase subunit I gene (COI) and cytochrome b (Cyt b). To amplify COI we used primer pair FishF1 and FishR1  and for Cyt b we used primer pair Glu and Thr . PCR conditions were as previously described  and for PCR cleanup, the NucleoSpin® Gel and PCR Clean-up kit (Macherey-Nagel, Düringen, Germany) was used according to the manufacturer’s instructions. PCR products were sequenced in both directions using the BigDye Terminator V3.1 Cycle Sequencing kit (Applied Biosystems, Foster City, CA, USA) and an ABI 3130 capillary sequencer (Applied Biosystems, Foster City, CA, USA). Sequences were edited with Geneious 6.0 (Biomatters, Auckland, New Zealand).
Genetic diversity indices
Observed heterozygosity (H O ), expected heterozygosity (H E ) and standardized allelic richness (AR) at microsatellite loci as well as pairwise F ST among catchments was calculated with FSTAT 2.9.4 . Differentiation among catchments was also estimated with Jost’s  differentiation estimator (D est ), using the R package DEMEtics . The software MicroChecker  was used to test for the presence of null alleles at the microsatellite loci. Missing data were negligible (four missing entries in the 167 individuals × 8 loci data matrix). By pooling all sites within a catchment we achieved sufficient sample sizes but thereby ignored potential substructure within these catchments, possibly violating assumptions of Hardy-Weinberg and linkage equilibrium. We used the tests implemented in FSTAT to assess if and to what extent these assumptions were indeed violated.
Microsatellite genotypes were subjected to a Bayesian Clustering analysis using the Markov chain Monte Carlo (MCMC) approach developed by Pritchard et al.  and implemented in Structure 2.3.4 . We used the admixture model with uninformative priors. Forty simulations for each number of genetic clusters (K) from K =1 to K = 10 were run with a burn-in of 50’000 iterations followed by 300’000 iterations. The most likely number of genetic clusters was inferred according to the method of Evanno et al. , using Structure Harvester . The 40 runs for the best-supported values of K were averaged using CLUMPP 1.1.2  for visualization with the software DiSTRUCT .
Microsatellite genotypes were also used for detection of potential hybrids among chub from the Po catchment. The algorithm for this MCMC approach was developed by Anderson and Thompson  and implemented in the program NewHybrids 1.1 Beta3 (available from http://ib.berkeley.edu/labs/slatkin/eriq/software/software.htm). The analysis was run without prior information as well as with priors, for the latter we assigned all individuals from Rhine, Rhone and Danube as pure parental species, because we had no evidence for introgression of Italian chub mitochondrial haplotypes in the Rhine, Rhone and Danube drainages. The simulation was run with a burn in of 50’000 and followed by 300’000 iterations.
Sequences of COI (618 bp) and Cyt b (1108 bp) were aligned using Clustal W in MEGA 6 . All sequences were indel-free and have been deposited in GenBank (accession nrs KU302616 - KU302623, KU302625 - KU302642). Because preliminary tree reconstruction from COI and Cyt b produced very similar results, the sequences of both genes were concatenated. Their evolutionary relationships were inferred by constructing a Maximum Likelihood tree based on the Tamura-Nei model  with 1000 bootstraps in MEGA 6 . The tree was rooted by including sequences from Petroleuciscus borysthenicus downloaded from Genbank (COI HM560281.1 and Cyt b HM560111.1 from the same voucher specimen and from Rutilus rutilus (accession nrs. KU302624, KU302643).
For comparability with earlier work inferring species status of fishes based on the DNA barcoding gene COI [45, 55], we also calculated pairwise haplotype differences just for COI using the Kimura 2-parameter model  which is the best model for sequences containing transitional an transversional substitutions.
Procrustes-based geometric morphometrics
Photographs of fish used in morphometric analyses were taken at a fixed distance and identical settings with a Nikon D5000 camera mounted on a tripod. Fish were placed on a white background with a fixed scale. Each fish was photographed twice but in the end one photograph was chosen at random because all were of sufficient quality. On these images we digitized 11 landmarks: tip of snout, posterior end of maxillary, posterior end of gills, posterior end of head, anterior insertion of fins (pectoral, pelvic, anal and dorsal) and superior and inferior insertion of caudal fin. Photographs were processed in a random order using TPSutil  for creating the input file to the landmark editor program TPS2 . Shape variation was analyzed in the software MorphoJ 2.03b  using a full Procrustes fit to remove variation in scale, position and orientation. To correct for allometric differences, size correction was done by using the residuals of a regression of the Procrustes coordinates on log-transformed centroid size. Size-corrected shape differences among chub from the different drainages were analyzed using a canonical variate analysis (CVA), expressed as Procrustes distances and compared with permutation tests as implemented in MorphoJ.
The measurements taken from each fish were a selection of morphological features previously presented . From the left side, the following measurements were taken: standard length, total length, predorsal length, postdorsal length, head length, dorsal head length, prepelvic length, preanal length, length of the fins (dorsal, pectoral, pelvic, anal), length of base of anal and dorsal fin, length of caudal peduncle, depth of caudal peduncle, snout length, eye diameter, postorbital length and interorbital width. All measurements were standardized as proportions of the standard length. Allometric changes were again corrected for by obtaining the residuals from a regression of the standardized measurements on size following the recommendations of Reist . We also counted the total number of rays of all fins, not discriminating between soft rays, branched rays and spiny rays. The last double ray of the dorsal and anal fin was counted as 1.5. The scales along the lateral line were also counted. Fin color was recorded as red (=0), black (=1) or mixed (=0.5). A Principal Components Analysis (PCA) executed in R  was used to reduce the dimensionality of the morphometric variables before comparisons. The PCs cumulatively explaining 80 % of the variance were retained for between-catchment comparisons of fish by MANOVAs on the scores of these PCs, carried out with the statistical software R .
We thank C. Elmiger, L. Wirtner-Bitterlin, M. Reyes and C. Baumgartner for help during field work and A. Giesen for support in the molecular lab. FORNAT AG provided electrofishing equipment. We also thank J. Van Buskirk, T. Torossi, D. Hefti, L. Keller and A. Peter for instructive discussions, and two anonymous reviewers for insightful criticism of the manuscript. Molecular data were generated in the Genetic Diversity Centre of ETH Zurich.
This study was supported by the Federal Office for the Environment (FOEN), the Swiss Federal Institute of Aquatic Science and Technology (Eawag), the Swiss Federal Institute of Technology Zurich (ETHZ) and the Fisheries Department of the canton Bern.
Availability of data and materials
Mitochondrial sequence data are deposited in GenBank under accession numbers KU302616 - KU302623 and KU302625 - KU302642.
Microsatellite genotypes and morphometric data are deposited in Dryad: http://dx.doi.org/10.5061/dryad.tj1jk
AG and CV conceived and designed the study, AG carried out the field sampling and laboratory work, AG and CV carried out the statistical analyses and jointly drafted the manuscript. Both authors read and approved the final manuscript.
The authors declare they have no competing interests.
Consent for publication
Ethics approval and consent to participate
This work was conducted with the ethical approval of the veterinarian office of canton Zurich, Switzerland (animal experimentation permit nr. 145/2010). Collection permits were provided by the fisheries authorities of the Swiss cantons Aargau, Bern, Geneva, Vaud, Solothurn, Thurgau and Zurich as well as the German state of Baden-Württemberg.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Hofreiter M, Stewart J. Ecological Change, Range Fluctuations and Population Dynamics during the Pleistocene. Curr Biol. 2009;19(14):R584–94.View ArticlePubMedGoogle Scholar
- Hewitt GM. Post-glacial re-colonization of European biota. Biol J Linnean Soc. 1999;68(1–2):87–112.View ArticleGoogle Scholar
- Reyjol Y, Hugueny B, Pont D, Bianco PG, Beier U, Caiola N, Casals F, Cowx I, Economou A, Ferreira T, et al. Patterns in species richness and endemism of European freshwater fish. Glob Ecol Biogeogr. 2007;16(1):65–75.View ArticleGoogle Scholar
- Ribeiro F, Leunda PM. Non-native fish impacts on Mediterranean freshwater ecosystems: current knowledge and research needs. Fisheries Manag Ecol. 2012;19(2):142–56.View ArticleGoogle Scholar
- Durand JD, Unlu E, Doadrio I, Pipoyan S, Templeton AR. Origin, radiation, dispersion and allopatric hybridization in the chub Leuciscus cephalus. Proc R Soc B Biol Sci. 2000;267(1453):1687–97.View ArticleGoogle Scholar
- Banarescu P. Zoogeography of Fresh Waters: Distribution and Dispersal of Freshwater Animals in North America and Eurasia. 2nd ed. Wiesbaden: AULA-Verlag; 1992.Google Scholar
- Durand JD, Persat H, Bouvet Y. Phylogeography and postglacial dispersion of the chub (Leuciscus cephalus) in Europe. Mol Ecol. 1999;8(6):989–97.View ArticlePubMedGoogle Scholar
- Seifertova M, Bryja J, Vyskocilova M, Martinkova N, Simkova A. Multiple Pleistocene refugia and post-glacial colonization in the European chub (Squalius cephalus) revealed by combined use of nuclear and mitochondrial markers. J Biogeogr. 2012;39(6):1024–40.View ArticleGoogle Scholar
- Perea S, Boehme M, Zupancic P, Freyhof J, Sanda R, Ozulug M, Abdoli A, Doadrio I. Phylogenetic relationships and biogeographical patterns in Circum-Mediterranean subfamily Leuciscinae (Teleostei, Cyprinidae) inferred from both mitochondrial and nuclear data. BMC Evol Biol. 2010;10:265.View ArticlePubMedPubMed CentralGoogle Scholar
- Kottelat M, Freyhof J. Handbook of European freshwater fishes. Cornol: Publications Kottelat; 2007: 646.Google Scholar
- Behrmann-Godel J, Gerlach G, Eckmann R. Postglacial colonization shows evidence for sympatric population splitting of Eurasian perch (Perca fluviatilis L.) in Lake Constance. Mol Ecol. 2004;13(2):491–7.View ArticlePubMedGoogle Scholar
- Keller O, Krayss E. Die Hydrographie des Bodenseeraums in Vergangenheit und Gegenwart. Berichte der St Gallischen Naturwissenschaftlichen Gesellschaft. 2000;89:39–56.Google Scholar
- Rutte E. Rhein Main Donau: Wie - wann - warum sie wurden. Eine geologische Geschichte: Sigmaringen: Jan Thorbecke Verlag Sigmaringen; 1987.Google Scholar
- Hötzl H. Origin of the Danube-Aach system. Environ Geol. 1996;27(2):87–96.Google Scholar
- Kimura M. A simple method of estimating evolutionary rate of base substitution through comparative studies of nucleotide sequences. J Mol Evol. 1980;16:111–20.View ArticlePubMedGoogle Scholar
- Jost L. G ST and its relatives do not measure differentiation. Mol Ecol. 2008;17(18):4015–26.View ArticlePubMedGoogle Scholar
- Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;155(2):945–59.PubMedPubMed CentralGoogle Scholar
- 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–20.View ArticlePubMedGoogle Scholar
- Anderson EC, Thompson EA. A model-based method for identifying species hybrids using multilocus genetic data. Genetics. 2002;160(3):1217–29.PubMedPubMed CentralGoogle Scholar
- Hebert PDN, Cywinska A, Ball SL, DeWaard JR. Biological identifications through DNA barcodes. Proc R Soc B Biol Sci. 2003;270(1512):313–21.View ArticleGoogle Scholar
- Nichols R. Gene trees and species trees are not the same. Trends Ecol Evol. 2001;16(7):358–64.View ArticlePubMedGoogle Scholar
- Balloux F. The worm in the fruit of the mitochondrial DNA tree. Heredity. 2010;104(5):419–20.View ArticlePubMedGoogle Scholar
- Doadrio I, Kottelat M, de Sostoa A. Squalius laietanus, a new species of cyprinid fish from north-eastern Spain and southern France (Teleostei : Cyprinidae). Ichthyol Explor Freshwat. 2007;18(3):247–56.Google Scholar
- Kottelat M, Economidis PS. Squalius orpheus, a new species of cyprinid fish from Evros drainage, Greece (Teleostei : Cyprinidae). Ichthyol Explor Freshwat. 2006;17(2):181–6.Google Scholar
- Kottelat M, Freyhof J. Handbook of European Freshwater Fishes. Switzerland and Berlin, Germany: Cornol; 2007.Google Scholar
- Pamilo P, Nei M. Relationships between gene trees and species trees. Mol Biol Evol. 1988;5(5):568–83.PubMedGoogle Scholar
- Gherardi F, Bertolino S, Bodon M, Casellato S, Cianfanelli S, Ferraguti M, Lori E, Mura G, Nocita A, Riccardi N, et al. Animal xenodiversity in Italian inland waters: distribution, modes of arrival, and pathways. Biol Invasions. 2008;10(4):435–54.View ArticleGoogle Scholar
- Vähä JP, Primmer CR. Efficiency of model-based Bayesian methods for detecting hybrid individuals under different hybridization scenarios and with different numbers of loci. Mol Ecol. 2006;15(1):63–72.View ArticlePubMedGoogle Scholar
- Denys GPJ, Dettai A, Persat H, Doadrio I, Cruaud C, Keith P. Status of the Catalan chub Squalius laietanus (Actinopterygii, Cyprinidae) in France: input from morphological and molecular data. Knowl Manage Aquat Ecosystems 2013(408)Google Scholar
- Freyhof, J. 2011. Squalius squalus. The IUCN Red List of Threatened Species 2011: e.T135619A4163501. http://dx.doi.org/10.2305/IUCN.UK.2008.RLTS.T135619A4163501.en.
- Allendorf FW. Genetic drift and the loss of alleles versus heterozygosity. Zoo Biol. 1986;5(2):181–90.View ArticleGoogle Scholar
- Largiadèr CR, Hefti D. Genetische Aspekte des Schutzes und der nachhaltigen Bewirtschaftung von Fischarten. In: BAFU, Mitteilung zur Fischerei, editor. Vollzug Umwelt. Bern: BAFU; 2002. p. 115.Google Scholar
- Vonlanthen P, Excoffier L, Bittner D, Persat H, Neuenschwander S, Largiader CR. Genetic analysis of potential postglacial watershed crossings in Central Europe by the bullhead (Cottus gobio L.). Mol Ecol. 2007;16(21):4572–84.View ArticlePubMedGoogle Scholar
- Riffel M, Schreiber A. Coarse-grained population structure in Central European sculpin (Cottus gobio L): Secondary contact or ongoing genetic drift? J Zool Syst Evol Res. 1995;33(4):173–84.Google Scholar
- Hantke R. Flussgeschichte Mitteleuropas. Stuttgart, Germany: Ferdinand Enke Verlag; 1993.Google Scholar
- Mader M. Die Flussgeschiche des Neckars und das Wandern des Albtraufs. Veröffentlichungen für Naturschutz und Landschaftspflege in Baden-Württemberg. 1978;47/48:443–507.Google Scholar
- Tancioni L, Russo T, Cataudella S, Milana V, Hett AK, Corsi E, Rossi AR. Testing species delimitations in four Italian sympatric Leuciscine fishes in the Tiber River: A combined morphological and molecular approach. Plos One. 2013;8(4):10.View ArticleGoogle Scholar
- Meraner A, Venturi A, Ficetola GF, Rossi S, Candiotto A, Gandolfi A. Massive invasion of exotic Barbus barbus and introgressive hybridization with endemic Barbus plebejus in Northern Italy: where, how and why? Mol Ecol. 2013;22(21):5295–312.View ArticlePubMedGoogle Scholar
- Meraner A, Gratton P, Baraldi F, Gandolfi A. Nothing but a trace left? Autochthony and conservation status of Northern Adriatic Salmo trutta inferred from PCR multiplexing, mtDNA control region sequencing and microsatellite analysis. Hydrobiologia. 2013;702(1):201–13.View ArticleGoogle Scholar
- Sunnucks P, Hales DF. Numerous transposed sequences of mitochondrial cytochrome oxidase I-II in aphids of the genus Sitobion (Hemiptera: Aphididae). Mol Biol Evol. 1996;13(3):510–24.View ArticlePubMedGoogle Scholar
- Larno V, Launey S, Devaux A, Laroche J. Isolation and characterization of microsatellite loci from chub Leuciscus cephalus (Pisces : Cyprinidae). Mol Ecol Notes. 2005;5(4):752–4.View ArticleGoogle Scholar
- Mesquita N, Cunha C, Hanfling B, Carvalho GR, Ze-Ze L, Tenreiro R, Coelho MM. Isolation and characterization of polymorphic microsatellite loci in the endangered Portuguese freshwater fish Squalius aradensis (Cyprinidae). Mol Ecol Notes. 2003;3(4):572–4.View ArticleGoogle Scholar
- Vyskocilova M, Simkova A, Martin JF. Isolation and characterization of microsatellites in Leuciscus cephalus (Cypriniformes, Cyprinidae) and cross-species amplification within the family Cyprinidae. Mol Ecol Notes. 2007;7(6):1150–4.View ArticleGoogle Scholar
- Gouskov A, Bitterlin L, Reyes M, Vorburger C. Fish Population genetic structure shaped by hydroelectric power plants in the upper Rhine catchment. Evol Appl. 2016;9:394–408.View ArticlePubMedPubMed CentralGoogle Scholar
- Ward RD, Zemlak TS, Innes BH, Last PR, Hebert PDN. DNA barcoding Australia's fish species. Philos Trans R Soc B Biol Sci. 2005;360(1462):1847–57.View ArticleGoogle Scholar
- Machordom A, Doadrio I. Evidence of a cenozoic Betic-Kabilian connection based on freshwater fish phylogeography (Luciobarbus, Cyprinidae). Mol Phylogenet Evol. 2001;18(2):252–63.View ArticlePubMedGoogle Scholar
- Goudet J: FSTAT 184.108.40.206, a program to estimate and test gene diversities and fixation indices. Available at: http://www2.unil.ch/popgen/softwares/fstat.htm (update from Goudet 1995). 2002.
- Gerlach G, Jueterbock A, Kraemer P, Deppermann J, Harmand P. Calculations of population differentiation based on G(ST) and D: forget G(ST) but not all of statistics! Mol Ecol. 2010;19(18):3845–52.View ArticlePubMedGoogle Scholar
- 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–8.View ArticleGoogle Scholar
- Earl DA, Vonholdt BM. STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv Genet Resour. 2012;4(2):359–61.View ArticleGoogle Scholar
- 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(14):1801–6.View ArticlePubMedGoogle Scholar
- Rosenberg NA. DISTRUCT: a program for the graphical display of population structure. Mol Ecol Notes. 2004;4(1):137–8.View ArticleGoogle Scholar
- Klingenberg CP. MorphoJ: an integrated software package for geometric morphometrics. Mol Ecol Resour. 2011;11:353–7.View ArticlePubMedGoogle Scholar
- Tamura K, Nei M. Estimation of the number of nucleotide substitutions in the control region of mitochondrial-DNA in humans and chimpanzees. Mol Biol Evol. 1993;10(3):512–26.PubMedGoogle Scholar
- Ward RD, Hanner R, Hebert PDN. The campaign to DNA barcode all fishes, FISH-BOL. J Fish Biol. 2009;74(2):329–56.View ArticlePubMedGoogle Scholar
- Rohlf FJ. tpsUtil, file utility program. version 1.40. Departement of Ecology and Evolution, State University of New York at Stony Brook, Stony Brook NY, USA; 2008. Available at: http://life.bio.sunysb.edu/morph.
- Rohlf FJ. tpsDig program, version 2.05. Departement of Ecology and Evolution, State University of New York at Stony Brook, Stony Brook NY, USA; 2005. Available at: http://life.bio.sunysb.edu/morph.
- Reist JD. An empirical-evaluation of several univariate methods that adjust for size variation in morphometric data. Can J Zool. 1985;63(6):1429–39.View ArticleGoogle Scholar
- R-Development-Core-Team. R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2007. Available at: http://www.R-project.org.Google Scholar