A range-wide postglacial history of Swiss stone pine based on molecular markers and palaeoecological evidence
Data files
Feb 20, 2023 version files 418.35 KB
-
Gugerli_etal_JoB2023_cpSSRs-BAPS_Pop_indices.csv
833 B
-
Gugerli_etal_JoB2023_cpSSRs-BAPS_Pop_names.csv
3.04 KB
-
Gugerli_etal_JoB2023_cpSSRs-BAPS_Pop_PopGroupsUnlinkedLoci.csv
57.64 KB
-
Gugerli_etal_JoB2023_nSSRs_STRUCTURE.csv
344.52 KB
-
Gugerli_etal_JoB2023_PopulationSampling.csv
8.55 KB
-
README.md
3.76 KB
Aug 17, 2026 version files 418.52 KB
-
Gugerli_etal_JoB2023_cpSSRs-BAPS_Pop_indices.csv
833 B
-
Gugerli_etal_JoB2023_cpSSRs-BAPS_Pop_names.csv
3.04 KB
-
Gugerli_etal_JoB2023_cpSSRs-BAPS_Pop_PopGroupsUnlinkedLoci.csv
57.64 KB
-
Gugerli_etal_JoB2023_nSSRs_STRUCTURE_rev08-2026.csv
344.52 KB
-
Gugerli_etal_JoB2023_PopulationSampling.csv
8.55 KB
-
README.md
3.92 KB
Abstract
Aim: Knowing a species’ response to historical climate shifts helps understanding its perspectives under global warming. We infer the hitherto unresolved postglacial history of Pinus cembra. Using independent evidence from genetic structure and demographic inference of extant populations, and from palaeoecological findings, we derive putative refugia and re-colonisation routes.
Location: European Alps and Carpathians.
Taxa: Pinus cembra.
Methods: We genotyped nuclear and chloroplast microsatellite markers in nearly 3,000 individuals from 147 locations across the entire natural range of P. cembra. Spatial genetic structure (Bayesian modelling) and demographic history (Approximate Bayesian Computation) were combined with palaeobotanical records (pollen, macrofossils) to infer putative refugial areas during the Last Glacial Maximum (LGM) and re-colonisation of the current range.
Results: We found distinct spatial genetic structure, despite low genetic differentiation even between the two disjunct mountain ranges. Nuclear markers revealed five genetic clusters aligned East–West across the range, while chloroplast haplotype distribution suggested nine clusters. Spatially congruent separation at both marker types highlighted two main genetic lineages in the East and West of the range. Demographic inference supported early separation of these lineages dating back to a previous interstadial or interglacial c. 210,000 years ago. Differentiation into five biologically meaningful genetic clusters likely established during post-glacial re-colonisation.
Main conclusions: Combining genetic and palaeoecological evidence suggests that P. cembra primarily survived the LGM in “cold period” refugia south of the Central European Alps and near the Carpathians, from where it expanded during the Late Glacial into its current Holocene “warm period” refugia. This colonisation history has led to the distinct East–West structure of five genetic clusters. The two main genetic lineages likely derived from ancient divergence during an interglacial or interstadial. The respective contact zone (Brenner line) matches a main biogeographic break in the European Alps also found in herbaceous alpine plant species.
Two data sets comprising (i) 1 file with individual genotypes of nuclear microsatellites, each individual represented on two lines, and (ii) 3 files formatted for input in BAPS software (Corander et al., Genetics 2003). A further file comprises coding and geographic information for each sampling location, including sample number for genotyping for the nuclear microsatellite data set (minor deviations for the chloroplast data set, indicated in Suppl. Material, Table S1.1 in Gugerli et al., Journal of Biogeography 2023). Note that the dataset was updated in August 2026 because a single genotype call was detected to be wrong: individual CH132_007: locus HJM is correctly homozygous for allele 180, whereas it was formerly called heterozygous 180/199.
Description of the data and file structure
Gugerli_etal_JoB2023_PopulationSampling.csv: Code: population code with 2 letters for country of origin + population number
Range: location of populations in the European Alps or the Carpathian Mountains Country: nationality of sampling site (cf. country code in Code)
Population: local name of sampling location LatNN: latitude North (WGS84, in decimal degrees)
LonE: longitude East (WGS84, in decimal degrees)
SampleNo: number of samples genotyped for nuclear microsatellites
Gugerli_etal_JoB2023_nSSRs_STRUCTURE_rev08-2026.csv
Code_individual: label of individual, composed of population code and individual number
Number_population: consecutive number of populations
Pc7: allele calls for nuclear microsatellite locus Pc7b
Pc1b: allele calls for nuclear microsatellite locus Pc1b
Pc23: allele calls for nuclear microsatellite locus Pc23
Pc18: allele calls for nuclear microsatellite locus Pc18
Pc22: allele calls for nuclear microsatellite locus Pc22
Pc35: allele calls for nuclear microsatellite locus Pc35
28Z: allele calls for nuclear microsatellite locus 28Z
BUG: allele calls for nuclear microsatellite locus BUG
CQG: allele calls for nuclear microsatellite locus CQG
HJM: allele calls for nuclear microsatellite locus HJM
YAU: allele calls for nuclear microsatellite locus YAU
(note: for running STRUCTURE, the first two column headings need be deleted and locus codes shifted to the left, so that the first row contains two cells less than all remaining rows)
Gugerli_etal_JoB2023_cpSSRs-BAPS_Pop_PopGroupsUnlinkedLoci.csv:
PT15169: allele call for chloroplast microsatellite PT15169
PT26081: allele call for chloroplast microsatellite PT26081
PT36480: allele call for chloroplast microsatellite PT36480
PT63718: allele call for chloroplast microsatellite PT63718
Pop_No: consecutive number of population
(note: for running BAPS, the first row with column labels needs to be removed)
Gugerli_etal_JoB2023_cpSSRs-BAPS_Pop_names.csv: Code: population code with 2 letters for country of origin + population number Class: ?population? for indicating population-wise sampling scheme
Pop_No: consecutive number of population
(note: for running BAPS, the first row with column labels needs to be removed)
Gugerli_etal_JoB2023_cpSSRs-BAPS_Pop_indices.csv
First_Ind: consecutive number of first individual of each population as listed in Gugerli_etal_JoB2023_cpSSRs-BAPS_Pop_PopGroupsUnlinkedLoci.csv
(note: for running BAPS, the first row with the column label needs to be removed)
Sharing/Access information
Links to other publicly accessible locations of the data:
none
Data were derived from the following source :
Original data, with parts of the dataset or samples previously included in previous publications (see respective list of DOIs)
Code/Software
Analyses were performed with the publicly available software STRUCTURE 2.3.4 (Pritchard et al., Genetics 2000) for the nuclear microsatellite dataset and BAPS 6.0 (Corander et al., Genetics 2008; currently under development as an R package).
All methodological approaches are described in the associated paper (Gugerli et al., Journal of Biogeography 2023).
Genotype file with individual genotypes from 11 nuclear microsatellites (2 lines per individual), prepared for input in the clustering software STRUCTURE (Pritchard et al., Genetics 2000), available at https://web.stanford.edu/group/pritchardlab/structure_software/release_versions/v2.3.4/html/structure.html:
Gugerli_etal_JoB2023_nSSRs_STRUCTURE.csv
Gugerli_etal_JoB2023_nSSRs_GenePop.csv
Three files with data from 4 chloroplast microsatellites, formatted for input in the clustering software BAPS (Corander et al., Genetics 2003), currently under development as an R package at https://github.com/ocbe-uio/rBAPS:
Gugerli_etal_JoB2023_cpSSRs-BAPS_Pop_PopGroupsUnlinkedLoci.csv
Gugerli_etal_JoB2023_cpSSRs-BAPS_Pop_names.csv
Gugerli_etal_JoB2023_cpSSRs-BAPS_Pop_indices.csv
Populations file with information on sampling sites (population codes, mountain range, country, population name, geographic coordinates in WGS84 format) and sample numbers.
Gugerli_etal_JoB2023_PopulationSampling.csv
Changes after Feb 20, 2023: One genotype call was found to be erroneously scored as heterozygote but is homozygote: individual CH132_007, locs HJM, correct genotype 180/180 instead of 180/199
- Lendvay, Bertalan et al. (2014), Genetic structure in Pinus cembra from the Carpathian Mountains inferred from nuclear and chloroplast microsatellites confirms post-glacial range contraction and identifies introduced individuals, Tree Genetics & Genomes, Journal-article, https://doi.org/10.1007/s11295-014-0770-9
- Höhn, Maria et al. (2009), Variation in the chloroplast DNA of Swiss stone pine (<i>Pinus cembra</i>L.) reflects contrasting post-glacial history of populations from the Carpathians and the Alps, Journal of Biogeography, Journal-article, https://doi.org/10.1111/j.1365-2699.2009.02122.x
- Gugerli, Felix; Rüegg, Martin; Vendramin, Giovanni G. (2009), Gradual decline in genetic diversity in Swiss stone pine populations (Pinus cembra) across Switzerland suggests postglacial re-colonization into the Alps from a common eastern glacial refugium, Botanica Helvetica, Journal-article, https://doi.org/10.1007/s00035-009-0052-6
- Gugerli, Felix et al. (2023), A range‐wide postglacial history of Swiss stone pine based on molecular markers and palaeoecological evidence, Journal of Biogeography, Journal-article, https://doi.org/10.1111/jbi.14586
