Genetic population structure of sympatric and allopatric populations of Baltic ciscoes (Coregonus albulacomplex, Teleostei, Coregonidae)
© Mehner et al; licensee BioMed Central Ltd. 2010
Received: 1 July 2009
Accepted: 29 March 2010
Published: 29 March 2010
Teleost fishes of the Coregonidae are good model systems for studying postglacial evolution, adaptive radiation and ecological speciation. Of particular interest is whether the repeated occurrence of sympatric species pairs results from in-situ divergence from a single lineage or from multiple invasions of one or more different lineages. Here, we analysed the genetic structure of Baltic ciscoes (Coregonus albula complex), examining 271 individuals from 8 lakes in northern Germany using 1244 polymorphic AFLP loci. Six lakes had only one population of C. albula while the remaining two lakes had C. albula as well as a sympatric species (C. lucinensis or C. fontanae).
AFLP demonstrated a significant population structure (Bayesian θB = 0.22). Lower differentiation between allopatric (θB = 0.028) than sympatric (0.063-0.083) populations contradicts the hypothesis of a sympatric origin of taxa, and there was little evidence for stocking or ongoing hybridization. Genome scans found only three loci that appeared to be under selection in both sympatric population pairs, suggesting a low probability of similar mechanisms of ecological segregation. However, removal of all non-neutral loci decreased the genetic distance between sympatric pairs, suggesting recent adaptive divergence at a few loci. Sympatric pairs in the two lakes were genetically distinct from the six other C. albula populations, suggesting introgression from another lineage may have influenced these two lakes. This was supported by an analysis of isolation-by-distance, where the drift-gene flow equilibrium observed among allopatric populations was disrupted when the sympatric pairs were included.
While the population genetic data alone can not unambiguously uncover the mode of speciation, our data indicate that multiple lineages may be responsible for the complex patterns typically observed in Coregonus. Relative differences within and among lakes raises the possibility that multiple lineages may be present in northern Germany, thus understanding the postglacial evolution and speciation in the C. albula complex requires a large-scale phylogenetic analysis of several potential founder lineages.
Fish in temperate lakes that have been formed since the last glaciation (12,000 - 15,000 yrs) provide good model systems for studying adaptive radiation and ecological speciation [1, 2]. According to the ecological theory of adaptive radiation, environmental gradients and intense competition for scarce resources lead to divergent natural selection that, in turn, produces phenotypic differentiation . Reproductive isolation is thought to evolve as a by-product of divergent natural selection because traits responsible for reproductive isolation are usually adaptive . Sympatric populations of North temperate freshwater fish with significant genetic differentiation and strong ecological divergence are well known in sticklebacks (Gasterosteus aculeatus L.) , rainbow smelt (Osmerus mordax (Mitchill)) , Arctic char (Salvelinus alpinus (L.))  or sockeye salmon (Oncorhynchus nerka (Walbaum)) . But no group of fishes is as well known for adaptive radiations and complex speciation patterns in postglacial lakes as coregonid fishes (Teleostei: Coregonidae) .
A recent review  identified six potential modes of speciation in Coregonus: sympatric speciation from a single founder lineage, repeated invasions from a single lineage, hybrid swarm radiation, and three variants of allopatric speciation; complete or following secondary contact with either subsequent genomic reinforcement or with subsequent ecological reinforcement, depending on the levels of reproductive isolation between the taxa . In order to distinguish from among these potential modes of speciation, population genetic information is required in addition to phylogenetic, biogeographical, morphological and ecological data. Specifically, lower genetic differentiation between sympatric than between allopatric populations provides evidence for hybrid swarm and sympatric speciation scenarios. Conversely, more pronounced reproductive isolation between sympatric populations than between allopatric populations provides evidence for allopatric speciation with reinforcement .
Ecological diversification leading to speciation has been repeatedly suggested in Coregonus [reviewed by ]; however, recent reviews [10, 11] have challenged the assumption that all cases of ecological segregation between sympatric populations provide valid evidence for ecological speciation. Nevertheless, particularly strong support for sympatric ecological speciation can be obtained from cases of parallel ecological divergence and genetic differentiation in similar systems [12, 13]. Further evidence for ecological speciation can be gained from genome scans that detect the same loci under selection in different systems with similar environments [14–16]. In a geographical context, the genetic distances among populations of the founder lineage of sympatric species should fit to an isolation-by-distance (IBD) model , whereas adaptive divergence between the sympatric populations may disrupt the IBD patterns . Subsequently, sympatric populations would be genetically different but with a geographical distance of zero, thereby reducing our ability to detect significant IBD.
Here, we studied genetic differentiation within the C. albula complex from eight lakes using amplified fragment length polymorphism (AFLP). Six lakes had single, autumn-spawning C. albula populations, whereas two lakes contained C. albula together with a sympatric species (spring-spawning C. fontanae in Lake Stechlin and C. lucinensis in Lake Breiter Luzin) (Fig. 1). We first examined the overall genetic structure of the 10 Coregonus populations, with particular emphasis on comparing differentiation between allopatric and sympatric populations. Under the scenario of sympatric speciation, we expected to find greater genetic similarity between the sympatric populations within lakes than among lakes . Species that have evolved recently may not have developed full genetic incompatibilities and thus may hybridize with little intrinsic fitness loss, therefore we also tested whether ongoing hybridization between sympatric populations may blur historical differentiation . Second, we fit genetic distances to a stepwise IBD model and hypothesized a disruption of IBD if sympatric populations are included . Third, we calculated the probability of hybridization between sympatric populations and fish from putative stocking source populations. Enhancement stocking in the commercially important Coregonus lake populations is a common management practice [36–38], and our earlier work has quantified how this can modify the genetic identity of the target populations and disrupt IBD . We hypothesized that hybridization with C. albula stocked from allopatric autumn-spawning populations would be more likely for the autumn-spawning C. albula from both lakes than for the sympatric spring-spawning C. fontanae or C. lucinensis. Finally, we tested for the presence of outlier loci that were shared in both sympatric pairs to find genetic evidence for parallel adaptive divergence [15, 16].
All Coregonus populations and all AFLP loci
The ten sampled populations of the Coregonus albula complex
% polymorphic Loci
95% credibility interval of hs
Lake Breiter Luzin - C. albula
Lake Breiter Luzin - C. lucinensis
Lake Stechlin - C. albula
Lake Stechlin - C. fontanae
Σ = 271
Total = 98.4
Hs = 0.247
There was significant genetic differentiation among the 10 populations, with 22.52% of variation found among populations (77.48% within populations; AMOVA, d.f. = 9.261, p < 0.0001). Global values of θ (0.225, p < 0.0001) and θB (0.222; credibility interval 0.216 - 0.228 95%) were nearly identical and also showed significant population structure. All pairwise θ-values were significantly different from zero (p < 0.0001) and ranged from 0.041 (C. albula from Lakes Stechlin and Breiter Luzin) to 0.405 (C. albula from Lakes Brueckentin and Behler) (see additional file 1: AdditionalFile_FST_Matrix).
Pairwise genetic and geographical distances were not correlated for the entire dataset (Mantel r = 0.26, n = 45, p = 0.056), providing evidence that an IBD model did not appropriately describe the pattern of genetic differentiation among all 10 populations (see additional file 2: AdditionalFile_IBD_Models). An IBD model was also not appropriate to describe the genetic differentiation between all eight C. albula populations including those from Lakes Stechlin and Breiter Luzin (Mantel r = 0.31, n = 28, p = 0.086). In contrast, when only the six lakes with only a C. albula population were included, genetic distances increased significantly with geographical distance as predicted by an IBD model (Mantel r = 0.53, n = 15, p = 0.042). However, the pairwise difference between the closely located Lakes Ploener and Behler (Fig. 1) appeared to be slightly higher (0.23) than predicted by strict linear pattern (see additional file 2: AdditionalFile_IBD_Models).
Pairwise genetic differentiation for sympatric populations of Lakes Stechlin and Breiter Luzin.
Lake Breiter Luzin C. albula
Lake Breiter Luzin C. lucinensis
Lake Stechlin C. albula
Lake Stechlin C. fontanae
Lake Breiter Luzin C. albula
Lake Breiter Luzin C. lucinensis
After removal of all 95 loci putatively under selection above the p = 0.975 quantile in either Lake Stechlin or Breiter Luzin, genetic differentiation based on the remaining loci declined considerably, in particular between the sympatric pairs (Table 2). θB-values were almost identical between C. albula and C. fontanae in Lake Stechlin (0.043) and C. albula and C. lucinensis in Lake Breiter Luzin (0.049), and did not differ from the differentiation between C. fontanae and C. lucinensis (0.043) (Table 2). Only differentiation between both C. albula populations was still lower (0.027), with its credibility interval not overlapping with those of any other pairwise comparisons (Table 2).
Genetic distances between the eight C. albula populations remained almost unaffected by the removal of outlier loci. The NJ-tree calculated from the putatively neutral AFLP loci had an identical topology to the NJ-tree based on all AFLP loci (Fig. 2), with only minimal changes in branch lengths (results not shown).
Conflicting evidence for sympatric speciation
Our results suggest that the post-glacial evolution of sympatric Coregonus pairs studied here is more complex than previously thought. Neither spring spawner occurred elsewhere in allopatry, supporting the hypothesis that they evolved by sympatric speciation from the autumn-spawning C. albula within each lake, as suggested previously . Nonetheless, the pattern of genetic differentiation that we observed contradicts this hypothesis. In cases of sympatric speciation, sympatric taxa should be more closely related , but here the autumn-spawning C. albula were more similar to each other than either was to its sympatric spring spawner. This was evident from the tests of pairwise population differentiation as well as from the NJ-tree based on genetic distance. These results were similar to an allozyme study of 19 Finnish lakes, where allopatric autumn spawners formed a separate group to winter- or spring spawners .
The complexity of the evolutionary history is evidenced by the somewhat conflicting results of previous genetic analyses. Six microsatellite loci revealed greater genetic similarity between sympatric spring- and autumn-spawners within lakes, thus suggesting ongoing gene flow and a sympatric origin . In contrast, mtDNA from the same study  qualitatively support the AFLP results here, that differentiation between sympatric taxa is greater than between allopatric taxa. In particular, the two most abundant of four recorded haplotypes were shared by C. albula in both lakes. One of these was shared with the Stechlin sympatric C. fontanae and the other was shared with the Lake Breiter Luzin sympatric C. lucinensis. C. fontanae and C. lucinensis shared no haplotypes .
The reason for the discrepancy between the microsatellite and AFLP data in one of the population pairs is not clear. The large number of AFLP loci employed here (>1200) should provide a better estimate of divergence than microsatellite loci (6), although we found evidence that not all loci were strictly neutral. The two marker systems produced similar estimates of genetic structure for populations of brown trout (Salmo trutta L.) in 11 Norwegian lakes, although admixture was greater based on AFLP . AFLP results were more congruent than mtDNA with morphological characters in mormyrid electric fishes, likely due to mitochondrial introgression among lineages . Accordingly, AFLP seem to be an appropriate tool to examine the genetic differentiation of incipient ecological speciation . But as noted above, our results suggest that not all AFLP markers were strictly neutral. Removing the loci putatively under selection reduced the within-lake genetic distances considerably, and differentiation based on the remaining neutral AFLP loci was more similar to the pattern obtained from the neutral microsatellite data. Interestingly, the genetic distance to other, allopatric C. albula populations remained largely unaffected. Previous studies similarly found that differentiation between ecologically or morphologically distinct populations can be maintained for a small portion of the genome, whereas gene exchange continues to prevent divergence at the majority of loci [16, 45].
Only three outlier loci occurred in both sympatric pairs at p > 0.975. It is possible that such a small number of loci positively exceeding neutral expectations might have occurred just by chance, thus providing little evidence in support of the hypothesis that the sympatric pairs result from parallel sympatric speciation . Alternatively, it may be that only a few genes are involved. The few shared loci deserve further attention since recent analyses in the C. clupeaformis complex in North America indicated that these θ outliers can be significantly clustered around quantitative trait loci (QTL) involved in ecologically based reproductive isolation [46–48]. Proportions of outlier loci substantially exceeding 1% were found in other studies applying a similar pairwise approach [13, 15, 16, 49, 50].
The presumed modes of ecological segregation have been well studied in Lake Stechlin, where sympatric taxa (C. albula and C. fontanae) perform diel vertical migrations with different amplitudes [51, 52] and they differ in their final temperature preferendum and their metabolic costs of swimming [31, 32]. The commonly observed benthic-limnetic habitat and diet segregations typical for whitefish do not occur in Stechlin . Such a pattern of physiological segregation but low trophic differentiation presumably exists for the species pair in Lake Breiter Luzin, since both lakes are deep, with steep temperature gradients and overall low productivity . Recent data suggest C. albula and C. lucinensis in Breiter Luzin have similar diets , as in Lake Stechlin . However, adaptive divergence may be weaker in Lake Breiter Luzin than in Lake Stechlin. For example, C. lucinensis and C. albula do not show the same strong vertical microhabitat segregation (S. Schiller, Leibniz-Institute of Freshwater Ecology and Inland Fisheries, unpubl. results) as found between C. fontanae and C. albula in Lake Stechlin . A weaker adaptive divergence is supported by the lower genetic distance and a substantially lower number of loci putatively under selection in Lake Breiter Luzin than in Lake Stechlin. Accordingly, these species pairs would resemble 'isolation by adaptation' (IBA), a process in which levels of genetic differentiation among population pairs will be positively correlated with degree of adaptive phenotypic divergence [13, 16].
Stocking and hybridization
An important part of our analysis was to determine whether the genetic similarity of Stechlin and Breiter Luzin C. albula was a result of modification by stocking and hybridization. Enhancement stocking of vendace populations by freshly hatched larvae (up to 6000 ind. ha-1 annually) obtained primarily from two source populations (Lakes Arend and Tollense) has been practised in Germany for at least 50 years, with the aim of stabilizing natural recruitment and fisheries yield of sympatric and allopatric C. albula populations . Because both putative stocking sources contain only autumn-spawning C. albula populations, hybridization with C. albula in Stechlin and Breiter Luzin was hypothesized to be more likely than hybridization with the spring-spawning C. fontanae and C. lucinensis. Our results support this hypothesis. Furthermore, our recent analysis of the genetic effects of stocking on 23 allopatric C. albula populations revealed a reduction of genetic distance by about 10% between lakes stocked intensively with a genetically similar source . Accordingly, a weak secondary homogenization between Lakes Stechlin and Breiter Luzin by stocking from Lakes Arend or Tollense cannot be ruled out, but the evidence is not strong enough to conclude that the similarity of the autumn spawners is an exclusive effect of current fisheries management practices.
Hybridization within the lakes could also contribute to the observed patterns, whereby introgression of C. albula genetic polymorphisms may obscure any historical similarities of the two spring spawners. Because species that have evolved since the most recent glaciation about 12,000 years ago may not have developed full genetic incompatibilities, they may hybridize with little fitness loss . The probability that single individuals have parents or grandparents from the other population was low, but greater than zero in all four populations. Therefore, gene flow between the sympatric populations may not be completely interrupted, despite the differing spawning times. Weak ongoing hybridization between the sympatric populations can also be inferred from the higher heterozygosity of sympatric compared with allopatric C. albula populations . In this analysis, we used a high prior probability of pure ancestry similarly applied to the analysis of sympatric whitefish populations from the Canadian Cliff Lake . This sympatric whitefish population pair was identified as being reproductively isolated by AFLP-markers (θB = 0.22) , but nevertheless encompassed seven out of 47 individuals with strong evidence for mixed ancestry . Accordingly, occasional hybridization within the two sympatric populations of the C. albula complex studied here cannot be excluded.
Large-scale genetic structure of Coregonus
An important finding was that the sympatric populations of Lakes Stechlin and Breiter Luzin were genetically distinct from the six other C. albula populations (Brueckentin, Glambeck, Ploener, Tollense, Ahrend, Behler). These patterns were confirmed by the IBD analyses, which demonstrated that increasing genetic distances corresponded to increasing geographical distances only for the six lakes with single C. albula populations. The interruption of the IBD became obvious when all eight C. albula populations were included, suggesting that the C. albula in sympatry with spring-spawners are genetically distinct from other C. albula populations. Even within the six allopatric C. albula populations, the IBD model was slightly biased by a higher-than-expected genetic distance between the neighboring Lakes Ploener and Behler. This is attributable to a higher genetic similarity of the population from Lake Behler to the sympatric populations and to C. albula from Lake Arend. The C. albula population in Lake Arend is not native, but has been introduced from Lake Enzig (Outer Pomerania, now Poland) about 100 years ago . Historical introduction of C. albula in Lake Behler cannot be excluded, whereas recent stocking can be ruled out . Altogether, the results from our analyses suggest a complex evolutionary scenario for the sympatric populations, including putative parallel sympatric speciation, introgression on secondary contact, divergent selection and within-lake hybridization.
Taken together, our results suggest the following evolutionary and phylogeographical hypothesis: All lakes in northern Germany including Lakes Stechlin and Breiter Luzin were invaded from a single refuge after the last glaciation about 12,000 years ago. This founder lineage of C. albula still exists in the majority of the lakes. The genetic differentiation between these populations can be appropriately described by an IBD model [this study] which is only slightly modified by more recent stocking . In two of the deepest lakes (Stechlin and Breiter Luzin), parallel sympatric speciation along steep vertical abiotic gradients may have occurred. The putative sympatric speciation is suggested by the low within-lake genetic segregation of sympatric pairs by neutral markers . Subsequent divergent selection may have induced a stronger genetic segregation of sympatric pairs, as suggested by the AFLP markers used here. Evidence for strong recent selection comes from the fact that the removal of AFLP loci putatively under selection substantially reduced the genetic differentiation of the sympatric pairs. Furthermore, physiological studies suggest that ecological divergence is driven by temperature-based adaptations of the sympatric pair at least in Lake Stechlin [31, 32]. Recent stocking from identical autumn-spawning sources  may have slightly reduced the genetic distance between the C. albula populations in Lakes Stechlin and Breiter Luzin, as indicated by a low frequency of hybridization between these populations and C. albula from Lake Tollense. Secondary contact with another Coregonus lineage and subsequent introgression may have further shaped the genetic signature of the sympatric populations. Introgression is more likely for the autumn-spawning C. albula because the majority of all other European populations of the C. albula complex are autumn-spawners . Ongoing within-lake hybridization between sympatric populations, as demonstrated here by the STRUCTURE analyses, may have resulted in the spreading of introgressed genetic material into both populations. Introgression after secondary contact may explain why both C. albula and the spring-spawning C. fontanae and C. lucinensis of Lakes Stechlin and Breiter Luzin were genetically distinct to the allopatric C. albula in the other German lakes.
The source of the secondary contact could be from populations living farther east. This is suggested by the similarity between Lakes Stechlin and Breiter Luzin and the population from Lake Arend, which has a historical origin in Poland. Brzuzan et al.  demonstrated that two distinct lineages exist within Polish C. albula populations, with one group occurring in the east, and one in central and west Poland. These authors speculated that these may have evolved by colonization from different refugia or by introgression from C. sardinella . Interestingly, introgression from C. sardinella has also been suspected for the Lake Breiter Luzin populations . An allozyme study  found low genetic distances between C. sardinella and C. albula populations in Russia. Accordingly, we cannot exclude that populations exist in East Europe which are part of a separate lineage involved into the secondary contact with the German sympatric populations, thus explaining the genetic similarity between Lakes Arend, Stechlin and Breiter Luzin. We would like to emphasize that this complex scenario is hypothetical at the moment. We cannot exclude that the sympatric populations had an allopatric divergence or originate from two different refugia , followed by reinforcement on secondary contact in Lakes Stechlin and Breiter Luzin. In this case, adaptive differentiation might have followed nonecological speciation . Based on all available evidence from this and other studies, further clarification of speciation mode will require samples of the C. albula complex from several populations in Poland and Russia to be included in the analysis, as well as additional sympatric populations existing in Russia and Scandinavia . Comparable evidence for multiple invasions as the main mechanism for adaptive divergence within single lakes was found for C. clupeaformis  and sticklebacks  in North America. Therefore, we cannot exclude that Russian or Scandinavian sympatric populations of the C. albula complex [25, 26] underwent an allopatric stage as well.
The effects of hybridization and stocking on the commercially important Coregonus can be quantified and accounted for, thus the group remains an excellent model for the study of speciation in sympatry. Despite having used more than 1200 marker loci in a study of 10 populations, our data failed to provide conclusive evidence for whether spring-spawning Coregonus stem from parallel sympatric speciation. Parallel trends in niche differentiation and phenotypic differentiation are evident in both lakes with sympatric populations, but neutral genetic differentiation suggests allopatric speciation and introgression on secondary contact. It is clear that population genetic criteria alone cannot resolve the mechanisms of speciation, and that our analysis benefited from a large spatial scale of sampling. Examination of relative differences within and among lakes has raised the possibility that multiple lineages may be present in northern Germany. As a result, phylogenetic analysis and sampling of additional lakes from potential source populations from a broad geographic range is required for an unequivocal identification of the evolutionary mechanisms and phylogeography which have led to the ecological and genetic diversity found in the C. albula complex.
Ten populations of Baltic ciscoes were sampled from eight lakes in Germany from 2000 to 2004 (Fig. 1). In two of the lakes, C. albula occurs sympatrically with a spring-spawner; C. fontanae in Lake Stechlin and C. lucinensis in Lake Breiter Luzin. The remaining six populations of C. albula occur in allopatry in northern Germany, at variable geographic distances to Lakes Stechlin and Breiter Luzin (Fig. 1, Table 1). Two lakes (Arend, Tollense) were documented source populations for enhancement stocking for a large number of lakes in northern Germany . The remaining four lakes (Glambeck, Brueckentin, Behler, Ploener) exhibit little or no genetic evidence of stocking , and are therefore considered as near-native populations. These lakes were included to determine whether the spring-spawning sympatric taxa (C. fontanae, C. lucinensis) occur elsewhere in allopatry.
In total, 271 individuals were analysed, with the number of analysed specimens per population varying between 12 and 38. To verify that all procedures and analyses produced consistent results, the Lake Stechlin populations were sampled twice, in 2000/2001 and 2007. These samples were later combined because the pairwise θ between the repeated samplings were not significantly different from zero (see results). All fish were caught with gill nets or trawls and stored on ice immediately after removal from the nets. Subsequently, clips from the pectoral fin were stored in pure ethanol.
Genetic (AFLP) analysis
Primers for the selective amplification step of AFLP
Number of loci per primer pair
Fragments were separated on an ABI Avant 3100 capillary sequencer with an internal size standard (GeneScan-500 LIZ; Applied Biosystems). Signal processing and binning were carried out using GeneMapper® v4.0 (Applied Biosystems). Optimization yielded the following settings: bin width 1 bp, maximal peak width of 1.5 bp, scoring to presence (1) or absence (0) of fragments between 100 and 400 bp. Normalization to the sum of signal versus a fixed cut off value did not cause any differences, thus we used peak detection levels of 50 relative fluorescent units. Correct fit of the size standard was checked manually for all electropherograms. Larger fragments were excluded from the analyses to avoid errors from signal-strength drop off resulting from PCR or capillary loading bias. Only electrograms were evaluated which met the following criteria: smallest detected fragments < 50 bp and largest > 450, at least 5 fragments > 400 bp, no visible drop outs (signal remaining flat at base line for > 20 bp) in the entire spectrum, a minimum of 100 fragments in each primer combination and reproducible standards. We used 2 samples as standards in each PCR and subsequent analyses. They were considered reproducible only if they yielded > 95% fragment reproducibility to previous runs. Analyses of individual samples were repeated from restriction/ligation on until meeting the above described criteria. A second complete AFLP analysis of a random 10% of the samples yielded 96.08% reproducibility.
Genetic diversity and structure of all Coregonus populations
Preparation of input data files for analyses was conducted by AFLPDAT. Data on population genetic structure were summarized using the frequency of polymorphic fragments as calculated in AFLP-SURV. The expected heterozygosity hs, the mean within-population expected heterozygosity Hs, and their 95% credibility values, were calculated using a Bayesian approach with 250,000 generations after a burn-in of 50,000 generations implemented in HICKORY 1.1 . To account for potential slight deviations from Hardy-Weinberg equilibrium, we ran the full model which calculates an inbreeding coefficient (FIS). Because this value can be unreliable , we report only the results on overall population differentiation here.
The genetic structure within and among all 10 populations was inferred from an analysis of molecular variance (AMOVA) and a global estimate of among-population differentiation (θ) , both calculated using ARLEQUIN 3.2 . Pairwise θ was calculated for all population pairs and tested for significant deviation from zero by a numerical resampling procedure in ARLEQUIN. All analyses used 50,175 permutations. An overall estimate of genetic differentiation using a Bayesian approach (θB)  was also calculated using HICKORY as described above (θII in the notation of the HICKORY manual). A modified Nei's genetic distance  was calculated between all population pairs in order to create an unrooted neighbor-joining (NJ) tree using PHYLIP 3.6 . Subsequently, the distance matrix was bootstrapped 100,000 times using AFLP-SURV to calculate a majority-rule consensus tree in PHYLIP 3.6 .
An isolation-by-distance model was calculated by correlating the matrices of log10 direct geographical distances (km) and pairwise genetic distance (θ/1-θ)  between the lakes by a Mantel test (n = 45 pairwise comparisons, PC-ORD 5.01 ). The geographical distances between the sympatric populations in either Lakes Stechlin or Breiter Luzin were arbitrarily set to 1 km (i.e., log10 = 0). Stepwise Mantel tests were subsequently separately performed by first removing the two spring-spawning populations C. lucinensis and C. fontanae (8 lakes, n = 28 pairwise comparisons) and then removing the two C. albula populations from Lakes Stechlin and Breiter Luzin (6 lakes, n = 15 pairwise comparisons). These tests aimed to determine whether IBD was disrupted by inclusion of the sympatric populations.
Differentiation and hybridization among sympatric pairs
We calculated pairwise θB using HICKORY in order to obtain 95% credibility intervals for differentiation among the four populations that occur as sympatric pairs (20,000 burn-in iterations, 100,000 final iterations). Furthermore, STRUCTURE 2.2 [40, 55] was used to investigate possible gene flow between the sympatric populations within the two lakes, and to detect potential hybridization with stocked fish from Lakes Arend and Tollense. We assumed that most individuals of both populations each had pure ancestry but that a small proportion of individuals may have some ancestry from alternative populations. We set GENSBACK = 2 thus assuming that the mixed ancestry may have originated from the two previous generations. To ensure that there is strong statistical support for any inference of mixed ancestry, we set the prior probability of migrants to a low value (MIGRPRIOR = 0.01). The calculations were performed on the final 50,000 runs, after discarding the first 10,000 runs as a burn-in.
To evaluate those AFLP loci that may be subject to selection, a genome scan was conducted by assuming that genetic differentiation between sympatric populations is higher for those loci that are under divergent selection [4, 72]. FST and heterozygosity per polymorphic locus were calculated for the two population pairs and compared with simulated null distributions (at p = 0.95 and p = 0.99) using DFDIST[73, 74]. The simulated quantiles and the locus-specific FST were plotted against heterozygosity to determine those loci for which the neutral model could be rejected [73, 74]. These outlier loci were compared between the two lakes. At p = 0.99, the probability that loci putatively under selection are shared in both lakes just by chance is 0.01 × 0.01 × 1264 = 0.1264. Accordingly, any locus with p > 0.995 or p < 0.005 found in both lakes would be a true outlier. At p = 0.95, up to three shared loci (0.05 × 0.05 × 1264 = 3.16) could be false positives.
Subsequently, we removed all loci putatively under positive selection in either Lakes Stechlin or Breiter Luzin at p = 0.95 (n = 95, see results), and again calculated pairwise θB using HICKORY in order to obtain 95% credibility intervals for differentiation among the four populations that occur as sympatric pairs (20,000 burn-in iterations, 100,000 final iterations). Furthermore, a NJ tree was calculated based on Nei's distance between all 10 populations, based on the reduced matrix of supposedly neutral AFLP loci.
Data evaluation benefited from discussions with F. Herder, M. Foll and P. Duchesne. Sampling was supported by R. Lemcke, W. Jansen and T. Schaarschmidt. M. Degebrodt greatly helped with AFLP genotyping. Two anonymous reviewers provided stimulating comments on an earlier version of this manuscript. The study was funded by grants from the German Bundesländer Mecklenburg-Vorpommern and Schleswig-Holstein and in part from a grant within the Aquashift priority program of the German Research Council (DFG, project numbers Me 1686/5-1+2).
- Schluter D: Ecological speciation in postglacial fishes. Philosophical Transactions of the Royal Society of London Series B-Biological Sciences. 1996, 351: 807-814. 10.1098/rstb.1996.0075.View ArticleGoogle Scholar
- Taylor EB: Species pairs of north temperate freshwater fishes: Evolution, taxonomy, and conservation. Rev Fish Biol Fish. 1999, 9: 299-324. 10.1023/A:1008955229420.View ArticleGoogle Scholar
- Schluter D: The ecology of adaptive radiation. 2000, Oxford: Oxford University PressGoogle Scholar
- Campbell D, Bernatchez L: Genomic scan using AFLP markers as a means to assess the role of directional selection in the divergence of sympatric whitefish ecotypes. Mol Biol Evol. 2004, 21: 945-956. 10.1093/molbev/msh101.View ArticlePubMedGoogle Scholar
- Schluter D, McPhail JD: Character displacement and replicate adaptive radiation. Trends Ecol Evol. 1993, 8: 197-200. 10.1016/0169-5347(93)90098-A.View ArticlePubMedGoogle Scholar
- Taylor EB, Bentzen P: Evidence for multiple origins and sympatric divergence of trophic ecotypes of smelt (Osmerus) in northeastern North-America. Evolution. 1993, 47: 813-832. 10.2307/2410186.View ArticleGoogle Scholar
- Jonsson B, Jonsson N: Polymorphism and speciation in Arctic charr. J Fish Biol. 2001, 58: 605-638. 10.1111/j.1095-8649.2001.tb00518.x.View ArticleGoogle Scholar
- Wood CC, Foote CJ: Evidence for sympatric genetic divergence of anadromous and nonanadromous morphs of sockeye salmon (Oncorhynchus nerka). Evolution. 1996, 50: 1265-1279. 10.2307/2410667.View ArticleGoogle Scholar
- Hudson AG, Vonlanthen P, Müller R, Seehausen O: Review: the geography of speciation and adaptive radiation in coregonines. Arch Hydrobiol Spec Issues Advanc Limnol. 2007, 60: 111-146.Google Scholar
- Hendry AP: Ecological speciation! Or the lack thereof?. Can J Fish Aquat Sci. 2009, 66: 1383-1398. 10.1139/F09-074.View ArticleGoogle Scholar
- Rundell RJ, Price TD: Adaptive radiation, nonadaptive radiation, ecological speciation and nonecological speciation. Trends Ecol Evol. 2009, 24: 394-399. 10.1016/j.tree.2009.02.007.View ArticlePubMedGoogle Scholar
- Johannesson K: Parallel speciation: a key to sympatric divergence. Trends Ecol Evol. 2001, 16: 148-153. 10.1016/S0169-5347(00)02078-4.View ArticlePubMedGoogle Scholar
- Nosil P, Funk DJ, Ortiz-Barrientos D: Divergent selection and heterogeneous genomic divergence. Mol Ecol. 2009, 18: 375-402. 10.1111/j.1365-294X.2008.03946.x.View ArticlePubMedGoogle Scholar
- Bonin A, Taberlet P, Miaud C, Pompanon F: Explorative genome scan to detect candidate loci for adaptation along a gradient of altitude in the common frog (Rana temporaria). Mol Biol Evol. 2006, 23: 773-783. 10.1093/molbev/msj087.View ArticlePubMedGoogle Scholar
- Herder F, Pfaender J, Schliewen UK: Adaptive sympatric speciation of polychromatic "roundfin" sailfin silverside fish in Lake Matano (Sulawesi). Evolution. 2008, 62: 2178-2195. 10.1111/j.1558-5646.2008.00447.x.View ArticlePubMedGoogle Scholar
- Nosil P, Egan SP, Funk DJ: Heterogeneous genomic differentiation between walking-stick ecotypes: "Isolation by adaptation" and multiple roles for divergent selection. Evolution. 2008, 62: 316-336. 10.1111/j.1558-5646.2007.00299.x.View ArticlePubMedGoogle Scholar
- Slatkin M: Isolation by distance in equilibrium and nonequilibrium populations. Evolution. 1993, 47: 264-279. 10.2307/2410134.View ArticleGoogle Scholar
- Wagner CE, Mccune AR: Contrasting patterns of spatial genetic structure in sympatric rock-dwelling cichlid fishes. Evolution. 2009, 63: 1312-1326. 10.1111/j.1558-5646.2009.00612.x.View ArticlePubMedGoogle Scholar
- Douglas MR, Brunner PC, Bernatchez L: Do assemblages of Coregonus (Teleostei : Salmoniformes) in the Central Alpine region of Europe represent species flocks?. Mol Ecol. 1999, 8: 589-603. 10.1046/j.1365-294x.1999.00581.x.View ArticleGoogle Scholar
- Turgeon J, Estoup A, Bernatchez L: Species flock in the North American Great Lakes: Molecular ecology of Lake Nipigon Ciscoes (Teleostei: Coregonidae: Coregonus). Evolution. 1999, 53: 1857-1871. 10.2307/2640446.View ArticleGoogle Scholar
- Bernatchez L, Chouinard A, Lu GQ: Integrating molecular genetics and ecology in studies of adaptive radiation: whitefish, Coregonus sp., as a case study. Biol J Linn Soc. 1999, 68: 173-194. 10.1111/j.1095-8312.1999.tb01165.x.View ArticleGoogle Scholar
- Lu G, Basley DJ, Bernatchez L: Contrasting patterns of mitochondrial DNA and microsatellite introgressive hybridization between lineages of lake whitefish (Coregonus clupeaformis); relevance for speciation. Mol Ecol. 2001, 10: 965-985. 10.1046/j.1365-294X.2001.01252.x.View ArticlePubMedGoogle Scholar
- Ostbye K, Bernatchez L, Naesje TF, Himberg KJM, Hindar K: Evolutionary history of the European whitefish Coregonus lavaretus (L.) species complex as inferred from mtDNA phylogeography and gill-raker numbers. Mol Ecol. 2005, 14: 4371-4387.View ArticlePubMedGoogle Scholar
- Kottelat M, Freyhof J: Handbook of European freshwater fishes. 2007, Cornot and Berlin: Kottelat and FreyhofGoogle Scholar
- Sendek DS: Electrophoretic studies of Coregonid fishes from across Russia. Arch Hydrobiol Spec Issues Advanc Limnol. 2002, 57: 35-55.Google Scholar
- Svärdson G: Speciation of the Scandinavian Coregonus. Reports of the Institute of Freshwater Research Drottningholm. 1979, 57: 1-95.Google Scholar
- Mehner T, Diekmann M, Brämick U, Lemcke R: Composition of fish communities in German lakes as related to lake morphology, trophic state, shore structure and human use intensity. Freshw Biol. 2005, 50: 70-85. 10.1111/j.1365-2427.2004.01294.x.View ArticleGoogle Scholar
- Schulz M, Freyhof J: Coregonus fontanae, a new spring-spawning cisco from Lake Stechlin, northern Germany (Salmoniformes: Coregonidae). Ichthyol Explor Freshw. 2003, 14: 209-216.Google Scholar
- Helland IP, Harrod C, Freyhof J, Mehner T: Coexistence of a pair of pelagic planktivorous coregonid fish. Evol Ecol Res. 2008, 10: 373-390.Google Scholar
- Ohlberger J, Mehner T, Staaks G, Hölker F: Is ecological segregation in a sympatric species pair of coregonines supported by divergent feeding efficiencies?. Can J Fish Aquat Sci. 2008, 65: 2105-2113. 10.1139/F08-120.View ArticleGoogle Scholar
- Ohlberger J, Mehner T, Staaks G, Hölker F: Temperature-related physiological adaptations promote ecological divergence in a sympatric species pair of temperate freshwater fish, Coregonus spp. Funct Ecol. 2008, 22: 501-508. 10.1111/j.1365-2435.2008.01391.x.View ArticleGoogle Scholar
- Ohlberger J, Staaks G, Petzoldt T, Mehner T, Hölker F: Physiological specialization by thermal adaptation drives ecological divergence in a sympatric fish species pair. Evol Ecol Res. 2008, 10: 1173-1185.Google Scholar
- Schulz M, Freyhof J, Saint-Laurent R, Ostbye K, Mehner T, Bernatchez L: Evidence for independent origin of two spring-spawning ciscoes in Germany (Salmoniformes: Coregonidae). J Fish Biol. 2006, 68 (Suppl A): 119-135. 10.1111/j.0022-1112.2006.01039.x.View ArticleGoogle Scholar
- Thienemann A: Coregonus albula lucinensis, eine Tiefenform der Kleinen Maräne aus einem norddeutschen See. Zugleich ein Beitrag zur Rassenbildung bei Coregonus albula L. Z Morph Ökol Tiere. 1933, 27: 654-683. 10.1007/BF00408521.View ArticleGoogle Scholar
- Grant PR, Grant BR, Petren K: Hybridization in the recent past. Am Nat. 2005, 166: 56-67. 10.1086/430331.View ArticlePubMedGoogle Scholar
- Eckmann R, Czerkies P, Helms C, Kleibs K: Evaluating the effectiveness of stocking vendace (Coregonus albula (L.)) eleutheroembryos by alizarin marking of otoliths. Arch Hydrobiol Spec Issues Advanc Limnol. 1998, 50: 457-463.Google Scholar
- Douglas MR, Brunner PC: Biodiversity of Central Alpine Coregonus (Salmoniformes): Impact of one-hundred years of management. Ecol Appl. 2002, 12: 154-172. 10.1890/1051-0761(2002)012[0154:BOCACS]2.0.CO;2.View ArticleGoogle Scholar
- Eckmann R, Kugler M, Ruhle C: Evaluating the success of large-scale whitefish stocking at Lake Constance. Arch Hydrobiol Spec Issues Advanc Limnol. 2006, 60: 361-368.Google Scholar
- Mehner T, Pohlmann K, Elkin C, Monaghan MT, Freyhof J: Genetic mixing from enhancement stocking in commercially exploited vendace populations. J Appl Ecol. 2009, 46: 1340-1349.Google Scholar
- Pritchard JK, Stephens M, Donnelly P: Inference of population structure using multilocus genotype data. Genetics. 2000, 155: 945-959.PubMed CentralPubMedGoogle Scholar
- Vuorinen J, Himberg MKJ, Lankinen P: Genetic differentiation in Coregonus albula (L) (Salmonidae) populations in Finland. Hereditas. 1981, 94: 113-121.View ArticleGoogle Scholar
- Sønstebø JH, Borgstrom R, Heun M: A comparison of AFLPs and microsatellites to identify the population structure of brown trout (Salmo trutta L.) populations from Hardangervidda, Norway. Mol Ecol. 2007, 16: 1427-1438. 10.1111/j.1365-294X.2007.03256.x.View ArticlePubMedGoogle Scholar
- Sullivan JP, Lavoue S, Arnegard ME, Hopkins CD: AFLPs resolve phylogeny and reveal mitochondrial introgression within a species flock of African electric fish (Mormyroidea: Teleostei). Evolution. 2004, 58: 825-841.View ArticlePubMedGoogle Scholar
- Bensch S, Akesson M: Ten years of AFLP in ecology and evolution: why so few animals?. Mol Ecol. 2005, 14: 2899-2914. 10.1111/j.1365-294X.2005.02655.x.View ArticlePubMedGoogle Scholar
- Wilding CS, Butlin RK, Grahame J: Differential gene exchange between parapatric morphs of Littorina saxatilis detected using AFLP markers. Journal of Evolutionary Biology. 2001, 14: 611-619. 10.1046/j.1420-9101.2001.00304.x.View ArticleGoogle Scholar
- Via S, West J: The genetic mosaic suggests a new role for hitchhiking in ecological speciation. Mol Ecol. 2008, 17: 4334-4345. 10.1111/j.1365-294X.2008.03921.x.View ArticlePubMedGoogle Scholar
- Rogers SM, Bernatchez L: The genetic architecture of ecological speciation and the association with signatures of selection in natural lake whitefish (Coregonus sp. Salmonidae) species pairs. Mol Biol Evol. 2007, 24: 1423-1438. 10.1093/molbev/msm066.View ArticlePubMedGoogle Scholar
- Rogers SM, Bernatchez L: Integrating QTL mapping and genome scans towards the characterization of candidate loci under parallel selection in the lake whitefish (Coregonus clupeaformis). Mol Ecol. 2005, 14: 351-361. 10.1111/j.1365-294X.2004.02396.x.View ArticlePubMedGoogle Scholar
- Egan SP, Nosil P, Funk DJ: Selection and genomic differentiation during ecological speciation: Isolating the contributions of host association via a comparative genome scan of Neochlamisus bebbianae leaf beetles. Evolution. 2008, 62: 1162-1181. 10.1111/j.1558-5646.2008.00352.x.View ArticlePubMedGoogle Scholar
- Gagnaire PA, Albert V, Jonsson B, Bernatchez L: Natural selection influences AFLP intraspecific genetic variability and introgression patterns in Atlantic eels. Mol Ecol. 2009, 18: 1678-1691. 10.1111/j.1365-294X.2009.04142.x.View ArticlePubMedGoogle Scholar
- Helland IP, Freyhof J, Kasprzak P, Mehner T: Temperature sensitivity of vertical distributions of zooplankton and planktivorous fish in a stratified lake. Oecologia. 2007, 151: 322-330. 10.1007/s00442-006-0541-x.View ArticlePubMedGoogle Scholar
- Mehner T, Kasprzak P, Hölker F: Exploring ultimate hypotheses to predict diel vertical migrations in coregonid fish. Can J Fish Aquat Sci. 2007, 64: 874-886. 10.1139/F07-067.View ArticleGoogle Scholar
- Scharf J, Krappe M, Koschel R, Waterstraat A: Feeding of European cisco (Coregonus albula and C. lucinensis) on the glacial relict crustacean Mysis relicta in Lake Breiter Luzin (Germany). Limnologica. 2008, 38: 147-158.View ArticleGoogle Scholar
- Araguas RM, Sanz N, Pla C, Garcia-Marin JL: Breakdown of the brown trout evolutionary history due to hybridization between native and cultivated fish. J Fish Biol. 2004, 65: 28-37. 10.1111/j.0022-1112.2004.00544.x.View ArticleGoogle Scholar
- Falush D, Stephens M, Pritchard JK: Inference of population structure using multilocus genotype data: dominant markers and null alleles. Mol Ecol Notes. 2007, 7: 574-578. 10.1111/j.1471-8286.2007.01758.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Kammerad B, Ellermann S, Mencke J, Wüstemann O, Zuppke U: Die Fischfauna von Sachsen-Anhalt - Verbreitungsatlas. 1997, Magdeburg: Ministerium für Raumordnung, Landwirtschaft und UmweltGoogle Scholar
- Brzuzan P, Kozlowski J, Fopp D: Genetic structure of Polish populations of vendace (Coregonus albula L.) inferred from mitochondrial DNA. Arch Hydrobiol Spec Issues Advanc Limnol. 2002, 57: 1-10.Google Scholar
- Bernatchez L, Dodson JJ: Allopatric origin of sympatric populations of lake whitefish (Coregonus clupeaformis) as revealed by mitochondrial-DNA restriction analysis. Evolution. 1990, 44: 1263-1271. 10.2307/2409287.View ArticleGoogle Scholar
- Taylor EB, McPhail JD: Historical contingency and ecological determinism interact to prime speciation in sticklebacks, Gasterosteus. Proc R Soc Lond B. 2000, 267: 2375-2384. 10.1098/rspb.2000.1167.View ArticleGoogle Scholar
- Meudt HM, Clarke AC: Almost forgotten or latest practice? AFLP applications, analyses and advances. Trends in Plant Science. 2007, 12: 106-117. 10.1016/j.tplants.2007.02.001.View ArticlePubMedGoogle Scholar
- Vos P, Hogers R, Bleeker M, Reijans M, Vandelee T, Hornes M, Frijters A, Pot J, Peleman J, Kuiper M, et al: AFLP - A new technique for DNA-fingerprinting. Nucleic Acids Research. 1995, 23: 4407-4414. 10.1093/nar/23.21.4407.PubMed CentralView ArticlePubMedGoogle Scholar
- Ehrich D: AFLPdat: a collection of R functions for convenient handling of AFLP data. Mol Ecol Notes. 2006, 6: 603-604. 10.1111/j.1471-8286.2006.01380.x.View ArticleGoogle Scholar
- Vekemans X: AFLP-SURV version 1.0. 2002, Distributed by the author: Laboratoire de Genetique et Ecologie Vegetale, Universite Libre de Bruxelles, BelgiumGoogle Scholar
- Holsinger KE, Lewis PO: Hickory: a package for analysis of population genetic data v 1.1. 2007, Distributed by the author: Department of Ecology and Evolutionary Biology, University of Conneticut, Storrs, CTGoogle Scholar
- Foll M, Beaumont MA, Gaggiotti O: An approximate Bayesian computation approach to overcome biases that arise when using amplified fragment length polymorphism markers to study population structure. Genetics. 2008, 179: 927-939. 10.1534/genetics.107.084541.PubMed CentralView ArticlePubMedGoogle Scholar
- Weir BS, Cockerham CC: Estimating F-statistics for the analysis of population structure. Evolution. 1984, 38: 1358-1370. 10.2307/2408641.View ArticleGoogle Scholar
- Excoffier L, Laval G, Schneider S: Arlequin ver 3.0: an integrated software package for population genetics data analysis. Evolutionary Bioinformatics Online. 2005, 1: 47-50.PubMed CentralGoogle Scholar
- Holsinger KE, Lewis PO, Dey DK: A Bayesian approach to inferring population structure from dominant markers. Mol Ecol. 2002, 11: 1157-1164. 10.1046/j.1365-294X.2002.01512.x.View ArticlePubMedGoogle Scholar
- Lynch M, Milligan BG: Analysis of population genetic-structure with RAPD markers. Mol Ecol. 1994, 3: 91-99. 10.1111/j.1365-294X.1994.tb00109.x.View ArticlePubMedGoogle Scholar
- Felsenstein J: PHYLIP (Phylogeny Inference Package), version 3.6. 2004, Distributed by the author: Department of Genome Sciences, University of Washington, Seattle, WashingtonGoogle Scholar
- McCune B, Mefford MM: PC-ORD - Multivariate analysis of ecological data, version 4. 1999, Gleneden Beach, Oregon, USA: MjM Software DesignGoogle Scholar
- Beaumont MA: Adaptation and speciation: what can F-st tell us?. Trends Ecol Evol. 2005, 20: 435-440. 10.1016/j.tree.2005.05.017.View ArticlePubMedGoogle Scholar
- Beaumont MA, Nichols RA: Evaluating loci for use in the genetic analysis of population structure. Proc R Soc Lond B. 1996, 263: 1619-1626. 10.1098/rspb.1996.0237.View ArticleGoogle Scholar
- Beaumont MA, Balding DJ: Identifying adaptive genetic divergence among populations from genome scans. Mol Ecol. 2004, 13: 969-980. 10.1111/j.1365-294X.2004.02125.x.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.