Winter conditions are rapidly changing in temperate ecosystems, particularly for those that experience periods of snow and ice cover. Relatively little is known of winter ecology in these systems, due to a historical research focus on summer ‘growing seasons’. We executed the first global quantitative synthesis on under‐ice lake ecology, including 36 abiotic and biotic variables from 42 research groups and 101 lakes, examining seasonal differences and connections as well as how seasonal differences vary with geophysical factors. Plankton were more abundant under ice than expected; mean winter values were 43.2% of summer values for chlorophyll a, 15.8% of summer phytoplankton biovolume and 25.3% of summer zooplankton density. Dissolved nitrogen concentrations were typically higher during winter, and these differences were exaggerated in smaller lakes. Lake size also influenced winter‐summer patterns for dissolved organic carbon (DOC), with higher winter DOC in smaller lakes. At coarse levels of taxonomic aggregation, phytoplankton and zooplankton community composition showed few systematic differences between seasons, although literature suggests that seasonal differences are frequently lake‐specific, species‐specific, or occur at the level of functional group. Within the subset of lakes that had longer time series, winter influenced the subsequent summer for some nutrient variables and zooplankton biomass.

Stable isotope analysis provides a powerful tool to identify the energy sources which fuel consumers, to understand trophic interactions and to infer consumer trophic position (TP), an important concept that describes the ecological role of consumers in food webs. However, current methods for estimating TP using stable isotopes are limited and do not fulfil the complete potential of the isotopic approach. For instance, researchers typically use point estimates for key parameters including trophic discrimination factors and isotopic baselines, and do not explicitly include variance associated with these parameters when calculating TP. We present “tRophicPosition,” anrpackage incorporating a Bayesian model for the calculation of consumer TP at the population level using stable isotopes, with one or two baselines. It combines Markov Chain Monte Carlo simulations through JAGS and statistical and graphical analyses using R. We model consumer and baseline observations using relevant statistical distributions, allowing them to be treated as random variables. The calculation of TP—a random parameter—for one baseline follows standard equations linking¹⁵N enrichment per trophic level and the trophic position of the baseline (e.g. a primary producer or primary consumer). In the case of two baselines, a simple mixing model incorporating δ¹³C allows for the differentiation between two distinct sources of nitrogen, thus including heterogeneity derived from alternatives sources of δ¹⁵N. Methods currently implemented in “tRophicPosition” include loading, plotting and summarizing stable isotope data either from multiple sites and/or communities or a local assemblage; loading trophic discrimination factors from an internal database or generating them; defining and initializing a Bayesian model of TP; sampling posterior parameters; analysing, comparing and plotting posterior estimates of TP and other parameters; and calculating a parametric (non‐Bayesian) TP estimate. Additionally, full documentation including examples, multiple vignettes and code are available for download.

Aim: Understanding the variation in community composition and species abundances (i.e., β‐diversity) is at the heart of community ecology. A common approach to examine β‐diversity is to evaluate directional variation in community composition by measuring the decay in the similarity among pairs of communities along spatial or environmental distance. We provide the first global synthesis of taxonomic and functional distance decay along spatial and environmental distance by analysing 148 datasets comprising different types of organisms and environments.Location: Global.Time period: 1990 to present.Major taxa studied: From diatoms to mammals.Method: We measured the strength of the decay using ranked Mantel tests (Mantel r) and the rate of distance decay as the slope of an exponential fit using generalized linear models. We used null models to test whether functional similarity decays faster or slower than expected given the taxonomic decay along the spatial and environmental distance. We also unveiled the factors driving the rate of decay across the datasets, including latitude, spatial extent, realm and organismal features.Results: Taxonomic distance decay was stronger than functional distance decay along both spatial and environmental distance. Functional distance decay was random given the taxonomic distance decay. The rate of taxonomic and functional spatial distance decay was fastest in the datasets from mid‐latitudes. Overall, datasets covering larger spatial extents showed a lower rate of decay along spatial distance but a higher rate of decay along environmental distance. Marine ecosystems had the slowest rate of decay along environmental distances.Main conclusions: In general, taxonomic distance decay is a useful tool for biogeographical research because it reflects dispersal‐related factors in addition to species responses to climatic and environmental variables. Moreover, functional distance decay might be a cost‐effective option for investigating community changes in heterogeneous environments.

Fish are an important source of eicosapentaenoic acid (EPA) and docosahexaenoic acid (DHA) for birds, mammals and humans. In aquatic food webs, these highly unsaturated fatty acids (HUFA) are essential for many physiological processes and mainly synthetized by distinct phytoplankton taxa. Consumers at different trophic levels obtain essential fatty acids from their diet because they cannot produce these sufficiently de novo. Here, we evaluated how the increase in phosphorus concentration (eutrophication) or terrestrial organic matter inputs (brownification) change EPA and DHA content in the phytoplankton. Then, we evaluated whether these changes can be seen in the EPA and DHA content of piscivorous European perch (Perca fluviatilis), which is a widely distributed species and commonly consumed by humans. Data from 713 lakes showed statistically significant differences in the abundance of EPA- and DHA-synthesizing phytoplankton as well as in the concentrations and content of these essential fatty acids among oligo-mesotrophic, eutrophic and dystrophic lakes. The EPA and DHA content of phytoplankton biomass (mgHUFAg-1) was significantly lower in the eutrophic lakes than in the oligo-mesotrophic or dystrophic lakes. We found a strong significant correlation between the DHA content in the muscle of piscivorous perch and phytoplankton DHA content (r=0.85) as well with the contribution of DHA-synthesizing phytoplankton taxa (r=0.83). Among all DHA-synthesizing phytoplankton this correlation was the strongest with the dinoflagellates (r=0.74) and chrysophytes (r=0.70). Accordingly, the EPA+DHA content of perch muscle decreased with increasing total phosphorus (r2=0.80) and dissolved organic carbon concentration (r2=0.83) in the lakes. Our results suggest that although eutrophication generally increase biomass production across different trophic levels, the high proportion of low-quality primary producers reduce EPA and DHA content in the food web up to predatory fish. Ultimately, it seems that lake eutrophication and brownification decrease the nutritional quality of fish for human consumers.

Summary: 1. The adaptive radiation of fishes into benthic (littoral) and pelagic (lentic) morphs in post‐glacial lakes has become an important model system for speciation. Although these systems are well studied, there is little evidence of the existence of morphs that have diverged to utilize resources in the remaining principal lake habitat, the profundal zone. 2. Here, we tested phenotype‐environment correlations of three whitefish (Coregonus lavaretus) morphs that have radiated into littoral, pelagic and profundal niches in northern Scandinavian lakes. We hypothesized that morphs in such trimorphic systems would have a morphology adapted to one of the principal lake habitats (littoral, pelagic or profundal niches). Most whitefish populations in the study area are formed by a single (monomorphic) whitefish morph, and we further hypothesized that these populations should display intermediate morphotypes and niche utilization. We used a combination of traditional (stomach content, habitat use, gill raker counts) and more recently developed (stable isotopes, geometric morphometrics) techniques to evaluate phenotype‐environment correlations in two lakes with trimorphic and two lakes with monomorphic whitefish. 3. Distinct phenotype‐environment correlations were evident for each principal niche in whitefish morphs inhabiting trimorphic lakes. Monomorphic whitefish exploited multiple habitats, had intermediate morphology, displayed increased variance in gillraker‐counts, and relied significantly on zooplankton, most likely due to relaxed resource competition. 4. We suggest that the ecological processes acting in the trimorphic lakes are similar to each other, and are driving the adaptive evolution of whitefish morphs, possibly leading to the formation of new species.

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