Ecological speciation in postglacial European whitefish: rapid adaptive radiations into the littoral, pelagic, and profundal lake habitats

Foraging-induced craniofacial plasticity is associated with an early, robust and dynamic transcriptional response

Phenotypic plasticity is the ability of a single genotype to vary its phenotype in response to the environment. Plasticity of the skeletal system in response to mechanical input is widely studied, but the timing of its transcriptional regulation is not well understood. Here, we used the cichlid feeding apparatus to examine the transcriptional dynamics of skeletal plasticity over time. Using three closely related species that vary in their ability to remodel bone and a panel of 11 genes, including well-studied skeletal differentiation markers and newly characterized environmentally sensitive genes, we examined plasticity at one, two, four and eight weeks following the onset of alternate foraging challenges. We found that the plastic species exhibited environment-specific bursts in gene expression beginning at one week, followed by a sharp decline in levels, while the species with more limited plasticity exhibited consistently low levels of gene expression. This trend held across nearly all genes, suggesting that it is a hallmark of the larger plasticity regulatory network. We conclude that plasticity of the cichlid feeding apparatus is not the result of slowly accumulating gene expression difference over time, but rather is stimulated by early bursts of environment-specific gene expression followed by a return to homeostatic levels.

Environmental correlates of adaptive diversification in postglacial freshwater fishes

Determining how environmental conditions contribute to divergence among populations and drive speciation is fundamental to resolving mechanisms and understanding outcomes in evolutionary biology. Postglacial freshwater fish species in the Northern Hemisphere are ideal biological systems to explore the effects of environment on diversification in morphology, ecology, and genetics (ecomorph divergences) within lakes. To date, various environmental factors have been implicated in the presence of multiple ecomorphs within particular lakes or regions. However, concerted evidence for generalizable patterns in environmental variables associated with speciation across geographical regions and across species and genera has been lacking. Here, we aimed to identify key biotic and abiotic factors associated with ecological divergence of postglacial freshwater fish species into multiple sympatric ecomorphs, focusing on species in the well-studied, widespread, and co-distributed genera Gasterosteus, Salvelinus, and Coregonus (stickleback, charr, and whitefish, respectively). We found that the presence of multiple sympatric ecomorphs tended to be associated with increasing lake surface area, maximum depth, and nutrient availability. In addition, predation, competition, and prey availability were suggested to play a role in divergence into multiple ecomorphs, but the effects of biotic factors require further study. Although we identified several environmental factors correlated with the presence of multiple ecomorphs, there were substantial data gaps across species and regions. An improved understanding of these systems may provide insight into both generalizable environmental factors involved in speciation in other systems, and potential ecological and evolutionary responses of species complexes when these variables are altered by environmental change.

Ecological Drivers of Mercury Bioaccumulation in Fish of a Subarctic Watercourse

Mercury (Hg) is a serious concern for aquatic ecosystems because it may biomagnify to harmful concentrations within food webs and consequently end up in humans that eat fish. However, the trophic transfer of mercury through the aquatic food web may be impacted by several factors related to network complexity and the ecology of the species present. The present study addresses the interplay between trophic ecology and mercury contamination in the fish communities of two lakes in a pollution-impacted subarctic watercourse, exploring the role of both horizontal (feeding habitat) and vertical (trophic position) food web characteristics as drivers for the Hg contamination in fish. The lakes are located in the upper and lower parts of the watercourse, with the lower site located closer to, and downstream from, the main pollution source. The lakes have complex fish communities dominated by coregonids (polymorphic whitefish and invasive vendace) and several piscivorous species. Analyses of habitat use, stomach contents, and stable isotope signatures (δ¹⁵ N, δ¹³ C) revealed similar food web structures in the two lakes except for a few differences chiefly related to ecological effects of the invasive vendace. The piscivores had higher Hg concentrations than invertebrate-feeding fish. Concentrations increased with size and age for the piscivores and vendace, whereas habitat differences were of minor importance. Most fish species showed significant differences in Hg concentrations between the lakes, the highest values typically found in the downstream site where the biomagnification rate also was higher. Mercury levels in piscivorous fish included concentrations that exceed health authorization limits, with possible negative implications for fishing and human consumption. Our findings accentuate the importance of acquiring detailed knowledge of the drivers that can magnify Hg concentrations in fish and how these may vary within and among aquatic systems, to provide a scientific basis for adequate management strategies. Environ Toxicol Chem 2023;42:873-887. © 2023 SETAC.

A taxonomic revision of ten whitefish species from the lakes Lucerne, Sarnen, Sempach and Zug, Switzerland, with descriptions of seven new species (Teleostei, Coregonidae)

The taxonomy of the endemic whitefish of the lakes of the Reuss River system (Lucerne, Sarnen, Zug) and Lake Sempach, Switzerland, is reviewed and revised. Lake Lucerne harbours five species. Coregonus intermundiasp. nov. and C. suspensussp. nov., are described. Coregonus nobilis Haack, 1882, C. suidteri Fatio, 1885, and C. zugensis Nüsslin, 1882, are redescribed. Genetic studies have shown that C. suidteri and C. zugensis are composed of several distinct species endemic to different lakes. The names C. suidteri and C. zugensis are restricted to the species of lakes Sempach and Zug, respectively. The whitefish populations previously referred to as C. suidteri and C. zugensis from Lake Lucerne are described as C. litoralissp. nov. and C. muellerisp. nov., respectively. Furthermore, the whitefish from Lake Zug that were previously referred to as C. suidteri are described as C. supersumsp. nov. A holotype is designated for C. supersum that was previously one of two syntypes of C. zugensis. The other syntype is retained for C. zugensis. Coregonus obliterussp. nov. is described from Lake Zug, and C. obliterus and C. zugensis from Lake Zug are extinct. Finally, we describe C. sarnensissp. nov. from lakes Sarnen and Alpnach. Coregonus suidteri from Lake Sempach shows strong signals of introgression from deliberately translocated non-native whitefish species, which questions if the extant population still carries a genetic legacy from the original species and thus may need to be considered extinct. Coregonus suspensus is genetically partially of allochthonous origin, closely related to the radiation of Lake Constance. It is therefore compared to all known and described species of Lake Constance: C. wartmanni Bloch, 1784, C. macrophthalmus Nüsslin, 1882, C. arenicolus Kottelat,1997, and C. gutturosus Gmelin, 1818.

The role of vital dietary biomolecules in eco-evo-devo dynamics

The physiological dependence of animals on dietary intake of vitamins, amino acids, and fatty acids is ubiquitous. Sharp differences in the availability of these vital dietary biomolecules among different resources mean that consumers must adopt a range of strategies to meet their physiological needs. We review the emerging work on omega-3 long-chain polyunsaturated fatty acids, focusing predominantly on predator-prey interactions, to illustrate that trade-off between capacities to consume resources rich in vital biomolecules and internal synthesis capacity drives differences in phenotype and fitness of consumers. This can then feedback to impact ecosystem functioning. We outline how focus on vital dietary biomolecules in eco-eco-devo dynamics can improve our understanding of anthropogenic changes across multiple levels of biological organization.

Genomic divergence, local adaptation, and complex demographic history may inform management of a popular sportfish species complex

The Neosho Bass (Micropterus velox), a former subspecies of the keystone top-predator and globally popular Smallmouth Bass (M. dolomieu), is endemic and narrowly restricted to small, clear streams of the Arkansas River Basin in the Central Interior Highlands (CIH) ecoregion, USA. Previous studies have detected some morphological, genetic, and genomic differentiation between the Neosho and Smallmouth Basses; however, the extent of neutral and adaptive divergence and patterns of intraspecific diversity are poorly understood. Furthermore, lineage diversification has likely been impacted by gene flow in some Neosho populations, which may be due to a combination of natural biogeographic processes and anthropogenic introductions. We assessed: (1) lineage divergence, (2) local directional selection (adaptive divergence), and (3) demographic history among Smallmouth Bass populations in the CIH using population genomic analyses of 50,828 single-nucleotide polymorphisms (SNPs) obtained through ddRAD-seq. Neosho and Smallmouth Bass formed monophyletic clades with 100% bootstrap support. We identified two major lineages within each species. We discovered six Neosho Bass populations (two nonadmixed and four admixed) and three nonadmixed Smallmouth Bass populations. We detected 29 SNPs putatively under directional selection in the Neosho range, suggesting populations may be locally adapted. Two populations were admixed via recent asymmetric secondary contact, perhaps after anthropogenic introduction. Two other populations were likely admixed via combinations of ancient and recent processes. These species comprise independently evolving lineages, some having experienced historical and natural admixture. These results may be critical for management of Neosho Bass as a distinct species and may aid in the conservation of other species with complex biogeographic histories.

Allopatric origin of sympatric whitefish morphs with insights on the genetic basis of their reproductive isolation

The European whitefish (Coregonus lavaretus) species complex is a classic example of recent adaptive radiation. Here, we examine a whitefish population introduced to northern Finnish Lake Tsahkal in the late 1960s, where three divergent morphs (viz. littoral, pelagic, and profundal feeders) were found 10 generations after. Using demographic modeling based on genomic data, we show that whitefish morphs evolved during a phase of strict isolation, refuting a rapid sympatric divergence scenario. The lake is now an artificial hybrid zone between morphs originated in allopatry. Despite their current syntopy, clear genetic differentiation remains between two of the three morphs. Using admixture mapping, we identify five SNPs associated with gonad weight variation, a proxy for sexual maturity and spawning time. We suggest that ecological adaptations in spawning time evolved in allopatry are currently maintaining partial reproductive isolation in the absence of other barriers to gene flow.

Complex and divergent histories gave rise to genome-wide divergence patterns amongst European whitefish (Coregonus lavaretus)

Pleistocene glaciations dramatically affected species distribution in regions that were impacted by ice cover and subsequent postglacial range expansion impacted contemporary biodiversity in complex ways. The European whitefish, Coregonus lavaretus, is a widely distributed salmonid fish species on mainland Europe, but in Britain it has only seven native populations, all of which are found on the western extremes of the island. The origins and colonization routes of the species into Britain are unknown but likely contributed to contemporary genetic patterns and regional uniqueness. Here, we used up to 25,751 genome‐wide polymorphic loci to reconstruct the history and to discern the demographic and evolutionary forces underpinning divergence between British populations. Overall, we found lower genetic diversity in Scottish populations but high differentiation (F_ST = 0.433–0.712) from the English/Welsh and other European populations. Differentiation was elevated genome‐wide rather than in particular genomic regions. Demographic modelling supported a postglacial colonization into western Scotland from northern refugia and a separate colonization route for the English/Welsh populations from southern refugia, with these two groups having been separated for more than ca. 50 Ky. We found cyto‐nuclear discordance at a European scale, with the Scottish populations clustering closely with Baltic population in the mtDNA analysis but not in the nuclear data, and with the Norwegian and Alpine populations displaying the same mtDNA haplotype but being distantly related in the nuclear tree. These findings suggest that neutral processes, primarily drift and regionally distinct pre‐glacial evolutionary histories, are important drivers of genomic divergence in British populations of European whitefish. This sheds new light on the establishment of the native British freshwater fauna after the last glacial maximum.

Rapid adaptive radiation in a hillstream cyprinid fish in the East African White Nile River basin

Adaptive radiation of freshwater fishes was long thought to be possible only in lacustrine environments. Recently, several studies have shown that riverine and stream environments also provide the ecological opportunity for adaptive radiation. In this study, we report on a riverine adaptive radiation of six ecomorphs of cyprinid hillstream fishes of the genus Garra in a river located in the Ethiopian Highlands in East Africa. Garra are predominantly highly specialized algae-scrapers with a wide distribution ranging from Southeast Asia to West Africa. However, adaptive phenotypic diversification in mouth type, sucking disc morphology, gut length and body shape have probably been found among these ecomorphs in a single Ethiopian river. Moreover, we found two novel phenotypes of Garra ("thick-lipped" and "predatory") that had not been discovered before in this species-rich genus (>160 species). Mitochondrial and genome-wide data suggest monophyletic, intrabasin evolution of Garra phenotypic diversity with signatures of gene flow from other local populations. Although sympatric ecomorphs are genetically distinct and can be considered to being young species as suggested by genome-wide single nucleotide polymorphism data, mitochondrial DNA was unable to identify any genetic structure suggesting recent and rapid speciation events. Some data suggest a hybrid origin of the novel "thick-lipped" ecomorph. Here we highlight how, driven by ecological opportunity, an ancestral trophically highly specialized lineage is likely to have rapidly radiated in a riverine environment promoted by the evolution of novel feeding strategies.

Contrasting Host-Parasite Population Structure: Morphology and Mitogenomics of a Parasitic Flatworm on Pelagic Deepwater Cichlid Fishes from Lake Tanganyika

Little phylogeographic structure is presumed for highly mobile species in pelagic zones. Lake Tanganyika is a unique ecosystem with a speciose and largely endemic fauna famous for its remarkable evolutionary history. In bathybatine cichlid fishes, the pattern of lake-wide population differentiation differs among species. We assessed the congruence between the phylogeographic structure of bathybatine cichlids and their parasitic flatworm Cichlidogyrus casuarinus to test the magnifying glass hypothesis. Additionally, we evaluated the use of a PoolSeq approach to study intraspecific variation in dactylogyrid monogeneans. The lake-wide population structure of C. casuarinus ex Hemibates stenosoma was assessed based on a portion of the cox1 gene combined with morphological characterisation. Additionally, intraspecific mitogenomic variation among 80 parasite samples from one spatially constrained metapopulation was assessed using shotgun NGS. While no clear geographic genetic structure was detected in parasites, both geographic and host-related phenotypic variation was apparent. The incongruence with the genetic north-south gradient observed in H. stenosoma may be explained by the broad host range of this flatworm including eupelagic bathybatine host species that form panmictic populations across the lake. In addition, we present the first parasite mitogenome from Lake Tanganyika and propose a methodological framework for studying the intraspecific mitogenomic variation of dactylogyrid monogeneans.

Among-individual diet variation within a lake trout ecotype: Lack of stability of niche use

In a polyphenic species, differences in resource use are expected among ecotypes, and homogeneity in resource use is expected within an ecotype. Yet, using a broad resource spectrum has been identified as a strategy for fishes living in unproductive northern environments, where food is patchily distributed and ephemeral. We investigated whether specialization of trophic resources by individuals occurred within the generalist piscivore ecotype of lake trout from Great Bear Lake, Canada, reflective of a form of diversity. Four distinct dietary patterns of resource use within this lake trout ecotype were detected from fatty acid composition, with some variation linked to spatial patterns within Great Bear Lake. Feeding habits of different groups within the ecotype were not associated with detectable morphological or genetic differentiation, suggesting that behavioral plasticity caused the trophic differences. A low level of genetic differentiation was detected between exceptionally large-sized individuals and other piscivore individuals. We demonstrated that individual trophic specialization can occur within an ecotype inhabiting a geologically young system (8,000-10,000 yr BP), a lake that sustains high levels of phenotypic diversity of lake trout overall. The characterization of niche use among individuals, as done in this study, is necessary to understand the role that individual variation can play at the beginning of differentiation processes.

Plasticity and genetic basis of cichlid gill arch anatomy reveal novel roles for Hedgehog signaling

Teleost gill arches are exquisitely evolved to maximize foraging efficiency, and include structures for the capture, filtering, and processing of prey. While both plasticity and a genetic basis for gill arch traits have been noted, the relative contributions of genetics and the environment in shaping these structures remains poorly understood. East African cichlids are particularly useful in this line of study due to their highly diverse and plastic feeding apparatus. Here we explore the gene-by-environmental effects on cichlid GRs by rearing pure bred species and their F₃ hybrids in different foraging environments. We find that anatomical differences between species are dependent on the environment. The genetic architecture of these traits is also largely distinct between foraging environments. We did, however, note a few genomic "hotspots" where multiple traits map to a common region. One of these, for GR number across multiple arches, maps to the ptch1 locus, a key component of the Hedgehog (Hh) pathway that has previously been implicated in cichlid oral jaw shape and plasticity. Since Hh signalling has not previously been implicated in GR development, we explored functional roles for this pathway. Using a small molecule inhibitor in cichlids, as well as zebrafish transgenic systems, we demonstrate that Hh levels negatively regulate GR number, and are both necessary and sufficient to maintain plasticity in this trait. In all these data underscore the critical importance of the environment in determining the relationship between genotype and phenotype, and provide a molecular inroad to better understand the origins of variation in this important foraging-related trait.

Evolutionary determinism and convergence associated with water-column transitions in marine fishes

Repeatable, convergent outcomes are prima facie evidence for determinism in evolutionary processes. Among fishes, well-known examples include microevolutionary habitat transitions into the water column, where freshwater populations (e.g., sticklebacks, cichlids, and whitefishes) recurrently diverge toward slender-bodied pelagic forms and deep-bodied benthic forms. However, the consequences of such processes at deeper macroevolutionary scales in the marine environment are less clear. We applied a phylogenomics-based integrative, comparative approach to test hypotheses about the scope and strength of convergence in a marine fish clade with a worldwide distribution (snappers and fusiliers, family Lutjanidae) featuring multiple water-column transitions over the past 45 million years. We collected genome-wide exon data for 110 (∼80%) species in the group and aggregated data layers for body shape, habitat occupancy, geographic distribution, and paleontological and geological information. We also implemented approaches using genomic subsets to account for phylogenetic uncertainty in comparative analyses. Our results show independent incursions into the water column by ancestral benthic lineages in all major oceanic basins. These evolutionary transitions are persistently associated with convergent phenotypes, where deep-bodied benthic forms with truncate caudal fins repeatedly evolve into slender midwater species with furcate caudal fins. Lineage diversification and transition dynamics vary asymmetrically between habitats, with benthic lineages diversifying faster and colonizing midwater habitats more often than the reverse. Convergent ecological and functional phenotypes along the benthic-pelagic axis are pervasive among different lineages and across vastly different evolutionary scales, achieving predictable high-fitness solutions for similar environmental challenges, ultimately demonstrating strong determinism in fish body-shape evolution.

A taxonomic revision of the whitefish of lakes Brienz and Thun, Switzerland, with descriptions of four new species (Teleostei, Coregonidae)

The alpha taxonomy of the endemic whitefish of lakes Brienz and Thun, Switzerland, is revised. We evaluate the status of seven known species: Coregonus steinmanni sp. nov., Coregonus profundus sp. nov. and Coregonus acrinasus sp. nov. are endemic to Lake Thun; Coregonus Brienzii sp. nov. is endemic to Lake Brienz; and C. alpinus, C. albellus, and C. fatioi from lakes Brienz and Thun are redescribed. One of these species, C. alpinus, is revised, since the lectotype for this species is incongruent with the species description given by Kottelat (1997) and Kottelat and Freyhof (2007). The name C. alpinus is thus retained for the lectotype designated by Kottelat (1997) and a new description of this taxon provided. For the species otherwise described by Kottelat (1997) and Kottelat and Freyhof (2007) as C. alpinus the new name C. profundus is designated. Coregonus acrinasus is genetically partially of allochthonous origin, closely related to the radiation of Lake Constance, and we therefore compare it to all recognized species of Lake Constance, C. wartmanni, C. macrophthalmus, C. arenicolus, and C. gutturosus.

Isotopic niche size of Coregonus artedi (sensu lato) Increases in the presence of Mysis diluviana, expanded habitat use and phenotypic diversity

Post-glacial colonization of lakes in Algonquin Park, Ontario, Canada resulted in food webs with cisco (Coregonus artedi sensu lato) and either Mysis diluviana or Chaoborus spp. as the dominant diel migrator. Mysis as prey, its diel movements and benthic occupancy, are hypothesized to be key elements of ecological opportunity for cisco diversity in the Laurentian Great Lakes. If correct, the hypothesis strongly implies that lakes with Mysis would have greater trophic niche size and drive greater adaptive radiation of cisco forms relative to lakes without Mysis. The dichotomy in diel migrator in Algonquin Park lakes was an opportunity to assess the isotopic niche size of cisco (δ15N and δ13C) and determine if niche size expands with Mysis presence. We found the presence of Mysis is necessary to expand isotopic niche size in our study lakes. The use of habitats not typically associated with the ancestral form of cisco (e.g., benthic habitats) and phenotypic diversity (blackfin and cisco) also continue to expand niche size in Mysis-based food webs. Partial ecological speciation based on a large niche space appears to be present in one lake (Cauchon Lake) where use of alternative habitats is the only real difference in cisco. The presence of blackfin expands niche space in Cedar and Radiant Lakes. This was not matched in Hogan Lake where niche space was relatively smaller with similar forms. Possible reasons for this discrepancy may be related to the asymmetric basin of Hogan Lake and whether the two forms overlap during cool and cold-water periods of the annual temperature cycle. By comparing trophic niche size among lakes with and without Mysis, we conclude that Mysis provides a key ecological opportunity for cisco diversity in our study lakes and likely more widely.

Geographic hierarchical population genetic structuring in British European whitefish (Coregonus lavaretus) and its implications for conservation

The European whitefish Coregonus lavaretus complex represents one of the most diverse radiations within salmonids, with extreme morphological and genetic differentiation across its range. Such variation has led to the assignment of many populations to separate species. In Great Britain, the seven native populations of C. lavaretus (two in Scotland, four in England, one in Wales) were previously classified into three species, and recent taxonomic revision resurrected the previous nomenclature. Here we used a dataset of 15 microsatellites to: (1) investigate the genetic diversity of British populations, (2) assess the level of population structure and the relationships between British populations. Genetic diversity was highest in Welsh (HO = 0.50, AR = 5.29), intermediate in English (HO = 0.41–0.50, AR = 2.83–3.88), and lowest in Scottish populations (HO = 0.28–0.35, AR = 2.56–3.04). Population structure analyses indicated high genetic differentiation (global FST = 0.388) between all populations but for the two Scottish populations (FST = 0.063) and two English populations (FST = 0.038). Principal component analysis and molecular ANOVA revealed separation between Scottish, English, and Welsh populations, with the Scottish populations being the most diverged. We argue that the data presented here are not sufficient to support a separation of the British European whitefish populations into three separate species, but support the delineation of different ESUs for these populations.

Pelagic-benthic resource polymorphism in Schizopygopsis thermalis Herzenstein 1891 (Pisces, Cyprinidae) in a headwater lake in the Salween River system on the Tibetan Plateau

Resource polymorphism is a ubiquitous phenomenon in vertebrates and may represent a critical intermediate stage in speciation. Freshwater lakes in high-altitude areas represent a natural system for understanding resource polymorphism in fishes benefiting from diverse lacustrine environments and species-poor fish assemblages. We report resource polymorphism in a cyprinid fish, Schizopygopsis thermalis, in Lake Amdo Tsonak Co, a headwater lake in the upper Salween River system. Two discrete intraspecific morphs, planktivorous and benthivorous, were identified according to geometric morphometrics and traditional univariate linear measures. The planktivorous morph exhibits a longer head, longer upper and lower jaw, larger asymptotic standard length (L∞ ), lower growth rate (k), and higher growth performance index (φ) than the benthivorous morph. With respect to descriptive traits, the planktivorous morph possesses a large, terminal mouth and obvious mucus pores on the cheek and chin, while the benthivorous morph is characterized by a more inferior mouth with a sharpen horny edge on the lower jaw and unconspicuous mucus pores. The discrete pelagic-benthic resources and low interspecific competition in the lake system might drive the initial differentiation of the two morphs, and partial spatial reproductive isolation in breeding further maintains and reinforces the differences between them.

Ecological speciation in European whitefish is driven by a large-gaped predator

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.

Genotyping-by-sequencing illuminates high levels of divergence among sympatric forms of coregonines in the Laurentian Great Lakes

Effective resource management depends on our ability to partition diversity into biologically meaningful units. Recent evolutionary divergence, however, can often lead to ambiguity in morphological and genetic differentiation, complicating the delineation of valid conservation units. Such is the case with the "coregonine problem," where recent postglacial radiations of coregonines into lacustrine habitats resulted in the evolution of numerous species flocks, often with ambiguous taxonomy. The application of genomics methods is beginning to shed light on this problem and the evolutionary mechanisms underlying divergence in these ecologically and economically important fishes. Here, we used restriction site-associated DNA (RAD) sequencing to examine genetic diversity and differentiation among sympatric forms in the Coregonus artedi complex in the Apostle Islands of Lake Superior, the largest lake in the Laurentian Great Lakes. Using 29,068 SNPs, we were able to clearly distinguish among the three most common forms for the first time, as well as identify putative hybrids and potentially misidentified specimens. Population assignment rates for these forms using our RAD data were 93%-100% with the only mis-assignments arising from putative hybrids, an improvement from 62% to 77% using microsatellites. Estimates of pairwise differentiation (F ST: 0.045-0.056) were large given the detection of hybrids, suggesting that reduced fitness of hybrid individuals may be a potential mechanism for the maintenance of differentiation. We also used a newly built C. artedi linkage map to look for islands of genetic divergence among forms and found widespread differentiation across the genome, a pattern indicative of long-term drift, suggesting that these forms have been reproductively isolated for a substantial amount of time. The results of this study provide valuable information that can be applied to develop well-informed management strategies and stress the importance of re-evaluating conservation units with genomic tools to ensure they accurately reflect species diversity.

Using mathematical modelling to investigate the adaptive divergence of whitefish in Fennoscandia

Modern speciation theory has greatly benefited from a variety of simple mathematical models focusing on the conditions and patterns of speciation and diversification in the presence of gene flow. Unfortunately the application of general theoretical concepts and tools to specific ecological systems remains a challenge. Here we apply modeling tools to better understand adaptive divergence of whitefish during the postglacial period in lakes of northern Fennoscandia. These lakes harbor up to three different morphs associated with the three major lake habitats: littoral, pelagic, and profundal. Using large-scale individual-based simulations, we aim to identify factors required for in situ emergence of the pelagic and profundal morphs in lakes initially colonized by the littoral morph. The importance of some of the factors we identify and study - sufficiently large levels of initial genetic variation, size- and habitat-specific mating, sufficiently large carrying capacity of the new niche - is already well recognized. In addition, our model also points to two other factors that have been largely disregarded in theoretical studies: fitness-dependent dispersal and strong predation in the ancestral niche coupled with the lack of it in the new niche(s). We use our theoretical results to speculate about the process of diversification of whitefish in Fennoscandia and to identify potentially profitable directions for future empirical research.

Long-Term Environmental Monitoring in an Arctic Lake Polluted by Metals under Climate Change

Lake Kuetsjarvi (in the lower reaches of the Pasvik River, Murmansk Region, Russia) in the border area between Russia and Norway, is one of the most polluted water reservoirs in the European Arctic. The operation of the Pechenganikel Smelter located on its shores has led to the extremely high concentrations of heavy metals observed in the waters and sediments of the lake. Long-term comprehensive studies of the ecosystem of Lake Kuetsjarvi have made it possible to identify the response of its components to the global and regional change in the environment and climate as a whole, resulting in increased water toxicity and eutrophication, reduction in the number of stenobiont species of aquatic organisms against the background of an increase in the number of eurybiontic and invasive species. Modern communities of Lake Kuetsjarvi are the result of a combination of long-term changes in the abiotic environment and biotic interactions. Heavy-metal pollution of Lake Kuetsjarvi, observed since the 1930s, has led to the formation of a community that is resistant to this type of impact and supports large populations of adapted species. Adaptations of communities to the dynamics of the environmental conditions that their members are exposed to include changes in the species composition, quantitative indicators, ratios between individual taxonomic groups, and the population structure. The development of sympatric forms that differ in the ecological niches they occupy, morphology, and life cycle strategies, including the transition to a short-cycle survival strategy, allows whitefish to remain the dominant species and maintain high population numbers. Unlike the organismal level, responses to medium-term environmental changes on the population and community level are less specific and characterized by stronger inertia.

Parallelism in eco-morphology and gene expression despite variable evolutionary and genomic backgrounds in a Holarctic fish

Understanding the extent to which ecological divergence is repeatable is essential for predicting responses of biodiversity to environmental change. Here we test the predictability of evolution, from genotype to phenotype, by studying parallel evolution in a salmonid fish, Arctic charr (Salvelinus alpinus), across eleven replicate sympatric ecotype pairs (benthivorous-planktivorous and planktivorous-piscivorous) and two evolutionary lineages. We found considerable variability in eco-morphological divergence, with several traits related to foraging (eye diameter, pectoral fin length) being highly parallel even across lineages. This suggests repeated and predictable adaptation to environment. Consistent with ancestral genetic variation, hundreds of loci were associated with ecotype divergence within lineages of which eight were shared across lineages. This shared genetic variation was maintained despite variation in evolutionary histories, ranging from postglacial divergence in sympatry (ca. 10-15kya) to pre-glacial divergence (ca. 20-40kya) with postglacial secondary contact. Transcriptome-wide gene expression (44,102 genes) was highly parallel across replicates, involved biological processes characteristic of ecotype morphology and physiology, and revealed parallelism at the level of regulatory networks. This expression divergence was not only plastic but in part genetically controlled by parallel cis-eQTL. Lastly, we found that the magnitude of phenotypic divergence was largely correlated with the genetic differentiation and gene expression divergence. In contrast, the direction of phenotypic change was mostly determined by the interplay of adaptive genetic variation, gene expression, and ecosystem size. Ecosystem size further explained variation in putatively adaptive, ecotype-associated genomic patterns within and across lineages, highlighting the role of environmental variation and stochasticity in parallel evolution. Together, our findings demonstrate the parallel evolution of eco-morphology and gene expression within and across evolutionary lineages, which is controlled by the interplay of environmental stochasticity and evolutionary contingencies, largely overcoming variable evolutionary histories and genomic backgrounds.

Population genetics and sympatric divergence of the freshwater gudgeon, Gobiobotia filifer, in the Yangtze River inferred from mitochondrial DNA

The ecosystem and Pleistocene glaciations play important roles in population demography. The freshwater gudgeon, Gobiobotia filifer, is an endemic benthic fish in the Yangtze River and is a good model for ecological and evolutionary studies. This study aimed to decode the population structure of G. filifer in the Yangtze River and reveal whether divergence occurred before or after population radiation. A total of 292 specimens from eight locations in the upper and middle reaches of the Yangtze River were collected from 2014 to 2016 and analyzed via mitochondrial DNA Cyt b gene sequencing. A moderately high level of genetic diversity was found without structures among the population. However, phylogenetic and network topology showed two distinct haplotype groups, and each group contained a similar proportion of individuals from all sampled sites. This suggested the existence of two genetically divergent source populations in G. filifer. We deduced that a secondary contact of distinct glacial refugia was the main factor creating sympatric populations of G. filifer, and climate improvement promoted population expansion and colonization.

Beyond trophic morphology: stable isotopes reveal ubiquitous versatility in marine turtle trophic ecology

The idea that interspecific variation in trophic morphology among closely related species effectively permits resource partitioning has driven research on ecological radiation since Darwin first described variation in beak morphology among Geospiza. Marine turtles comprise an ecological radiation in which interspecific differences in trophic morphology have similarly been implicated as a pathway to ecopartition the marine realm, in both extant and extinct species. Because marine turtles are charismatic flagship species of conservation concern, their trophic ecology has been studied intensively using stable isotope analyses to gain insights into habitat use and diet, principally to inform conservation management. This legion of studies provides an unparalleled opportunity to examine ecological partitioning across numerous hierarchical levels that heretofore has not been applied to any other ecological radiation. Our contribution aims to provide a quantitative analysis of interspecific variation and a comprehensive review of intraspecific variation in trophic ecology across different hierarchical levels marshalling insights about realised trophic ecology derived from stable isotopes. We reviewed 113 stable isotope studies, mostly involving single species, and conducted a meta-analysis of data from adults to elucidate differences in trophic ecology among species. Our study reveals a more intricate hierarchy of ecopartitioning by marine turtles than previously recognised based on trophic morphology and dietary analyses. We found strong statistical support for interspecific partitioning, as well as a continuum of intraspecific trophic sub-specialisation in most species across several hierarchical levels. This ubiquity of trophic specialisation across many hierarchical levels exposes a far more complex view of trophic ecology and resource-axis exploitation than suggested by species diversity alone. Not only do species segregate along many widely understood axes such as body size, macrohabitat, and trophic morphology but the general pattern revealed by isotopic studies is one of microhabitat segregation and variation in foraging behaviour within species, within populations, and among individuals. These findings are highly relevant to conservation management because they imply ecological non-exchangeability, which introduces a new dimension beyond that of genetic stocks which drives current conservation planning. Perhaps the most remarkable finding from our data synthesis is that four of six marine turtle species forage across several trophic levels. This pattern is unlike that seen in other large marine predators, which forage at a single trophic level according to stable isotopes. This finding affirms suggestions that marine turtles are robust sentinels of ocean health and likely stabilise marine food webs. This insight has broader significance for studies of marine food webs and trophic ecology of large marine predators. Beyond insights concerning marine turtle ecology and conservation, our findings also have broader implications for the study of ecological radiations. Particularly, the unrecognised complexity of ecopartitioning beyond that predicted by trophic morphology suggests that this dominant approach in adaptive radiation research likely underestimates the degree of resource overlap and that interspecific disparities in trophic morphology may often over-predict the degree of realised ecopartitioning. Hence, our findings suggest that stable isotopes can profitably be applied to study other ecological radiations and may reveal trophic variation beyond that reflected by trophic morphology.

Adaptive Radiation from a Chromosomal Perspective: Evidence of Chromosome Set Stability in Cichlid Fishes (Cichlidae: Teleostei) from the Barombi Mbo Lake, Cameroon

Cichlid fishes are the subject of scientific interest because of their rapid adaptive radiation, resulting in extensive ecological and taxonomic diversity. In this study, we examined 11 morphologically distinct cichlid species endemic to Barombi Mbo, the largest crater lake in western Cameroon, namely Konia eisentrauti, Konia dikume, Myaka myaka, Pungu maclareni, Sarotherodon steinbachi, Sarotherodon lohbergeri, Sarotherodon linnellii, Sarotherodon caroli, Stomatepia mariae, Stomatepia pindu, and Stomatepia mongo. These species supposedly evolved via sympatric ecological speciation from a common ancestor, which colonized the lake no earlier than one million years ago. Here we present the first comparative cytogenetic analysis of cichlid species from Barombi Mbo Lake using both conventional (Giemsa staining, C-banding, and CMA₃/DAPI staining) and molecular (fluorescence in situ hybridization with telomeric, 5S, and 28S rDNA probes) methods. We observed stability on both macro and micro-chromosomal levels. The diploid chromosome number was 2n = 44, and the karyotype was invariably composed of three pairs of meta/submetacentric and 19 pairs of subtelo/acrocentric chromosomes in all analysed species, with the same numbers of rDNA clusters and distribution of heterochromatin. The results suggest the evolutionary stability of chromosomal set; therefore, the large-scale chromosomal rearrangements seem to be unlikely associated with the sympatric speciation in Barombi Mbo.

A way forward with eco evo devo: an extended theory of resource polymorphism with postglacial fishes as model systems

A major goal of evolutionary science is to understand how biological diversity is generated and altered. Despite considerable advances, we still have limited insight into how phenotypic variation arises and is sorted by natural selection. Here we argue that an integrated view, which merges ecology, evolution and developmental biology (eco evo devo) on an equal footing, is needed to understand the multifaceted role of the environment in simultaneously determining the development of the phenotype and the nature of the selective environment, and how organisms in turn affect the environment through eco evo and eco devo feedbacks. To illustrate the usefulness of an integrated eco evo devo perspective, we connect it with the theory of resource polymorphism (i.e. the phenotypic and genetic diversification that occurs in response to variation in available resources). In so doing, we highlight fishes from recently glaciated freshwater systems as exceptionally well-suited model systems for testing predictions of an eco evo devo framework in studies of diversification. Studies on these fishes show that intraspecific diversity can evolve rapidly, and that this process is jointly facilitated by (i) the availability of diverse environments promoting divergent natural selection; (ii) dynamic developmental processes sensitive to environmental and genetic signals; and (iii) eco evo and eco devo feedbacks influencing the selective and developmental environments of the phenotype. We highlight empirical examples and present a conceptual model for the generation of resource polymorphism - emphasizing eco evo devo, and identify current gaps in knowledge.

Genome-Wide Assessment of Diversity and Divergence Among Extant Galapagos Giant Tortoise Species

Genome-wide assessments allow for fuller characterization of genetic diversity, finer-scale population delineation, and better detection of demographically significant units to guide conservation compared with those based on "traditional" markers. Galapagos giant tortoises (Chelonoidis spp.) have long provided a case study for how evolutionary genetics may be applied to advance species conservation. Ongoing efforts to bolster tortoise populations, which have declined by 90%, have been informed by analyses of mitochondrial DNA sequence and microsatellite genotypic data, but could benefit from genome-wide markers. Taking this next step, we used double-digest restriction-site associated DNA sequencing to collect genotypic data at >26000 single nucleotide polymorphisms (SNPs) for 117 individuals representing all recognized extant Galapagos giant tortoise species. We then quantified genetic diversity, population structure, and compared results to estimates from mitochondrial DNA and microsatellite loci. Our analyses detected 12 genetic lineages concordant with the 11 named species as well as previously described structure within one species, C. becki. Furthermore, the SNPs provided increased resolution, detecting admixture in 4 individuals. SNP-based estimates of diversity and differentiation were significantly correlated with those derived from nuclear microsatellite loci and mitochondrial DNA sequences. The SNP toolkit presented here will serve as a resource for advancing efforts to understand tortoise evolution, species radiations, and aid conservation of the Galapagos tortoise species complex.

Evolution and diversity of two cisco forms in an outlet of glacial Lake Algonquin

The diversity of Laurentian Great Lakes ciscoes (Coregonus artedi, sensu lato) arose via repeated local adaptive divergence including deepwater ciscoes that are now extirpated or threatened. The nigripinnis form, or Blackfin Cisco, is extirpated from the Great Lakes and remains only in Lake Nipigon. Putative nigripinnis populations were recently discovered in sympatry with artedi in a historical drainage system of glacial Lake Algonquin, the precursor of lakes Michigan and Huron. Given the apparent convergence on Great Lakes form, we labeled this form blackfin. Here, we test the hypothesis that nigripinnis may have colonized this area from the Great Lakes as a distinct lineage. It would then represent a relict occurrence of the historical diversity of Great Lakes ciscoes. Alternatively, blackfin could have evolved in situ in several lakes. We captured more than 600 individuals in the benthic or pelagic habitat in 14 lakes in or near Algonquin Provincial Park (Ontario, Canada). Fish were compared based on habitat, morphology, and genetic variation at 6,676 SNPs. Contrary to our expectations, both cisco and blackfin belonged to an Atlantic lineage that colonized the area from the east, not from the Great Lakes. Sympatric cisco and blackfin were closely related while fish from different lakes were genetically differentiated, strongly suggesting the repeated in situ origin of each form. Across lakes, there was a continuum of ecological, morphological, and genetic differentiation that could be associated with alternative resources and lake characteristics. This study uncovers a new component of cisco diversity in inland lakes of Canada that evolved independently from ciscoes of the Laurentian Great lakes. The diversity of cisco revealed in this study and across their Canadian range presents a challenge for designating conservation units at the intraspecific level within the framework of the Committee on the Status of Endangered Wildlife in Canada (COSEWIC).

Resource polymorphism in European whitefish: Analysis of fatty acid profiles provides more detailed evidence than traditional methods alone

Resource polymorphism—whereby ancestral generalist populations give rise to several specialised morphs along a resource gradient—is common where species colonise newly formed ecosystems. This phenomenon is particularly well documented in freshwater fish populations inhabiting postglacial lakes formed at the end of the last ice age. However, knowledge on how such differential exploitation of resources across contrasting habitats might be reflected in the biochemical compositions of diverging populations is still limited, though such patterns might be expected. Here, we aimed to assess how fatty acids (FA)—an important biochemical component of animal tissues—diverged across a polymorphic complex of European whitefish (Coregonus lavaretus) and their closely related monomorphic specialist congener vendace (Coregonus albula) inhabiting a series of six subarctic lakes in northern Fennoscandia. We also explored patterns of FA composition in whitefish’s predators and invertebrate prey to assess how divergence in trophic ecology between whitefish morphs would relate to biochemical profiles of their key food web associates. Lastly, we assessed how information on trophic divergence provided by differential FA composition compared to evidence of resource polymorphism retrieved from more classical stomach content and stable isotopic (δ¹³C, δ¹⁵N) information. Examination of stomach contents provided high-resolution information on recently consumed prey, whereas stable isotopes indicated broad-scale patterns of benthic-pelagic resource use differentiation at different trophic levels. Linear discriminant analysis based on FA composition was substantially more successful in identifying whitefish morphs and their congener vendace as distinct groupings when compared to the other two methods. Three major FA (myristic acid, stearic acid, and eicosadienoic acid) proved particularly informative, both in delineating coregonid groups, and identifying patterns of pelagic-benthic feeding throughout the wider food web. Myristic acid (14:0) content and δ¹³C ratios in muscle tissue were positively correlated across fish taxa, and together provided the clearest segregation of fishes exploiting contrasting pelagic and benthic niches. In general, our findings highlight the potential of FA analysis for identifying resource polymorphism in animal populations where this phenomenon occurs, and suggest that this technique may provide greater resolution than more traditional methods typically used for this purpose.

Cisco diversity in a historical drainage of glacial Lake Algonquin

Cisco (Coregonus artedi (sensu lato) Lesueur, 1818) forms matching in appearance to Blackfin Cisco from the Laurentian Great Lakes occur in four lakes in Algonquin Park, Ontario, Canada, a historical drainage of glacial Lake Algonquin (precursor of lakes Michigan and Huron). Their occurrence may represent colonization from glacial Lake Algonquin drainage patterns 13 000 calibrated years BP or independent evolution within each lake. Gill-raker numbers, temperature at capture depth during lake stratification, and hurdle models of habitat distribution are summarized. Blackfin (nigripinnis-like) in the four lakes had higher gill-raker numbers than artedi-like cisco captured in nearby lakes or within the same lake. Two lakes have a bimodal gill-raker distribution that indicate co-occurrence of two forms. Blackfin occupied the hypolimnion with a peak depth distribution at 20–25 m. Maximum depth for blackfin was 35–40 m. The presence of the opossum shrimp (Mysis diluviana Audzijonyte and Väinölä, 2005) appears necessary for the occurrence of cisco diversity in lakes but not sufficient in all cases. The presence of two forms of cisco in at least two lakes points to the possibility of the colonization hypothesis or the ecological speciation hypothesis as accounting for this phenomenon. Genetic analysis is needed to determine which of these hypotheses best accounts for the occurrence of blackfin in Algonquin Park.

Model-based demographic inference of introgression history in European whitefish species pairs'

Parallel phenotypic differentiation is generally attributed to parallel adaptive divergence as an evolutionary response to similar environmental contrasts. Such parallelism may actually originate from several evolutionary scenarios ranging from repeated parallel divergence caused by divergent selection to a unique divergence event followed by gene flow. Reconstructing the evolutionary history underlying parallel phenotypic differentiation is thus fundamental to understand the relative contribution of demography and selection on genomic divergence during speciation. In this study, we investigate the divergence history of replicate European whitefish (Coregonus lavaretus), limnetic and benthic species pairs from two lakes in Norway and two lakes in Switzerland. Demographic models accounting for semi-permeability and linked selection were fitted to the unfolded joint allele frequency spectrum built from genome-wide SNPs and compared to each other in each species pair. We found strong support for a model of asymmetrical post-glacial secondary contact between glacial lineages in all four lakes. Moreover, our results suggest that heterogeneous genomic differentiation has been shaped by the joint action of linked selection accelerating lineage sorting during allopatry, and heterogeneous migration eroding divergence at different rates along the genome following secondary contact. Our analyses reveal how the interplay between demography, selection and historical contingency has influenced the levels of diversity observed in previous whitefish phylogeographic studies. This study thus provides new insights into the historical demographic and selective processes that shaped the divergence associated with ecological speciation in European whitefish.

The Effects of Demography and Genetics on the Neutral Distribution of Quantitative Traits

Neutral models for quantitative trait evolution are useful for identifying phenotypes under selection. These models often assume normally distributed phenotypes. This assumption may be violated when a trait is affected by relatively few variants or when the effects of those variants arise from skewed or heavy tailed distributions. Molecular phenotypes such as gene expression levels may have these properties. To accommodate deviations from normality, models making fewer assumptions about the underlying genetics and patterns of variation are needed. Here, we develop a general neutral model for quantitative trait variation using a coalescent approach. This model allows interpretation of trait distributions in terms of familiar population genetic parameters because it is based on the coalescent. We show how the normal distribution resulting from the infinitesimal limit, where the number of loci grows large as the effect size per mutation becomes small, depends only on expected pairwise coalescent times. We then demonstrate how deviations from normality depend on demography through the distribution of coalescence times as well as through genetic parameters. In particular, population growth events exacerbate deviations while bottlenecks reduce them. We demonstrate the practical applications of this model by showing how to sample from the neutral distribution of [Formula: see text], the ratio of the variance between subpopulations to that in the overall population. We further show it is likely impossible to distinguish sparsity from skewed or heavy tailed mutational effects using only sampled trait values. The model analyzed here greatly expands the parameter space for neutral trait models.

RADseq and mate choice assays reveal unidirectional gene flow among three lamprey ecotypes despite weak assortative mating: Insights into the formation and stability of multiple ecotypes in sympatry

Adaptive divergence with gene flow often results in complex patterns of variation within taxa exhibiting substantial ecological differences among populations. One example where this may have occurred is the parallel evolution of freshwater-resident nonparasitic lampreys from anadromous-parasitic ancestors. Previous studies have focused on transitions between these two phenotypic extremes, but here, we considered more complex evolutionary scenarios where an intermediate freshwater form that remains parasitic is found sympatrically with the other two ecotypes. Using population genomic analysis (restriction-associated DNA sequencing), we found that a freshwater-parasitic ecotype was highly distinct from an anadromous-parasitic form (Q_lake-P  = 96.8%, F_st  = 0.154), but that a freshwater-nonparasitic form was almost completely admixed in Loch Lomond, Scotland. Demographic reconstructions indicated that both freshwater populations likely derived from a common freshwater ancestor. However, while the nonparasitic ecotype has experienced high levels of introgression from the anadromous-parasitic ecotype (Q_anad-P  = 37.7%), there is no evidence of introgression into the freshwater-parasitic ecotype. Paradoxically, mate choice experiments predicted high potential for gene flow: Males from all ecotypes were stimulated to spawn with freshwater-parasitic females, which released gametes in response to all ecotypes. Differentially fixed single nucleotide polymorphisms identified genes associated with growth and development, which could possibly influence the timing of metamorphosis, resulting in significant ecological differences between forms. This suggests that multiple lamprey ecotypes can persist in sympatry following shifts in adaptive peaks, due to environmental change during their repeated colonization of post-glacial regions, followed by periods of extensive gene flow among such diverging populations.

Rapid buildup of sympatric species diversity in Alpine whitefish

Adaptive radiations in postglacial fish offer excellent settings to study the evolutionary mechanisms involved in the rapid buildup of sympatric species diversity from a single lineage. Here, we address this by exploring the genetic and ecological structure of the largest Alpine whitefish radiation known, that of Lakes Brienz and Thun, using microsatellite data of more than 2000 whitefish caught during extensive species-targeted and habitat-randomized fishing campaigns. We find six strongly genetically and ecologically differentiated species, four of which occur in both lakes, and one of which was previously unknown. These four exhibit clines of genetic differentiation that are paralleled in clines of eco-morphological and reproductive niche differentiation, consistent with models of sympatric ecological speciation along environmental gradients. In Lake Thun, we find two additional species, a profundal specialist and a species introduced in the 1930s from another Alpine whitefish radiation. Strong genetic differentiation between this introduced species and all native species of Lake Thun suggests that reproductive isolation can evolve among allopatric whitefish species within 15,000 years and persist in secondary sympatry. Consistent with speciation theory, we find stronger correlations between genetic and ecological differentiation for sympatrically than for allopatrically evolved species.

The Genetic Architecture Underlying the Evolution of a Rare Piscivorous Life History Form in Brown Trout after Secondary Contact and Strong Introgression

Identifying the genetic basis underlying phenotypic divergence and reproductive isolation is a longstanding problem in evolutionary biology. Genetic signals of adaptation and reproductive isolation are often confounded by a wide range of factors, such as variation in demographic history or genomic features. Brown trout (Salmo trutta) in the Loch Maree catchment, Scotland, exhibit reproductively isolated divergent life history morphs, including a rare piscivorous (ferox) life history form displaying larger body size, greater longevity and delayed maturation compared to sympatric benthivorous brown trout. Using a dataset of 16,066 SNPs, we analyzed the evolutionary history and genetic architecture underlying this divergence. We found that ferox trout and benthivorous brown trout most likely evolved after recent secondary contact of two distinct glacial lineages, and identified 33 genomic outlier windows across the genome, of which several have most likely formed through selection. We further identified twelve candidate genes and biological pathways related to growth, development and immune response potentially underpinning the observed phenotypic differences. The identification of clear genomic signals divergent between life history phenotypes and potentially linked to reproductive isolation, through size assortative mating, as well as the identification of the underlying demographic history, highlights the power of genomic studies of young species pairs for understanding the factors shaping genetic differentiation.

From top to bottom: Do Lake Trout diversify along a depth gradient in Great Bear Lake, NT, Canada?

Depth is usually considered the main driver of Lake Trout intraspecific diversity across lakes in North America. Given that Great Bear Lake is one of the largest and deepest freshwater systems in North America, we predicted that Lake Trout intraspecific diversity to be organized along a depth axis within this system. Thus, we investigated whether a deep-water morph of Lake Trout co-existed with four shallow-water morphs previously described in Great Bear Lake. Morphology, neutral genetic variation, isotopic niches, and life-history traits of Lake Trout across depths (0-150 m) were compared among morphs. Due to the propensity of Lake Trout with high levels of morphological diversity to occupy multiple habitat niches, a novel multivariate grouping method using a suite of composite variables was applied in addition to two other commonly used grouping methods to classify individuals. Depth alone did not explain Lake Trout diversity in Great Bear Lake; a distinct fifth deep-water morph was not found. Rather, Lake Trout diversity followed an ecological continuum, with some evidence for adaptation to local conditions in deep-water habitat. Overall, trout caught from deep-water showed low levels of genetic and phenotypic differentiation from shallow-water trout, and displayed higher lipid content (C:N ratio) and occupied a higher trophic level that suggested an potential increase of piscivory (including cannibalism) than the previously described four morphs. Why phenotypic divergence between shallow- and deep-water Lake Trout was low is unknown, especially when the potential for phenotypic variation should be high in deep and large Great Bear Lake. Given that variation in complexity of freshwater environments has dramatic consequences for divergence, variation in the complexity in Great Bear Lake (i.e., shallow being more complex than deep), may explain the observed dichotomy in the expression of intraspecific phenotypic diversity between shallow- vs. deep-water habitats. The ambiguity surrounding mechanisms driving divergence of Lake Trout in Great Bear Lake should be seen as reflective of the highly variable nature of ecological opportunity and divergent natural selection itself.

Diversifying selection drives parallel evolution of gill raker number and body size along the speciation continuum of European whitefish

Adaptive radiation is the evolution of ecological and phenotypical diversity. It arises via ecological opportunity that promotes the exploration of underutilized or novel niches mediating specialization and reproductive isolation. The assumed precondition for rapid local adaptation is diversifying natural selection, but random genetic drift could also be a major driver of this process. We used 27 populations of European whitefish (Coregonus lavaretus) from nine lakes distributed in three neighboring subarctic watercourses in northern Fennoscandia as a model to test the importance of random drift versus diversifying natural selection for parallel evolution of adaptive phenotypic traits. We contrasted variation for two key adaptive phenotypic traits correlated with resource utilization of polymorphic fish; the number of gill rakers and the total length of fish, with the posterior distribution of neutral genetic differentiation from 13 microsatellite loci, to test whether the observed phenotypic divergence could be achieved by random genetic drift alone. Our results show that both traits have been under diversifying selection and that the evolution of these morphs has been driven by isolation through habitat adaptations. We conclude that diversifying selection acting on gill raker number and body size has played a significant role in the ongoing adaptive radiation of European whitefish morphs in this region.

DNA Barcoding Marine Biodiversity: Steps from Mere Cataloguing to Giving Reasons for Biological Differences

DNA barcoding has become a useful tool in many contexts and has opened up a completely new avenue for taxonomy. DNA barcoding has its widest application in biodiversity and ecological research to detect and describe diversity whenever morphological discrimination is difficult or impossible (e.g., in the case of species lacking diagnostic characters, early life stages, or cryptic species). In this chapter, we outline the utility of including physiological parameters as part of species description in publicly available databases that catalog taxonomic information resulting from barcoding projects. Cryptic species or different life stages of a species often differ in their physiological traits. Thus, if physiological aspects were included in species definitions, the presently cryptic species could be distinguished. We furthermore give suggestions for physiological information that should be included in a species description and describe potential applications of DNA barcoding for research with physiological components.

Total mercury concentrations in liver and muscle of European whitefish (Coregonus lavaretus (L.)) in a subarctic lake - Assessing the factors driving year-round variation

Subarctic lakes are characterised by extreme seasonal variation in light and temperature which influences growth, maturation, condition and resource use of fishes. However, our understanding of how seasonal changes affect mercury concentrations of fishes is limited. We conducted a year-round study (3 ice-covered months, 3 open-water months) with open-water inter-annual aspect (3 years: samples from August/September), focusing on total mercury (THg) concentrations and ecological characteristics of a common freshwater fish, European whitefish (Coregonus lavaretus (L.)) from a subarctic lake. We measured THg concentrations from tissues with fast (liver, n = 164) and moderate (muscle, n = 225) turnover rates, providing information on THg dynamics over different temporal scales. In both tissues, lipid-corrected THg concentrations were highest in winter (liver: 1.70 ± 0.88 μg/g, muscle: 0.24 ± 0.05 μg/g) and lowest in summer (liver: 0.87 ± 0.72 μg/g, muscle: 0.19 ± 0.04 μg/g). THg concentrations increased in winter following the summer-autumn dietary shift to pelagic zooplankton and starvation after spawning. Whitefish THg concentrations decreased towards summer, and were associated with consumption of benthic macroinvertebrates and subsequent growth dilution. Mercury bioaccumulated in both tissues with age, both showing the strongest regression slopes in winter and lowest in summer. THg concentrations in liver and muscle tissue were correlated throughout the year, however the correlation was lowest in summer, indicating high metabolism during somatic growing season in summer and growth dilution. Multiple linear regression models explained 50% and 55% of the THg variation in liver and muscle both models dominated by seasonally-variable factors i.e. sexual maturity, δ¹³C, and condition factor. Seasonally varying bioaccumulation slopes and the higher level of intra-annual variation (21%) in whitefish THg concentration in muscle than the inter-annual accumulation (8%) highlight the importance of including seasonal factors in future THg studies.

Ecomorphological divergence drives differential mercury bioaccumulation in polymorphic European whitefish (Coregonus lavaretus) populations of subarctic lakes

Resource polymorphism, whereby ancestral trophic generalists undergo divergence into multiple specialist morphs, is common in salmonid fish populations inhabiting subarctic lakes. However, the extent to which such resource specialization into the three principal lake habitats (littoral, profundal, and pelagic) affects patterns of contaminant bioaccumulation remains largely unexplored. We assessed total mercury concentrations (THg) of European whitefish (Coregonus lavaretus (L.)) and their invertebrate prey in relation to potential explanatory variables across 6 subarctic lakes, of which three are inhabited by polymorphic (comprised of four morphs) and three by monomorphic populations. Among invertebrate prey, the highest THg concentrations were observed in profundal benthic macroinvertebrates, followed by pelagic zooplankton, with concentrations lowest in littoral benthic macroinvertebrates in both lake types. Broadly similar patterns were apparent in whitefish in polymorphic systems, where average age-corrected THg concentrations and bioaccumulation rates were the highest in pelagic morphs, intermediate in the profundal morph, and the lowest in the littoral morph. In monomorphic systems, age-corrected THg concentrations were generally lower, and showed pronounced lake-specific variation. In the polymorphic systems, we found significant relationships between whitefish muscle tissue THg concentration and gill raker count, resource use, lipid content and maximum length, whilst no such relationships were apparent in the monomorphic systems. Across all polymorphic lakes, the major variables explaining THg in whitefish were gill raker count and age, whereas in monomorphic systems, the factors were lake-specific. Whitefish resource polymorphism across the three main lake habitats therefore appears to have profound impacts on THg concentration and bioaccumulation rate. This highlights the importance of recognizing such intraspecific diversity in both future scientific studies and mercury monitoring programs.

Ecological speciation in a generalist consumer expands the trophic niche of a dominant predator

Ecological speciation – whereby an ancestral founder species diversifies to fill vacant niches – is a phenomenon characteristic of newly formed ecosystems. Despite such ubiquity, ecosystem-level effects of such divergence remain poorly understood. Here, we compared the trophic niche of European whitefish (Coregonus lavaretus) and their predators in a series of contrasting subarctic lakes where this species had either diversified into four ecomorphologically distinct morphs or instead formed monomorphic populations. We found that the trophic niche of whitefish was almost three times larger in the polymorphic than in the monomorphic lakes, due to an increase in intraspecific specialisation. This trophic niche expansion was mirrored in brown trout (Salmo trutta), a major predator of whitefish. This represents amongst the first evidence for ecological speciation directly altering the trophic niche of a predator. We suggest such mechanisms may be a common and important – though presently overlooked – factor regulating trophic interactions in diverse ecosystems globally.

Allometric trajectories of body and head morphology in three sympatric Arctic charr (Salvelinus alpinus (L.)) morphs

A study of body and head development in three sympatric reproductively isolated Arctic charr (Salvelinus alpinus (L.)) morphs from a subarctic lake (Skogsfjordvatn, northern Norway) revealed allometric trajectories that resulted in morphological differences. The three morphs were ecologically assigned to a littoral omnivore, a profundal benthivore and a profundal piscivore, and this was confirmed by genetic analyses (microsatellites). Principal component analysis was used to identify the variables responsible for most of the morphological variation of the body and head shape. The littoral omnivore and the profundal piscivore morph had convergent allometric trajectories for the most important head shape variables, developing bigger mouths and relatively smaller eyes with increasing head size. The two profundal morphs shared common trajectories for the variables explaining most of the body and head shape variation, namely head size relative to body size, placement of the dorsal and pelvic fins, eye size and mouth size. In contrast, the littoral omnivore and the profundal benthivore morphs were not on common allometric trajectories for any of the examined variables. The findings suggest that different selective pressures could have been working on traits related to their trophic niche such as habitat and diet utilization of the three morphs, with the two profundal morphs experiencing almost identical environmental conditions.

Morphological divergence of lake and stream Phoxinus of Northern Italy and the Danube basin based on geometric morphometric analysis

Minnows of the genus Phoxinus are promising candidates to investigate adaptive divergence, as they inhabit both still and running waters of a variety of altitudes and climatic zones in Europe. We used landmark-based geometric morphometric methods to quantify the level of morphological variability in Phoxinus populations from streams and lakes of Northern Italy and the Danube basin. We analyzed body shape differences of populations in the dorsal, lateral, and ventral planes, using a large array of landmarks and semilandmarks. As the species identification of Phoxinus on morphological characters is ambiguous, we used two mitochondrial genes to determine the genetic background of the samples and to ensure we are comparing homogenous groups. We have found significant body shape differences between habitats: Minnow populations inhabiting streams had a deeper body and caudal peduncle and more laterally inserted pectoral fins than minnows inhabiting lakes. We have also found significant body shape differences between genetic groups: Italian minnows had deeper bodies, deeper and shorter caudal peduncles, and a shorter and wider gape than both groups from the Danube. Our results show that the morphology of Phoxinus is highly influenced by habitat and that body shape variation between habitats was within the same range as between genetic groups. These morphological differences are possibly linked to different modes of swimming and foraging in the respective habitats and are likely results of phenotypic plasticity. However, differences in shape and interlandmark distances between the groups suggest that some (though few) morphometric characters might be useful for separating Phoxinus species.

Managing cryptic biodiversity: Fine-scale intralacustrine speciation along a benthic gradient in Alpine whitefish (Coregonus spp.)

Whitefish (Coregonus spp.) are an important catch for many freshwater fisheries, particularly in Switzerland. In support of this, supplemental stocking of whitefish species is carried out, despite lacking complete knowledge of the extent, distribution and origin of whitefish diversity in these lakes, potentially threatening local endemics via artificial gene flow. Here, we investigate phenotypic and genetic differentiation among coexisting whitefish species spawning along a depth gradient in a subalpine Swiss lake to better delineate intralacustrine whitefish biodiversity. We find depth‐related clines in adaptive morphology and in neutral genetic markers. This individual variation is structured in three distinct clusters with spatial overlap. Individual genetic distances correlate strongly with differences in growth rate and gill‐raker number, consistent with predictions of isolation‐by‐adaptation and ecological speciation. Genetic differentiation between species suggests reproductive isolation, despite demographic admixture on spawning grounds. Our results are consistent with clinal speciation resulting in three species coexisting in close ecological parapatry, one (C. sp. “benthic intermediate”) being previously unknown. A second unknown species spawning in close proximity was found to be of potential allochthonous origin. This study highlights the importance of taxonomically unbiased sampling strategies to both understand evolutionary mechanisms structuring biodiversity and to better inform conservation and fisheries management.