DOI: 10.5061/dryad.ck6q52n
title: Data from: Spatiotemporally explicit demographic modelling supports a joint effect of historical barriers to dispersal and contemporary landscape composition on structuring genomic variation in a red-listed grasshopper

description: Inferring the processes underlying spatial patterns of genomic variation
 is fundamental to understand how organisms interact with landscape
 heterogeneity and to identify the factors determining species
 distributional shifts. Here, we employ genomic data (ddRADSeq) to test
 biologically-informed models representing historical and contemporary
 demographic scenarios of population connectivity for the Iberian
 cross-backed grasshopper Dociostaurus hispanicus, a species with a narrow
 distribution that currently forms highly fragmented populations. All
 models incorporated biological aspects of the focal taxon that could
 hypothetically impact its geographical patterns of genomic variation,
 including (a) spatial configuration of impassable barriers to dispersal
 defined by topographic landscapes not occupied by the species, (b)
 distributional shifts resulted from the interaction between the species
 bioclimatic envelope and Pleistocene glacial cycles, and (c) contemporary
 distribution of suitable habitats after extensive land clearing for
 agriculture. Spatiotemporally-explicit simulations under different
 scenarios considering these aspects and statistical evaluation of
 competing models within an Approximate Bayesian Computation (ABC)
 framework supported spatial configuration of topographic barriers to
 dispersal and human-driven habitat fragmentation as the main factors
 explaining the geographical distribution of genomic variation in the
 species, with no apparent impact of hypothetical distributional shifts
 linked to Pleistocene climatic oscillations. Collectively, this study
 supports that both historical (i.e., topographic barriers) and
 contemporary (i.e., anthropogenic habitat fragmentation) aspects of
 landscape composition have shaped major axes of genomic variation in the
 studied species and emphasizes the potential of model-based approaches to
 gain insights into the temporal scale at which different processes impact
 the demography of natural populations.
description: Samples_Dociostaurus_hispanicusSampling localities for all genotyped individuals and codes used in the different analyses and genomic datasets.Dociostaurus_hispanicusGenomic dataset (VCF format) obtained with STACKS for all genotyped populations and individuals.Dociostaurus_hispanicus_10000lociInput file (STR format) with 10,000 unlinked SNPs used for performing Bayesian clustering (STRUCTURE) and principal component (ADEGENET) analyses.Dociostaurus_hispanicus_1250lociInput file (STR format) with 1,250 unlinked SNPs used for obtaining observed summary statistics (ARLSUMSTAT) and performing Bayesian clustering (STRUCTURE) and principal component (ADEGENET) analyses.Dociostaurus_hispanicus_MIROCRaster map (TIF format) obtained by projecting the environmental niche model (ENM) from MAXENT to the last glacial maximum (LGM, 21 ka) bioclimatic conditions derived from the MIROC-ESM general atmospheric circulation model.Dociostaurus_hispanicus_CCSMRaster map (TIF format) obtained by projecting the environmental niche model (ENM) from MAXENT to the last glacial maximum (LGM, 21 ka) bioclimatic conditions derived from the CCSM4 general atmospheric circulation model.Dociostaurus_hispanicus_CURRENTRaster map (TIF format) obtained by projecting the environmental niche model (ENM) from MAXENT to current bioclimatic conditions (WORLDCLIM).StairwayPlotGenomic datasets (VCF format, exported from STACKS) and site frequency spectra (SFS) for each studied population and used for demographic reconstructions in STAIRWAY PLOT.LAYERS_Splatche2Raster files (ASCII format) used for running demographic and genetic simulations in SPLATCHE2 under each tested model. For model description, see METHODS section and Table 2.