Data from: Ecology and biogeography of sexual size dimorphism in squamates
Abstract
Aim: Sexual size dimorphism (SSD) is widespread in the animal kingdom. The direction and magnitude of SSD differ considerably across taxa, potentially due to different selective forces acting on female and male sizes. We assembled a comprehensive database of mean body sizes of females and males within squamate species. We then tested for associations between the degree and direction of sexual size dimorphism and environmental factors, clutch sizes, reproductive modes, substrate types, and species richness (a common measure of interspecific competition).
Location: Global
Time period: Present
Major taxa studied: Squamata (Reptilia)
Methods: We studied SSD patterns and their correlates for 11,792 squamate species. We also tested the effect of the number of putative competitors on SSD within (~9915 km2) grid cells. We applied phylogenetic path analysis and phylogenetic generalized least squares regression (PGLS) at the species, and applied spatial auto-regressive (SAR) multiple regressions at the assemblage levels.
Results: In general, snake females are larger than males, whereas male lizards are larger, on average, than females. Female squamates in general are larger than males in cold regions, while in warm regions, particularly in deserts, males were usually larger than females. SSD became more female-biased (i.e., larger females) as clutch size increased, and viviparous taxa had more female-biased SSD. There was little relationship between the magnitude of SSD and species richness. Sexual size dimorphism did not vary significantly across substrate types.
Main conclusion: We suspect that the mechanisms driving squamate SSD differ between oviparous and viviparous taxa. The more female-biased SSD in colder regions is likely driven by fecundity selection, while a male bias in warmer regions may be associated with sexual selection. However, we found little evidence to suggest natural selection for substrates, or resource-based competition, affects squamate sexual size dimorphism and suspect the underlying hypotheses may be flawed, and/or that species richness is a poor measure of the intensity of interspecific competition.
https://doi.org/10.5061/dryad.r7sqv9sp5
Description of the data and file structure
The data were collected from published papers and books, museum collections, and personal communications, primarily by S.M. and G.M.
Files and variables
File: Data.zip
Description: Prior to sharing the files, we wish to confirm that the supporting data contains no precise species occurrence or distribution details. We have ensured that the provided datasets will not adversely affect conservation efforts or pose any risk to vulnerable, endangered, or otherwise threatened species.
This dataset includes 2 types:
1: scripts, these are the R scripts that were used in this paper, including data imputation, making a hexagon grid, and analyses.
CodeMkeGrid.R This script generates a hexagonal spatial grid derived from squamate distribution polygons (hexagonal binning was performed using squamate range maps sourced from the GARD dataset (updated from Roll et al., 2017 and Caetano et al., 2022). Note: High-resolution polygons are maintained by the GARD team and are available upon request.)
Imputation.R This script generates the missing traits for suqamates.
SSD.AnalysisJune25.R This script executes all regression analyses presented in the study.
2: data, includes: 1 csv data of species level data; 2 R.data of assemblage level (based on 1 csv data); 3 imputation data (for missing species); 4 tree data.
Assemblage folder
Assem.107KM.Dec26.RData This RData is the data of assemble level (hexagon)
- "grid_id" grid id of each hexagon (no unit)
- "uniqueComm" the community (i.e., species composition) ID of each hexagon (no unit)
- "Long" the longitude of each hexagon (°)
- "Lat" the latitude of each hexagon (°)
- "Tem" median temperature of each hexahon (℃)
- "TemSea" median temperature seasonality of each hexagon (Coefficient of Variation, no unit)
- "Pre" median precipitation of each hexagon(mm/year)
- "PreSea" median precipitation seasonality of each hexagon (Coefficient of Variation, no unit)
- "NPP" median Net primary productivety of each hexagon (
g
⋅m−2
⋅yr−1)
- "cell_area_proportion" The proportion of each hexagonal cell covered by land (no unit)
- "spRichness" species richness of each hexagon (no unit)
- "massSSD" median sexual size dimorphism of squamates (g)
- "ABSmassSSD" median absolute sexual size dimorphism of squamates (no unit)
- "CS" cluth size of squamates of each hexagon (no unit)
- "FemaleMass" median female body mass of each hexagon (g)
- "MaleMass" median male body mass of each hexagon (g)
- "ViviNumb" richness of viviparous squamates of each hexagon (no unit)
- "VivPro.all" proportion of viviparous squamates of each hexagon (no unit)
- "RawRichness" species richness for all squamates in the hexagon with observed trait data (no unit)
- "RawmassSSD" median sexual size dimorphism for all squamates in the hexagon with observed trait data (no unit)
- "RawABSmassSSD" median sexual size dimorphism for all squamates in the hexagon with observed trait data (no unit)
- "RawCS" median clutch size for all squamates in the hexagon with observed trait data (no unit)
- "RawFemaleMass" median female mass for all squamates in the hexagon with observed trait data (g)
- "RawMaleMass" median male mass for all squamates in the hexagon with observed trait data (g)
- "RawViviNumb" species richness of viviparous squamates in the hexagon with observed body mass data (no unit)
- "VivPro.Raw" proportion of viviparous squamates in the hexagon with observed body mass data (no unit)
- "ORichness" species richness of oviparous squamates in the hexagon (no unit)
- "OmassSSD" median sexual size dimorphism of oviparous squamates in the hexagon (no unit)
- "OFemaleMass" media mass of female squamates in the hexagon (g)
- "OMaleMass" median mass of male squamates in the hexagon (g)
- "OCS" median clutch size of oviparous squamates in the hexagon (no unit)
- "VRichness" richness of viviparous squamates in the hexagon (no unit)
- "VmassSSD" median sexual size dimorphism of viviparous squamates in the hexagon (no unit)
- "VFemaleMass" median female mass of viviparious squamates in the hexagon (g)
- "VMaleMass" median male mass of viviparious squamates in the hexagon (g)
- "VCS" median clutch size of viviparous squamates in the hexagon (no unit)
- "LizardRichness" richness of lizards in the hexagon (no unit)
- "LizardmassSSD" median sexual size dimorphism of lizards in the hexagon (no unit)
- "LizardABSmassSSD" median absolute sexual size dimorphism of lizards in the hexagon (no unit)
- "LizardFemaleMass" median female mass of lizards in the hexagon (g)
- "LizardMaleMass" median male mass of lizards in the hexagon (g)
- "LizardCS" median clutch size of lizards in the hexagon (no unit)
- "LizardViviNumb" richness of viviparous lizard in the hexagon (no unit)
- "VivPro.Lizard" proportion of viviparous lizards in the hexagon (no unit)
- "SnakeRichness" richness of snakes in the hexagon (no unit)
- "SnakemassSSD"median sexual size dimorphism of snakes in the hexagon (no unit)
- "SnakeABSmassSSD"median absolute sexual size dimorphism of snakes in the hexagon (no unit)
- "SnakeFemaleMass"median female mass of snakes in the hexagon (g)
- "SnakeMaleMass"median male mass of snakes in the hexagon (g)
- "SnakeCS" median clutch size of snakes in the hexagon (no unit)
- "SnakeViviNumb" richness of viviparous snakes in the hexagon (no unit)
- "VivPro.Snake" proportion of viviparous snakes in the hexagon (no unit)
- "squamateNRI" Net relatedness index of squamates in the hexagon (no unit)
- "lizardNRI" Net relatedness index of lizards in the hexagon (no unit)
- "snakeNRI" Net relatedness index of snakes in the hexagon (no unit)
- "geometry" geometry of each hexagon (sf object) (no unit)
Imputation folder
LizardDataJuly21.csv imputed date for lizards. NA=missing traits
- "species": binomial (no unit)
- "ReproMode": reproductive mode (no unit)
- "CS": clutch size (no unit)
- "FemaleMass": female mass (g)
- "MaleMass": male mass (g)
- "HabitatType": habitat type of each species (no unit)
SnakeDataJuly21.csv imputed date for snakes. NA=missing traits
- "species": binomial (no unit)
- "ReproMode": reproductive mode (no unit)
- "CS": clutch size (no unit)
- "FemaleMass": female mass (g)
- "MaleMass": male mass (g)
- "HabitatType": habitat type of each species (no unit)
Tree Folder
- 10_10792tree.nex 10 trees of squamates that were used for missing data imputation
- phylizards.txt 1 of the 10 trees that were used for the analyses in the text.
Code/software
The scripts are based on R software, with the necessary packages and data to be read described in the code.
