Abstract
Northern ecosystems are experiencing rapid and large-scale changes driven by accelerated warming, which have profound effects on the terrestrial and freshwater biodiversity. A comprehensive understanding of the distribution of aquatic biodiversity of subarctic ecosystems is therefore needed to better predict future trajectories of their unique biodiversity. In this study, we examined the functional diversity of chironomid communities in subarctic lakes across a 1000 m-elevation gradient, reflecting gradual changes in temperature and landscape characteristics. Using fuzzy correspondence analyses, we investigated spatial variability in trait composition of chironomid communities from 100 lakes in northern Sweden, and tested the hypotheses that (1) climate directly and indirectly shapes chironomid trait composition across the studied gradient, and (2) that generalist taxa with smaller body size and broader food preferences are more able to persist in cold environments. Our results showed that complex interplays between direct (e.g. temperature) and indirect climate processes (e.g. elevation-driven changes in vegetation/habitats) affect the functional diversity of chironomid communities. Specifically, traits such as larval size, food preference and feeding habits were well separated along the gradient, and this pattern revealed that low elevation lakes with forested catchments tended to have more sediment-feeding taxa and larger larvae than those above the tree line. As expected, food resource availability in lakes is strongly linked to vegetation composition/cover, and traits related to resource exploitation in chironomid communities are therefore well constrained by landscape characteristics. Furthermore, our findings suggested that short life cycles could facilitate the development of viable population in northern and high-elevation lakes where the short ice-free period is a limiting factor, thus contradicting patterns showing smaller organisms in warmer environments reported for other invertebrates. As a consequence of climate warming, the highest elevation lakes in subarctic landscapes will likely lose their typical cold-adapted chironomid taxa along with their functional attributes leading to potential impacts on the food web structure and the overall functioning of northern lake ecosystems.
Similar content being viewed by others
Introduction
Patterns in biodiversity (i.e. taxonomic composition and functional attributes) are determined by a combination of large-scale landscape attributes and habitat properties (Townsend et al. 2003), among which climate plays a central role. Northern ecosystems are experiencing rapid and large-scale changes driven by accelerated warming (AMAP 2011), and climate change is already responsible for ecological shifts observed on freshwater biodiversity (Heino et al. 2009; Lento et al. 2019; Lau et al. 2020). Also, the recent Arctic Biodiversity Assessment by the Conservation of Arctic Flora and Fauna (CAFF) has identified climate change as the most significant threat to Arctic biodiversity and pointed out subarctic freshwaters as hotspots for ecological change (Lento et al. 2019). A comprehensive understanding of the effects of climate warming on structural and functional diversity of subarctic ecosystems is therefore needed to better predict future trajectories of their unique biodiversity.
Functional traits reflect how species interact with their environment and with co-occurring species (Dı́az and Cabido 2001) and can be used to estimate the functional diversity of communities. The effects of environmental changes on ecosystem functions are partly set by the functional diversity in a community (Hulot et al. 2000; Reynolds et al. 2002; Hooper et al. 2005), and trait composition of communities therefore helps to define ecological resilience (McLean et al. 2019). Although functional diversity responses to changing environmental conditions have been well documented for temperate lakes (see also Schmera et al. 2017), we are still lacking information about their effects on ecosystem properties of subarctic lakes.
Reponses of species distribution to environmental changes are mainly controlled by the ability of species to disperse and to persist in suitable environments (Lavergne et al. 2010). For example, short generation times (represented by traits such as multivoltinism and small body size) may facilitate dispersal, while the ability to exploit available food resources (e.g. feeding preference, food type) influences the establishment and persistence of viable populations (Urban et al. 2012). Accordingly, generalist species with smaller body size and broader food preferences should be more likely to find suitable habitats in extreme environments (Angert et al. 2011). In freshwater ecosystems, invertebrate trait composition responses to climate change are still debated, and most studies have considered only a few traits (see also Gardner et al. 2011). According to Bergmann’s rule, smaller invertebrates should be found in warmer environments (see also Daufresne et al. 2009; Hayden et al. 2017), but this rule might not be valid for all taxonomic groups of invertebrates and/or ecosystem types (see also MacLean and Beissinger 2017). Furthermore, traits related to resource exploitation have been largely neglected in this framework but should be considered as they could supplement our knowledge of community trait responses to climate change.
Chironomidae (Arthropoda; Diptera; Nematocera) are non-biting midges whose larvae constitute a predominant group in northern freshwater ecosystems (Northington et al. 2010; Lento et al. 2019), thus providing an ideal model for studying biodiversity patterns in subarctic environments. Chironomid species are well studied and have been good indicators of environmental change, as they have well defined ecological requirements across environmental gradients, such as temperature, nutrient and oxygen concentrations (Brodersen and Quinlan 2006; Verbruggen et al. 2011; Heiri et al. 2011). Furthermore, chironomid species exhibit distinct functional traits related to their biological and physiological properties (Usseglio-Polatera et al. 2000; Serra et al. 2016). As a single taxon may have several traits, minor changes in taxonomic composition of chironomid assemblage can shift in assemblage trait composition and ecosystem functioning (c.f. van Kleef et al. 2015; Serra et al. 2017; Desrosiers et al. 2019). Hence, novel functional (rather than taxonomic) approaches must be developed to explore functional diversity changes of chironomid communities across environmental gradients (Motta and Massaferro 2019).
Due to the lack of long-term monitoring in most ecosystems, a comprehensive understanding of the effects of climate change on the functional diversity of chironomid communities of northern ecosystems is lacking, but space-for-time substitutions could help to overcome this problem. Relationships between the functional attributes of chironomid communities and environmental conditions can be studied using calibration training data sets that describe the spatial distribution and abundance of taxa across a temperature gradient (Eggermont and Heiri 2012). Training sets are produced by collecting sediment samples from a large number of lakes across latitudinal or elevation gradients and by identifying the chironomid exoskeleton remains (i.e. head capsules) in these sediments, and provide insight into the composition of the contemporary chironomid communities of these lakes (Heiri et al. 2011). For this purpose, surface sediments are usually retrieved from the largest depth of lakes where chironomid remains from the entire lake accumulate (Heiri 2004; van Hardenbroek et al. 2010). Training sets could thus give valuable insights into biodiversity patterns in subarctic environments and their relationships with environmental conditions.
This study used samples collected across an elevation gradient (169–1183 m.a.s.l.) of subarctic lakes in Sweden (Larocque et al. 2001) to evaluate the effects of temperature and catchment characteristics on chironomid trait diversity. We tested the hypotheses (1) that climate directly and indirectly shapes trait composition of chironomid communities across the studied 1000 m-elevation gradient, and consequently (2) that generalist taxa with smaller body size and large ranges in food preferences are more able to persist in cold environments than specialist taxa.
Materials and methods
Study sites and chironomid training set
We used the previously published training set (Larocque et al. 2001) describing chironomid assemblage composition of 100 lakes in the northern boreal and Arctic/alpine ecoregions of Sweden (see Table S1 for limnological and geographic characteristics). Due to influences of elevation (ranging 169–1183 m a.s.l.), the area is characterised by a strong gradient in temperature (estimated mean July air temperatures ranging 7.0–14.7 °C, Larocque et al. 2001), and the ice-free season ranges from ca. 6 months at low elevation lakes to < 2 months at the highest elevations (Table S1). Moreover, catchment vegetation ranges from conifer-dominated (Pinus sylvestris, Picea abies) boreal forests at low elevations, through mountain-birch woodland (Betula pubescens ssp. tortuosa) and alpine meadows, to areas with sparse alpine-tundra vegetation (lichens and mosses) or barren ground at the highest elevations. The lakes are generally small (< 20 ha), shallow (average maximum depth of 6 m) and oligotrophic (total phosphorus concentrations ranging 8.3 ± 7 µg L−1).
Surface sediment samples (0–1 cm) were retrieved from the deepest part of each lake using a modified Kajak-gravity corer. In the lab, chironomid head capsules were hand-sorted from samples of wet sediment (3–5 mL) following recommendations from Walker (2001) and identified under a microscope using Wiederholm (1983). Chironomid counts were then pooled to the genus level to fit the taxonomic resolution of the European Chironomidae Database (Serra et al. 2016). Data are expressed as relative abundance, and only taxa occurring in at least two lakes were used.
Functional trait database
The functional diversity of chironomid communities was characterized using four life-history traits (i.e. 4th instar larval size, emergence season, voltinism and flight period), as well as traits describing food preference, and feeding habits obtained from the European Chironomidae Database published by Serra et al. (2016; Table 1). These six traits are considered responsive to climate change (Bonada et al. 2007) and represent biological features that directly influence ecosystem functions (Hulot et al. 2000; Hevia et al. 2017). Trait categories were standardized to sum 1 for each trait (depicted as the standardized trait profile following a fuzzy coding procedure, Chevenet et al. 1994), ensuring that all taxa were weighted equally in further analyses. A trait-by-lake dataset was computed as the cross-product between the relative abundance of chironomid genera and the standardized trait composition, thus containing the relative abundance of each trait category in each lake.
Data analysis
An overview of the data analysis workflow adopted in this study is provided in Fig. S1. A principal component analysis (PCA) was performed on lake characteristics to provide an overview of environmental gradients in the training set (step 1 in Fig. S1). Functional diversity of chironomid communities was analysed following a procedure similar to that by Desrosiers et al. (2019). In this procedure, correlation between functional traits (Table 1) and chironomid taxa found in the training set was analysed using a fuzzy correspondence analysis (FCA) applied to the taxa-by-trait category dataset (18 taxa × 33 categories of 6 traits; (step 2.1 in Fig. S1). The contribution of each trait to total variation was quantified by correlation ratios representing the amount of variance explained by each trait along the first two FCA axes (Picazo et al. 2012). Cluster analysis, based on a matrix of Euclidean distances among the taxa calculated on taxon coordinates along the first four FCA axes, was applied to define functional groups of taxa with similar trait combinations (step 2.2 in Fig. S1). Chironomid genera were then grouped into successively larger groups using the Ward’s method (Ward 1963) and the significance of the functional groups was assessed using the broken-stick model (Bennett 1996).
An FCA was also performed to compare lakes along the elevation gradient based on the trait composition of their chironomid communities (i.e. trait-by-lake dataset: 100 lakes × 33 categories of 6 traits; (step 3.1 in Fig. S1). Cluster analysis, based on a matrix of Euclidean distances calculated on lake coordinates along the first four FCA axes, was applied to identify lake clusters with similar trait composition, and lakes were then grouped into successively larger groups using the Ward’s method and the broken-stick model for selecting significant clusters (step 3.2 in Fig. S1). In addition, the “IndVal” method (Dufrêne and Legendre 1997) was applied to identify indicator categories of traits associated with each lake cluster. Only significant categories (p < 0.05) were selected as indicator trait category in the given group of lakes.
Statistical comparisons of functional group relative abundances and site elevation among lake clusters were also made by Kruskal–Wallis tests followed by multiple comparison tests (step 4.1 in Fig. S1). Finally, a redundancy analysis (RDA) using Monte Carlo unrestricted 999 permutation tests (Legendre and Legendre 2012) was performed to establish a model associating chironomid trait compositions (100 lakes × 33 categories of 6 traits) to environmental characteristics of lakes (Table S1; step 4.2 in Fig. S1). All statistical tests and graphical displays were done using R 3.5.2 statistical software (R Core Team 2018).
Results
Distribution of subfossil chironomid along the elevation gradient
The first two PCA-axes of lake characteristics accounted for 34.5% and 21.2% of the total variance, respectively (Fig. 1). PCA revealed the existence of two distinct gradients: the first axis explained elevation, total organic carbon concentration in water (TOC), organic matter concentration in sediments (as loss on ignition; LOI) and climate parameters (mean air temperature of January and July: TJan and TJul, respectively), whereas the second axis was strongly correlated with specific conductance, pH, and concentrations of major cations (Ca, Mg, K; Fig. 1). The first gradient was interpreted as an elevation/humic pattern, where high elevation and clear water induced negative PCA1 scores and positive PCA2 scores, along which most lakes were spread (Fig. 1). The second gradient consisted of a conductivity/alkalinity pattern, where low-alkaline water and low conductivity induced negative PCA1 and PCA2 scores.
As expected, chironomid communities showed large changes across the elevation gradient (Fig. 2a), similar to that in the original taxonomic resolution applied by Larocque et al. (2001). Several genera including Pseudodiamesa, Abiskomyia and Corynoneura were found only in high-elevation lakes while others such as Heterotanytarsus, Dicrotendipes, Microtendipes, Cladotanytarsus and Polypedilum were associated with warmer lakes at low-elevation (Fig. 2a). A third group of genera including Tanytarsus, Paratanytarsus, Micropsectra, Corynocera, Sergentia and Procladius were instead evenly distributed across the elevation gradient (Fig. 2a). In total, 18 chironomid genera had sufficient information regarding their ecological requirements and biological traits in the European Chironomidae Database (plotted in black; Fig. 2a), and were kept for the second part of the study. For the remaining 11 taxa (plotted in grey), the trait database contained very little information, and these taxa were therefore excluded from further analyses. The exclusion of these taxa should, however, not have biased subsequent analyses as they consisted of rare taxa (accounting for less than ca. 20% of the total chironomid counts on average, Fig. 2b), and were evenly distributed along the elevation gradient (Fig. 2b).
Relationships among traits and trait-based functional groups
The first two axes of the FCA (F1–F2) applied to chironomid taxa and their functional traits (taxa-by-trait category dataset: 18 taxa × 33 categories of 6 traits) explained 15.4% and 12.8%, respectively, and accounted for 28.2% of the total variance (Fig. 3a). The traits “4th instar larval size” (SIZE1 and SIZE2 vs. SIZE4 and SIZE5, see Table 1 for full names of trait categories), and to a lesser extent “Food preference” (primary producers vs. other food resources) and “Feeding habits” showed the highest correlation ratios along the F1-axis (0.502, 0.152 and 0.134, respectively). These results revealed that some categories of these traits were well separated along this axis (Fig. 3a). Along the F2-axis, some categories of the traits “Feeding habits” (0.372, predators vs. other trophic guilds) and “Food preference” showed the highest correlation ratios (0.281; animals vs. other food resources; Fig. 3a), i.e. higher for the F2-axis than those for the F1-axis (Fig. 3a).
Relationships between chironomid taxa and their functional traits were analysed by quantifying the proximity of the taxa in the FCA factorial space. Using the Ward’s linkage method. This analysis revealed four functional groups (FG1–4), each sharing the same trait-categories (Fig. 3b). Functional group 1 (FG1; consisting of Chironomus, Microtendipes and Sergentia) was mainly composed of large, sediment-feeding larvae (according to trait distribution along F1–F2 axes; Fig. 3b), while group 2 (FG2; consisting of Pseudodiamesa, Cladopelma and Procladius) mainly consisted of intermediate-size predators (Fig. 3b). Functional group 3 (FG3; 11 taxa including among others Paratanytarsus, Heterotrissocladius and Micropsectra) was composed of moderate-size detritivorous and algae-feeding larvae (SHR, SCR), while group 4 (FG4) consisted only of Corynoneura and characterized by small larvae (SIZE1 and SIZE2, Fig. 3b).
Trait-based site clustering and relationships with environmental gradients
The first four axes of the FCA applied to the chironomid trait compositions of lakes (trait-by-lake dataset: 100 lakes × 33 categories of 6 traits) explained 22.5% of the total variability among lakes (Fig. 4). The Ward’s linkage method revealed four significantly distinct groups of lakes (i.e. cluster 1–4 in Fig. 4). Cluster 1 and 4 consisted of 32 and 11 lakes, respectively, and their locations were opposed along the second FCA axis (Fig. 4). Cluster 2 was composed of 54 lakes and located near the centre of F1-F2 biplot. Cluster 3 was composed of only 3 lakes and was opposed to other clusters along the first FCA axis (Fig. 4). The “IndVal” method highlighted significant indicator trait categories within each lake cluster, categories of traits characterising each cluster were listed in Fig. 4.
Relative abundances of genera from all the functional groups (FG1–4) significantly varied among all four lake clusters (Fig. 5a and Table 2). Chironomid larvae from FG3 were the most abundant in all lake clusters. Chironomid communities in clusters 1 and 2 were strongly predominated by taxa of FG3, while other functional groups were poorly represented (FG1, FG2 and FG4). The genera within FG1 were mostly represented in lakes of cluster 4 (Fig. 5a). Chironomid communities in cluster 3 were instead predominated by genera belonging to FG2 and FG3, while FG1 was poorly represented and FG4 not found (Fig. 5a and Table 2). Elevation also significantly varied among lake clusters, being higher for clusters 1 and 3 than for lake clusters 2 and 4 (Fig. 5b).
Overall, the RDA model linking trait categories to environmental variables explained 16.9% of the total variation in chironomid trait compositions of lakes. Most of the variance contained in both first and second RDA axes is described by elevation and TOC (e.g. low elevation and humic sites inducing positive RDA1 and negative RDA2 scores; Fig. 6), suggesting a strong elevation/humic pattern in chironomid trait compositions among lakes. The traits “4th instar larval size”, “Food preference” and “Feeding habits” were separated along RDA axes, and chironomid trait compositions of lakes were well distributed across the elevation gradient. Some indicator traits previously highlighted by the “IndVal” approach, such as large larvae (SIZE5), detritivore feeding habits (DEFEE, SHR) and food preference (SEDMIC and WOOD; Fig. 6), were typical for the chironomid communities of the low-elevation lakes (Clusters 2 and 4; Fig. 6). In contrast, several trait categories, including small larvae (SIZE1), predator and herbivore feeding habits (PRED, FFEEDT and SCR) and associated food preference (MICPHY, MACINVERT and MICINV), were characteristic for lakes high elevations (clusters 1 and 3). The remaining trait categories were judged indifferent, i.e. not clearly associated with any environmental descriptor (Fig. 6).
Discussion
We examined how chironomid biodiversity of subarctic lakes is distributed across an elevation gradient in temperature and landscape characteristics. Based on spatial patterns of functional diversity of chironomids in subarctic lakes, we show that some categories of the traits larval size, food preference, and feeding habits distinctly separate along the studied gradients. Our study therefore suggests complex interplays between direct (e.g. temperature, ice-free period) and indirect climate-related processes (e.g. elevation-driven changes in vegetation/habitats) on the functional diversity of chironomid communities. Resource availability in lakes is strongly linked to vegetation composition/cover, and predominant traits related to resource exploitation in chironomid communities are therefore strongly constrained by landscape characteristics. Finally, our results also show that short life cycles could facilitate the development of viable population in northern and high-elevation lakes where the ice-free period is a limiting factor, thus contradicting patterns showing smaller organisms in warmer environments reported for other invertebrates.
Trait-based functional groups of chironomid taxa
In our training-set, four different groups of chironomid genera were defined by three functional traits, i.e. larval size, feeding habits, and food preferences. These results support previous findings that traits such as food types and larval size are the main contributors to trait variability among aquatic invertebrates, (see also Usseglio-Polatera et al. 2000; Poff et al. 2006; Desrosiers et al. 2019), and also for chironomid larvae (Cummins 1973; Berg 1995; Luoto and Nevalainen 2015). As expected, these trait-based groupings corresponded poorly with known phylogenetic links among main chironomid subfamilies (i.e. Diamesinae, Orthocladiinae, Tanypodinae and Chironominae; Serra et al. 2016). Finally, these traits represent biological features that directly influence ecosystem functions (Hulot et al. 2000), such as energy flow within food webs, making these functional groups excellent indicators of the effects of climate change on the structure and function of subarctic lakes.
Direct and indirect climate processes shape assemblage structure
Overall, our results suggest complex interplays between direct (e.g. temperature) and indirect climate processes (e.g. elevation-driven changes in vegetation/habitats) on the functional diversity of chironomid communities. The observed strong effects of elevation in structuring trait composition of chironomid communities demonstrate that both temperature and large-scale landscape attributes (e.g. land use and vegetation composition/cover) strongly influence the functional diversity of northern aquatic invertebrate assemblages, confirming similar findings for temperate lakes (Poff 1997; Richards et al. 1997; Statzner et al. 2004; Townsend et al. 2003; Motta and Massaferro 2019). However, our results also show a strong influence on trait composition of habitat characteristics, such as water depth, suggesting that some functional features of chironomid communities also showed strong relationships with ecosystem types and local environmental conditions (c.f. Richards et al. 1996, Motta and Massaferro 2019). This is corroborated by Johnson and Goedkoop (2002), who found that local-scale habitat descriptors had a strong unique effect on the functional composition of lake’s macroinvertebrate communities.
Effects of elevation gradient on trait composition
Specifically, our study shows that some categories of the traits larval size, food preference, and feeding habits are significantly separated along an elevation/humic pattern, and trait composition of chironomid communities were therefore distributed across the elevation gradient in temperature and landscape characteristics. At low elevation sites, communities were mainly composed of large bodied larvae belonging to the shredder and sediment feeding trophic guilds, whereas larvae in higher elevation lakes tended to be smaller and commonly had predatory feeding mode. Our results thereby support studies that have shown that temperature might play both a direct and/or indirect role in selection for smaller chironomid taxa in high-elevation and cold-water lakes (Hamerlík et al. 2017; Motta and Massaferro 2019). Furthermore, low-elevation lakes with forested catchments typically have more allochthonous organic material that largely covaries with the forest cover (Fig. 1), and therefore tend to have more sediment-feeding taxa than lakes above the tree line (Luoto and Nevalainen 2015). As expected, resource availability in lakes is strongly linked to catchment characteristics, and traits related to resource exploitation in chironomid communities are therefore well constrained by landscape characteristics.
Our results, however, conflict with patterns previously reported for other aquatic invertebrate groups showing that organisms tend to be smaller and better dispersers in warmer environments (c.f. Chown and Gaston 2010; Hayden et al. 2017). Indeed, these studies suggested that increasing body size in colder environments is the result of prolonged lifecycles (slower growth) in response to low temperatures (Angilletta et al. 2004). However, short life cycles could also facilitate the development of viable populations in subarctic and high-elevation lakes where the ice-free period is a dispersal-limiting factor (Nolte and Hoffmann 1992). In addition, the ability to feed on available food resources, indirectly linked to catchment characteristics (Heino 2008; Kivilä et al. 2019) might also influence the persistence of taxa. Interestingly, our analysis also show that lake water depth contributes to shape structure and trait composition of chironomid communities, but additional studies are needed to better understand how lake type (including deep vs. shallow lakes) could affect the trait composition of chironomid communities and their responses to climate change.
Implications for climate change effects, and future challenges
Predicted climate change and concurrent climate-induced changes in catchment characteristics will affect Arctic and subarctic lakes and their unique biodiversity (Heino 2008; Lento et al. 2019). Using a space-for-time approach, our study can help to predict future trajectories of freshwater biodiversity in these subarctic lakes. As a consequence of climate warming, migration of chironomid species towards higher elevation sites will be induced by both temperature tolerance of taxa and changes in catchment characteristics (i.e. vegetation composition/cover and associated patterns in organic matter run-off). Hence, the highest elevation lakes in subarctic landscapes might lose their typical cold-adapted chironomid taxa along with their functional attributes, when these are dispersal-limited and replaced by better adapted species from lower elevations or southern latitudes (see Lento et al. 2019).
In addition to space-for-time substitutions, paleolimnological studies are powerful approaches to reconstruct the effects of past environmental changes archived in lake sediments (Smol 1992) and provide excellent opportunities to explore the responses of functional diversity to past climate change. In fact, these paleo-ecological studies cover far more than the variation in climate that we have seen so far since the 1990s, and analysis of subfossil chironomid communities has provided valuable insights on past environmental and climatic impacts on lake’s biodiversity (see also Heiri et al. 2003; Luoto and Nevalainen 2018). Most of these studies have, however, applied a taxonomy-based approach, and impacts on the trait composition have been so far largely neglected. Our study, combining trait-based and subfossil communities, takes this one step further by quantifying effects on trait compositions, and contributes to the growing hierarchical understanding of the roles of environmental and climatic changes in structuring functional diversity of chironomid communities.
Trait plasticity is well documented for insects (Cummins 1973; Nylin and Gotthard 1998) leading to differences between actual and theoretical functional traits of species (Bennett et al. 2015) and a potential limitation to the use of chironomid traits as indicators of environmental change. Moreover, most aquatic insects are assumed omnivorous in their early instars (Cummins and Klug 1979), and classification based on feeding strategies of later instars may therefore mismatch with those of all larval stages. Serra et al. (2017) also reported large differences in chironomid traits between North America and Europe and suggested alternative adaptation to similar ecological environments. Functional plasticity implies that many combinations of traits appear for a single species depending on environmental conditions (Jourdan et al. 2019). A better understanding of functional trait plasticity and how it is affected by environmental changes could strengthen the reliability of ecological inferences of trait-based chironomid properties.
Conclusion
Our study shows that complex interplays between direct (e.g. temperature) and indirect climate processes (e.g. elevation-driven changes in vegetation/habitats) affect the functional diversity of chironomid communities. As expected, traits related to resource exploitation (i.e. food preference and feeding habits) in chironomid communities are indirectly constrained by landscape characteristics as these strongly impact on resource availability. Furthermore, our findings suggest that short generation times could facilitate the development of viable populations in subarctic and high elevation lakes where the short ice-free period is a limiting factor, thus contradicting patterns reported for other invertebrates. These results also highlighted how species traits can contribute to a mechanistic understanding of functional changes along climate gradients and provide new knowledge about functional changes in subarctic lakes under climate-warming scenarios. We envision that functional approaches in paleolimnological studies have the inherent strength to forecast future trajectories of freshwater biodiversity in high-latitude lakes by learning from the past changes archived in their sediments.
Data availability statement
The datasets utilized in this study are available upon request and will be archived in Pangea (Larocque et al. 2001).
References
AMAP (2011) Snow, water, ice and permafrost in the ARCTIC (SWIPA): climate change and the cryosphere. Arctic Monitoring and Assessment Programme (AMAP), Oslo
Angert AL, Crozier LG, Rissler LJ, Gilman SE, Tewksbury JJ, Chunco AJ (2011) Do species’ traits predict recent shifts at expanding range edges? Ecol Lett 14(7):677–689
Angilletta MJ, Steury TD, Sears MW (2004) Temperature, growth rate, and body size in ectotherms: fitting pieces of a life-history puzzle. Integr Comp Biol 44(6):498–509
Bennett KD (1996) Determination of the number of zones in a biostratigraphical sequence. New Phytol 132(1):155–170
Bennett NL, Severns PM, Parmesan C, Singer MC (2015) Geographic mosaics of phenology, host preference, adult size and microhabitat choice predict butterfly resilience to climate warming. Oikos 124(1):41–53
Berg MB (1995) Larval food and feeding behaviour. In: Armitage PD, Cranston PS, Pinder LCV (eds) The Chironomidae. Chapman and Hall, New York, pp 136–168
Bonada N, Dolédec S, Statzner B (2007) Taxonomic and biological trait differences of stream macroinvertebrate communities between mediterranean and temperate regions: implications for future climatic scenarios. Glob Change Biol 13(8):1658–1671
Brodersen KP, Quinlan R (2006) Midges as palaeoindicators of lake productivity, eutrophication and hypolimnetic oxygen. Quat Sci Rev 25(15–16):1995–2012
Chevenet F, Doléadec S, Chessel D (1994) A fuzzy coding approach for the analysis of long-term ecological data. Freshwat Biol 31(3):295–309
Chown SL, Gaston KJ (2010) Body size variation in insects: a macroecological perspective. Biol Rev 85(1):139–169
Cummins KW (1973) Trophic relations of aquatic insects. Annu Rev Entomol 18(1):183–206
Cummins KW, Klug MJ (1979) Feeding ecology of stream invertebrates. Annu Rev Ecol Syst 10(1):147–172
Daufresne M, Lengfellner K, Sommer U (2009) Global warming benefits the small in aquatic ecosystems. Proc Natl A Sci 106(31):12788–12793
Desrosiers M, Usseglio-Polatera P, Archaimbault V, Larras F, Méthot G, Pinel-Alloul B (2019) Assessing anthropogenic pressure in the St. Lawrence River using traits of benthic macroinvertebrates. Sci Total Environ 649:233–246
Dı́az S, Cabido M (2001) Vive la différence: plant functional diversity matters to ecosystem processes. Trends Ecol Evol 16(11):646–655
Dufrêne M, Legendre P (1997) Species assemblages and indicator species: the need for a flexible asymmetrical approach. Ecol Monogr 67(3):345–366
Eggermont H, Heiri O (2012) The chironomid-temperature relationship: expression in nature and palaeoenvironmental implications. Biol Rev 87(2):430–456
Gardner JL, Peters A, Kearney MR, Joseph L, Heinsohn R (2011) Declining body size: a third universal response to warming? Trends Ecol Evol 26(6):285–291
Hamerlík L, Svitok M, Novikmec M, Veselská M, Bitušík P (2017) Weak altitudinal pattern of overall chironomid richness is a result of contrasting trends of subfamilies in high-altitude ponds. Hydrobiologia 793(1):67–81
Hayden B, Myllykangas JP, Rolls RJ, Kahilainen KK (2017) Climate and productivity shape fish and invertebrate community structure in subarctic lakes. Freshwat Biol 62(6):990–1003
Heino J (2008) Patterns of functional biodiversity and function-environment relationships in lake littoral macroinvertebrates. Limnol Oceanogr 53(4):1446–1455
Heino J, Virkkala R, Toivonen H (2009) Climate change and freshwater biodiversity: detected patterns, future trends and adaptations in northern regions. Biol Rev 84(1):39–54
Heiri O (2004) Within-lake variability of subfossil chironomid assemblages in shallow Norwegian lakes. J Paleolimn 32(1):67–84
Heiri O, Wick L, van Leeuwen JFN, van der Knaap WO, Lotter AF (2003) Holocene tree immigration and the chironomid fauna of a small Swiss subalpine lake (Hinterburgsee, 1515 m asl). Palaeogeogr Palaeoecolo 189(1):35–53
Heiri O, Brooks SJ, Birks HJB, Lotter AF (2011) A 274-lake calibration data-set and inference model for chironomid-based summer air temperature reconstruction in Europe. Quat Sci Rev 30(23–24):3445–3456
Hevia V, Martín-López B, Palomo S, García-Llorente M, de Bello F, González JA (2017) Trait-based approaches to analyze links between the drivers of change and ecosystem services: synthesizing existing evidence and future challenges. Ecol Evol 7(3):831–844
Hooper DU, Chapin FS, Ewel JJ, Hector A, Inchausti P, Lavorel S, Wardle DA (2005) Effects of Biodiversity on ecosystem functioning: a consensus of current knowledge. Ecol Monogr 75(1):3–35
Hulot FD, Lacroix G, Lescher-Moutoué F, Loreau M (2000) Functional diversity governs ecosystem response to nutrient enrichment. Nature 405(6784):340
Johnson RK, Goedkoop W (2002) Littoral macroinvertebrate communities: spatial scale and ecological relationships. Freshwat Biol 47:1840–1854
Jourdan J, Baranov V, Wagner R, Plath M, Haase P (2019) Elevated temperatures translate into reduced dispersal abilities in a natural population of an aquatic insect. J Anim Ecol 88(10):1498–1509
Kivilä EH, Luoto TP, Rantala MV, Kiljunen M, Rautio M, Nevalainen L (2019) Environmental controls on benthic food web functions and carbon resource use in subarctic lakes. Freshwat Biol 64(4):643–658
Larocque I, Hall RI, Grahn E (2001) Chironomids as indicators of climate change: a 100-lake training set from a subarctic region of northern Sweden (Lapland). J Paleolimnol 26(3):307–322
Lau DCP, Christoffersen KS, Erkinaro J, Hayden B, Heino J, Hellsten S, Holmgren K, Kahilainen KK, Kahlert M, Karjalainen SM, Karlsson J, Forsström L, Lento J, Mjelde M, Ruuhijärvi J, Sandøy S, Schartau AK, Svenning M-A, Vrede T, Goedkoop W (2020) Multitrophic biodiversity patterns and environmental descriptors of sub-Arctic lakes in northern Europe. Freshwat Biol. https://onlinelibrary.wiley.com/doi/abs/10.1111/fwb.13477(in press)
Lavergne S, Mouquet N, Thuiller W, Ronce O (2010) Biodiversity and climate change: integrating evolutionary and ecological responses of species and communities. Annu Rev Ecology Evol S 41:321–350
Legendre P, Legendre LFJ (2012) Numerical ecology. Elsevier, Amsterdam
Lento J, Goedkoop W, Culp J et al (2019) Report state of the Arctic Freshwater Biodiversity. Conservation of Arctic Flora and Fauna International Secretariat, Akureyri. ISBN 978-9935-431-77-6
Luoto TP, Nevalainen L (2015) Climate-forced patterns in midge feeding guilds. Hydrobiologia 742(1):141–152
Luoto TP, Nevalainen L (2018) Temperature-precipitation relationship of the Common Era in northern Europe. Theor Appl Climatol 132:933–938
MacLean SA, Beissinger SR (2017) Species’ traits as predictors of range shifts under contemporary climate change: a review and meta-analysis. Glob Change Biol 23(10):4094–4105
McLean M, Auber A, Graham NAJ, Houk P, Villéger S, Violle C, Mouillot D (2019) Trait structure and redundancy determine sensitivity to disturbance in marine fish communities. Glob Change Biol 25:3424–3437
Motta L, Massaferro J (2019) Climate and site-specific factors shape chironomid taxonomic and functional diversity patterns in northern Patagonia. Hydrobiologia 839:131–143
Nolte U, Hoffmann T (1992) Fast life in cold water: Diamesa incallida (Chironomidae). Ecography 15(1):25–30
Northington RM, Keyse MD, Beaty SR, Whalen SC, Sokol ER, Hershey AE (2010) Benthic secondary production in eight oligotrophic arctic Alaskan lakes. J N Am Benthol Soc 29(2):465–479
Nylin S, Gotthard K (1998) Plasticity in life-history traits. Annu Rev Entomol 43(1):63–83
Picazo F, Millán A, Dolédec S (2012) Are patterns in the taxonomic, biological and ecological traits of water beetles congruent in Mediterranean ecosystems? Freshwat Biol 57(11):2192–2210
Poff NL (1997) Landscape filters and species traits: towards mechanistic understanding and prediction in stream ecology. J N Am Benthol Soc 16(2):391–409
Poff NL, Olden JD, Vieira NKM, Finn DS, Simmons MP, Kondratieff BC (2006) Functional trait niches of North American lotic insects: traits-based ecological applications in light of phylogenetic relationships. J N Am Benthol Soc 25(4):730–755
R Core Team (2018) R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. ISBN 3-900051-07-0. http://www.R-project.org
Reynolds CS, Huszar V, Kruk C, Naselli-Flores L, Melo S (2002) Towards a functional classification of the freshwater phytoplankton. J Plankton Res 24(5):417–428
Richards C, Johnson LB, Host GE (1996) Landscape-scale influences on stream habitats and biota. Can J Fish Aquat Sci 53(S1):295–311
Richards C, Haro R, Johnson L, Host G (1997) Catchment and reach-scale properties as indicators of macroinvertebrate species traits. Freshwat Biol 37(1):219–230
Schmera D, Heino J Podani, Erős T, Dolédec S (2017) Functional diversity: a review of methodology and current knowledge in freshwater macroinvertebrate research. Hydrobiologia 787:27–44
Serra SRQ, Cobo F, Graça MAS, Dolédec S, Feio MJ (2016) Synthesising the trait information of European Chironomidae (Insecta: Diptera): towards a new database. Ecol Indic 61:282–292
Serra SRQ, Graça MAS, Dolédec S, Feio MJ (2017) Chironomidae of the Holarctic region: a comparison of ecological and functional traits between North America and Europe. Hydrobiol 794(1):273–285
Smol JP (1992) Paleolimnology: an important tool for effective ecosystem management. J Aquat Ecosyst Health 1(1):49–58
Statzner B, Dolédec S, Hugueny B (2004) Biological trait composition of European stream invertebrate communities: assessing the effects of various trait filter types. Ecography 27(4):470–488
Townsend CR, Dolédec S, Norris R, Peacock K, Arbuckle C (2003) The influence of scale and geography on relationships between stream community composition and landscape variables: description and prediction. Freshw Biol 48(5):768–785
Urban MC, Tewksbury JJ, Sheldon KS (2012) On a collision course: competition and dispersal differences create no-analogue communities and cause extinctions during climate change. Proc R Soc B-Biol Sci 279:2072–2080
Usseglio-Polatera P, Bournaud M, Richoux P, Tachet H (2000) Biological and ecological traits of benthic freshwater macroinvertebrates: relationships and definition of groups with similar traits. Freshwat Biol 43(2):175–205
van Hardenbroek M, Heiri O, Wilhelm MF, Lotter AF (2010) How representative are subfossil assemblages of Chironomidae and common benthic invertebrates for the living fauna of Lake De Waay, the Netherlands? Aquat Sci 73(2):247–259
van Kleef H, Verberk WCEP, Kimenai FFP, van der Velde G, Leuven RSEW (2015) Natural recovery and restoration of acidified shallow soft-water lakes: Successes and bottlenecks revealed by assessing life-history strategies of chironomid larvae. Basic Appl Ecol 16(4):325–334
Verbruggen F, Heiri O, MeriläInen JJ, Lotter AF (2011) Subfossil chironomid assemblages in deep, stratified European lakes: relationships with temperature, trophic state and oxygen. Freshwat Biol 56(3):407–423
Walker IR (2001) Midges: chironomidae and related diptera. In: Smol JP, Birks HJB, Last WM (eds) Tracking environmental change using lake sediments. Springer, Dordrecht, pp 43–66
Ward J (1963) Hierarchical grouping to optimize an objective function. J Am Stat Assoc 58(301):236–244
Wiederholm T (ed) (1983) Chironomidae of the Holarctic region: keys and diagnoses. Part 1 Larvae Entomologica Scandinavica Supplement No. 19:1–457
Acknowledgements
Open access funding provided by Swedish University of Agricultural Sciences. We acknowledge Isabelle Larocque for sharing this 100-lake training set with us, and Oliver Heiri for valuable comments on earlier versions of this manuscript.
Author information
Authors and Affiliations
Contributions
All authors contributed to the development of theory; SB analysed data; SB wrote the initial draft of the paper and the co-author contributed substantially to revisions.
Corresponding author
Ethics declarations
Conflict of interest
The authors have no conflict of interest to declare.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Handling Editor: Jamie Kneitel.
Electronic supplementary material
Below is the link to the electronic supplementary material.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Belle, S., Goedkoop, W. Functional diversity of chironomid communities in subarctic lakes across gradients in temperature and catchment characteristics. Limnology 22, 5–16 (2021). https://doi.org/10.1007/s10201-020-00624-0
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10201-020-00624-0