Figure legends
Figure 1. A) Known distribution range of Psammodromus algirus (based on Bons and Geniez 1996 [north-west Africa] and Pleguezuelos 1997 [Iberian Peninsula]). The gaps within the range that correspond to artificial land uses according to Corine database wereconsidered as presences when embedded within the range. A question mark is placed where the species is known to inhabit but there is no accurate information about its distribution (27 presences scattered all over North Africa; grey squares in the figure). The discontinuous line defines the assumed southern edge of the distribution range according to IUCN. B) Bioclimatic gradient defined by the environmental PCA for all the western Mediterranean region (potential distribution range ofP. algirus ). Black circles mark the location of the sampled populations (L = Lerma, B = Brihuega, N = Navacerrada, P = El Pardo, and A = Aranjuez). On the right, the location of these populations and range edges (C/W = coldest/wettest, W/D = warmest/driest) within the environmental gradient from a single PC axis that was defined by using the scree plot criterion.
Figure 2. Distribution of adjusted R2 values for regression models obtained with four different randomized datasets (see text for details).
Figure 3. Inferred distribution range. In dark green, predicted range #1 (grid cells with the same environmental scores than the sampled populations); and in clear green, predicted range #2 (inferred by extrapolating GEAM to any possible combination of alleles at the loci included in the final model). All cells within predicted range #1 were also included in range #2. Areas beyond 8 km of the continuum inferred range were removed (see text).
Figure 4. Inferred distribution ranges #1 (cells with same environmental values than sampling locations; A) and #2 (cells that were inferred to be suitable for any possible combination of alleles at the loci included in our GEAM; B) at the resolution scale that was comparable with previously published ranges.
Figure 5. Distribution of the proportion of the actual species’ range that was inferred by completely randomized datasets incorporating the outlier selection step (see text for details). The proportion of range inferred by extrapolation of GEAM is marked with an arrow.
Figure 6. Heatmap representing the number of adapted genotypes per grid cell according to our GEAM.