Data and code from: Using landscape genomics to define species distributions, delineate seed zones, and predict genomic offset to future climate for the interior spruce hybrid complex (Picea glauca, Picea engelmannii, and their hybrids
Data files
Jun 09, 2026 version files 39.79 MB
-
Data_and_code.zip
39.78 MB
-
README.md
8.58 KB
Abstract
Understanding how tree species adapt to climate is crucial for forest management under climate change. This study employs landscape genomics to investigate climate adaptation in the interior spruce complex (Picea glauca, P. engelmannii, and their hybrids) across western Canada. Using gradient forest modeling with 41,253 SNPs genotyped in 1692 natural interior spruce individuals from 252 populations, we identified winter temperature and moisture-related variables as key drivers of genomic variation. Both adaptive and neutral genetic variation showed similar patterns along climatic gradients, suggesting that population structure largely follows environmental clines. We delineated 11 seed zones based on genomic variation and climate relationships, achieving 88.2% concordance with previous genetic studies in defining species boundaries. Analysis of genetic offsets under future climate scenarios revealed potential risks of maladaptation, particularly in northern and eastern British Columbia. These predictions were validated using fitness-related measurements from common-garden experiments, which showed negative correlations between genetic offsets and traits such as height and diameter at breast height (DBH). Maladaptation predictions based on genetic offsets showed some differences compared with common-garden-based assessments in several regions, offering new perspectives on population vulnerability. Our results demonstrate the effectiveness of landscape genomics as a complementary approach to traditional methods for assessing climate adaptation in tree species. Our findings provide a scientific foundation for climate-smart reforestation strategies and offer practical guidance for forest management in the context of climate change.
Dataset DOI: 10.5061/dryad.xsj3tx9tp
Description of the data and file structure
This data and code .zip package contains the mandatory materials to reproduce the analyses and results of Ye et al. (2025): Using landscape genomics to define species distributions, delineate seed zones, and predict genomic offset to future climate for the interior spruce hybrid complex (Picea glauca, Picea engelmannii, and their hybrids)
Any questions or errors should be brought to the attention of Zhengyang (Lambert) Ye (yezhyang@mail.ubc.ca)
Files and variables
File: Data_and_code.zip
Description:
CODE
There are three R scripts with numerical prefixes that represent the order they should be executed in because some scripts generate files used in subsequent scripts. Their brief descriptions are:
- 1_Ye_et_al_SNP_formatting_calculation.R
This script formats the raw SNP data and calculates the SNP allele frequencies.
- 2_Ye_et_al_GF_modeling_seed_zone_delineation.R
This script develops gradient forest models for interior spruce and reveals its genomic variation.
Based on the models, it then defines species distribution and delineates seed zones.
- 3_Ye_et_al_Genetic_offset_analysis.R
This script does genetic offset calculations, including forward and reverse offset, and genetic offset predicted by taxon-specific models. It also contains content of validation using phenotypic data.
DATA FILES
3 data folders accompany Ye et al. (2025). Their brief descriptions are:
- 1. Spruce_SDM
This directory contains species distribution rasters for white spruce, Engelmann spruce, Sitka spruce (not used), and hybrid for the normal period 1961-1990 (Degner, 2023). The data was used for the study area map and validation of species distribution by the gradient forest model.
- 2. Spruce_SNP_data
This directory contains the complete SNP dataset along with the subset of putatively neutral SNPs, along with the information about the samples.
Spruce_NatSeedlotGeoData_21stOct2014.csv
This file contains seedlot information.
- Working_Seedlot: Seedlot number
- Prov_Seedlot: Province seedlot code
- lat: Latitude of the seedlot origin
- long: Longitude of the seedlot origin
- el: Elevation of the seedlot origin
Spruce_NaturalOrchard_IndivSeedlingInfo_July28th2014.xlsx
This file contains seedlings' information collected from the seedlots.
- Species: Picea spp. Spruces, the genus
- Subsp: Taxon. White spruce, Engelmann spruce, or hybrid
- WorkingSeedlot: Seedlot number
- Internal_ID: Seedlot number and individual number, separated by "D."
- SeedlingNo: Individual number
- No.Individuals_Seedlot: number of individuals in this seedlot
- FTGSC_Barcode: Sample barcode identifier. No values were available for this dataset
- SampleType: Seedlings. How the materials were sampled.
- ExperimentType: Outdoor raised bed. Type of the experiment.
- ExperimentalClimate: Experimental condition. Ambient in UBC,
- Province: seedlot origin is British Columbia or Alberta
- ProvSeedlotCode: Province seedlot code
- SeedlotType: Natural or selected. Only natural individuals were used in this study.
- SeedlotClass: Natural/Natural comparative/Orchard. Natural comparative set is a subset of Natural seedlings, which was used to compare with selected seedlings. N and NC were not distinguished in this study.
- Actual_BZ: Breeding zone.
- Comparative_BZ: Not used in this study.
- Block & LocationInBlock: Block number and location in the block. Block design.
Spruce_Totem_SNPtableAll_July13th.txt
This file contains the complete SNP data in IUPAC call letters. Each individual is displayed in "Seedlot number"D"individual number" format.
- probeset_id: The identifier assigned to a SNP probe on the genotyping array
- gcontig__pos: The genomic contig and the physical base-pair position of the SNP on the reference genome
- pos_gcontig: The physical position of the corresponding contig of the SNP
Spr_NaturalComp_Neutral85SNPfilter.map
This file contains the position of intergenic Non-Coding regions included in the Spruce array design. Used as the control group.
- 3. Phenotypic_data_for_validation
This directory contains the complete dataset of the Interior spruce climate change/genecology provenance trial (EP 670.71.12) of the BC Ministry of Forests (MoF).
masterfile16.csv
This file contains the recordings of fitness measurements of interior spruce individuals from different seed sources and growing sites. The measurements were taken at ages 3, 6, 10, and 16. The numbers following specific variables are the age of the tree.
- Tree_ID: Tree number across site
- Meas_seq: Tree number in each site (17 sites in extent)
- Site: Name of the sites
- Col, Row, Rep, Blk: location of the tree in the site.
- Prov: Provenance (seedlot, seed source)
- Tree: For each combination of site and provenance, there are four trees planted.
- HT: Height
- Dam: Damage
- DBH: Diameter at Breast Height (1.3m), only taken at 16 years old.
- Weev: Weevil damage
- Sunscold: sunscold damage
- Numbers following: xxx (e.g., HT) at age xx.
- Thin_age: Stand age when the last thinning treatment occurred
- "NA" or "." means there was no such damage, or the measurement is not applicable (e.g, 0 DBH for trees under breast height)
Sx CC seedlots ClimBC v5.4_Normal_1961_1990Y.csv
This file contains the seedlot information of the trials.
seedlot_code: seedlot code (different from those in SNP dataset)
Province: province or state of the origin
- lat: Latitude of the seedlot origin
- long: Longitude of the seedlot origin
- el: Elevation of the seedlot origin
- "_P" represents the variable from the provenance (seedlot)
- The following climate data were from ClimateBC ver 5.40 (about 66 provinces) and ClimateWNA 5.4 (the remainder, except one Ontario province, which is from ClimateNA 5.40) in the reference period 1961-1990
- Climate variable abbreviations:
- MAT mean annual temperature (°C)
- MWMT means the warmest month temperature (°C)
- MCMT means the coldest month temperature (°C)
- TD temperature difference between MWMT and MCMT, or continentality (°C)
- MAP mean annual precipitation (mm)
- MSP mean annual summer (May to Sept.) precipitation (mm)
- AHM annual heat-moisture index (MAT+10)/(MAP/1000))
- SHM summer heat-moisture index ((MWMT)/(MSP/1000))
- DD<0 (or DD_0) degree-days below 0°C, chilling degree-days
- DD>5 (or DD5) degree-days above 5°C, growing degree-days
- DD<18 (or DD_18) degree-days below 18°C, heating degree-days
- DD>18 (or DD18) degree-days above 18°C, cooling degree-days
- NFFD: the number of frost-free days, FFP: frost-free period
- bFFP the day of the year on which FFP begins
- eFFP the day of the year on which FFP ends
- PAS precipitation as snow (mm). For individual years, it covers the period between August in the previous year and July in the current year
- EMT extreme minimum temperature over 30 years (°C)
- EXT extreme maximum temperature over 30 years (°C)
- Eref Hargreaves reference evaporation (mm)
- CMD Hargreaves climatic moisture deficit (mm)
- MAR mean annual solar radiation (MJ m‐2 d‐1)
- RH mean annual relative humidity (%)
- CMI Hogg’s climate moisture index (mm)
- DD1040 degree-days above 10°C and below 40°C
Sx CC site ClimBC v5.4_Normal_1961_1990Y.csv
This file contains the site information of the trials.
- Site: site name, all in British Columbia. 18 test sites while Ftne was destroyed by wildfire, so 17 in extent.
- lat: Latitude of the seedlot origin
- long: Longitude of the seedlot origin
- el: Elevation of the seedlot origin
- "_S" represents the variable is from the site. The following climate data were from ClimateBC ver 5.40, except the Alberta and YK sites, which used Climate WNA v5.40 in the reference period 1961-1990
- Climate variable abbreviations are the same as the seedlot table
Access information
Other publicly accessible locations of the data:
- None
Data was derived from the following sources:
- Degner, J. (2023). Genomics of adaptation in interior spruce to past, present, and future climates of western Canada [Doctoral dissertation]. University of British Columbia.
