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

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.

Geographic destiny trumps taxonomy in the Roundtail Chub, Gila robusta species complex (Teleostei, Leuciscidae)

The Gila robusta species complex in the lower reaches of the Colorado River includes three nominal and contested species (G. robusta, G. intermedia, and G. nigra) originally defined by morphological and meristic characters. In subsequent investigations, none of these characters proved diagnostic, and species assignments were based on capture location. Two recent studies applied conservation genomics to assess species boundaries and reached contrasting conclusions: an ezRAD phylogenetic study resolved 5 lineages with poor alignment to species categories and proposed a single species with multiple population partitions. In contrast, a dd-RAD coalescent study concluded that the three nominal species are well-supported evolutionarily lineages. Here we developed a draft genome (~ 1.229 Gbp) to apply genome-wide coverage (10,246 SNPs) with nearly range-wide sampling of specimens (G. robusta N = 266, G. intermedia N = 241, and G. nigra N = 117) to resolve this debate. All three nominal species were polyphyletic, whereas 5 of 8 watersheds were monophyletic. AMOVA partitioned 23.1% of genetic variance among nominal species, 30.9% among watersheds, and the Little Colorado River was highly distinct (FST ranged from 0.79 to 0.88 across analyses). Likewise, DAPC identified watersheds as more distinct than species, with the Little Colorado River having 297 fixed nucleotide differences compared to zero fixed differences among the three nominal species. In every analysis, geography explains more of the observed variance than putative taxonomy, and there are no diagnostic molecular or morphological characters to justify species designation. Our analysis reconciles previous work by showing that species identities based on type location are supported by significant divergence, but natural geographic partitions show consistently greater divergence. Thus, our data confirm Gila robusta as a single polytypic species with roughly a dozen highly isolated geographic populations, providing a strong scientific basis for watershed-based future conservation.

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.

Competition among small individuals hinders adaptive radiation despite ecological opportunity

Ontogenetic diet shifts, where individuals change their resource use during development, are the rule rather than the exception in the animal world. Here, we aim to understand how such changes in diet during development affect the conditions for an adaptive radiation in the presence of ecological opportunity. We use a size-structured consumer–resource model and the adaptive dynamics approach to study the ecological conditions for speciation. We assume that small individuals all feed on a shared resource. Large individuals, on the other hand, have access to multiple food sources on which they can specialize. We find that competition among small individuals can hinder an adaptive radiation to unfold, despite plenty of ecological opportunity for large individuals. When small individuals experience strong competition for food, they grow slowly and only a few individuals are recruited to the larger size classes. Hence, competition for food among large individuals is weak and there is therefore no disruptive selection. In addition, initial conditions determine if an adaptive radiation occurs or not. A consumer population initially dominated by small individuals will not radiate. On the other hand, a population initially dominated by large individuals may undergo adaptive radiation and diversify into multiple species.

Intraspecific structure of the Coregonus lavaretus complex in water bodies of Siberia: a case of postglacial allopatric origin of Yukagirian whitefish

The results of morphological and genetic analyses of forms/species of the Coregonus lavaretus pidschian (Gmelin, 1789) complex from the Indigirka and Kolyma river basins are presented in the context of there being recent postglacial speciation events. It has been found that the studied whitefishes belong to the sparsely rakered and low lateral-line forms and have previously been described as Coregonus lavaretus pidschian n. jucagiricus Drjagin (Berg), 1932. Based on these characters, this whitefish does not differ from most Arctic whitefish populations (in particular from Coregonus lavaretus glacialis Kirillov, 1972). Analysis of variability of the ND1 gene (mtDNA) showed that whitefishes from the Indigirka and Kolyma basins belong to a distant phylogenetic lineage, which is significantly different from all previously studied whitefish lineages from the Ob, Yenisei, Lena, Anadyr, and Amur river basins. Analysis of variability of the ITS1 fragment (nDNA) showed that all studied forms/species (from Ob River to Amur River basins), including C. l. pidschian n. jucagiricus, have a tandem arrangement of two identical nucleotide fragments and very similar nucleotide composition of the ITS1 region. Based on contemporary data, this phylogenetic lineage of the C. pidschian complex could be considered a young postglacial allopatric species.

Varying Intensities of Introgression Obscure Incipient Venom-Associated Speciation in the Timber Rattlesnake (Crotalus horridus)

Ecologically divergent selection can lead to the evolution of reproductive isolation through the process of ecological speciation, but the balance of responsible evolutionary forces is often obscured by an inadequate assessment of demographic history and the genetics of traits under selection. Snake venoms have emerged as a system for studying the genetic basis of adaptation because of their genetic tractability and contributions to fitness, and speciation in venomous snakes can be associated with ecological diversification such as dietary shifts and corresponding venom changes. Here, we explored the neurotoxic (type A)–hemotoxic (type B) venom dichotomy and the potential for ecological speciation among Timber Rattlesnake (Crotalus horridus) populations. Previous work identified the genetic basis of this phenotypic difference, enabling us to characterize the roles geography, history, ecology, selection, and chance play in determining when and why new species emerge or are absorbed. We identified significant genetic, proteomic, morphological, and ecological/environmental differences at smaller spatial scales, suggestive of incipient ecological speciation between type A and type B C. horridus. Range-wide analyses, however, rejected the reciprocal monophyly of venom type, indicative of varying intensities of introgression and a lack of reproductive isolation across the range. Given that we have now established the phenotypic distributions and ecological niche models of type A and B populations, genome-wide data are needed and capable of determining whether type A and type B C. horridus represent distinct, reproductively isolated lineages due to incipient ecological speciation or differentiated populations within a single species.

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.

Ambient temperature as a factor contributing to the developmental divergence in sympatric salmonids

Factors and mechanisms promoting resource-based radiation in animals still represent a main challenge to evolutionary biology. The modifications of phenotype tied with adaptive diversification may result from an environmentally related shift having occurred at the early stage of development. Here, we study the role of temperature dynamics on the reproductive sites in the early-life divergence and adaptive radiation of the salmonid fish Salvelinus malma dwelling in the Lake Kronotskoe basin (North-East Asia). Local sympatric charr ecomorphs demonstrate strict homing behaviour guiding the preordained distribution along tributaries and, hence, further development under different temperatures. We thoroughly assessed the annual temperature dynamics at the spawning grounds of each morph as compared to an ancestral anadromous morph. Then we carried out an experimental rearing of both under naturally diverging and uniformed temperatures. To compare the morphs' development under the dynamically changing temperatures, we have designed a method based on calculating the accumulated heat by the Arrhenius equation. The proposed equation shows a strong predictive power and, at the same time, is not bias-susceptible when the developmental temperature approximates 0°C. The temperature was found to significantly affect the charrs' early ontogeny, which underlies the divergence of developmental and growth rates between the morphs, as well as morph-specific ontogenetic adaptations to the spawning site's temperatures. As opposed to the endemic morphs from Lake Kronotskoe, the anadromous S. malma, being unexposed to selection оn highly specific reproduction conditions, showed a wide temperature tolerance, Our findings demonstrate that the hatch, onset timing of external feeding, and size dissimilarities between the sympatric morphs reveal themselves during the development under contrast temperatures. As a result of the observed developmental disparities, the morphs occupy specific definitive foraging niches in the lake.

Genetic relationships between sympatric and allopatric Coregonus ciscoes in North and Central Europe

BACKGROUND

Sympatric speciation along ecological gradients has been studied repeatedly, in particular in freshwater fishes. Rapid post-glacial ecological divergence has resulted in numerous endemic species or ecologically distinct populations in lakes of the temperate zones. Here, we focus on the Baltic cisco (Coregonus albula) complex, to study the genetic similarity among two pairs of sympatric autumn- and spring-spawning populations from post-glacial German Lakes Stechlin and Breiter Luzin. For comparison, we included a similar pair of sympatric populations from the Swedish Lake Fegen. We wanted to explore potential genetic similarities between the three sympatric cisco population pairs in the three lakes, to evaluate whether the pairs may have emerged independently in the three lakes, or whether two different species may have colonized all three lakes independently. Furthermore, we considered allopatric C. albula populations from three Polish, three Finnish, and four Swedish locations, added one Siberian population of the sister species C. sardinella and a Swedish C. maraena (whitefish) population. By genotyping nine microsatellite markers in 655 individuals from these 18 populations, we wanted to elucidate how strongly the cisco populations differ across a larger geographical area within Europe. Finally, we compared the genetic differences between the spring- and autumn-spawning populations of ciscoes in the two German lakes to infer the potentially deteriorating effect of strong anthropogenic pressure on the lakes.

RESULTS

Dendrogram, Principal Coordinate Analysis and admixture analysis all indicated strong correspondence between population differentiation and geographical location for most cisco populations in Europe, including the Siberian population of C. sardinella. However, populations from some Swedish lakes deviated from this general pattern, by showing a distinct genetic structure. We found evidence for independent evolution of the three sympatric population pairs, because the populations co-occurring in the same lake were always most closely related. However, genetic differentiation was weak in the two German population pairs, but strong in the Swedish Lake Fegen, indicating that the weak differentiation in the German pairs reported earlier has eroded further.

CONCLUSIONS

Our results suggest that the genetic differentiation at neutral genetic markers among populations of the Baltic cisco complex has evolved (and is maintained) by random genetic drift in isolated populations. However, earlier studies on the Swedish populations combining mitochondrial DNA and microsatellite data indicate that also post-glacial immigration from separate glacial refugia has shaped the present genetic population structure. The low neutral differentiation of the German sympatric pairs in contrast to the Swedish pair suggests that recent anthropogenic effects on the lakes in Germany may put the endemic spring-spawners at risk to extinction.

Multidimensional divergent selection, local adaptation, and speciation

Divergent selection applied to one or more traits drives local adaptation and may lead to ecological speciation. Divergent selection on many traits might be termed "multidimensional" divergent selection. There is a commonly held view that multidimensional divergent selection is likely to promote local adaptation and speciation to a greater extent than unidimensional divergent selection. We disentangle the core concepts underlying dimensionality as a property of the environment, phenotypes, and genome. In particular, we identify a need to separate the overall strength of selection and the number of loci affected from dimensionality per se, and to distinguish divergence dimensionality from dimensionality of stabilizing selection. We then critically scrutinize this commonly held view that multidimensional selection promotes speciation, re-examining the evidence base from theory, experiments, and nature. We conclude that the evidence base is currently weak and generally suffers from confounding of possible causal effects. Finally, we propose several mechanisms by which multidimensional divergent selection and related processes might influence divergence, both as a driver and as a barrier.

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.

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.