Shared ancestry, sexual size dimorphism, and reproductive ecology jointly shape sex differences in adult lifespan in insects
Data files
Sep 22, 2026 version files 286.57 KB
-
Dataset_longevity.csv
279.10 KB
-
README.md
7.48 KB
Abstract
Understanding the factors underlying sex differences in lifespan can provide valuable insights into lifespan evolution. Despite abundant species-level data, comparative studies on the drivers of sex differences in lifespan in invertebrates remain scarce compared with vertebrates. Here, we present a broad comparative analysis of sex differences in adult lifespan (SDL) in insects, using literature-derived data for 395 species spanning 17 orders. Across the dataset, female adults outlived males in 65 % of species, with a median lifespan 9 % longer than that of males. Mean lifespan was consistently female-biased across orders and functional groups, but the direction of SDL varied extensively at all taxonomic levels, even among congeneric species, consistent with the low phylogenetic signal observed for SDL (Pagel’s λ = 0.16). SDL was positively, albeit weakly, associated with sexual size dimorphism. Lifespan differences were generally more female-biased in functional groups where extended lifespan is expected to more strongly enhance female reproductive success, either because of time-consuming resource acquisition for egg production (e.g., predators) or prolonged search for oviposition substrates (e.g., parasitoids). Overall, our findings indicate that shared ancestry, sexual size dimorphism, and sex-specific selection pressures associated with reproductive ecology each contribute to variation in SDL, with no single factor clearly predominating.
Data set creators:
Tiit Teder & Barbora Žabova
Date created:
28 January 2026
Licensing
This data set is covered by CC0 1.0 Universal (CC0 1.0).
We encourage researchers interested in the re-use of this dataset and code to contact the authors of the dataset.
Data Files Overview
- File count: 1
- Total file size: 279 KB
- File formats: .csv
#1 'Dataset_longevity.csv':
csv-formatted file containing literature-derived data to analyze sex differences in adult lifespan in insects; corresponding R-script 'R_longevity.R'
Variables
The dataset contains the following variables:
- ID - species identifier
- genus, species - genus and species names as reported in the source article
- fullspecies - species name as used in the source article
- newspecies - species name as used in the analysis
- 1) a few species (5) names were updated to their currently accepted names
- 2) a few species (5) not included in the Open Tree Taxonomy were replaced by a species from the same genus
- order, family - taxonomic affiliation of the species
- hemi_holo - species classified according to developmental mode (hemimetabolous/holometabolous)
- larval_food - species classified according to larval diet (herbivores sensu stricto, granivores and frugivores were combined into a single herbivore category in the analyses); "." = data not available or not classified unambiguously
- adult_food - species classified according to the adult diet (herbivores sensu stricto and granivores were combined into a single herbivore category in the analyses); "." = data not available or not classified unambiguously
- mating - mating opportunities provided (yes/no); "." = not reported
- food - food provided to adults (yes/no); "." = not reported
- longevity_measure - reported measure of adult lifespan (mean/median/marginal_mean); "." = not reported
- M_longevity_N, F_longevity_N - sample sizes of males and females, respectively, for adult lifespan measurements; "." = not reported
- M_longevity, F_longevity - male and female adult lifespan, respectively
- M_long_SE, F_long_SE - standard errors of male and female lifespan estimates; "." = not reported
- size_stage - life stage at which body size was measured (pupal/adult)
- size_measure - type of body size metric used (DM = dry mass, FM = fresh mass, BL - body length, HW - head width, HL - head length, PNW - pronotum width, TL - thorax length, ABL - abdomen length, FWL - forewing length, HFL - hind femur length, HTL - hind tibia length)
- M_size, F_size - male and female body size, respectively
- SSD - sexual size dimorphism (see main text for calculation details)
- SSDsource - origin of SSD values, specifying whether they were obtained from the same study as lifespan data or from another study (ibidem/elsewhere)
- lnRR - sex difference in adult lifespan, calculated as the natural logarithm of the ratio of mean female to mean male adult lifespan
- author_year - authors and publication year of the source study
- journal_volume_pages - journal, volume and pages of the source study
- title - title of the source study
- full_reference - full bibliographic reference of the study from which data on sex differences in adult lifespan (and, when available, on sexual size dimorphism) were retrieved
- body_size_reference - full bibliographic reference of the study from which sexual size dimorphism data were obtained when these were unavailable in the study reporting sex differences in adult lifespan
Missing data code: '.'
Session Information = sessionInfo()
R version 4.5.1 (2025-06-13 ucrt)
Platform: x86_64-w64-mingw32/x64
Running under: Windows 10 x64 (build 19045)
Matrix products: default
LAPACK version 3.12.1
locale:
[1] LC_COLLATE=English_United Kingdom.utf8 LC_CTYPE=English_United Kingdom.utf8 LC_MONETARY=English_United Kingdom.utf8 LC_NUMERIC=C
[5] LC_TIME=English_United Kingdom.utf8
time zone: Europe/Tallinn
tzcode source: internal
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] phyloch_1.5-5 XML_3.99-0.19 colorspace_2.1-2 flextable_0.9.10 emmeans_2.0.1 orchaRd_2.2.1 R.rsp_0.46.0 patchwork_1.3.2 lubridate_1.9.4
[10] forcats_1.0.1 stringr_1.6.0 dplyr_1.1.4 purrr_1.1.0 readr_2.1.6 tidyr_1.3.2 tibble_3.3.1 tidyverse_2.0.0 devtools_2.4.6
[19] usethis_3.2.1 phytools_2.5-2 maps_3.4.3 ggplot2_4.0.1 data.table_1.17.8 nlme_3.1-168 metafor_4.8-0 numDeriv_2016.8-1.1 metadat_1.4-0
[28] Matrix_1.7-3 ape_5.8-1 rotl_3.1.0 openxlsx_4.2.8.1 readxl_1.4.5
loaded via a namespace (and not attached):
[1] mathjaxr_2.0-0 RColorBrewer_1.1-3 jsonlite_2.0.0 magrittr_2.0.4 estimability_1.5.1 ggbeeswarm_0.7.3 farver_2.1.2 rmarkdown_2.30
[9] fs_1.6.6 ragg_1.5.0 vctrs_0.7.0 memoise_2.0.1 askpass_1.2.1 htmltools_0.5.9 progress_1.2.3 curl_7.0.0
[17] DEoptim_2.2-8 cellranger_1.1.0 desc_1.4.3 cachem_1.1.0 uuid_1.2-1 igraph_2.2.0 lifecycle_1.0.5 iterators_1.0.14
[25] pkgconfig_2.0.3 R6_2.6.1 fastmap_1.2.0 digest_0.6.37 ps_1.9.1 pkgload_1.4.1 textshaping_1.0.4 labeling_0.4.3
[33] latex2exp_0.9.8 clusterGeneration_1.3.8 timechange_0.3.0 mgcv_1.9-3 httr_1.4.7 compiler_4.5.1 remotes_2.5.0 fontquiver_0.2.1
[41] withr_3.0.2 doParallel_1.0.17 S7_0.2.1 optimParallel_1.0-2 pkgbuild_1.4.8 R.utils_2.13.0 MASS_7.3-65 openssl_2.3.4
[49] sessioninfo_1.2.3 scatterplot3d_0.3-44 tools_4.5.1 vipor_0.4.7 rncl_0.8.7 otel_0.2.0 beeswarm_0.4.0 zip_2.3.3
[57] rentrez_1.2.4 R.oo_1.27.1 glue_1.8.0 quadprog_1.5-8 callr_3.7.6 R.cache_0.17.0 grid_4.5.1 generics_0.1.4
[65] gtable_0.3.6 tzdb_0.5.0 R.methodsS3_1.8.2 hms_1.1.4 xml2_1.5.1 foreach_1.5.2 pillar_1.11.1 splines_4.5.1
[73] lattice_0.22-7 tidyselect_1.2.1 fontLiberation_0.1.0 pbapply_1.7-4 knitr_1.51 fontBitstreamVera_0.1.1 xfun_0.55 expm_1.0-0
[81] stringi_1.8.7 pacman_0.5.1 evaluate_1.0.5 codetools_0.2-20 officer_0.7.0 gdtools_0.4.4 cli_3.6.5 xtable_1.8-4
[89] systemfonts_1.3.1 processx_3.8.6 Rcpp_1.1.0 coda_0.19-4.1 parallel_4.5.1 ellipsis_0.3.2 prettyunits_1.2.0 phangorn_2.12.1
[97] mvtnorm_1.3-3 scales_1.4.0 crayon_1.5.3 combinat_0.0-8 rlang_1.1.7 fastmatch_1.1-6 mnormt_2.1.1
Data collection
We aimed to collect data on sex differences in adult lifespan representing a taxonomically and ecologically broad range of insect species. Data were gathered from published case studies using two complementary approaches. The larger portion, roughly two-thirds of the data, was obtained through targeted searches in Google Scholar and Clarivate Web of Science conducted in 2023 and 2025. Since lifespan data are often reported as secondary outcomes and may not be mentioned in titles or abstracts, Google Scholar’s full-text search capability was particularly useful for locating such data. Our search queries included combinations of four generic terms that were neutral with respect to the magnitude and direction of our focal variable (i.e., SDL): one of several synonyms for adult lifespan (longevity, life span, lifespan), male*, female*, and insect* (or the name of a specific insect order). We used longevity as the primary synonym throughout, whereas other synonyms of adult lifespan were used only when searches with longevity and routine screening (see below) yielded insufficient data for a particular taxonomic group.
The rest of the data was collected by one of the authors through systematic screening of the tables of contents of a broad list of journals in entomology, ecology, and evolutionary biology. This effort involved both a one-time retrospective screening (conducted in 2003) and ongoing screening for articles published thereafter (until 2025). The screening process included regularly examining journals' tables of contents and reviewing the full-text of papers whose titles indicated they might contain relevant data. As lifespan data are typically presented in tables and figures, their identification within articles was generally straightforward. Because studies reporting sex-specific lifespan estimates in insects vary widely in their contexts, aims, and terminology for adult lifespan, assembling a complete set of published information for the more speciose and intensively studied orders (e.g., Lepidoptera, Coleoptera, Hymenoptera, Diptera) was not considered feasible from the outset. Although extensive, our dataset still should be regarded as a representative sample rather than complete record of the available data for these taxonomic groups. For smaller and less-studied orders, our targeted searches likely allowed us to assemble most of the published information on sex-specific adult lifespan. Inevitably, taxonomic and functional groups are not evenly represented in our dataset, reflecting differences in research focus as well as the practical challenges inherent to studying particular taxa. In total, we assembled sex-specific adult lifespan data for 395 insect species from 17 insect orders.
Where reported, sex-specific data on adult (or pupal) body size (preferably body mass) were extracted from the same studies reporting adult lifespan. Additional sex-specific body size data were derived from several previously published syntheses. Further body size data on specific species were obtained through targeted literature searches in Google Scholar using the species’ scientific name combined with relevant body mass-related terms (‘adult weight’, ‘adult mass’, ‘pupal weight’, ‘pupal mass’, ‘body size’).
Data extraction and criteria for eligibility
With the scope of the synthesis clearly defined, the eligibility of the identified studies was evaluated during the data collection process. The primary criterion for inclusion was that a study had to report mean values of both male and female adult lifespan of an insect species. Six studies presenting median values of lifespan were also accepted. Only laboratory-based data were considered, where major external sources of mortality, predators and parasitoids in particular, were explicitly excluded.
For a substantial proportion of species, lifespan data were available from multiple studies and/or multiple experimental treatments. To maximize comparability among species, we selected the dataset obtained under conditions that most closely approximated those considered optimal, rather than using all available datasets. Specifically, selection was based on the extent to which the experimental conditions met the following criteria: 1) ad libitum access to water and food, 2) the opportunity to mate, 3) the availability of an oviposition substrate (e.g., host plant), 4) the highest food quality (defined as the treatment associated with the highest larval survival, shortest development time, or largest body size, depending on the data available), and 5) an optimal temperature regime (assessed based on survival rates or conditions typical for the species’ biogeographic realm). When multiple datasets met these criteria equally well, priority was given to the dataset with the largest sample size. We explicitly excluded data from stress treatments, such as exposure to chill coma, heat shock, or biological and artificial toxins.
While we generally did not impose any taxonomic restrictions on the insects included in this study, we did not consider those with complex caste systems, such as ants, social bees and wasps. Also, we did not include a few studies in which the species' adult lifespan included diapause.
Volume-based measures (i.e., fresh or dry body mass) were the preferred metric for sex-specific body size, although linear measurements were also accepted when body mass data were unavailable. When body size and lifespan data were reported together, both were extracted from the same study and treatment. If size data came from a different study, we used data from treatments where the conditions for juvenile development were likely closest to optimal. The choice of treatment was guided primarily by measures of pre-adult survival (egg-to-adult, larval-pupal, or larval, depending on availability) or, when survival data were unavailable, by juvenile development time. Both metrics are commonly used indicators of developmental environment in insects, with higher survival and shorter development time generally reflecting more benign (laboratory) conditions. Although SSD can show some population-level variation within species, it is typically limited relative to the variation observed among insect species.
Graphically presented lifespan and body size data were extracted using WebPlotDigitizer 4.3.
Species were categorized into functional groups, where possible, based on the feeding ecology of their adult and larval stages.
Deriving SDL and SSD values from collected data
For each species, we calculated the natural logarithm of the ratio of mean female to mean male adult lifespan. This metric served as the effect size quantifying sex differences in adult lifespan (SDL) in our meta-analytic models. The resulting species-specific SDL values are symmetrically distributed around zero, with SDL = 0 indicating no sex difference in adult lifespan, and positive and negative SDL values indicating female-biased and male-biased lifespan, respectively. This choice of effect size permits direct calculation of the standard error of SDL when the standard errors of the sex-specific means are available, as was the case for most species in our dataset.
SSD was quantified using the widely used algorithm of Lovich and Gibbons (1992): [(size of the larger sex) / (size of the smaller sex)] – 1, with a positive sign assigned to the resulting values when females are larger and a negative sign when males are larger. Because body size was reported either as body mass (85 % of the 266 species included in the meta-analysis) or as various linear measurements, we standardized all size data to ensure comparability of SSD estimates across species. When body mass was available, it was used directly in the calculations. When only linear measurements (e.g., body length, head width) were available, these were converted to mass equivalents prior to estimating SSD. In particular, when body length was provided, we used published taxon-specific allometric equations to estimate body mass. For other linear size measures lacking appropriate allometric relationships, we applied cube transformation to approximate mass under isometric scaling assumptions before calculating SSD.
