Data and code from: Landscape context affects physiological but not genetic proxies of butterfly movement
Data files
May 27, 2026 version files 20.03 MB
-
Data.zip
20 MB
-
R_scripts.zip
19.52 KB
-
README.md
7.62 KB
Abstract
Land-use change is a major driver of biodiversity loss, yet its effects on intraspecific variation in dispersal remain understudied. Dispersal influences species’ capacity to persist in fragmented landscapes, but whether agricultural landscapes filter individuals based on their dispersal is unclear. We quantified intra- and interspecific variation in dispersal proxies (flight metabolic rate, body mass, and wing size) in three butterfly species differing in dispersal ability, sampled from agricultural and meadow landscapes. To distinguish landscape filtering (spatial sorting), phenotypic plasticity, and genetic adaptation, we combined phenotypic measurements with population genomic analyses based on whole-genome sequencing. In Maniola jurtina, females from agricultural landscapes exhibited higher peak flight metabolic rates, body mass, and an increased ability to sustain flight. More sustained flight in agricultural landscapes was also observed for Pieris napi. No landscape-related differences were detected in Thymelicus lineola. Despite phenotypic divergence, population genetic structure was weak across landscapes, although elevated runs of homozygosity in two species suggested limited dispersal and an increase in inbreeding. Our results indicate that agricultural landscapes can act as filters on dispersal-related traits at small spatial scales. Variation in dispersal is present within and among species and is likely mediated by plasticity or spatial sorting rather than local genetic adaptation.
Dataset DOI: 10.5061/dryad.sn02v6xkb
Description of the data and file structure
Data on three species of butterflies were collected during the summer of 2023, and it includes morphological measures (forewing length, body mass) and flight metabolic rate (FMR) measurements collected from two different landscape types (agricultural vs. meadow). A genetic sample was taken from each individual for whole-genome sequencing. Files include morphological and flight metabolic rate measurements and VCFs, heterozygosity, and ROH results that can be used to re-create the genetic analyses in the manuscript.
Files and variables
File: R_scripts.zip
Description: This folder includes all the R scripts used in the analysis and plots. The script needed to extract the FMR measures from the raw FMR data is numbered 1, and this creates the files used in the script numbered 2 (these files are also available in the data folder). Other scripts can be run in any order; the files necessary are in the data folder. File names and the first line of each script describe which analyses/plots from the manuscript can be re-created with that script.
File: Data.zip
Description: This folder includes all the different data files required to re-create different analyses. Files include both separate files for each species for raw FMR files (folder FMR_files23_species), background information, and FMR results (with species name in the file name), and combined dataframes including all species that were used for creating the plots in the manuscript. A genetic sample was taken from each individual for whole-genome sequencing. For genetic data files, include separate VCFs for each species and heterozygosity and ROH results for all species combined in a single file.
Files and variables
Folder: FMR_files23_species
Description: Includes rawtab files for each species required to extract the FMR values used in the analyses. Each species has its own folder named FMR_files23_speciesname. 23 refers to the year the data was collected (2023), and species are written with the first initial of the species name and genus (e.g., Mjurtina). "2nd" at the end of the Pnapi folder refers to the second generation (this species is multivoltine in Finland). Files are named following the pattern ID = species field ID, same as in other data files; f = indicating female; and the measurement temperature either a temperature range during the temperature tests (e.g. 23-24) or the final measurement temperature (e.g. 31), all in Celsius.
File: Measurements_all_species.csv
Description: Different dispersal proxies measured from each individual: flight metabolic rate (FMR), body mass, and forewing length
Variables
- ID: individual field ID, running number
- Cap_date: date of capturing the individual
- Species: species of the individual
- Location: sampling location
- Landscape: A=agriculture, M=meadow, the surrounding landscape type of the sampling site
- M_date: date of the FMR measurement
- Temp: set measurement temperature (Celsius)
- Mass: body mass of an individual in milligrams
- FW_length: forewing length in mm, measured to the closest 0.5 mm
- Measurer: person doing the FMR measurement
- M_time: time at which the FMR measurement was started
- Time_min: time in minutes after 8 am
- File.Name: file name for the raw .exp files with individual ID, f for female, and measurement temperature
- RMR: average resting metabolic rate extracted from the raw FMR measurement
- RMR.Temperature: average temperature during resting metabolic rate
- MMR: maximum metabolic rate
- MMR.Temperature: temperature during maximum metabolic rate
- Total.CO2: the total amount of CO2 in ml produced during the measurement period
- Total.End.CO2: the total amount of CO2 in ml produced during the last three minutes of the measurement
- Corrected.MMR: MMR corrected for body mass
- Corrected.CO2: total CO2 corrected for body mass
File: FMR_results_species.csv
Description: Different dispersal proxies measured from each individual: flight metabolic rate (FMR), body mass, and forewing length
Variables
- Same as for Measurements_all_species.csv variable,s except Mass is in grams
- Habitat: A=agriculture, M=meadow, the surrounding landscape type of the sampling site
File: Species_background_2023.csv
Description: Background information and morphological measures (body mass in g and forewing length in mm)
- Same as for FMR_results_species.csv variables
File: roh_all_species.csv
Description: Runs of homozygosity for each individual
Variables
- param_config: parameters used to assess runs of homozygosity
- id: individual field ID, running number
- Froh_genome: representing the proportion of the genome in ROH
- mean_lengthBps: the mean length of ROH in an individual
- Landscape: A=agriculture, M=meadow, the surrounding landscape type of the sampling site
- Species: species of the individual
File: Heterozygosity_all_species.csv
Description: heterozygosity estimates for each individual
Variables
- Field.ID: individual field ID, running number
- Species: species of the individual
- Location: sampling location
- Landscape: A=agriculture, M=meadow, the surrounding landscape type of the sampling site
- Specimen_ID: unique ID with a species and year prefix (e.g., PN23 = Pieris napi) and the field ID
- Heterozygosity: estimated heterozygosity of an individual
File: FMR_raw_all_species.csv
Description: Raw FMR data combined into a single file to run the curve analysis part of the FMR_curve_analysis.R script
Variables
- ID: individual field ID, running number
- Cap_date: date of capturing the individual
- Species: species of the individual
- Location: sampling location
- Landscape: A=agriculture, M=meadow, the surrounding landscape type of the sampling site
- M_date: date of the FMR measurement
- Temp: set measurement temperature (Celsius)
- Mass: body mass of an individual in grams
- FW_length: forewing length in mm, measured to the closest 0.5 mm
- Measurer: person doing the FMR measurement
- File.Name: file name for the raw .exp files with individual ID, f for female, and measurement temperature
- Seconds: seconds since start of the measurement
- Temperature: temperature at that second
- Smooth.Temperature: smoothed temperature to remove noise
- Flow.Volts: flow rate
- Smooth.Flow.Volts: smoothed flow volts to remove noise
- Flow.Lmin: flow rate in litres per minute
- CO2.Volts: CO2 flow rate
- CO2.PPM: CO2 in ppm
- CO2.MLhour: CO2 in ml/hour
- UID: unique ID with the field ID and species combined
File: Maniola_jurtina_sorted_5bp_depth_0.8_missing.pruned_0.1.vcf
Description: The variant calling file for Maniola jurtina used in the admixture and PCA analysis.
File: Pieris_napi_sorted_5bp_depth_0.8_missing.pruned_0.1.vcf
Description: The variant calling file for Pieris napi used in the admixture and PCA analysis.
File: Thymelicus_lineola_sorted_5bp_depth_0.8_missing.pruned_0.1.vcf
Description: The variant calling file for Thymelicus lineola used in the admixture and PCA analysis.
Code/software
Scripts are written in the R language and provided in R Markdown file format.
Access information
Other publicly accessible locations of the data:
- Raw sequencing data available at ENA with accession number PRJEB105221
Data was derived from the following sources:
- NA
