An increasing number of studies are using landscape genomics to investigate local adaptation in wild and domestic populations. Implementation of this approach requires the sampling phase to consider the complexity of environmental settings and the burden of logistical constraints. These important aspects are often underestimated in the literature dedicated to sampling strategies. In this study, we computed simulated genomic data sets to run against actual environmental data in order to trial landscape genomics experiments under distinct sampling strategies. These strategies differed by design approach (to enhance environmental and/or geographical representativeness at study sites), number of sampling locations and sample sizes. We then evaluated how these elements affected statistical performances (power and false discoveries) under two antithetical demographic scenarios. Our results highlight the importance of selecting an appropriate sample size, which should be modified based on the demographic characteristics of the studied population. For species with limited dispersal, sample sizes above 200 units are generally sufficient to detect most adaptive signals, while in random mating populations this threshold should be increased to 400 units. Furthermore, we describe a design approach that maximizes both environmental and geographical representativeness of sampling sites and show how it systematically outperforms random or regular sampling schemes. Finally, we show that although having more sampling locations (between 40 and 50 sites) increase statistical power and reduce false discovery rate, similar results can be achieved with a moderate number of sites (20 sites). Overall, this study provides valuable guidelines for optimizing sampling strategies for landscape genomics experiments.
The vulnerability of alpine environments to climate change presses an urgent need to accurately model and understand these ecosystems. Popularity in the use of digital elevation models (DEMs) to derive proxy environmental variables has increased over the past decade, particularly as DEMs are relatively cheaply acquired at very high resolutions (VHR; <1 m spatial resolution). Here, we implement a multiscale framework and compare DEM-derived variables produced by Light Detection and Ranging (LiDAR) and stereo-photogrammetry (PHOTO) methods, with the aim of assessing their relevance and utility in species distribution modelling (SDM). Using a case study on the arctic-alpine plant, Arabis alpina, in two valleys in the western Swiss Alps, we show that both LiDAR and PHOTO technologies can be relevant for producing DEM-derived variables for use in SDMs. We demonstrate that PHOTO DEMs, up to a spatial resolution of at least 1 m, rivalled the accuracy of LiDAR DEMs, largely owing to the customizability of PHOTO DEMs to the study sites compared to commercially available LiDAR DEMs. We obtained DEMs at spatial resolutions of 6.25 cm–8 m for PHOTO and 50 cm–32 m for LiDAR, where we determined that the optimal spatial resolutions of DEM-derived variables in SDM were between 1 and 32 m, depending on the variable and site characteristics. We found that the reduced extent of PHOTO DEMs altered the calculations of all derived variables, which had particular consequences on their relevance at the site with heterogenous terrain. However, for the homogenous site, SDMs based on PHOTO-derived variables generally had higher predictive powers than those derived from LiDAR at matching resolutions. From our results, we recommend carefully considering the required DEM extent to produce relevant derived variables. We also advocate implementing a multiscale framework to appropriately assess the ecological relevance of derived variables, where we caution against the use of VHR-DEMs finer than 50 cm in such studies.
Environmental change is intensifying the biodiversity crisis and threatening species across the tree of life. Conservation genomics can help inform conservation actions and slow biodiversity loss. However, more training, appropriate use of novel genomic methods and communication with managers are needed. Here, we review practical guidance to improve applied conservation genomics. We share insights aimed at ensuring effectiveness of conservation actions around three themes: (1) improving pedagogy and training in conservation genomics including for online global audiences, (2) conducting rigorous population genomic analyses properly considering theory, marker types and data interpretation and (3) facilitating communication and collaboration between managers and researchers. We aim to update students and professionals and expand their conservation toolkit with genomic principles and recent approaches for conserving and managing biodiversity. The biodiversity crisis is a global problem and, as such, requires international involvement, training, collaboration and frequent reviews of the literature and workshops as we do here.
Landscape genomic analyses associating genetic variation with environmental variables are powerful tools for studying molecular signatures of species' local adaptation and for detecting candidate genes under selection. The development of landscape genomics over the past decade has been spurred by improvements in resolutions of genomic and environmental datasets, allegedly increasing the power to identify putative genes underlying local adaptation in non‐model organisms. Although these associations have been successfully applied to numerous species across a diverse array of taxa, the spatial scale of environmental predictor variables has been largely overlooked, potentially limiting conclusions to be reached with these methods. To address this knowledge gap, we systematically evaluated performances of genotype–environment association (GEA) models using predictor variables at multiple spatial resolutions. Specifically, we used multivariate redundancy analyses to associate whole‐genome sequence data from the plant Arabis alpina L. collected across four neighboring valleys in the western Swiss Alps, with very high‐resolution topographic variables derived from digital elevation models of grain sizes between 0.5 m and 16 m. These comparisons highlight the sensitivity of landscape genomic models to spatial resolution, where the optimal grain sizes were specific to variable type, terrain characteristics, and study extent. To assist in selecting variables at appropriate spatial resolutions, we demonstrate a practical approach to produce, select, and integrate multiscale variables into GEA models. After generalizing fine‐grained variables to multiple spatial resolutions, a forward selection procedure is applied to retain only the most relevant variables for a particular context. Depending on the spatial resolution, the relevance for topographic variables in GEA studies calls for integrating multiple spatial scales into landscape genomic models. By carefully considering spatial resolutions, candidate genes under selection by a more realistic range of pressures can be detected for downstream analyses, with important applied implications for experimental research and conservation management of natural populations.