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

Phylogenetic analysis of museum specimens of houting Coregonus oxyrinchus shows the need for a revision of its extinct status

According to the IUCN Red List the anadromous houting Coregonus oxyrinchus is categorized as 'extinct'. However, this extinct status might be incorrect because taxonomic difference between C. oxyrinchus and the closely related C. lavaretus is based on a disputable morphological comparison. Also, phylogenetic studies on mtDNA only focused on recent obtained coregonids. We are the first to perform a mtDNA analysis on both historic and recent specimens, including the syntype specimen which was used for species description by Linnaeus originally. Two primer pairs for mitochondrial CytB and ND3 were used to extract sequences for phylogenetic analysis. Sequences from 14 out of 21 C. oxyrinchus museum specimens were successfully obtained and compared with sequences from recent obtained C. lavaretus. The sequences were combined with GenBank data from a previous phylogenetic study on houting to create a phylogenetic tree and two minimum spanning haplotype networks. Results show that C. oxyrinchus and C. lavaretus form a clade with limited genetic variation. Low bootstrap values also show weak support for geographical patterns in distribution of mitochondrial haplotypes. Statistical analysis of the haplotype networks also shows that historic and recent specimens are similar species. Our results suggest that C. oxyrinchus is a junior synonym of C. lavaretus. A definitive taxonomic revision could not be made because only CytB sequencing was successful for the syntype specimen. We discuss taxonomic consequences and the species-specific focus in nature conservation. We propose a shift in nature conservation to a more functional approach based on traits rather than species.

Parallel shifts in trout feeding morphology suggest rapid adaptation to alpine lake environments

Eco-evolutionary interactions following ecosystem change provide critical insight into the ability of organisms to adapt to shifting resource landscapes. Here we explore evidence for the rapid parallel evolution of trout feeding morphology following eco-evolutionary interactions with zooplankton in alpine lakes stocked at different points in time in the Wind River Range (Wyoming, USA). In this system, trout predation has altered the zooplankton species community and driven a decrease in average zooplankton size. In some lakes that were stocked decades ago, we find shifts in gill raker traits consistent with the hypothesis that trout have rapidly adapted to exploit available smaller-bodied zooplankton more effectively. We explore this morphological response in multiple lake populations across two species of trout (cutthroat trout, Oncorhynchus clarkii, and golden trout Oncorhynchus aguabonita) and examine the impact of resource availability on morphological variation in gill raker number among lakes. Furthermore, we present genetic data to provide evidence that historically stocked cutthroat trout populations likely derive from multiple population sources, and incorporate variation from genomic relatedness in our exploration of environmental predictors of feeding morphology. These findings describe rapid adaptation and eco-evolutionary interactions in trout and document an evolutionary response to novel, contemporary ecosystem change.

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.

Winter ecology of specialist and generalist morphs of European whitefish, Coregonus lavaretus, in subarctic northern Europe

European whitefish is a model species for adaptive radiation of fishes in temperate and subarctic lakes. In northern Europe the most commonly observed morphotypes are a generalist (LSR) morph and a pelagic specialist (DR) morph. The evolution of a pelagic specialist morph is something of an enigma, however, as this region is characterized by long, dark winters with pelagic primary production limited to a brief window in late summer. We conducted the first winter-based study of polymorphic whitefish populations to determine the winter ecology of both morphs, and we combined seasonal diet and stable isotope analysis with several proxies of condition in three polymorphic whitefish populations. The generalist LSR morph fed on benthic and pelagic prey in summer but was solely reliant on benthic prey in winter. This was associated with a noticeable but moderate reduction in condition, lipid content and stomach fullness in winter relative to summer. In contrast, the DR whitefish occupied a strict pelagic niche in both seasons. A significant reduction in pelagic prey during winter resulted in severe decrease in condition, lipid content and stomach fullness in DR whitefish in winter relative to summer, with the pelagic morph apparently approaching starvation in winter. We suggest that this divergent approach to seasonal foraging is associated with the divergent life-history traits of both morphs.

Genotyping-by-Sequencing Reveals the Impact of Restocking on Wild Common Carp Populations of the Southern Caspian Basin

Understanding the population structure and level of genetic diversity of wild populations is fundamental for appropriate stock management and species conservation. The common carp (Cyprinus carpio) is one of the most important bony fish throughout the Southern coastline of the Caspian Sea, but captures of this species have seen a dramatic reduction during the last decade. As a consequence, a restocking program has been put in place to maintain C. carpio populations, but its impact is not clear. In the present study, the population structure and genetic diversity of C. carpio in the Southern Caspian basin was determined using 17,828 single-nucleotide polymorphism (SNP) markers. A total of 117 individuals collected from four different locations in the southern Caspian basin and a farm were genotyped by genotyping-by-sequencing. The overall Fst obtained was 0.04, indicating a low level of differentiation between populations, and most genetic diversity was attributed to within population variation (97%). The low Fst values suggest that frequent migration events between different locations occur, and three migration events were inferred in the present study. However, each population still showed a distinct genetic profile, which allowed distinguishing the origin of the fish. This indicates that the ongoing restocking program is maintaining the differences between populations to some extent. Nonetheless, high inbreeding and low heterozygosity were detected in all populations, suggesting that additional conservation efforts are required to protect C. carpio populations in the Southern coast of the Caspian Sea.

Species interactions, environmental gradients and body size shape population niche width

Competition for shared resources is commonly assumed to restrict population-level niche width of coexisting species. However, the identity and abundance of coexisting species, the prevailing environmental conditions, and the individual body size may shape the effects of interspecific interactions on species' niche width. Here we study the effects of interspecific and intraspecific interactions, lake area and altitude, and fish body size on the trophic niche width and resource use of a generalist predator, the littoral-dwelling large, sparsely rakered morph of European whitefish (Coregonus lavaretus; hereafter LSR whitefish). We use stable isotope, diet and survey fishing data from 14 subarctic lakes along an environmental gradient in northern Norway. The isotopic niche width of LSR whitefish showed a humped-shaped relationship with increasing relative abundance of sympatric competitors, suggesting widest population niche at intermediate intensity of interspecific interactions. The isotopic niche width of LSR whitefish tended to decrease with increasing altitude, suggesting reduced niche in colder, less productive lakes. LSR whitefish typically shifted to a higher trophic position and increased reliance on littoral food resources with increasing body size, although between-lake differences in ontogenetic niche shifts were evident. In most lakes, LSR whitefish relied less on littoral food resources than coexisting fishes and the niche overlap between sympatric competitors was most evident among relatively large individuals (>250 mm). Individual niche variation was highest among >200 mm long LSR whitefish, which likely have escaped the predation window of sympatric predators. We demonstrate that intermediate intensity of interspecific interactions may broaden species' niche width, whereas strong competition for limited resources and high predation risk may suppress niche width in less productive environments. Acknowledging potential humped-shaped relationships between population niche width and interspecific interactions can help us understand species' responses to environmental disturbance (e.g. climate change and species invasions) as well as the driving forces of niche specialization.

The interaction of resource use and gene flow on the phenotypic divergence of benthic and pelagic morphs of Icelandic Arctic charr (Salvelinus alpinus)

Conceptual models of adaptive divergence and ecological speciation in sympatry predict differential resource use, phenotype-environment correlations, and reduced gene flow among diverging phenotypes. While these predictions have been assessed in past studies, connections among them have rarely been assessed collectively. We examined relationships among phenotypic, ecological, and genetic variation in Arctic charr (Salvelinus alpinus) from six Icelandic localities that have undergone varying degrees of divergence into sympatric benthic and pelagic morphs. We characterized morphological variation with geometric morphometrics, tested for differential resource use between morphs using stable isotopes, and inferred the amount of gene flow from single nucleotide polymorphisms. Analysis of stable isotopic signatures indicated that sympatric morphs showed similar difference in resource use across populations, likely arising from the common utilization of niche space within each population. Carbon isotopic signature was also a significant predictor of individual variation in body shape and size, suggesting that variation in benthic and pelagic resource use is associated with phenotypic variation. The estimated percentage of hybrids between sympatric morphs varied across populations (from 0% to 15.6%) but the majority of fish had genotypes (ancestry coefficients) characteristic of pure morphs. Despite evidence of reduced gene flow between sympatric morphs, we did not detect the expected negative relationship between divergence in resource use and gene flow. Three lakes showed the expected pattern, but morphs in the fourth showed no detectable hybridization and had relatively low differences in resource use between them. This coupled with the finding that resource use and genetic differentiation had differential effects on body shape variation across populations suggests that reproductive isolation maintains phenotypic divergence between benthic and pelagic morphs when the effects of resource use are relatively low. Our ability to assess relationships between phenotype, ecology, and genetics deepens our understanding of the processes underlying adaptive divergence in sympatry.

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.

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.

"And if you gaze long into an abyss, the abyss gazes also into thee": four morphs of Arctic charr adapting to a depth gradient in Lake Tinnsjøen

The origin of species is a central topic in biology. Ecological speciation might be a driver in adaptive radiation, providing a framework for understanding mechanisms, level, and rate of diversification. The Arctic charr Salvelinus alpinus L. is a polymorphic species with huge morphological and life‐history diversity in Holarctic water systems. We studied adaptive radiation of Arctic charr in the 460‐m‐deep Lake Tinnsjøen to (a) document eco‐morphology and life‐history traits of morphs, (b) estimate reproductive isolation of morphs, and (c) illuminate Holarctic phylogeography and lineages colonizing Lake Tinnsjøen. We compared Lake Tinnsjøen with four Norwegian outgroup populations. Four field‐assigned morphs were identified in Lake Tinnsjøen: the planktivore morph in all habitats except deep profundal, the dwarf morph in shallow‐moderate profundal, the piscivore morph mainly in shallow‐moderate profundal, and a new undescribed abyssal morph in the deep profundal. Morphs displayed extensive life‐history variation in age and size. A moderate‐to‐high concordance was observed among morphs and four genetic clusters from microsatellites. mtDNA suggested two minor endemic clades in Lake Tinnsjøen originating from one widespread colonizing clade in the Holarctic. All morphs were genetically differentiated at microsatellites (F_ST: 0.12–0.20), associated with different mtDNA clade frequencies. Analyses of outgroup lakes implied colonization from a river below Lake Tinnsjøen. Our findings suggest postglacial adaptive radiation of one colonizing mtDNA lineage with niche specialization along a depth–temperature–productivity–pressure gradient. Concordance between reproductive isolation and habitats of morphs implies ecological speciation as a mechanism. Particularly novel is the extensive morph diversification with depth into the often unexplored deepwater profundal habitat, suggesting we may have systematically underestimated biodiversity in lakes. In a biological conservation framework, it is imperative to protect endemic below‐species‐level biodiversity, particularly so since within‐species variation comprises an extremely important component of the generally low total biodiversity observed in the northern freshwater systems.

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.

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.

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).