Metabarcoding is by now a well‐established method for biodiversity assessment in terrestrial, freshwater, and marine environments. Metabarcoding data sets are usually used for α‐ and β‐diversity estimates, that is, interspecies (or inter‐MOTU [molecular operational taxonomic unit]) patterns. However, the use of hypervariable metabarcoding markers may provide an enormous amount of intraspecies (intra‐MOTU) information—mostly untapped so far. The use of cytochrome oxidase (COI) amplicons is gaining momentum in metabarcoding studies targeting eukaryote richness. COI has been for a long time the marker of choice in population genetics and phylogeographic studies. Therefore, COI metabarcoding data sets may be used to study intraspecies patterns and phylogeographic features for hundreds of species simultaneously, opening a new field that we suggest to name metaphylogeography. The main challenge for the implementation of this approach is the separation of erroneous sequences from true intra‐MOTU variation. Here, we develop a cleaning protocol based on changes in entropy of the different codon positions of the COI sequence, together with co‐occurrence patterns of sequences. Using a data set of community DNA from several benthic littoral communities in the Mediterranean and Atlantic seas, we first tested by simulation on a subset of sequences a two‐step cleaning approach consisting of a denoising step followed by a minimal abundance filtering. The procedure was then applied to the whole data set. We obtained a total of 563 MOTUs that were usable for phylogeographic inference. We used semiquantitative rank data instead of read abundances to perform AMOVAs and haplotype networks. Genetic variability was mainly concentrated within samples, but with an important between seas component as well. There were intergroup differences in the amount of variability between and within communities in each sea. For two species, the results could be compared with traditional Sanger sequence data available for the same zones, giving similar patterns. Our study shows that metabarcoding data can be used to infer intra‐ and interpopulation genetic variability of many species at a time, providing a new method with great potential for basic biogeography, connectivity and dispersal studies, and for the more applied fields of conservation genetics, invasion genetics, and design of protected areas.

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 sympatric European whitefish morphs from three northern Fennoscandian 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.

Natural populations often vary in their degree of ecological, morphological and genetic divergence. This variation can be arranged along an ecological speciation continuum of increasingly discrete variation, with high inter-individual variation at one end and well defined species in the other. In postglacial fishes, evolutionary divergence has commonly resulted in the co-occurrence of a pelagic and a benthic specialist. We studied three replicate lakes supporting sympatric pelagic and benthic European whitefish (Coregonus lavaretus (L.)) morphs in search for early signs of possible further divergence into more specialized niches. Using stomach content data (recent diet) and stable isotope analyses (time-integrated measure of trophic niche use), we observed a split in the trophic niche within the benthic whitefish morph, with individuals specializing on either littoral or profundal resources. This divergence in resource use was accompanied by small but significant differences in an adaptive morphological trait (gill raker number) and significant genetic differences between fish exploiting littoral and profundal habitats and foraging resources. The same pattern of parallel divergence was found in all three lakes, suggesting similar natural selection pressures driving and/or maintaining the divergence. The two levels of divergence (a clear and robust benthic -pelagic and a more subtle littoralprofundal divergence) observed in this study apparently represent different stages in the process of ecological speciation.Electronic supplementary material The online version of this article

Climate change and the intensification of land use practices are causing widespread eutrophication of subarctic lakes. The implications of this rapid change for lake ecosystem function remain poorly understood. To assess how freshwater communities respond to such profound changes in their habitat and resource availability, we conducted a space‐for‐time analysis of food‐web structure in 30 lakes situated across a temperature‐productivity gradient equivalent to the predicted future climate of subarctic Europe (temperature +3°C, precipitation +30% and nutrient +45 μg L−1 total phosphorus). Along this gradient, we observed an increase in the assimilation of pelagic‐derived carbon from 25 to 75% throughout primary, secondary and tertiary consumers. This shift was overwhelmingly driven by the consumption of pelagic detritus by benthic primary consumers and was not accompanied by increased pelagic foraging by higher trophic level consumers. Our data also revealed a convergence of the carbon isotope ratios of pelagic and benthic food web endmembers in the warmest, most productive lakes indicating that the incorporation of terrestrial derived carbon into aquatic food webs increases as land use intensifies. These results, reflecting changes along a gradient characteristic of the predicted future environment throughout the subarctic, indicate that climate and land use driven eutrophication and browning are radically altering the function and fuelling of aquatic food webs in this biome.

The biodiverse Neotropical ecoregion remains insufficiently assessed, poorly managed, and threatened by unregulated human activities. Novel, rapid and cost-effective DNA-based approaches are valuable to improve understanding of the biological communities and for biomonitoring in remote areas. Here, we evaluate the potential of environmental DNA (eDNA) metabarcoding for assessing the structure and distribution of fish communities by analysing sediments and water from 11 locations along the Jequitinhonha River catchment (Brazil). Each site was sampled twice, before and after a major rain event in a five-week period and fish diversity was estimated using high-through-put sequencing of 12S rRNA amplicons. In total, 252 Molecular Operational Taxonomic Units (MOTUs) and 34 fish species were recovered, including endemic, introduced, and previously unrecorded species for this basin. Spatiotemporal variation of fish assemblages was detected, richness during the first campaign was nearly twice as high as in the second sampling round; though peaks of diversity were primarily associated with only four locations. No correlation between β-diversity and longitudinal distance or presence of dams was detected, but low species richness observed at sites located near dams indicates that these anthropogenic barriers might have an impact on local fish diversity. Unexpectedly high α-diversity levels recorded at the river mouth suggest that these sections should be further evaluated as putative "eDNA reservoirs" for rapid monitoring. By uncovering spatio-temporal changes, unrecorded biodiversity components, and putative anthropogenic impacts on fish assemblages, we further strengthen the potential of eDNA metabarcoding as a biomonitoring tool, especially in regions often neglected or difficult to access.

In contrast to the plethora of studies focusing on the genomic basis of adaptive phenotypic divergence, the role of gene expression during speciation has been much less investigated and consequently less understood. Yet, the convergence of differential gene expression patterns between closely related species‐pairs might reflect the role of natural selection during the process of ecological speciation. Here, we test for intercontinental convergence in differential transcriptional signatures between limnetic and benthic sympatric species‐pairs of Lake Whitefish (Coregonus clupeaformis) and its sister lineage, the European Whitefish (Coregonus lavaretus), using six replicated sympatric species‐pairs (two in North America, two in Norway and two in Switzerland). We characterized both sequence variation in transcribed regions and differential gene expression between sympatric limnetic and benthic species across regions and continents. Our first finding was that differentially expressed genes (DEG) between limnetic and benthic whitefish tend to be enriched in shared polymorphism among sister lineages. We then used both genotypes and covariation in expression in order to infer polygenic selection at the gene level. We identified parallel outliers and DEG involving genes primarily overexpressed in limnetic species relative to the benthic species. Our analysis finally revealed the existence of shared genomic bases underlying parallel differential expression across replicated species‐pairs from both continents, such as a cis‐eQTL affecting the pyruvate kinase expression level involved in glycolysis. Our results are consistent with a long‐standing role of natural selection in maintaining trans‐continental diversity at phenotypic traits involved in ecological speciation between limnetic and benthic whitefishes.

Gene flow has tremendous importance for local adaptation, by influencing the fate of de novo mutations, maintaining standing genetic variation and driving adaptive introgression. Furthermore, structural variation as chromosomal rearrangements may facilitate adaptation despite high gene flow. However, our understanding of the evolutionary mechanisms impending or favouring local adaptation in the presence of gene flow is still limited to a restricted number of study systems. In this study, we examined how demographic history, shared ancestral polymorphism, and gene flow among glacial lineages contribute to local adaptation to sea conditions in a marine fish, the capelin (Mallotus villosus). We first assembled a 490‐Mbp draft genome of M. villosus to map our RAD sequence reads. Then, we used a large data set of genome‐wide single nucleotide polymorphisms (25,904 filtered SNPs) genotyped in 1,310 individuals collected from 31 spawning sites in the northwest Atlantic. We reconstructed the history of divergence among three glacial lineages and showed that they probably diverged from 3.8 to 1.8 million years ago and experienced secondary contacts. Within each lineage, our analyses provided evidence for large Ne and high gene flow among spawning sites. Within the Northwest Atlantic lineage, we detected a polymorphic chromosomal rearrangement leading to the occurrence of three haplogroups. Genotype–environment associations revealed molecular signatures of local adaptation to environmental conditions prevailing at spawning sites. Our study also suggests that both shared polymorphisms among lineages, resulting from standing genetic variation or introgression, and chromosomal rearrangements may contribute to local adaptation in the presence of high gene flow.

Summary Many species are expanding their distribution towards higher latitudes and altitudes in response to climate change. These range shifts are expected to change fish community structure and alter food‐web dynamics in subarctic lakes. However, the impacts of invading species on native fish and invertebrate prey communities remain understudied. The trophic ecology of invasive species determines the likelihood of direct resource competition with native taxa. In Northern Europe, perch (Perca fluviatilis), a trophic generalist, and ruffe (Gymnocephalus cernuus), a benthic specialist, are expanding their distribution ranges northwards, colonising lakes inhabited by a native generalist, whitefish (Coregonus lavaretus). We predicted that increased fish diversity and density would deplete the invertebrate community and increase resource competition between native and invasive species. To assess the degree of resource competition between native and invasive species, we compared (i) fish and invertebrate community structure; (ii) diet and stable carbon and nitrogen ratios of whitefish, ruffe and perch and (iii) growth, condition and relative population size of whitefish in two non‐invaded lakes with two lakes containing one and two lakes containing both invasive species. Each lake was sampled on a single occasion between August and September. Benthic macroinvertebrate density and community structure were unaffected by increased consumer diversity, while top‐down control of pelagic zooplankton density and size was evident in lakes with increased fish diversity. Differences in diet and stable isotope ratios were evident between all whitefish populations, although these were not directly related to the presence of invasive species. Specialised adaptations of invasive species may confer a competitive advantage in invaded lakes; ruffe dominated the profundal niche, while perch displayed an ontogenetic shift to piscivory, reducing niche overlap with native whitefish. Growth rate and population density of whitefish were largely independent of fish community structure and were governed by local variations in lake productivity. However, there was a sign of lowered condition of whitefish in invaded systems. Shallow and more productive lakes with higher food availability supported populations of native and invasive species. Our findings indicate that trophic specialisations of invasive species play a key role in determining their impacts on the systems they invade. This study focussed on early stages of invasion, and the outcome of species interactions may change following the establishment of new species. In addition, these impacts will not be uniform across the invaded landscape as lake‐specific variations in morphometry and resource availability will alter the competitive balance between native and invasive species.