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.
Within the Coregonidae, the postglacial evolution of whitefish (Coregonus clupeaformis (Mitchill) and C. lavaretus (L.) complexes) and the North America ciscoes (C. artedi Lesueur complex) are comparatively well studied [4, 19–23]. Several speciation modes have been suggested for these groups, but evidence for parallel sympatric speciation from a single founder lineage is weak in any of the systems . An alternative system to study putative parallel ecological speciation can be found within populations of the Baltic ciscoes (C. albula (L.) complex), a group that is widespread in northern Europe . Sympatric populations of the C. albula complex occur in Scandinavia and Russia [25, 26] and have been reported from two deep, stratified clear-water lakes  situated geographically close to one another in northern Germany (Fig. 1). In Lake Stechlin (Fig. 1), the autumn-spawning C. albula coexists with the smaller, spring-spawning C. fontanae . Both species are pelagic zooplanktivores [29, 30], but differ substantially with respect to temperature-dependent metabolic adaptations [31, 32]. Analysis of six microsatellite loci revealed very low levels of genetic differentiation between the sympatric populations in Lake Stechlin . In the nearby Lake Breiter Luzin (Fig. 1), two populations also coexist (C. albula and spring-spawning C. lucinensis) , between which the microsatellite differentiation was also very low . While these results suggest a sympatric origin of each spring spawner, shared mtDNA haplotypes revealed lower genetic relatedness between sympatric taxa . An important advance to our understanding of the origin of the sympatric pairs could come from a genetic comparison with other allopatric populations from the same geographical area. Furthermore, the few neutral microsatellite markers employed did not allow for the study of potentially parallel adaptive divergence in both lakes. Accordingly, the specific speciation mode of the sympatric populations has yet to be unambiguously determined.
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
Using a multiplex AFLP method (see Methods) we scored 1264 loci from 271 individuals, of which 1244 (98.4%) were polymorphic over the 10 populations. Differentiation (θ) between the repeated samplings of populations in Lake Stechlin were not significantly different from zero (C. albula: θ = -0.004, p = 0.999; C. fontanae, θ = -0.007, p = 0.99; 50175 iterations in ARLEQUIN), thus samples were pooled for each population in subsequent analyses. The percentage of polymorphic loci varied from 64 - 92% within the 10 populations (Table 1). Estimated with hs, the expected heterozygosity ranged from 0.219 to 0.300, with the highest values being found in the two C. albula populations that occur with the sympatric species (Table 1).
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
Lake names, geographical position and genetic diversity of the 10 Coregonus albula complex populations sampled. Allopatric populations consist exclusively of C. albula. Sympatric populations are named by lake origin and the lake-specific species names. n = number of individuals analysed, % polymorphic loci refers to the total of 1264 loci analysed, hs = gene diversity with credibility intervals calculated by a Bayesian approach as implemented in HICKORY, Hs = the mean within-population expected heterozygosity (= Nei's gene diversity within populations).
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).
All nodes of the NJ tree of Nei's genetic distances were supported by bootstrap values of at least 95% (Fig. 2). The tree grouped four lakes with single populations (Brueckentin, Glambeck, Tollense, Ploener) as distinct from the two lakes with sympatric population pairs (Stechlin, Breiter Luzin) (Fig. 2). Lakes Arend and Behler appeared more closely related to Stechlin and Breiter Luzin. Interestingly, the allopatric spring-spawners, C. fontanae and C. lucinensis, grouped more closely to one another than either did to its sympatric C. albula population (Fig. 2).
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 differentiation (θB) within and among the sympatric population pairs in Lakes Stechlin and Breiter Luzin ranged from 0.028 - 0.113 (Table 2). C. fontanae appeared most distinct, while C. lucinensis was equally differentiated from the three other populations (Table 2). The lowest differentiation observed was between allopatric C. albula in Stechlin and Breiter Luzin (Table 2). This was significantly lower (i.e., outside the 95% credibility interval) than the differentiation with either sympatric taxon (C. fontanae in Lake Stechlin; C. lucinensis in Lake Breiter Luzin.
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
Genetic differentiation (θB with 95% credibility intervals in parentheses) for the two sympatric population pairs of the Coregonus albula complex. Values were calculated with HICKORY based on all 1264 AFLP loci (above diagonals) and after removal of n = 95 loci putatively under selection at p = 0.95 (below diagonals).
The analyses of potential gene flow within the two lakes by STRUCTURE found some evidence for the existence of hybrids between sympatric populations (Fig. 3). The mean proportion (P) of any individual's genome that originated from its own population  was generally ≥ 87% (Lake Stechlin: C. albula, n = 38, P = 96%, C. fontanae, n = 32, P = 94%, Lake Breiter Luzin: C. albula, n = 23, P = 87%, C. lucinensis, n = 27, P = 96%). Three individuals (one C. fontanae and two C. albula from Breiter Luzin) showed elevated probabilities of misclassification (alternative assignment) or hybridization in the F1-generation (Fig. 3). For 23 individuals, there was a high probability of F2- or older hybridizations with the sympatric population. In contrast, the probability of F1-hybrids formed with stocked fish from Lakes Arend and Tollense was zero for all populations, and only six individuals showed a high probability of F2- and older hybrids with stocked fish (Fig. 3).
The explorative genome scans revealed very few loci to be positively exceeding neutral expectations (at p = 0.995) in the two lakes with sympatric population pairs. More loci under putative selection were found in Lake Stechlin (n = 12, ca. 0.95%; Fig. 4a) than in Lake Breiter Luzin (n = 3, ca. 0.24%; Fig. 4b). No single locus was non-neutral in both lakes. At the p = 0.975 quantile, 5.30% of loci (n = 67) were positively non-neutral in Lake Stechlin, and 2.21% (n = 28) in Breiter Luzin. At this quantile, three loci (numbers 233, 879, 1100) were shared from both lakes (Fig. 4a, b).
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
We examined population genetic diversity using amplified fragment length polymorphism (AFLP) in order to generate the large number of characters needed for genomic scans of non-model organisms [4, 60]. We developed a multiplex method of fragment generation in order to increase the number of scorable bands in a given reaction. Total genomic DNA was extracted from pectoral fin tissue using the DNeasy Tissue Kit (Qiagen) according to the manufacturer's protocol. DNA-quality was visually inspected under UV-light on a 1% agarose-gel stained with ethidium bromide. For subsequent analysis, only samples with a clearly visible high molecular weight band were used. AFLP analysis was carried out following a slightly modified protocol of the AFLP-method . Restriction and ligation were carried out in a single step in a thermocycler (2 h at 37°C, 8 h at 16°C). Genomic DNA (300-400 ng DNA per sample) was digested with restriction enzymes MseI (1 unit) and EcoRI (5 units). Polymerase chain reaction (PCR) adaptors specific to the cutting sites were ligated using DNA-Ligase (1 unit; all enzymes from New England Biolabs). Preselective PCR was performed with one selective base on each primer (MseI+C and EcoRI+A; all oligonucleotides from Metabion). PCR product quality was visually confirmed on an agarose gel after preselective amplification. For the selective amplifications, two additional bases were added to the 3' end of each primer (MseI/EcoRIDYE) und the utilization of two different fluorescent dyes permitted multiplexing, where two EcoRI primers were combined with one MseI primer (Table 3).
Primers for the selective amplification step of AFLP
Number of loci per primer pair
All primers had 3 additional bases on the 3' end. EcoRI primers A and B were labelled with different fluorescent dyes (fam and joe) and 3 different multiplex reactions were run with two forward and one backward primer each.
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).
Leibniz-Institute of Freshwater Ecology and Inland Fisheries
ETH Zurich, Forest Ecology
Zoological State Collection
Institute of Epidemiology, Helmholtz Center Munich, German Research Center for Environmental Health
Schluter D: Ecological speciation in postglacial fishes.Philosophical Transactions of the Royal Society of London Series B-Biological Sciences 1996, 351:807–814.View Article
Taylor EB: Species pairs of north temperate freshwater fishes: Evolution, taxonomy, and conservation.Rev Fish Biol Fish 1999, 9:299–324.View Article
Schluter D: The ecology of adaptive radiation. Oxford: Oxford University Press; 2000.
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.PubMedView Article
Schluter D, McPhail JD: Character displacement and replicate adaptive radiation.Trends Ecol Evol 1993, 8:197–200.PubMedView Article
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.View Article
Jonsson B, Jonsson N: Polymorphism and speciation in Arctic charr.J Fish Biol 2001, 58:605–638.View Article
Wood CC, Foote CJ: Evidence for sympatric genetic divergence of anadromous and nonanadromous morphs of sockeye salmon (Oncorhynchus nerka).Evolution 1996, 50:1265–1279.View Article
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.
Hendry AP: Ecological speciation! Or the lack thereof?Can J Fish Aquat Sci 2009, 66:1383–1398.View Article
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.PubMedView Article
Herder F, Pfaender J, Schliewen UK: Adaptive sympatric speciation of polychromatic "roundfin" sailfin silverside fish in Lake Matano (Sulawesi).Evolution 2008, 62:2178–2195.PubMedView Article
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.PubMedView Article
Slatkin M: Isolation by distance in equilibrium and nonequilibrium populations.Evolution 1993, 47:264–279.View Article
Wagner CE, Mccune AR: Contrasting patterns of spatial genetic structure in sympatric rock-dwelling cichlid fishes.Evolution 2009, 63:1312–1326.PubMedView Article
Douglas MR, Brunner PC, Bernatchez L: Do assemblages ofCoregonus(Teleostei : Salmoniformes) in the Central Alpine region of Europe represent species flocks?Mol Ecol 1999, 8:589–603.View Article
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.View Article
Bernatchez L, Chouinard A, Lu GQ: Integrating molecular genetics and ecology in studies of adaptive radiation: whitefish,Coregonussp., as a case study.Biol J Linn Soc 1999, 68:173–194.View Article
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.PubMedView Article
Ostbye K, Bernatchez L, Naesje TF, Himberg KJM, Hindar K: Evolutionary history of the European whitefishCoregonus lavaretus(L.) species complex as inferred from mtDNA phylogeography and gill-raker numbers.Mol Ecol 2005, 14:4371–4387.PubMedView Article
Kottelat M, Freyhof J: Handbook of European freshwater fishes. Cornot and Berlin: Kottelat and Freyhof; 2007.
Sendek DS: Electrophoretic studies of Coregonid fishes from across Russia.Arch Hydrobiol Spec Issues Advanc Limnol 2002, 57:35–55.
Svärdson G: Speciation of the ScandinavianCoregonus.Reports of the Institute of Freshwater Research Drottningholm 1979, 57:1–95.
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.View Article
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.
Helland IP, Harrod C, Freyhof J, Mehner T: Coexistence of a pair of pelagic planktivorous coregonid fish.Evol Ecol Res 2008, 10:373–390.
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.View Article
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,Coregonusspp.Funct Ecol 2008, 22:501–508.View Article
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.
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.View Article
Thienemann A: Coregonus albula lucinensis, eine Tiefenform der Kleinen Maräne aus einem norddeutschen See. Zugleich ein Beitrag zur Rassenbildung beiCoregonus albulaL.Z Morph Ökol Tiere 1933, 27:654–683.View Article
Grant PR, Grant BR, Petren K: Hybridization in the recent past.Am Nat 2005, 166:56–67.PubMedView Article
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.
Douglas MR, Brunner PC: Biodiversity of Central AlpineCoregonus(Salmoniformes): Impact of one-hundred years of management.Ecol Appl 2002, 12:154–172.View Article
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.
Sønstebø JH, Borgstrom R, Heun M: A comparison of AFLPs and microsatellites to identify the population structure of brown trout (Salmo truttaL.) populations from Hardangervidda, Norway.Mol Ecol 2007, 16:1427–1438.PubMedView Article
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.PubMed
Bensch S, Akesson M: Ten years of AFLP in ecology and evolution: why so few animals?Mol Ecol 2005, 14:2899–2914.PubMedView Article
Wilding CS, Butlin RK, Grahame J: Differential gene exchange between parapatric morphs ofLittorina saxatilisdetected using AFLP markers.Journal of Evolutionary Biology 2001, 14:611–619.View Article
Via S, West J: The genetic mosaic suggests a new role for hitchhiking in ecological speciation.Mol Ecol 2008, 17:4334–4345.PubMedView Article
Rogers SM, Bernatchez L: The genetic architecture of ecological speciation and the association with signatures of selection in natural lake whitefish (Coregonussp. Salmonidae) species pairs.Mol Biol Evol 2007, 24:1423–1438.PubMedView Article
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.PubMedView Article
Egan SP, Nosil P, Funk DJ: Selection and genomic differentiation during ecological speciation: Isolating the contributions of host association via a comparative genome scan ofNeochlamisus bebbianaeleaf beetles.Evolution 2008, 62:1162–1181.PubMedView Article
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.PubMedView Article
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.PubMedView Article
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.View Article
Scharf J, Krappe M, Koschel R, Waterstraat A: Feeding of European cisco (Coregonus albulaandC. lucinensis) on the glacial relict crustaceanMysis relictain Lake Breiter Luzin (Germany).Limnologica 2008, 38:147–158.
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.View Article
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.PubMedView Article
Kammerad B, Ellermann S, Mencke J, Wüstemann O, Zuppke U: Die Fischfauna von Sachsen-Anhalt - Verbreitungsatlas. Magdeburg: Ministerium für Raumordnung, Landwirtschaft und Umwelt; 1997.
Brzuzan P, Kozlowski J, Fopp D: Genetic structure of Polish populations of vendace (Coregonus albulaL.) inferred from mitochondrial DNA.Arch Hydrobiol Spec Issues Advanc Limnol 2002, 57:1–10.
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.View Article
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.View Article
Meudt HM, Clarke AC: Almost forgotten or latest practice? AFLP applications, analyses and advances.Trends in Plant Science 2007, 12:106–117.PubMedView Article
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.PubMedView Article
Ehrich D: AFLPdat: a collection of R functions for convenient handling of AFLP data.Mol Ecol Notes 2006, 6:603–604.View Article
Vekemans X: AFLP-SURV version 1.0. Distributed by the author: Laboratoire de Genetique et Ecologie Vegetale, Universite Libre de Bruxelles, Belgium; 2002.
Holsinger KE, Lewis PO: Hickory: a package for analysis of population genetic data v 1.1. Distributed by the author: Department of Ecology and Evolutionary Biology, University of Conneticut, Storrs, CT; 2007.
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.PubMedView Article
Weir BS, Cockerham CC: Estimating F-statistics for the analysis of population structure.Evolution 1984, 38:1358–1370.View Article
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
Holsinger KE, Lewis PO, Dey DK: A Bayesian approach to inferring population structure from dominant markers.Mol Ecol 2002, 11:1157–1164.PubMedView Article
Lynch M, Milligan BG: Analysis of population genetic-structure with RAPD markers.Mol Ecol 1994, 3:91–99.PubMedView Article
Felsenstein J: PHYLIP (Phylogeny Inference Package), version 3.6. Distributed by the author: Department of Genome Sciences, University of Washington, Seattle, Washington; 2004.
McCune B, Mefford MM: PC-ORD - Multivariate analysis of ecological data, version 4. Gleneden Beach, Oregon, USA: MjM Software Design; 1999.
Beaumont MA: Adaptation and speciation: what can F-st tell us?Trends Ecol Evol 2005, 20:435–440.PubMedView Article
Beaumont MA, Nichols RA: Evaluating loci for use in the genetic analysis of population structure.Proc R Soc Lond B 1996, 263:1619–1626.View Article
Beaumont MA, Balding DJ: Identifying adaptive genetic divergence among populations from genome scans.Mol Ecol 2004, 13:969–980.PubMedView Article
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.