Data and code from: Climatic origin and plasticity shape emergence timing and fitness in bees and wasps under experimental climate regimes
Data files
Mar 12, 2026 version files 22.15 MB
-
Detailed_temperature_experiment.csv
707.62 KB
-
Emergence_mass.csv
1.11 MB
-
Emergence_time.csv
4.68 MB
-
Emergence-mass_script.html
3.04 MB
-
Emergence-time_script.html
12.55 MB
-
README.md
8 KB
-
Temp_comparison.csv
51.70 KB
Abstract
Climate warming threatens the persistence of insect populations by forcing them to adjust their phenology—responses that may be constrained by local adaptation and involve physiological trade-offs. We investigated how hymenopterans with contrasting seasonal phenology and from different climates adjust their emergence timing under current and projected future climates. We also examined the consequences of these adjustments on body mass at emergence as a key fitness trait. We analysed the emergence timing of 14,921 individuals from five cavity-nesting bee and wasp species exposed to cold, warm, and hot post-winter temperature treatments in a common garden experiment to identify potential genotype-environment interactions. Insects had developed at 161 sites of origin across southern Germany, covering multi-annual mean temperatures (MAT) of 5.9 to 10 ºC and differing in the temperature deviation (ΔT) that insects experienced during development in the pre-emergence year relative to MAT. Emergence timing was highly plastic to post-winter temperatures, with insects emerging earlier in warmer treatments. However, emergence was modulated by MAT and ΔT, suggesting genetic adaptation to long-term climatic conditions and adjustments to short-term temperature deviations. For spring-emerging species, individuals from sites with higher MAT and warmer treatments emerged the earliest (cogradient variation). In contrast, the latest summer-emerging species exhibited countergradient variation: in the cold treatment, individuals from higher MAT emerged later than those from lower MAT. In spring species, mass loss was higher in warmer post-winter treatments, with the strongest reductions observed in cool-adapted individuals. Mass loss was particularly rapid for summer females in warmer treatments, with individuals emerging later losing up to 34 % of their mass. However, body mass of summer insects was independent of MAT and ΔT. Our results demonstrate high plasticity of cavity-nesting Hymenoptera to post-winter temperatures but also suggest that local adaptation and responses to early-life temperature can compromise fitness under rapid climate changes. This large-scale experimental study highlights the complex drivers of insect emergence phenology and fitness and suggests that cool-adapted, spring-emerging species may be most vulnerable to ongoing climate warming.
Dataset DOI: 10.5061/dryad.3ffbg7b0b
Description of the data and file structure
This folder contains four data frames and two Markdown files with the code to obtain the main results included in the manuscript.
Files and variables
File: Detailed_temperature_experiment.csv
Description: data to produce Fig. S1.2 (code in Emergence-time_script)
Variables
- Date: date since insects were brought from the field in format dd.mm.yyyy.
- Week: week number, from week 1 (the first week that insects spent at the university of Wuerzburg, after being brought from the field) to week 50 (the end of the climate chamber experiment).
- Hour: hour of the day, from 0 to 1, with 0 starting at 00:00.
- Temp_CT: temperature programmed for the cool treatment on that day, at that hour.
- Temp_WT: temperature programmed for the warm treatment on that day, at that hour.
- Temp_HT: temperature programmed for the hot treatment on that day, at that hour.
- Temp_control: temperature recorded with iButtons under outdoor conditions on that day, at that hour.
- Moment: pre-experiment (before March 1st, when the climate chamber experiment started) or post-experiment (after March 1st).
File: Temp_comparison.csv
Description: data to produce Fig. S1.1 (code in Emergence-time_script)
Variables
- PlotID: unique identifier of each study site.
- time_period: six different periods. Either spring or summer temperatures of the periods 1981-2010, 2019 or 2020.
- mean_temperature: mean temperature of the time period considered.
- season: spring or summer.
File: Emergence-mass_script.html
Description: Code to produce the main results of the manuscript related to emergence mass.
File: Emergence-time_script.html
Description: Code to produce the main results of the manuscript related to emergence time.
File: Emergence_mass.csv
Description:
Variables
- Bee_No: individual identifier for each weighed bee.
- NestID: identifier of the nest from which the insect emerged.
- Treatment: climate chamber treatment (factor with four levels: HT, WT, CT, control).
- BeeSpecies: bee species name.
- sex: the sex of the emerged bee (male/female).
- weight_mg: the dry weight of the bee in mg.
- head_size_mm: wider part of the head of the bees in mm, measured with a digital calliper to the nearest 0.01 mm.
- PlotID: unique identifier of each study site.
- POINT_X: longitude of the site point of this region, coordinate system as "ETRS89 / UTM zone 32N".
- POINT_Y: latitude of the site point of this region, coordinate system as "ETRS89 / UTM zone 32N".
- Region: unique identifier of each study region.
- MAT: 30 year means of the monthly averaged mean daily air temperature in °C in 2 m height above ground from 01.01.1981 - 31.12.2010.
- AprilSept_PT: means of the temperature recorded with iButtons from April to September 2019.
- AprilAug_PT: means of the temperature recorded with iButtons from April to August 2019.
- AprilJuly_PT: means of the temperature recorded with iButtons from April to July 2019.
- AprilJune_PT: means of the temperature recorded with iButtons from April to June 2019.
- AprilMay_PT: means of the temperature recorded with iButtons from April to May 2019.
- MaySept_PT: means of the temperature recorded with iButtons from May to September 2019.
- MayAug_PT: means of the temperature recorded with iButtons from May to August 2019.
- MayJuly_PT: means of the temperature recorded with iButtons from May to July 2019.
- MayJune_PT: means of the temperature recorded with iButtons from May to June 2019.
- JuneSept_PT: means of the temperature recorded with iButtons from June to September 2019.
- JuneAug_PT: means of the temperature recorded with iButtons from June to August 2019.
- JuneJuly_PT: means of the temperature recorded with iButtons from June to July 2019.
- JulySept_PT: means of the temperature recorded with iButtons from July to September 2019.
- JulyAug_PT: means of the temperature recorded with iButtons from July to August 2019.
- AugSept_PT: means of the temperature recorded with iButtons from August to September 2019.
- Date_emergence: date of emergence in format m/d/yyyy.
- Day_emergence: day of emergence. Day 1 is the first day a weighed-bee emerged.
- DateOpening: the date the nest was processed in the laboratory in format m/d/yyyy.
- DayOpening: the day the nest was processed in the laboratory (day 1 is the first day of processing).
File: Emergence_time.csv
Description:
Variables
- IndividualNo: individual identifier for each emerged insect.
- NestID: identifier of the nest from which the insect emerged.
- Treatment: climate chamber treatment (factor with four levels: HT, WT, CT, control).
- Species: bee or wasp species name.
- Spring_Summer: if the species is considered as spring- or summer-emerging.
- Sex: the sex of the emerged insect (Male/Female).
- Date_emergence: date of emergence in format m/d/yyyy.
- Day_emergence: day of emergence. Day 1 is the first day of the experiment. Insects emerged between days 5 and 196.
- PlotID: unique identifier of each study site.
- Region: unique identifier of each study region.
- POINT_X: longitude of the site point of this region, coordinate system as "ETRS89 / UTM zone 32N".
- POINT_Y: latitude of the site point of this region, coordinate system as "ETRS89 / UTM zone 32N".
- Lat: latitude of the site point of this region, coordinate system as "GCS_WGS_1984".
- Lon: longitude of the site point of this region, coordinate system as "GCS_WGS_1984".
- MAT: 30 year means of the monthly averaged mean daily air temperature in °C in 2 m height above ground from 01.01.1981 - 31.12.2010.
- AprilSept_PT: means of the temperature recorded with iButtons from April to September 2019.
- AprilAug_PT: means of the temperature recorded with iButtons from April to August 2019.
- AprilJuly_PT: means of the temperature recorded with iButtons from April to July 2019.
- AprilJune_PT: means of the temperature recorded with iButtons from April to June 2019.
- AprilMay_PT: means of the temperature recorded with iButtons from April to May 2019.
- MaySept_PT: means of the temperature recorded with iButtons from May to September 2019.
- MayAug_PT: means of the temperature recorded with iButtons from May to August 2019.
- MayJuly_PT: means of the temperature recorded with iButtons from May to July 2019.
- MayJune_PT: means of the temperature recorded with iButtons from May to June 2019.
- JuneSept_PT: means of the temperature recorded with iButtons from June to September 2019.
- JuneAug_PT: means of the temperature recorded with iButtons from June to August 2019.
- JuneJuly_PT: means of the temperature recorded with iButtons from June to July 2019.
- JulySept_PT: means of the temperature recorded with iButtons from July to September 2019.
- JulyAug_PT: means of the temperature recorded with iButtons from July to August 2019.
- AugSept_PT: means of the temperature recorded with iButtons from August to September 2019.
- DateOpening: the date the nest was processed in the laboratory in format m/d/yyyy.
- DayOpening: the day the nest was processed in the laboratory (day 1 is the first day of processing).
Code/software
Analyses were conducted with R version 4.5.0. The code can be found in the html files 'Emergence-mass_script' and 'Emergence-time_script'. The versions of the packages used were:
plyr_1.8.9
dplyr_1.1.4
tidyverse_2.0.0
DHARMa_0.4.7
emmeans_1.11.1
performance_0.15.3
MuMIn_1.48.11
glmmTMB_1.1.11
mgcv_1.9-3
car_3.1-3
ggeffects_2.3.0
ggplot2_4.0.1
ggpubr_0.6.0
patchwork_1.3.2
ggtext_0.1.2
dagitty_0.3-4
lavaan_0.6-19
sf_1.0-21
terra_1.8-70
ggnewscale_0.5.2
rnaturalearth_1.1.0
rnaturalearthhires_1.0.0.9000
Note: "NA" is used in csv files to indicate "not applicable"
