Fabienne Maihoff*
Simone Sahler
Simon Schoger
Kristof BrenzingerKatharina KallnikNikki SauerLukas Bofinger
Thomas SchmittSabine S. Nooten†
Alice Classen†- Department of Animal Ecology and Tropical Biology, University of Würzburg, Am Hubland, Würzburg, Germany
Alpine bumble bees are the most important pollinators in temperate mountain ecosystems. Although they are used to encounter small-scale successions of very different climates in the mountains, many species respond sensitively to climatic changes, reflected in spatial range shifts and declining populations worldwide. Cuticular hydrocarbons (CHCs) mediate climate adaptation in some insects. However, whether they predict the elevational niche of bumble bees or their responses to climatic changes remains poorly understood. Here, we used three different approaches to study the role of bumble bees’ CHCs in the context of climate adaptation: using a 1,300 m elevational gradient, we first investigated whether the overall composition of CHCs, and two potentially climate-associated chemical traits (proportion of saturated components, mean chain length) on the cuticle of six bumble bee species were linked to the species’ elevational niches. We then analyzed intraspecific variation in CHCs of Bombus pascuorum along the elevational gradient and tested whether these traits respond to temperature. Finally, we used a field translocation experiment to test whether CHCs of Bombus lucorum workers change, when translocated from the foothill of a cool and wet mountain region to (a) higher elevations, and (b) a warm and dry region. Overall, the six species showed distinctive, species-specific CHC profiles. We found inter- and intraspecific variation in the composition of CHCs and in chemical traits along the elevational gradient, but no link to the elevational distribution of species and individuals. According to our expectations, bumble bees translocated to a warm and dry region tended to express longer CHC chains than bumble bees translocated to cool and wet foothills, which could reflect an acclimatization to regional climate. However, chain lengths did not further decrease systematically along the elevational gradient, suggesting that other factors than temperature also shape chain lengths in CHC profiles. We conclude that in alpine bumble bees, CHC profiles and traits respond at best secondarily to the climate conditions tested in this study. While the functional role of species-specific CHC profiles in bumble bees remains elusive, limited plasticity in this trait could restrict species’ ability to adapt to climatic changes.
Introduction
Mountains are hotspots of biodiversity with many endemic and cold-adapted species (Rahbek et al., 2019; Trew and Maclean, 2021). Yet climate change affects biodiversity by changing spatial ranges, phenology and species interaction (Bellard et al., 2012; Halsch et al., 2021). Since the temperature increase associated with climate change is likely to be more pronounced at higher elevations than at lower elevations (Pepin et al., 2015), our urgent concern should be to understand the responses of these particularly threatened ecosystems, which can also serve as early warning system.
Alpine bumble bees (Apidae: Bombus) are the most efficient and widespread pollinators in temperate mountains (Bingham and Orthner, 1998; Gorenflo et al., 2017). Within the genus Bombus, species differ substantially in their preferred elevational niche (Rasmont et al., 2015; Minachilis et al., 2020; Sponsler et al., 2022a), with few species being restricted to the cool and exposed conditions above the tree line. Within a species, populations and individuals can face substantial differences in their climatic environments along mountain slopes or during foraging flights under the typically fluctuating weather conditions of the mountains (Sponsler et al., 2022a), occasionally, when the cloud cover shifts abruptly, within minutes. Despite their broad temperature tolerance (Peters et al., 2016), many bumble bee species are sensitive to climate change. Many species are globally declining (Potts et al., 2010; Arbetman et al., 2017) and shift their range to higher elevations or toward the poles to track their preferred temperature niche and to avoid heat and desiccation stress (Kerr et al., 2015; Pyke et al., 2016; Marshall et al., 2020). The magnitude of such shifts varies greatly among species and is likely associated with species capacity to withstand new environmental conditions (Kerr et al., 2015; Pyke et al., 2016; Marshall et al., 2020; Maihoff et al., 2022b). One factor known to influence adaptation to climate change is desiccation resistance. Although studies focusing on gene expression and heat resistance sugguest that the ability to tolerate desiccation, caused by high temperatures and lower percipitation, is relevant for determining bumble bees’ response to climate change (Jackson et al., 2020; Maebe et al., 2021; Martinet et al., 2021) and, consequently, the maintenance of their pollination service, mechanisms that determine desiccation resistance can be versatile. Desiccation resistance in insects depend on the ability of insects cuticle to prevent water loss (Chown et al., 2011). Thus, differences in physiological traits modulating water loss are likely to contribute to their uneven ability to cope with environmental conditions along the elevational gradient (Williams et al., 2010; Aguirre-Gutiérrez et al., 2016; Wong et al., 2019; Stemkovski et al., 2020; Maebe et al., 2021). However, such traits have not yet been identified in bumble bees.
The cuticular hydrocarbons (CHCs) composition on the insects’ cuticle is such a trait. CHCs provide waterproofing, and comprise a complex mixture of n-alkanes, methyl-branched alkanes and unsaturated hydrocarbons that cover the cuticle (Gibbs et al., 1991; Gibbs, 1998; Blomquist and Bagnères, 2010). Two features are particularly relevant in the context of waterproofing (Gibbs and Pomonis, 1995): the proportion of saturated components and the chain length of hydrocarbons. Saturated hydrocarbons (n-alkanes and mono-methyl-branched alkanes) aggregate more tightly than cis-configurated unsaturated hydrocarbons (n-alkenes) due to increased van-der Waals forces. Thus, a higher proportion of saturated components in the CHC profile enhances waterproofing (Gibbs and Pomonis, 1995). Aggregating forces between hydrocarbons also increase with chain length; CHC profiles composed of on average longer hydrocarbons thus provide a better protection against desiccation than profiles with shorter CHCs (Menzel et al., 2017b).
Various studies on insects confirm that the composition of CHC profiles can reflect climate adaptation on an inter- and intraspecific level (Rajpurohit et al., 2017; Menzel et al., 2017a; Michelutti et al., 2018; Sprenger et al., 2019; Mayr et al., 2021). For example ant species from habitats with high rainfall produce various alkenes, alkadienes and methyl-branched alkenes, i.e., substance classes with reduced protection against desiccation stress (Menzel et al., 2017a), while Drosophila populations from warmer regions produce longer CHC chains than population from colder regions (Rajpurohit et al., 2017). Thus, it can be hypothesized that bumble bee species that differ in their preferred environmental niche have CHC profiles that reflect the degree of desiccation stress. Importantly, in some insect species CHC profiles alter plastically in a short-term under changing climatic conditions (Menzel et al., 2018; Sprenger et al., 2018). This intraspecific acclimatization of CHC changes can be elicited even within hours (Stinziano et al., 2015) and can effectively reduce heat stress (Menzel et al., 2018). Thus, the capacity for a short-term acclimatization response, alongside long-term adaptation, may be crucial in determining the climatic range along the elevational gradient of a species and for its survival under changing climatic conditions.
In addition to the function of water balance regulation, CHCs are involved in insect communication. In social insects communication is crucial for colony maintenance. CHCs are important for nestmate recognition and task signalization, in, e.g., wasps (Polistes) honey bees (Apis), stingless bees and ants (Akino et al., 2004; Dani et al., 2005; Nunes et al., 2008; Leonhardt et al., 2016; Maihoff et al., 2022a) and also in caste recognition and signalizing of health status (Leonhardt et al., 2016; Beani et al., 2019). Thus, CHCs that are involved in species communication and recognition of nest affiliations should be under strong pressure to be maintained, as any deviations may have negative effects on fitness and colony survival. Even though the role of CHCs in inter- and intraspecific bumble bee communication remains unclear, the potential dual function of methyl-alkanes, alkanes and alkenes in communication and waterproofing could potentially limit species capacity to adapt or acclimatize to new climatic conditions (Dani et al., 2005; Colazza et al., 2007; Lacey et al., 2008; Sprenger et al., 2019; Awater-Salendo et al., 2020). Furthermore, like many other functional traits, CHCs may be constrained by phylogeny, which should be considered when analyzing their variation within a clade (De Oliveira et al., 2011; Flynn et al., 2011; Kellermann et al., 2012; Kather and Martin, 2015; Menzel et al., 2017a). In general, less is known about CHC profile variation within the clade of bumble bees, but in the context of studies addressing climate change, genetic variation at cuticle formation suggests a likely role in climatic adaptation (Jackson et al., 2020; Straub et al., 2022).
In this study, we used three different approaches to study the role of CHCs in climatic adaptation across and within bumble bee species. We hypothesized that the CHC profiles of bumble bees exposed to higher desiccation stress at higher temperatures have longer chains and a higher proportion of saturated components. A 1.2 km elevational gradient serves as a natural model system that provides variation in temperatures and associated relative humidity patterns also expected under climate change (up to +5°C by the end of the century; Körner, 2007; IPCC et al., 2018) as temperature declines linearly with elevation (∼0.5°C per 100 m increase in elevation; Körner, 2007). Furthermore, the tree line, i.e., the transition from forest to shrub or grassland, represents a critical threshold, as outside the protection of the tree canopy the various abiotic conditions, to which bumble bees are exposed, become more extreme (Beals, 1969; Slatyer and Noble, 1992). We first compared the CHC profiles of six bumble bee species that differed in their average elevational distribution and hypothesized that CHC composition and assumed to be climate-associated chemical traits (i.e., mean chain length and proportion of saturated components) correlate with the preferred climate niche of the respective species (Mayr et al., 2021). For simplicity, we refer to these traits climate-associated chemical traits in the following. We also consider phylogenetic constraints, assuming that traits are more similarly expressed in closely related species than in distant ones (Flynn et al., 2011).
In a second approach, we investigated intraspecific variation in the CHC profile of Bombus pascuorum (Scopoli, 1763)—a species, which was previously restricted to lower elevation and which has expanded its range to higher elevation during recent climate change (Marshall et al., 2020; Maihoff et al., 2022b). Specifically, we hypothesized that variation in climate-associated chemical traits on the cuticle correlates systematically with environmental change along the elevational gradient and becomes most pronounced above the tree line where abiotic conditions are more extreme.
The third approach entails a field translocation experiment, to test the acclimatization capacity of CHCs in Bombus lucorum (Linnaeus, 1761), a species known to prefer cold and humid forest habitats (Geue and Thomassen, 2020). Field translocation experiments are powerful tools to assess the extent of species’ acclimatization and genetic determination, as they elucidate the factors that limit species distributions and have the potential to predict species range shifts in a changing climate (Nooten and Hughes, 2017). We transferred young bumble bee colonies derived from queens, that were collected in the foothills of a comparably cold and wet mountainous region (Berchtesgadener Land, Germany) to (a) higher elevations within the mountainous region and (b) a warmer and drier region in Bavaria’s lowlands (Würzburg, Germany). We hypothesized that after 4 weeks differences in CHC profiles and climate-associated chemical traits are indicative for new nesting sites in the way that in the warmer and dryer regions CHC profiles are characterized by a high proportion of saturated components and on average longer mean chain lengths.
Materials and methods
Study area and bumble bee sampling
Approach 1: CHC differences between species occupying different elevational niches
To test whether differences in CHC profiles between species can be explained by species’ elevational niche, bumble bee workers (foraging at flowers) of six species were collected within the National Park Berchtesgaden and its close vicinity (Lat: 47.5477, Lon: 12.9247) for one month (29.07.-31.08.) in two years (2019 and 2020). We selected six coexisting species in the study region to represent differences in elevational distribution. Bombus mucidus Gerstäcker (1869) and Bombus monticola Smith (1849) mainly occur above the tree line at the highest sites of our study region, Bombus soroeensis Fabricius (1777), and Bombus wurflenii Radoszkowski (1859) across the entire studied elevational gradient, and B. pascuorum and B. lucorum at lower elevations below the tree line on average at an elevation of 1,200 m.a.s.l. [monitored in Maihoff et al., 2022b; see Figure 1 and further information in Supplementary material S1]. The National Park is located within the limestone Alps in southern east Germany—a region characterized by coniferous forest, alpine meadows, and mountain pastures (Konnert and Siegrist, 2000). The tree line in the study area is at an elevation of about 1,500 m.a.s.l. (Köstler and Mayer, 1970). We collected individuals from sites (60 m × 60 m) representative of the species’ range and also aimed for different species at the same sites. Therefore, mainly individuals above 1,300 m were collected, except for B. lucorum and B. pascuorum, where we included individuals from lower elevations according to their distribution. In total we collected 65 individuals [B. soroeensis (n = 9), B. mucidus (n = 9), B. lucorum (n = 12), B. monticola (n = 12), B. wurflenii (n = 11), and B. pascuorum (n = 12)] from 18 sites covering a gradient from 641 to 2,114 m.a.s.l. (Figure 1). Note that the species were selected not only according to their preferred distribution, but also according to their relatedness, which does not correlate with the distribution (see Figure 2A for phylogeny).
Figure 1. Study design. (A) Study region and sites in Berchtesgaden (Bavaria, Germany). Elevation level is indicated in shades (at 250 m intervals), with darker shades signaling increasing elevation. The National Park border is represented with a green line. Points show study sites. Red sites were used for bumble bee collection in approach 1 and 2, while yellow, green, and blue sites represent sites which were additionally used for the translocation approach (approach 3). (B) Schematic description of the colony translocation (approach 3). Queens collected from the valleys in spring in our study region were reared in the lab and then settled with their first offspring in nest boxes. Numbers indicate new colony locations. In the cold and wet region of origin: mountain foothills (1), mountain mid (2), mountain high (3) (see also Panel A, where mountain sites are shown in respective colors), and in the warm and dry region (4), where the mountain sites are shown according to their color in Panel A. The table lists the temperature regimes to which the colonies were exposed either over the course of the translocation (=Texp) or within the week before sampling along the elevational gradient (=Tacclim). Color code in the map refers to multi annual means of air temperature conditions. Maps were produced in QGIS. Data was obtained from https://www.lfu.bayern.de/umweltdaten, https://search.earthdata.nasa.gov, and http://www.dwd.de.
Figure 2. Interspecific differences and phylogenetic signals in climate-associated chemical traits of six bumble bee species. (A) Phylogenetic tree extracted from phylogeny of bumble bees estimated from Bayesian analysis of combined sequence data from five gene fragments (16S, opsin, ArgK, EF-1α, and PEPCK) by Cameron et al. (2007). (B) Individual-based Cluster dendrogram showing the similarity of the total CHC profiles composition based on Bray-Curtis distances. (C) Proportion of saturated components and (D) mean chain lengths within the CHC profiles of six bumble bee species. Letters indicate significant differences between species. (E,F) Phylogenetic correlograms for the two traits are presented: (E) proportion of saturated components, (F) mean chain length. The solid bold black line shows the Moran’s I index of autocorrelation, and the dashed black lines the lower and upper bounds of the 95% confidence interval. The horizontal black line indicates the expected value of Moran’s I under the null hypothesis of no phylogenetic autocorrelation. The colored bar shows whether the autocorrelation is significant (based on the confidence interval): red for significant positive autocorrelation, black for nonsignificant autocorrelation, and blue for significant negative autocorrelation. The vertical gray lines indicate the predefined intraspecific distance (0.02) and the interspecific distance (≥0.5). Thus, the interpretation of modeled Moran’s I index is limited between a phylogenetic distance of 0.02 and 0.5.
Approach 2: Intraspecific CHC variation of Bombus pascuorum along an elevational gradient
For analyzing intraspecific variation in CHCs, we collected B. pascuorum workers foraging at flowers (n = 58) from 12 sites (60 m × 60 m) across an elevational gradient from 641 to 2,032 m.a.s.l. within the National Park Berchtesgaden and its close vicinity (Lat: 47.5477, Lon: 12.9247). We selected B. pascuorum because of its high abundance in the study region, and its characteristics as a species expanding its range toward higher elevations under recent climate change (Maihoff et al., 2022b).
Approach 3: Intraspecific CHC variation in translocated workers of Bombus lucorum
To test whether bumble bees can adjust their CHC profile to new climatic conditions, we translocated young self-reared colonies, originally from cold and wet climate to a warm and dry climate and to distinct elevations. For this, we collected queens, emerging from diapause, of B. lucorum in the alpine study region in spring and reared them in the laboratory under constant climatic conditions [30°C, 60% humidity following Requier et al. (2020)]. Rearing continued until at least 11 workers per queen were hatched. By the end of May (28th of May +/− 3 days), we settled 13 young colonies (queens + workers) in nest boxes and placed them in two climate regions [warm and dry in the Franconian lowlands (264 m.a.s.l.) around Würzburg (n = 5) and the cold and rather wet mountain region around Berchtesgaden, where queens originated from (n = 8); Figure 1B]. The regions differ strongly in their multi-annual mean temperature (Figure 1B) and precipitation (Supplementary material S2): The warm and dry region has an average annual temperature of 9°C and annual precipitation of about 650 mm, while the alpine region of origin has an average annual temperature of 7°C and > 1,500 mm annual precipitation. Within the warm and dry lowlands colonies were placed at one site. Within the alpine region colonies were placed at three sites [mountain foothill: 752 m.a.s.l. (n = 3); mountain mid: 1,100 m.a.s.l. (n = 3); and mountain high: 1,933 m.a.s.l. (n = 2) for temperature differences per site see Figure 1]. Because temperature increase and precipitation decrease are changing in the same direction within this approach (Supplementary material S2) we consider parallel changing conditions between regions. After at least 4 weeks in the field, we collected 3 individuals per colony (individuals leaving and entering the nest entrance). We cannot prove that the collected individuals really hatched in the field, because bumble bee workers’ lifespan averages from 22 to 69 days and depends on the species and environmental conditions (Goldblatt and Fell, 1987; Smeets and Duchateau, 2003; Kelemen et al., 2019). Thus, we here speak of individuals exposed to different conditions and not of individuals that have developed under different climatic conditions.
Collected individuals were flash frozen on dry ice in the field and stored at −20°C for further analyses.
Temperature analyses
In the study region Berchtesgaden, temperature data were modeled at the site level from daily temperature data averaged from neighboring climate stations (18 in total). The statistical details of temperature predictions are provided in Supplementary material S3. Temperature data were used to calculate short term acclimatization temperature (mean temperature one week before sampling = Tacclim) and mean annual temperature (MAT). In the warm and dry study region near Würzburg, temperature data originated from the closest German Meteorological Service (Deutscher Wetterdienst) climate station, which was at the same elevation and 1 km away from the colonies. Within the translocation experiment we calculated mean exposure temperature (i.e., mean temperature experienced by colonies from the establishment in the field to the collection of the workers = Texp) and the short-term acclimatization temperature one week before sampling (=Tacclim). Coordinates for all climate stations are given in the Supplementary material S3.
Chemical analyses
CHC profiles were extracted from legs and wings of all workers, as these body parts were shown to have the highest volume to surface ratio and best reflect climate-associated chemical traits in ants (Sprenger et al., 2021) while being representative for a species and revealing individual differences at the same time (Young et al., 2000; Wang et al., 2016; Mayr et al., 2021; Sprenger et al., 2021). Extraction was performed by immersing pooled body parts in n-hexane for 10 min per individual. The extracts were concentrated under gentle CO2 stream to approximately 20 μL and transferred to a micro insert. CHC extracts were analyzed with an Agilent 6890 gas chromatograph coupled with an Agilent 5975 Mass Selective Detector (GC–MS, Agilent, Waldbronn, Germany): The GC (split/splitless injector in splitless mode for 1 min, injected volume 1 μL at 300°C) was equipped with a DB-5 Fused Silica capillary column (30 m × 0.25 mm ID, df = 0.25 μm; J&W Scientific, Folsom, United States). Helium served as carrier gas at a constant flow of 1 mL/min. The following temperature program was used: Start temperature 60°C, temperature increase by 5°C per min up to 300°C, isotherm at 300°C for 10 min. The electron ionization mass spectra (EI-MS) were acquired at an ionization voltage of 70 eV (source temperature: 230°C). Chromatograms and mass spectra were recorded and quantified via integrated peak areas with the software HP Enhanced ChemStation G1701AA (version A.03.00; Hewlett Packard). CHC compounds were identified by the compound specific retention indices and their detected diagnostic ions (Carlson et al., 1998). Alongside CHC samples, we run an analytical alkane standard (C8–C20 and C21–C40; Sygma Aldrich) for the calculation of the retention index and to check for the sensitivity of the GCMS. Double-bond position in alkenes were identified by DMDS-derivatization following Dunkelblum et al. (1985). The GC settings were the same as described before, but oven temperature increased by 5°C per min up to 325°C and isotherm at 325°C for 10 min.
Statistical analyses
All analyses were performed using the software R version 4.1.3. (R Core Team, 2021).
Approach 1: Differences in CHC composition and climate-associated chemical traits between species occupying different elevational niches
In all six bumble bee species, we compared the relative abundances of compounds in the CHC profile of workers. Only CHC compounds which contributed to at least 0.1% of the total abundance of compounds in the CHC profile were analyzed. CHC profiles were assessed by non-metric multidimensional scaling (NMDS), a two-dimensional ordination method to visualize similarity, and agglomerative hierarchical cluster analysis (Figure 2B). We used the “plot_ordination” function from the “phyloseq” package and, respectively, the “hclust” function from the “vegan” package and “as.dendrogram” function from the “dendextend” package. Dissimilarities between worker profiles were calculated using Bray-Curtis distances. We assessed CHC profile composition differences between bumble bee species by using permutational multivariate analysis of variance (PERMANOVA) in the packages vegan (Oksanen et al., 2020) and pairwiseAdonis (Martinez Arbizu, 2020). Permutations were set on 10,000 or 999, respectively.
Overall differences in CHC profiles between species, however, do not inform per se about the protection capacity of a cuticle against desiccation. We therefore calculated more informative climate-associated chemical traits, i.e., (a) the proportion of saturated components, and (b) the abundance weighted mean chain length. Both values are predicted to increase with evaporation protection capacity (Gibbs and Pomonis, 1995). We used a generalized linear mixed model (function “glmmTMB” Brooks et al., 2017) with mean chain length and the proportion of saturated components, respectively, as response variables and species as the predictive variable. Since individuals were collected in two different years, we included year as a random effect in the model. We run a gaussian-distributed model in each case. Models were determined and checked with the DHARMa package. Model selection (null vs. full model) was performed using the dredge function from the MuMIN package (Barton, 2020) based on Akaike information criterion corrected for small sample sizes (AICc) where the lowest AICc relative to other models (here null model) indicates the preferred model (ΔAICc > 2 indicates statistically relevant differences; Burnham and Anderson, 2004). Post hoc comparison among species were performed with the emmeans function with a Tukey adjustment.
We attempted to elucidate the mechanisms that underly CHC responses to temperature, by testing and comparing CHC responses to two different temperature regimes (MAT vs. Tacclim): We assume that CHC traits that are better explained by MAT than by Tacclim, indicate species adaptation to the respective temperature niche along the elevational gradient. In contrast, we assume that CHC traits that are better explained by Tacclim, suggest a rather spontaneous and likely reversible response of CHCs to the current temperature conditions. We addressed these different temperature effects with separate generalized linear mixed models, as MAT and Tacclim were moderately correlated (corr = 0.76; see Supplementary material S4). In linear mixed models we included temperature as fixed effect and species nested in temperature as random effect (because certain temperatures exist only for certain species). We run a gaussian-distributed model in each case.
To address phylogenetic constraints within CHC profile compositions we performed a Mantel test, which compares phylogenetic distances of species with chemical distances (=calculated Bray-Curtis dissimilarities) between averaged species profiles following Buellesbach et al. (2013) and Martin et al. (2013). We obtained phylogenetic distances from the comprehensive phylogeny provided by Cameron et al. (2007). Furthermore, we analyzed phylogenetic signals (Morans’I) in the climate-associated chemical traits by phylogenetic correlograms (Gittleman and Kot, 1990), which measure the correlation between phylogenetic distances and trait distances, using the R package phylosignal (Keck et al., 2016). Higher correlation indicates that CHC variation is stronger driven by phylogeny. As an initial tree we used again the comprehensive phylogeny of Cameron et al. (2007), but this time we assigned individuals to species (Figure 2A). Hereby we defined the distance between individuals of a given species to 0.02. By doing so, we standardized the unknown genetic distances between individuals to conduct our analyses at the individual level. Our analysis should therefore only be considered as an approximation for the relatedness of individuals.
Approach 2: Intraspecific CHC variation of Bombus pascuorum along an elevational gradient
To test for intraspecific CHC differences in B. pascuorum along the elevational gradient, we first assessed profile similarity by NMDS and performed a PERMANOVA to test for an elevational effect on CHC profile composition. We then conducted a generalized linear mixed model with proportion of saturated components and mean chain length as response variables and elevation as predictor variable. We included site as random effect, because individuals sampled on the same site might per se show a higher CHC similarity than individuals sampled on different sites (e.g., due to relatedness or response to other, non-considered site factors). We run a beta-distributed model for the proportion of saturated components and a gaussian-distributed model for the mean chain length, both were determined and checked with the DHARMa package. Due to high correlation between elevation, MAT and Tacclim (Supplementary material S5) we cannot separate the mechanisms underlying temperature responses here. We thus analyzed elevational patterns, which can be interpreted as response to temperature changes. We used a likelihood ratio test (LRT) to test whether the model explains more than a random change. Test was performed with the “anova” function from the “stats” package.
Further, we assessed differences in the coefficient of variation (CV) of the two climate-associated chemical traits along the elevational gradient. The CV, which equals the standard deviation of the respective trait divided by the trait mean, is assumed to decrease with increasing elevation, as increased environmental filtering processes, especially above the tree line may constrain phenotypic variability. We tested for an elevational effect on the coefficient of variation for each site with generalized additive models (GAMs) allowing the flexible detection of linear and non-linear relationships (Wood, 2011). The basis dimension of smoothing term was set to five to have enough degrees of freedom to represent the underlying pattern well, but small enough to maintain reasonable computational efficiency. GAMs were computed with the ´gam´ function from the mgcv package (Wood, 2006).
Approach 3: CHC differences in workers of Bombus lucorum after translocation
To test for a translocation effect on CHC profile composition in B. lucorum, we assessed profile similarity using NMDS and PERMANOVA. We further, tested for differences in climate-associated chemical traits in translocated workers with generalized linear mixed models. We included colony identity as a random effect in the generalized mixed effect models. See Supplementary material S6 for NMDS displaying CHC-profile similarity from individuals out of the same colonies. We run a gaussian-distributed model in each case. Models were determined and checked with the DHARMa package. To elucidate the prevailing mechanisms behind responses, we again tested two different temperature measures against each other: the average exposed temperature throughout the study period (=Texp), whose explanatory power should be high if workers respond with a change in CHC composition to temperature exposure over several weeks vs. Tacclim, i.e., mean temperature one week before sampling, which suggest a short-time response. Due to high correlation of both Texp and Tacclim with colony location (warm and dry region; mountain foothill, mountain mid, mountain high; see Supplementary material S6) we cannot separate the mechanisms underlying temperature responses here and only present the results of the model with colony location as fixed factor.
Results
CHC difference between species occupying different elevational niches
A high proportion of CHC profile variation was explained by species (F = 56.90, df = 5, p < 0.001, R2 = 0.83; Figure 3), resulting in species-specific CHC profiles (Figure 3; Supplementary material S8, PERMANOVA between all species pairs: p = 0.001). The species-specific compositions stayed consistent within species when individuals are depicted in a cluster dendrogram (Figure 2B) with minor exceptions in B. lucorum and B. mucidus. The averaged CHC profile distances between species were not correlated with species phylogenetic distance (Mantel: R2 = 0.08, p = 0.361). See Supplementary material S9 for mean proportion of components for each species.
Figure 3. Similarity of the CHC profiles of six bumble bee species and distribution of these species along the elevational gradient. (A) Similarity of CHC profiles of bumble bees displayed in a two-dimensional graph by non-metric multidimensional scaling (NMDS) based on Bray-Curtis distances. Dots represent CHC profiles of individuals, dot colors the respective species. The closer the dots, the more similar the CHC profile. (B) Bumble bee distribution along the elevational gradient recorded within an intensive monitoring conducted between June and September 2019 (see Supplementary material S1). The green line refers to the tree line in the study region. CHC profiles are clearly determined by species however, species profile similarity does not reflect similarity in the elevational distribution.
Interspecific comparison of the mean proportions of saturated components in the CHC profiles revealed that B. pascuorum differed from the five other species, by having a lower mean proportion of saturated components (average: 45% ± 1.7%) than the other species (~56%; Figure 2C, glmm: χ2 = 36.246, df = 5, p < 0.001). Moran’s I correlation was significant at the phylogenetic distance of 0.02, representing the intraspecific level (indicated in red in Figure 2E). This suggest that proportions of saturated components are most similar among individuals within a species. The phylogenetic signal of this trait did not diverge from random expectations at the inter-species level at larger phylogenetic distances. This indicates that the proportion of saturated components show no pattern with phylogenetic distance of the considered species pool. Across all samples, neither MAT nor Tacclim affected the proportion of saturated components (MAT: χ2 = 0.372, p = 0.541; Tacclim: χ2 = 0.056, p = 0.815).
Within the CHC profiles, mean chain lengths differed between species (glmm: χ2 = 352.320, df = 5, p < 0.001). Bombus soroeensis had the longest mean chain length (mean: 27.9 ± 0.13 carbons; Figure 2D). Bombus lucorum (mean: 26.9 ± 0.1) and B. mucidus (mean: 27.1 ± 0.14) had a longer mean chain length than B. pascuorum (mean: 26.4 ± 0.07; Figure 2D). B.monticola and B. wurflenii had the shortest averaged chain length (mean: 25.6 ± 0.12 and 25.5 ± 0.12 respectively). Moran’s I correlation was significant at a phylogenetic distance of 0.02 (intraspecific level, Figure 2F), i.e., mean chain length is similar among individuals of the same species. At the interspecific level, Moran’s I correlation was negative and increased with phylogenetic distances. This indicates that within the here considered species pool, closely related species express stronger differences in mean chain length than less related ones. Across all samples, neither MAT nor Tacclim affected mean chain length (MAT: χ2 = 0.292, p = 0.589; Tacclim: χ2 = 0.187, p = 0.674).
Intraspecific CHC variation along an elevational gradient in Bombus pascuorum
In B. pascuorum (n = 58) the CHC composition varied along the elevational gradient, but the degree of variance explained by elevation was very low (PERMANOVA: F = 3.293; p = 0.015, R2 = 0.06; Figure 4A). Both climate-associated chemical traits, i.e., the proportion of saturated components (Figure 4B) and the mean chain length (Figure 4C), did not systematically change with elevation (glmm: χ2 = 1.46; p = 0.226 and χ2 = 0.39; p = 0.532). The coefficient of variation (CV) in the proportion of saturated components between individuals from the same site decreased with elevation (gam: ~elevation: edf = 1.94; DE = 85.6%; p < 0.001; Figure 4D)—a pattern driven by the highest sites above the tree line. The CV in mean chain length between individuals from the same site did not change with elevation (gam: ~ elevation: edf = 1; DE = 0.26%; p = 0.87; Figure 4E).
Figure 4. Intraspecific CHC variation of Bombus pascuorum workers along the elevational gradient. (A) Similarity of CHC profiles of individuals displayed in a two-dimensional graph by non-metric multidimensional scaling (NMDS) based on Bray-Curtis distances. Dots represent CHC profiles of individuals, dot colors the respective elevation. The closer the dots, the more similar the CHC profile. NMDS revealed a significant elevational effect on CHC profile composition but with a high degree of unexplained variance between samples. (B) The proportion of saturated components along the elevational gradient, dot colors represent mean annual temperature (=MAT), which was highly correlated with elevation (Person’s R = 0.99). (D) Coefficient of variation of the proportion of saturated components, green line represents the tree line. (C) Mean chain length along the elevational gradient, dot colors represent mean annual temperature (=MAT) and respective coefficient of variation (E).
CHC variation in translocated workers
CHC profile composition of translocated B. lucorum workers did not change with colony location (PERMANOVA: F = 1.494, R2 = 0.11 p = 0.143, lowlands: n = 16, mountain foothill: n = 9, mountain mid: n = 9, mountain high: n = 6; Figure 5A). Colony location did not affect the proportion of saturated components (glmm: df = 3; χ2 = 0.694; p = 0.875; Figure 5B) but had a marginal effect on mean chain length (glmm: df = 3; χ2 = 7.85; p = 0.049, LRT: χ2 = 6.09; p = 0.10; Figure 5C). The mean chain length was significant longer in the warm and dry region (lowlands) compared to the foothills of the cool and wet mountain region.
Figure 5. Differences in CHC profiles of Bombus lucorum translocated to four climatic different locations: warm and dry lowland region (orange), and mountain foothill (light green), mountain mid (green) and mountain high (blue; see Figure 1) in a cool and wet region. (A) Similarity of CHC profiles of individuals displayed in a two-dimensional graph by non-metric multidimensional scaling (NMDS) based on Bray-Curtis distances. Dots represent CHC profiles of individuals, dot colors the respective locations. The closer the dots, the more similar the CHC profile. Ellipse lines assuming a multivariate t-distribution, are based on the variance observed among each group of samples are presented. Non-metric multidimensional scaling revealed no differences in CHC profile compositions. (B) Proportion of saturated components. (C) Mean chain length.
Discussion
CHC profiles are species-specific
In this study, we explored the role of cuticular hydrocarbons in mediating bumble bees’ capacity to cope with different climates in mountainous regions. Bumble bee species preferring different elevational niches revealed distinct and species-specific CHC profiles. Similarities between CHC profiles were not explained by similarities in species elevational niche preference, indicating that CHCs are at least not exclusively shaped by environmental factors correlating with elevation (e.g., mean annual temperature which correlates to 98% with elevation in this study). Rather, species specificity of the profiles suggests a strong genetic component in bumble bee CHCs. Genetic heritability of CHC patterns is known from Drosophila (Ferveur and Jallon, 1996; Holze et al., 2021; Ward and Moehring, 2021) and social insect species like ants and termites (Dronnet et al., 2006; Guillem et al., 2016). Notwithstanding the apparently strong genetic component in the CHCs, the similarity of CHC profiles did not reflect phylogenetic distances among species, suggesting that CHC profiles are under diverging selection within the pool of considered species.
Higher proportion of saturated components within CHCs have been shown to increase the desiccation resistance in some insects (Menzel et al., 2018). In alpine bumble bees, five out of six investigated species did not differ in this trait, making it unlikely that the proportion of saturated components is of major importance in structuring the occurrence of bumble bee species along the elevational gradient. Only B. pascuorum, which tends to occur in the warmer foothills, expressed a lower proportion of saturated components than all other species. This contradicts the expectation that higher temperatures select for higher proportion of saturated components as improved desiccation barrier. Yet, it is partly in line with findings from other studies analyzing insect desiccation capacity along elevational gradients, which assume higher desiccation stress in higher elevations despite lower mean annual temperatures (Parkash et al., 2008; Mayr et al., 2021). Assuming that the proportion of saturated components really increases desiccation resistance in bumble bees, we suggest that low elevations may not select per se for a higher desiccation barrier, despite higher mean annual temperatures. Instead, less solar radiation, higher oxygen partial pressure, diverse microclimates and sufficient water supply (Körner, 2007) may release species like B. pascuorum from desiccation stress in the foothills. In particular, the canopy cover below the tree line can provide refugia (De Frenne et al., 2019), which allow bumble bees to actively seek places where desiccation stress is low.
Mean chain lengths differed between bumble bee species. Given that each additional carbon atom increases the melting temperature of CHC compounds by 1°C–3°C (Gibbs, 2002), differences in desiccation resistances between, e.g., B. soroeensis (on average 27.9 ± 0.13 carbons) and B. wurflenii (on average 25.5 ± 0.12 carbons) should be striking and require further examination—at best, under consideration of the absolute abundances of the different substances. Along the considered elevational gradient, however, differences in mean chain lengths did not predict the preferred elevational niche of bumble bee species. However, habitat preferences in a broader geographical context (Europe; Rasmont et al., 2015) which established 3.5 Ma ago (Hines, 2008) could explain differences in traits mediating desiccation resistance. Interestingly, species differences in the mean chain lengths were more pronounced in closely related than in less related species (Figure 2F). CHCs and mean chain lengths are known to be involved in the context of chemical communication in various insects (Nunes et al., 2009; Menzel and Schmitt, 2012; Leonhardt et al., 2016; Maihoff et al., 2022a). The negative correlation between trait similarity and phylogenetic distance suggests that mean chain length could play a role in the recognition of conspecifics, as a character displacement from closely related species might confer an evolutionary advantage in this context of communication. However, it should be noted, that the strength of a phylogenetic signal always depends on the taxonomic level studied and that the power of these analyses depends on the number of considered species (Blomberg et al., 2003; Menzel et al., 2017b), which restricts the interpretability of the phylogenetic signal in our study.
High intraspecific CHC variation and potential environmental filtering at high elevations in Bombus pascuorum
The CHC profile of B. pascuorum showed a high intraspecific variation along the 1,300 m elevational gradient. Profile variation between individuals collected at the same site was comparable to the variation between individuals collected at different elevations (Figure 4A). We also found no evidence that either of the two climate-associated chemical traits was shaped by factors linked to elevation. We are aware that bumble bee workers of B. pascuorum can perform long vertical flights to exploit resources in different elevational belts (Knight et al., 2005). Such movements could impede the detection of intraspecific acclimatization in the CHCs of this species if the bumble bees are exposed to the changing environmental conditions for a shorter period of time than a change in CHCs can happen. Although the relevant time frame for changes in CHCs in bumble bees is still unknown, laboratory studies with flies indicate, that changes even within hours are conceivable (Stinziano et al., 2015). An experimental approach under controlled climatic conditions, testing short time acclimatization within hours, however, could shed light on this. Meanwhile, our findings suggest that the chemical composition of B. pascuorums’ hydrocarbons does not seem to restrict individuals in the performance of such vertical flights in the considered elevational range.
Albeit CHC profiles and climate-associated chemical traits did not change systematically along the elevational gradient in B. pascuorum, we detected a decline in trait variance: Individuals collected in high elevations—especially above the tree line—differed less in the proportion of saturated components than individuals from low- or mid-elevations. A similar reduction of trait variation with elevation, and with the tree line as potential threshold, was detected in other functional traits, including bumble bee tongue length (Sponsler et al., 2022b) and wild bee body size (Classen et al., 2017). Pronounced environmental filtering at habitats with less variable microclimatic conditions (Hoiss et al., 2012) might cause such declines. Here, the phenotype that withstands the potential filtering process, shows a comparable low proportion of saturated components, i.e., on the intraspecific level the desiccation barrier seems rather reduced in individuals that occurred in high elevations. The environmental filtering suggests that the maintenance of phenotypes expressing fewer saturated components leads to an adaptive advantage, even though the mechanisms behind remain elusive. Genetic bottleneck effects at higher elevations (i.e., fewer colonies at high elevation equals less genetic variation) may further constrain trait variations in the highlands (Jackson et al., 2018).
CHC variation in translocated workers of Bombus lucorum
We found no effect of colony location on the CHC profile patterns of translocated B. lucorum workers. The translocation experiment showed that a 4-weeks exposure to different climatic regions, reflecting a 3-fold difference in annual precipitation (Supplementary material S2) and a 25% difference in mean temperature (Figure 1B), did not trigger a systematic plasticity in the overall CHC composition as a response of warmer and dryer nest location. This is surprising, as precipitation/humidity and temperature differences can affect CHC composition in insects (e.g., ants) earlier—even within weeks—in both field and laboratory experiments (Buellesbach et al., 2018; Sprenger et al., 2018).
Colony location influenced chain length in translocated workers: workers translocated to warm and dry lowlands tended to have longer chains than those from the cool and wet mountainous foothills. This can be interpreted as a response to differences in regional climate, as longer chains are needed in warmer/drier conditions to prevent water loss through desiccation (Rajpurohit et al., 2017; Menzel et al., 2017a). Similarly, ants subjected to laboratory experiments increased their average chain length in warm conditions (Menzel et al., 2018), and so did flies (Rajpurohit et al., 2017). However, in our study, mean chain lengths was not found to be shorter in even cooler habitats along the elevational gradient, suggesting again that other factors than temperature predict this trait. Further, chain length increased on average by only 0.3carbons, challenging the question about the physiological relevance of this effect. While chain length can also be influenced by age (Jackson and Bartelt, 1986; Nunes et al., 2009) we found this to be an unlikely effect here, as the colonies for the translocation experiment were reared at the same time, and thus workers were of comparable age. Nevertheless, including life history conditions and other thermal performance related functional traits, e.g., body size, would be informative in future studies as desiccation resistance is a complex trait, likely driven by multiple factors (Davis and Moyle, 2019). Other traits that may influence climate tolerance and hence the distribution pattern of bumble bee species along elevational gradients (Martinet et al., 2021) are genetic regulatory mechanisms like, e.g., the upregulation of heat shock proteins (Jackson et al., 2018; Pimsler et al., 2020) and pile color (Williams, 2007). Bumble bees—unlike ants—are characterized by exhaustive hairiness. The dense cover of hair might add insulation that reduces water loss in general (Kühsel et al., 2017). Thus, bumble bees might be less susceptible to desiccation stress (Parsons, 2019), which might reduce the pressure to alter CHC profile in response to desiccation.
Conclusion
In this study we highlight the high variability of cuticular chemical phenotypes in six species in the genus Bombus. Our results suggest that the proportion of saturated components and mean chain length are not decisive traits explaining bumble bees’ distribution along elevation in our study area. Profiles and chain lengths were rather species-specific, indicating that species recognition or long-established evolutionary divergence may limit trait variation. Yet, differences in profiles and climate-associated chemical traits may explain differences in species response to climate change whereas plastic adjustments within species seem to be unlikely. At present, species range shifts to higher elevations do not appear to be hindered by any particular CHC phenotype—however, a link to species fitness is urgently needed to clarify this claim. With increasing temperatures or under extreme events, the differences in CHC composition of species might become relevant and determine species’ resistance toward warming. Extending this study to more species, including those that already inhabit much warmer habitats, and under the consideration of multiple traits and related fitness consequences may shed light on bumble bee responses to climate change. Our study provides first insights that conserved species-specific signals are maintained under different environmental conditions. This is an important step toward deciphering how potentially conflicting functions in insect CHC profiles can be unified.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.
Author contributions
FM, AC, SN, and TS contributed to conception and design of the study. FM collected the specimens with help of KB, LB, SSc, and AC, performed the statistical analysis with input from SN and AC, and wrote the first draft of the manuscript. FM, SSa, SSc, KK, NS, and AC reared and maintained the bumble bee colonies. All authors contributed to the article and approved the submitted version.
Funding
The GC/MS was granted by the DFG to TS. This project was sponsored by the Bavarian State Ministry of Science and the Arts in the context of the Bavarian Climate Research Network (bayklif).
Acknowledgments
We thank the National Park Berchtesgaden for logistic support and for sharing climatic data, in consultation with the avalanche warning center Bavaria. We also thank the owners of the study sites outside and inside the National Park for their permission and support. Special thanks to Jie Zhang for providing the Bavaria map with temperature and precipitation data and to Julia Zetzsche for the beautiful bumble bee illustrations. Mountain illustrations were obtained from vecteezy.com and klebefux.de and modified.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fevo.2023.1082559/full#supplementary-material
References
Aguirre-Gutiérrez, J., Kissling, W. D., Carvalheiro, L. G., WallisDeVries, M. F., Franzén, M., and Biesmeijer, J. C. (2016). Functional traits help to explain half-century long shifts in pollinator distributions. Sci. Rep. 6:24451. doi: 10.1038/srep24451
Akino, T., Yamamura, K., Wakamura, S., and Yamaoka, R. (2004). Direct behavioral evidence for hydrocarbons as nestmate recognition cues in Formica japonica (hymenoptera: Formicidae). Appl. Entomol. Zool. 39, 381–387. doi: 10.1303/aez.2004.381
Arbetman, M. P., Gleiser, G., Morales, C. L., Williams, P., and Aizen, M. A. (2017). Global decline of bumblebees is phylogenetically structured and inversely related to species range size and pathogen incidence. Proc. R. Soc. B Biol. Sci. 284:20170204. doi: 10.1098/rspb.2017.0204
Awater-Salendo, S., Schulz, H., Hilker, M., and Fürstenau, B. (2020). The importance of methyl-branched Cuticular hydrocarbons for successful host recognition by the larval Ectoparasitoid Holepyris sylvanidis. J. Chem. Ecol. 46, 1032–1046. doi: 10.1007/s10886-020-01227-w
Barton, K. (2020). MuMIn: Multi-Model Inference. Available at: https://cran.r-project.org/web/packages/MuMIn/MuMIn.pdf (Accessed 15 December 2021).
Beals, E. W. (1969). Vegetational change along altitudinal gradients. Science 165, 981–985. doi: 10.1126/science.165.3897.981
Beani, L., Bagnères, A.-G., Elia, M., Petrocelli, I., Cappa, F., and Lorenzi, M. C. (2019). Cuticular hydrocarbons as cues of sex and health condition in Polistes dominula wasps. Insect. Soc. 66, 543–553. doi: 10.1007/s00040-019-00721-z
Bellard, C., Bertelsmeier, C., Leadley, P., Thuiller, W., and Courchamp, F. (2012). Impacts of climate change on the future of biodiversity. Ecol. Lett. 15, 365–377. doi: 10.1111/j.1461-0248.2011.01736.x
Bingham, R. A., and Orthner, A. R. (1998). Efficient pollination of alpine plants. Nature 391, 238–239. doi: 10.1038/34564
Blomberg, S. P., Garland, J. R. T., and Ives, A. R. (2003). Testing for phylogenetic signal in comparative data: behavioral traits are more labile. Evolution 57, 717–745. doi: 10.1111/j.0014-3820.2003.tb00285.x
Blomquist, G. J., and Bagnères, A.-G. (2010). Insect hydrocarbons: Biology, biochemistry, and chemical ecology. Cambridge. New York: Cambridge University Press.
Brooks, M. E., Kristensen, K., Van Benthem, K. J., Magnusson, A., Berg, C. W., Nielsen, A., et al. (2017). glmmTMB balances speed and flexibility among packages for zero-inflated generalized linear mixed modeling. R J. 9:378. doi: 10.32614/RJ-2017-066
Buellesbach, J., Gadau, J., Beukeboom, L. W., Echinger, F., Raychoudhury, R., Werren, J. H., et al. (2013). Cuticular hydrocarbon divergence in the jewel wasp Nasonia: evolutionary shifts in chemical communication channels? J. Evol. Biol. 26, 2467–2478. doi: 10.1111/jeb.12242
Buellesbach, J., Whyte, B. A., Cash, E., Gibson, J. D., Scheckel, K. J., Sandidge, R., et al. (2018). Desiccation resistance and micro-climate adaptation: Cuticular hydrocarbon signatures of different argentine ant Supercolonies across California. J. Chem. Ecol. 44, 1101–1114. doi: 10.1007/s10886-018-1029-y
Burnham, K. P., and Anderson, D. R. (2004). Multimodel inference: understanding AIC and BIC in model selection. Sociol. Methods Res. 33, 261–304. doi: 10.1177/0049124104268644
Cameron, S. A., Hines, H. M., and Williams, P. H. (2007). A comprehensive phylogeny of the bumble bees (Bombus). Biol. J. Linn. Soc. 91, 161–188. doi: 10.1111/j.1095-8312.2007.00784.x
Carlson, D. A., Bernier, U. R., and Sutton, B. D. (1998). Elution patterns from capillary GC for methyl-branched alkanes. J. Chem. Ecol. 24, 1845–1865. doi: 10.1023/A:1022311701355
Chown, S. L., Sørensen, J. G., and Terblanche, J. S. (2011). Water loss in insects: an environmental change perspective. J. Insect Physiol. 57, 1070–1084. doi: 10.1016/j.jinsphys.2011.05.004
Classen, A., Steffan-Dewenter, I., Kindeketa, W. J., and Peters, M. K. (2017). Integrating intraspecific variation in community ecology unifies theories on body size shifts along climatic gradients. Funct. Ecol. 31, 768–777. doi: 10.1111/1365-2435.12786
Colazza, S., Aquila, G., De Pasquale, C., Peri, E., and Millar, J. G. (2007). The egg parasitoid Trissolcus basalis uses n-nonadecane, a Cuticular hydrocarbon from its stink bug host Nezara viridula, to discriminate between female and male hosts. J. Chem. Ecol. 33, 1405–1420. doi: 10.1007/s10886-007-9300-7
Dani, F. R., Jones, G. R., Corsi, S., Beard, R., Pradella, D., and Turillazzi, S. (2005). Nestmate recognition cues in the honey bee: differential importance of Cuticular alkanes and alkenes. Chem. Senses 30, 477–489. doi: 10.1093/chemse/bji040
Davis, J. S., and Moyle, L. C. (2019). Desiccation resistance and pigmentation variation reflects bioclimatic differences in the Drosophila Americana species complex. BMC Evol. Biol. 19:204. doi: 10.1186/s12862-019-1536-7
De Frenne, P., Zellweger, F., Rodríguez-Sánchez, F., Scheffers, B. R., Hylander, K., Luoto, M., et al. (2019). Global buffering of temperatures under forest canopies. Nat. Ecol. Evol. 3, 744–749. doi: 10.1038/s41559-019-0842-1
De Oliveira, C. C., Manfrin, M. H., Sene, F. D. M., Jackson, L. L., and Etges, W. J. (2011). Variations on a theme: diversification of cuticular hydrocarbons in a clade of cactophilic drosophila. BMC Evol. Biol. 11:179. doi: 10.1186/1471-2148-11-179
Dronnet, S., Lohou, C., Christides, J.-P., and Bagnères, A.-G. (2006). Cuticular hydrocarbon composition reflects genetic relationship among colonies of the introduced termite Reticulitermes santonensis Feytaud. J. Chem. Ecol. 32, 1027–1042. doi: 10.1007/s10886-006-9043-x
Dunkelblum, E., Tan, S. H., and Silk, P. J. (1985). Double-bond location in monounsaturated fatty acids by dimethyl disulfide derivatization and mass spectrometry: application to analysis of fatty acids in pheromone glands of four lepidoptera. J. Chem. Ecol. 11, 265–277. doi: 10.1007/BF01411414
Ferveur, J.-F., and Jallon, J.-M. (1996). Genetic control of male cuticular hydrocarbons in Drosophila melanogaster. Genet. Res. 67, 211–218. doi: 10.1017/S0016672300033693
Flynn, D. F. B., Mirotchnick, N., Jain, M., Palmer, M. I., and Naeem, S. (2011). Functional and phylogenetic diversity as predictors of biodiversity–ecosystem-function relationships. Ecology 92, 1573–1581. doi: 10.1890/10-1245.1
Geue, J. C., and Thomassen, H. A. (2020). Unraveling the habitat preferences of two closely related bumble bee species in Eastern Europe. Ecol. Evol. 10, 4773–4790. doi: 10.1002/ece3.6232
Gibbs, A. G. (1998). Water-proofing properties of Cuticular lipids. Am. Zool. 38, 471–482. doi: 10.1093/icb/38.3.471
Gibbs, A. G. (2002). Lipid melting and cuticular permeability: new insights into an old problem. J. Insect Physiol. 48, 391–400. doi: 10.1016/S0022-1910(02)00059-8
Gibbs, A., Mousseau, T. A., and Crowe, J. H. (1991). Genetic and acclimatory variation in biophysical properties of insect cuticle lipids. Proc. Natl. Acad. Sci. 88, 7257–7260. doi: 10.1073/pnas.88.16.7257
Gibbs, A., and Pomonis, J. G. (1995). Physical properties of insect cuticular hydrocarbons: the effects of chain length, methyl-branching and unsaturation. Comp. Biochem. Physiol. B Biochem. Mol. Biol. 112, 243–249. doi: 10.1016/0305-0491(95)00081-X
Gittleman, J. L., and Kot, M. (1990). Adaptation: statistics and a null model for estimating phylogenetic effects. Syst. Biol. 39, 227–241. doi: 10.2307/2992183
Goldblatt, J. W., and Fell, R. D. (1987). Adult longevity of workers of the bumble bees Bombus fervidus (F.) and Bombus pennsylvanicus (De Geer) (hymenoptera: Apidae). Can. J. Zool. 65, 2349–2353. doi: 10.1139/z87-354
Gorenflo, A., Diekötter, T., Van Kleunen, M., Wolters, V., and Jauker, F. (2017). Contrasting pollination efficiency and effectiveness among flower visitors of Malva sylvestris, Borago officinalis and Onobrychis viciifolia. J. Pollinat. Ecol. 21, 62–70. doi: 10.26786/1920-7603(2017)421
Guillem, R. M., Drijfhout, F. P., and Martin, S. J. (2016). Species-specific Cuticular hydrocarbon stability within European Myrmica ants. J. Chem. Ecol. 42, 1052–1062. doi: 10.1007/s10886-016-0784-x
Halsch, C. A., Shapiro, A. M., Fordyce, J. A., Nice, C. C., Thorne, J. H., Waetjen, D. P., et al. (2021). Insects and recent climate change. Proc. Natl. Acad. Sci. 118:e2002543117. doi: 10.1073/pnas.2002543117
Hines, H. M. (2008). Historical biogeography, divergence times, and diversification patterns of bumble bees (hymenoptera: Apidae: Bombus). Syst. Biol. 57, 58–75. doi: 10.1080/10635150801898912
Hoiss, B., Krauss, J., Potts, S. G., Roberts, S., and Steffan-Dewenter, I. (2012). Altitude acts as an environmental filter on phylogenetic composition, traits and diversity in bee communities. Proc. R. Soc. B Biol. Sci. 279, 4447–4456. doi: 10.1098/rspb.2012.1581
Holze, H., Schrader, L., and Buellesbach, J. (2021). Advances in deciphering the genetic basis of insect cuticular hydrocarbon biosynthesis and variation. Heredity 126, 219–234. doi: 10.1038/s41437-020-00380-y
IPCC (2018). “Global Warming of 1.5°C,” in An IPCC Special Report on the Impacts of Global Warming of 1.5°C Above Pre-Industrial Levels and Related Global Greenhouse Gas Emission Pathways, in the Context of Strengthening the Global Response to the Threat of Climate Change, Sustainable Development, and Efforts to Eradicate Poverty. eds. V. Masson-Delmotte, P. Zhai, H.-O. Pörtner, D. Roberts, J. Skea, and P. R. Shukla, et al. (Cambridge, UK and New York, NY, USA: Cambridge University Press), 616.
Jackson, L. L., and Bartelt, R. J. (1986). Cuticular hydrocarbons of Drosophila virilis: comparison by age and sex. Insect Biochem. 16, 433–439. doi: 10.1016/0020-1790(86)90056-9
Jackson, J. M., Pimsler, M. L., Oyen, K. J., Koch-Uhuad, J. B., Herndon, J. D., Strange, J. P., et al. (2018). Distance, elevation and environment as drivers of diversity and divergence in bumble bees across latitude and altitude. Mol. Ecol. 27, 2926–2942. doi: 10.1111/mec.14735
Jackson, J. M., Pimsler, M. L., Oyen, K. J., Strange, J. P., Dillon, M. E., and Lozier, J. D. (2020). Local adaptation across a complex bioclimatic landscape in two montane bumble bee species. Mol. Ecol. 29, 920–939. doi: 10.1111/mec.15376
Kather, R., and Martin, S. J. (2015). Evolution of Cuticular hydrocarbons in the hymenoptera: a meta-analysis. J. Chem. Ecol. 41, 871–883. doi: 10.1007/s10886-015-0631-5
Keck, F., Rimet, F., Bouchez, A., and Franc, A. (2016). Phylosignal: an R package to measure, test, and explore the phylogenetic signal. Ecol. Evol. 6, 2774–2780. doi: 10.1002/ece3.2051
Kelemen, E. P., Cao, N., Cao, T., Davidowitz, G., and Dornhaus, A. (2019). Metabolic rate predicts the lifespan of workers in the bumble bee Bombus impatiens. Apidologie 50, 195–203. doi: 10.1007/s13592-018-0630-y
Kellermann, V., Loeschcke, V., Hoffmann, A. A., Kristensen, T. N., Fløjgaard, C., David, J. R., et al. (2012). Phylogenetic constraints in key functional traits behind species’ climate niches: patterns of desiccation and cold resistance across 95 drosophila species. Evolution 66, 3377–3389. doi: 10.1111/j.1558-5646.2012.01685.x
Kerr, J. T., Pindar, A., Galpern, P., Packer, L., Potts, S. G., Roberts, S. M., et al. (2015). Climate change impacts on bumblebees converge across continents. Science 349, 177–180. doi: 10.1126/science.aaa7031
Knight, M. E., Martin, A. P., Bishop, S., Osborne, J. L., Hale, R. J., Sanderson, R. A., et al. (2005). An interspecific comparison of foraging range and nest density of four bumblebee (Bombus) species. Mol. Ecol. 14, 1811–1820. doi: 10.1111/j.1365-294X.2005.02540.x
Konnert, V., and Siegrist, J. (2000). Waldentwicklung im Nationalpark Berchtesgaden von 1983 bis 1997: gemeinsame Auswertung der 1. und 2. permanenten Stichproben-Inventur: natürliche Wiederbewaldung eines ungeräumten Fichten-Windwurfs auf potentiellemBergmischwaldstandort. Berchtesgaden: Nationalparkverwaltung Berchtesgaden.
Körner, C. (2007). The use of ‘altitude’ in ecological research. Trends Ecol. Evol. 22, 569–574. doi: 10.1016/j.tree.2007.09.006
Köstler, J. N., and Mayer, H. (1970). “Waldgrenzen Im Berchtesgadener Land.,” in Jubiläumsjahrbuch 1900–1970. Available at: https://www.zobodat.at/pdf/Jb-Verein-Schutz-Alpenpfl-Tiere_35_1970_0121-0153.pdf (Accessed 20 April 2022).
Kühsel, S., Brückner, A., Schmelzle, S., Heethoff, M., and Blüthgen, N. (2017). Surface area–volume ratios in insects. Insect Sci. 24, 829–841. doi: 10.1111/1744-7917.12362
Lacey, E. S., Ginzel, M. D., Millar, J. G., and Hanks, L. M. (2008). 7-methylheptacosane is a major component of the contact sex pheromone of the cerambycid beetle Neoclytus acuminatus acuminatus. Physiol. Entomol. 33, 209–216. doi: 10.1111/j.1365-3032.2008.00624.x
Leonhardt, S. D., Menzel, F., Nehring, V., and Schmitt, T. (2016). Ecology and evolution of communication in social insects. Cells 164, 1277–1287. doi: 10.1016/j.cell.2016.01.035
Maebe, K., Hart, A. F., Marshall, L., Vandamme, P., Vereecken, N. J., Michez, D., et al. (2021). Bumblebee resilience to climate change, through plastic and adaptive responses. Glob. Chang. Biol. 27, 4223–4237. doi: 10.1111/gcb.15751
Maihoff, F., Bohlke, K., Brockmann, A., and Schmitt, T. (2022a). Increased complexity of worker CHC profiles in Apis dorsata correlates with nesting ecology. PLoS One 17:e0271745. doi: 10.1371/journal.pone.0271745
Maihoff, F., Friess, N., Hoiss, B., Schmid-Egger, C., Kerner, J., Neumayer, J., et al. (2022b). Smaller, more diverse and on the way to the top: rapid community shifts of montane wild bees within an extraordinary hot decade. Divers. Distrib. 29, 272–288. doi: 10.1111/ddi.13658
Marshall, L., Perdijk, F., Dendoncker, N., Kunin, W., Roberts, S., and Biesmeijer, J. C. (2020). Bumblebees moving up: shifts in elevation ranges in the Pyrenees over 115 years. Proc. R. Soc. B Biol. Sci. 287:20202201. doi: 10.1098/rspb.2020.2201
Martin, S. J., Vitikainen, E., Shemilt, S., Drijfhout, F. P., and Sundström, L. (2013). Sources of variation in Cuticular hydrocarbons in the ant Formica exsecta. J. Chem. Ecol. 39, 1415–1423. doi: 10.1007/s10886-013-0366-0
Martinet, B., Dellicour, S., Ghisbain, G., Przybyla, K., Zambra, E., Lecocq, T., et al. (2021). Global effects of extreme temperatures on wild bumblebees. Conserv. Biol. 35, 1507–1518. doi: 10.1111/cobi.13685
Martinez Arbizu, P. (2020). pairwiseAdonis: Pairwise multilevel comparison using Adonis. R package. Available at: https://github.com/pmartinezarbizu/pairwiseAdonis.
Mayr, A. V., Keller, A., Peters, M. K., Grimmer, G., Krischke, B., Geyer, M., et al. (2021). Cryptic species and hidden ecological interactions of halictine bees along an elevational gradient. Ecol. Evol. 11, 7700–7712. doi: 10.1002/ece3.7605
Menzel, F., Blaimer, B. B., and Schmitt, T. (2017a). How do cuticular hydrocarbons evolve? Physiological constraints and climatic and biotic selection pressures act on a complex functional trait. Proc. R. Soc. B Biol. Sci. 284:20161727. doi: 10.1098/rspb.2016.1727
Menzel, F., and Schmitt, T. (2012). Tolerance requires the right smell: first evidence for interspecific selection on chemical recognition cues. Evolution 66, 896–904. doi: 10.1111/j.1558-5646.2011.01489.x
Menzel, F., Schmitt, T., and Blaimer, B. B. (2017b). The evolution of a complex trait: cuticular hydrocarbons in ants evolve independent from phylogenetic constraints. J. Evol. Biol. 30, 1372–1385. doi: 10.1111/jeb.13115
Menzel, F., Zumbusch, M., and Feldmeyer, B. (2018). How ants acclimate: impact of climatic conditions on the cuticular hydrocarbon profile. Funct. Ecol. 32, 657–666. doi: 10.1111/1365-2435.13008
Michelutti, K. B., Soares, E. R. P., Sguarizi-Antonio, D., Piva, R. C., Súarez, Y. R., Cardoso, C. A. L., et al. (2018). Influence of temperature on survival and cuticular chemical profile of social wasps. J. Therm. Biol. 71, 221–231. doi: 10.1016/j.jtherbio.2017.11.019
Minachilis, K., Kantsa, A., Devalez, J., Trigas, P., Tscheulin, T., and Petanidou, T. (2020). Bumblebee diversity and pollination networks along the elevation gradient of Mount Olympus. Greece Divers. Distrib. 26, 1566–1581. doi: 10.1111/ddi.13138
Nooten, S. S., and Hughes, L. (2017). The power of the transplant: direct assessment of climate change impacts. Clim. Change 144, 237–255. doi: 10.1007/s10584-017-2037-6
Nunes, T. M., Nascimento, F. S., Turatti, I. C., Lopes, N. P., and Zucchi, R. (2008). Nestmate recognition in a stingless bee: does the similarity of chemical cues determine guard acceptance? Anim. Behav. 75, 1165–1171. doi: 10.1016/j.anbehav.2007.08.028
Nunes, T. M., Turatti, I. C. C., Lopes, N. P., and Zucchi, R. (2009). Chemical signals in the stingless bee, Frieseomelitta varia, indicate caste, gender, age, and reproductive status. J. Chem. Ecol. 35, 1172–1180. doi: 10.1007/s10886-009-9691-8
Oksanen, J., Blanchet, G. F., Friendly, M., Kindt, R., Legendre, P., McGlinn, D., et al. (2020). Vegan: Community ecology package. Available at: https://CRAN.R-project.org/package=vegan.
Parkash, R., Rajpurohit, S., and Ramniwas, S. (2008). Changes in body melanisation and desiccation resistance in highland vs. lowland populations of D. melanogaster. J. Insect Physiol. 54, 1050–1056. doi: 10.1016/j.jinsphys.2008.04.008
Parsons, Z. (2019). Influence of pile length and body size on rates of heat loss in the bumble bee, Bombus vosnesenskii, The. Environ. Sci. doi: 10.15786/13700854.v3
Pepin, N., Bradley, R. S., Diaz, H. F., Baraer, M., Caceres, E. B., Forsythe, N., et al. (2015). Elevation-dependent warming in mountain regions of the world. Nat. Clim. Chang. 5, 424–430. doi: 10.1038/nclimate2563
Peters, M. K., Peisker, J., Steffan-Dewenter, I., and Hoiss, B. (2016). Morphological traits are linked to the cold performance and distribution of bees along elevational gradients. J. Biogeogr. 43, 2040–2049. doi: 10.1111/jbi.12768
Pimsler, M. L., Oyen, K. J., Herndon, J. D., Jackson, J. M., Strange, J. P., Dillon, M. E., et al. (2020). Biogeographic parallels in thermal tolerance and gene expression variation under temperature stress in a widespread bumble bee. Sci. Rep. 10:17063. doi: 10.1038/s41598-020-73391-8
Potts, S. G., Biesmeijer, J. C., Kremen, C., Neumann, P., Schweiger, O., and Kunin, W. E. (2010). Global pollinator declines: trends, impacts and drivers. Trends Ecol. Evol. 25, 345–353. doi: 10.1016/j.tree.2010.01.007
Pyke, G. H., Thomson, J. D., Inouye, D. W., and Miller, T. J. (2016). Effects of climate change on phenologies and distributions of bumble bees and the plants they visit. Ecosphere 7:e01267. doi: 10.1002/ecs2.1267
R Core Team. (2021). R: A language and environment for statistical computing. Vienne, Austria Available at: https://www.R-project.org/.
Rahbek, C., Borregaard, M. K., Colwell, R. K., Dalsgaard, B., Holt, B. G., Morueta-Holme, N., et al. (2019). Humboldt’s enigma: what causes global patterns of mountain biodiversity? Science 365, 1108–1113. doi: 10.1126/science.aax0149
Rajpurohit, S., Hanus, R., Vrkoslav, V., Behrman, E. L., Bergland, A. O., Petrov, D., et al. (2017). Adaptive dynamics of cuticular hydrocarbons in drosophila. J. Evol. Biol. 30, 66–80. doi: 10.1111/jeb.12988
Rasmont, P., Franzen, M., Lecocq, T., Harpke, A., Roberts, S., Biesmeijer, K., et al. (2015). Climatic risk and distribution atlas of European bumblebees. BioRisk 10, 1–236. doi: 10.3897/biorisk.10.4749
Requier, F., Jowanowitsch, K. K., Kallnik, K., and Steffan-Dewenter, I. (2020). Limitation of complementary resources affects colony growth, foraging behavior, and reproduction in bumble bees. Ecology 101:e02946. doi: 10.1002/ecy.2946
Slatyer, R. O., and Noble, I. R. (1992). Dynamics of montane Treelines. Landsc. Boundaries Consequences Biot. Divers. Ecol. Flows 92, 346–359. doi: 10.1007/978-1-4612-2804-2_17
Smeets, P., and Duchateau, M. J. (2003). Longevity of Bombus terrestris workers (hymenoptera: Apidae) in relation to pollen availability, in the absence of foraging. Apidologie 34, 333–337. doi: 10.1051/apido:2003026
Sponsler, D. B., Kallnik, K., Requier, F., Classen, A., Maihoff, A. F., Sieger, J., et al. (2022b). Floral preferences of mountain bumble bees are constrained by functional traits but flexible through elevation and season. Oikos 2022:e08902. doi: 10.1111/oik.08902
Sponsler, D. B., Requier, F., Kallnik, K., Classen, A., Maihoff, F., Sieger, J., et al. (2022a). Contrasting patterns of richness, abundance, and turnover in mountain bumble bees and their floral hosts. Ecology 103:e3712. doi: 10.1002/ecy.3712
Sprenger, P. P., Burkert, L. H., Abou, B., Federle, W., and Menzel, F. (2018). Coping with the climate: cuticular hydrocarbon acclimation of ants under constant and fluctuating conditions. J. Exp. Biol. 221:jeb171488. doi: 10.1242/jeb.171488
Sprenger, P. P., Gerbes, L. J., Sahm, J., and Menzel, F. (2021). Cuticular hydrocarbon profiles differ between ant body parts: implications for communication and our understanding of CHC diffusion. Curr. Zool. 67, 531–540. doi: 10.1093/cz/zoab012
Sprenger, P. P., Hartke, J., Feldmeyer, B., Orivel, J., Schmitt, T., and Menzel, F. (2019). Influence of mutualistic lifestyle, mutualistic partner, and climate on Cuticular hydrocarbon profiles in Parabiotic ants. J. Chem. Ecol. 45, 741–754. doi: 10.1007/s10886-019-01099-9
Stemkovski, M., Pearse, W. D., Griffin, S. R., Pardee, G. L., Gibbs, J., Griswold, T., et al. (2020). Bee phenology is predicted by climatic variation and functional traits. Ecol. Lett. 23, 1589–1598. doi: 10.1111/ele.13583
Stinziano, J. R., Sové, R. J., Rundle, H. D., and Sinclair, B. J. (2015). Rapid desiccation hardening changes the cuticular hydrocarbon profile of Drosophila melanogaster. Comp. Biochem. Physiol. A Mol. Integr. Physiol. 180, 38–42. doi: 10.1016/j.cbpa.2014.11.004
Straub, F., Kuppler, J., Fellendorf, M., Teuscher, M., Vogt, J., and Ayasse, M. (2022). Land-use stress alters cuticular chemical surface profile and morphology in the bumble bee Bombus lapidarius. PLoS One 17:e0268474. doi: 10.1371/journal.pone.0268474
Trew, B. T., and Maclean, I. M. D. (2021). Vulnerability of global biodiversity hotspots to climate change. Glob. Ecol. Biogeogr. 30, 768–783. doi: 10.1111/geb.13272
Wang, Q., Goodger, J. Q. D., Woodrow, I. E., and Elgar, M. A. (2016). Location-specific cuticular hydrocarbon signals in a social insect. Proc. Biol. Sci. 283:20160310. doi: 10.1098/rspb.2016.0310
Ward, H. K. E., and Moehring, A. J. (2021). Genes underlying species differences in cuticular hydrocarbon production between Drosophila melanogaster and D. simulans. Genome 64, 87–95. doi: 10.1139/gen-2019-0224
Williams, P. (2007). The distribution of bumblebee colour patterns worldwide: possible significance for thermoregulation, crypsis, and warning mimicry: BUMBLEBEE COLOUR-PATTERN GROUPS. Biol. J. Linn. Soc. 92, 97–118. doi: 10.1111/j.1095-8312.2007.00878.x
Williams, N. M., Crone, E. E., Roulston, T. H., Minckley, R. L., Packer, L., and Potts, S. G. (2010). Ecological and life-history traits predict bee species responses to environmental disturbances. Biol. Conserv. 143, 2280–2291. doi: 10.1016/j.biocon.2010.03.024
Wong, M. K. L., Guénard, B., and Lewis, O. T. (2019). Trait-based ecology of terrestrial arthropods. Biol. Rev. 94, 999–1022. doi: 10.1111/brv.12488
Wood, S. N. (2006). Generalized additive models: An introduction with R. New York: Chapman and Hall/CRC.
Wood, S. N. (2011). Fast stable restricted maximum likelihood and marginal likelihood estimation of semiparametric generalized linear models. J. R. Stat. Soc. Ser. B Stat. Methodol. 73, 3–36. doi: 10.1111/j.1467-9868.2010.00749.x
Keywords: pollinators, altitudinal gradient, cuticular hydrocarbon, desiccation, mountain, global change, translocation experiment, drought stress
Citation: Maihoff F, Sahler S, Schoger S, Brenzinger K, Kallnik K, Sauer N, Bofinger L, Schmitt T, Nooten SS and Classen A (2023) Cuticular hydrocarbons of alpine bumble bees (Hymenoptera: Bombus) are species-specific, but show little evidence of elevation-related climate adaptation. Front. Ecol. Evol. 11:1082559. doi: 10.3389/fevo.2023.1082559
Edited by:
Cesar Rodriguez-Saona, Rutgers, The State University of New Jersey, United StatesReviewed by:
Corinne Hertaeg, ETH Zürich, SwitzerlandFederico Cappa, University of Florence, Italy
Todd Johnson, Louisiana State University, United States
Copyright © 2023 Maihoff, Sahler, Schoger, Brenzinger, Kallnik, Sauer, Bofinger, Schmitt, Nooten and Classen. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Fabienne Maihoff, 
YW5uZV9mYWJpZW5uZS5tYWlob2ZmQHVuaS13dWVyemJ1cmcuZGU=
†These authors share last authorship