Understanding how a monophyletic lineage of a species diverges into several adaptive forms has received increased attention in recent years, but the underlying mechanisms in this process are still under debate. Postglacial fishes are excellent model organisms for exploring this process, especially the initial stages of ecological speciation, as postglacial lakes represent replicated discrete environments with variation in available niches. Here, we combine data of niche utilization, trophic morphology, and 17 microsatellite loci to investigate the diversification process of three sympatricEuropean whitefish morphs from three northernFennoscandian lakes. The morphological divergence in the gill raker number among the whitefish morphs was related to the utilization of different trophic niches and was associated with reproductive isolation within and across lakes. The intralacustrine comparison of whitefish morphs showed that these systems represent two levels of adaptive divergence: (1) a consistent littoral–pelagic resource axis; and (2) a more variable littoral–profundal resource axis. The results also indicate that the profundal whitefish morph has diverged repeatedly from the ancestral littoral whitefish morph in sympatry in two different watercourses. In contrast, all the analyses performed revealed clustering of the pelagic whitefish morphs across lakes suggesting parallel postglacial immigration with the littoral whitefish morph into each lake. Finally, the analyses strongly suggested that the trophic adaptive trait, number of gill rakers, was under diversifying selection in the different whitefish morphs. Together, the results support a complex evolutionary scenario where ecological speciation acts, but where both allopatric (colonization history) and sympatric (within watercourse divergence) processes are involved.

Environmental DNA (eDNA) metabarcoding is a relatively new monitoring tool featuring in an increasing number of applications such as the facilitation of the accurate and cost effective detection of species in environmental samples. eDNA monitoring is likely to have a major impact on the ability of salmonid aquaculture industry producers and their regulators to detect the presence and abundance of pathogens and other biological threats in the surrounding environment. However, for eDNA metabarcoding to develop into a useful bio-monitoring tool it is necessary to (a) validate that sequence datasets derived from amplification of metabarcoding markers reflect the true species’ identity, (b) test the sensitivity under different abundance levels and environmental noise and (c) establish a low-cost sequencing method to enable the bulk processing of field samples. In this study, we employed an elaborate experimental design whereby different combinations of five biological agents were crossed at three abundance levels and exposed to sterile pre-filtered and unfiltered seawater, prior to coarse filtering and then eDNA ultrafiltration of the resultant material. We then benchmarked the low-cost, scalable, Ion Torrent sequencing method against the current gold-standard Illumina platform for eDNA surveys in aquaculture. Based on amplicon-seq of the 18S SSU rDNA v9 region, we were able to identify two parasites (Lepeophtheirus salmonis and Paramoeba perurans) to species level, whereas the microalgae species Prymnesium parvum, Pseudo-nitzschia seriata, and P. delicatissima could be assigned correctly only to the genus level. Illumina and Ion Torrent provided near identical results in terms of community composition in our samples, whereas Ion Torrent was more sensitive in detecting species richness when the medium was unfiltered seawater. Both methods were able to reflect the difference in relative abundance between treatments in 4 out of 5 species when samples were exposed to the unfiltered seawater, despite the significant amount of background noise from both bacteria and eukaryotes. Our findings indicate that eDNA metabarcoding offers significant potential in the monitoring of species harmful to aquaculture and for this purpose, the low-cost Ion Torrent sequencing is as accurate as Illumina in determining differences in their relative abundance between samples.

Lake‐dwelling fish that form species pairs/flocks characterized by body size divergence are important model systems for speciation research. Although several sources of divergent selection have been identified in these systems, their importance for driving the speciation process remains elusive. A major problem is that in retrospect, we cannot distinguish selection pressures that initiated divergence from those acting later in the process. To address this issue, we studied the initial stages of speciation in European whitefish (Coregonus lavaretus) using data from 358 populations of varying age (26–10,000 years). We find that whitefish speciation is driven by a large‐growing predator, the northern pike (Esox lucius). Pike initiates divergence by causing a largely plastic differentiation into benthic giants and pelagic dwarfs: ecotypes that will subsequently develop partial reproductive isolation and heritable differences in gill raker number. Using an eco‐evolutionary model, we demonstrate how pike's habitat specificity and large gape size are critical for imposing a between‐habitat trade‐off, causing prey to mature in a safer place or at a safer size. Thereby, we propose a novel mechanism for how predators may cause dwarf/giant speciation in lake‐dwelling fish species.

Invasion of exotic species has caused the loss of biodiversity and imparts evolutionary and ecological changes in the introduced systems. In northern Fennoscandia, European whitefish (Coregonus lavaretus (L.)) is a highly polymorphic species displaying adaptive radiations into partially reproductively isolated and thus genetically differentiated sympatric morphs utilizing the planktivorous and benthivorous food niche in many lakes. In 1993, Lake Skrukkebukta was invaded by vendace (Coregonus albula (L.)) which is a zooplanktivorous specialist. The vendace displaced the densely rakered whitefish from its preferred pelagic niche to the benthic habitat harbouring the large sparsely rakered whitefish. In this study, we investigate the potential influence of the vendace invasion on the breakdown of reproductive isolation between the two whitefish morphs. We inferred the genotypic and phenotypic differentiation between the two morphs collected at the arrival (1993) and 15 years after (2008) the vendace invasion using 16 microsatellite loci and gill raker numbers, the most distinctive adaptive phenotypic trait between them. The comparison of gill raker number distributions revealed two modes growing closer over 15 years following the invasion. Bayesian analyses of genotypes revealed that the two genetically distinct whitefish morphs that existed in 1993 had collapsed into a single population in 2008. The decline in association between the gill raker numbers and admixture values over 15 years corroborates the findings from the Bayesian analysis. Our study thus suggests an apparent decrease of reproductive isolation in a morph-pair of European whitefish within 15 years (≃ 3 generations) following the invasion of a superior trophic competitor (vendace) in a subarctic lake, reflecting a situation of “speciation in reverse”.