skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: Optimising Sampling Design for Landscape Genomics
ABSTRACT Landscape genomic approaches for detecting genotype‐environment associations (GEA), isolation by distance (IBD) and isolation by environment (IBE) have seen a dramatic increase in use, but there have been few thorough analyses of the influence of sampling strategy on their performance under realistic genomic and environmental conditions. We simulated 24,000 datasets across a range of scenarios with complex population dynamics and realistic landscape structure to evaluate the effects of the spatial distribution and number of samples on common landscape genomics methods. Our results show that common analyses are relatively robust to sampling scheme as long as sampling covers enough environmental and geographic space. We found that for detecting adaptive loci and estimatingIBE, sampling schemes that were explicitly designed to increase coverage of available environmental space matched or outperformed sampling schemes that only considered geographic space. When sampling does not cover adequate geographic and environmental space, such as with transect‐based sampling, we detected fewer adaptive loci and had higher error when estimatingIBDandIBE. We found thatIBDcould be detected with as few as nine sampling sites, while large sample sizes (e.g., greater than 100 individuals) were crucial for detecting adaptive loci andIBE. We also demonstrate that, even with optimal sampling strategies, landscape genomic analyses are highly sensitive to landscape structure and migration—when spatial autocorrelation and migration are weak, commonGEAmethods fail to detect adaptive loci.  more » « less
Award ID(s):
1845682
PAR ID:
10601209
Author(s) / Creator(s):
; ;
Publisher / Repository:
Wiley
Date Published:
Journal Name:
Molecular Ecology Resources
Volume:
25
Issue:
3
ISSN:
1755-098X
Subject(s) / Keyword(s):
adaptation ecological genetics genotype-environment association isolation by environment population genetics simulation
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; (, Molecular Ecology)
    ABSTRACT Coalescent modelling of hybrid zones can provide novel insights into the historical demography of populations, including divergence times, population sizes, introgression proportions, migration rates and the timing of hybrid zone formation. We used coalescent analysis to determine whether the hybrid zone between phylogeographic lineages of the Plateau Fence Lizard (Sceloporus tristichus) in Arizona formed recently due to human‐induced landscape changes, or if it originated during Pleistocene climatic shifts. Given the presence of mitochondrial DNA from another species in the hybrid zone (Southwestern Fence Lizard,S. cowlesi), we tested for the presence ofS. cowlesinuclear DNA in the hybrid zone as well as reassessed the species boundary betweenS. tristichusandS. cowlesi. No evidence ofS. cowlesinuclear DNA is found in the hybrid zone, and the paraphyly of both species raises concerns about their taxonomic validity. Introgression analysis placed the divergence time between the parental hybrid zone populations at approximately 140 kya and their secondary contact and hybridization at approximately 11 kya at the end of the Pleistocene. Introgression proportions estimated for hybrid populations are correlated with their geographic distance from parental populations. The multispecies coalescent with migration provided significant support for unidirectional migration moving from south to north, which is consistent with spatial cline analyses that suggest a slow but steady northward shift of the centre of the hybrid zone over the last two decades. When analysing hybrid populations sampled along a linear transect, coalescent methods can provide novel insights into hybrid zone dynamics. 
    more » « less
  2. ; ; ; ; ; (, Molecular Ecology)
    ABSTRACT In integrative distributional, demographic and coalescent (iDDC) modelling, a critical component is the statistical relationship between habitat suitability and local population sizes. This study explores this relationship in twoEnyaliuslizard species from the Brazilian Atlantic Forest: the high‐elevationE. iheringiiand low‐elevationE. catenatusand how this transformation affects spatiotemporal demographic inference. Most previous iDDC studies assumed a linear relationship, but this study hypothesises that the relationship may be nonlinear, especially for high‐elevation species with broader environmental tolerances. We test two key hypotheses: (1) The habitat suitability to population size relationship is nonlinear forE. iheringii(high‐elevation) and linear forE. catenatus(low‐elevation); and (2)E. iheringiiexhibits higher effective migration across populations thanE. catenatus. Our findings provide clear support for hypothesis (2), but mixed support for hypothesis (1), with strong model support for a nonlinear transformation in the high‐elevationE. iheringiiand some (albeit weak) support for a nonlinear transformation also inE. catenatus. The iDDC models allow us to generate landscape‐wide maps of predicted genetic diversity for both species, revealing that genetic diversity predictions for the high‐elevationE. iheringiialign with estimated patterns of historical range stability, whereas predictions for low‐elevationE. catenatusare distinct from range‐wide stability predictions. This research highlights the importance of accurately modelling the habitat suitability to population size relationship in iDDC studies, contributing to our understanding of species' demographic responses to environmental changes. 
    more » « less
  3. ; ; ; ; ; ; (, BMC Evolutionary Biology)
    Abstract Background Disentangling the drivers of genetic differentiation is one of the cornerstones in evolution. This is because genetic diversity, and the way in which it is partitioned within and among populations across space, is an important asset for the ability of populations to adapt and persist in changing environments. We tested three major hypotheses accounting for genetic differentiation—isolation-by-distance (IBD), isolation-by-environment (IBE) and isolation-by-resistance (IBR)—in the annual plant Arabidopsis thaliana across the Iberian Peninsula, the region with the largest genomic diversity. To that end, we sampled, genotyped with genome-wide SNPs, and analyzed 1772 individuals from 278 populations distributed across the Iberian Peninsula. Results IBD, and to a lesser extent IBE, were the most important drivers of genetic differentiation in A. thaliana . In other words, dispersal limitation, genetic drift, and to a lesser extent local adaptation to environmental gradients, accounted for the within- and among-population distribution of genetic diversity. Analyses applied to the four Iberian genetic clusters, which represent the joint outcome of the long demographic and adaptive history of the species in the region, showed similar results except for one cluster, in which IBR (a function of landscape heterogeneity) was the most important driver of genetic differentiation. Using spatial hierarchical Bayesian models, we found that precipitation seasonality and topsoil pH chiefly accounted for the geographic distribution of genetic diversity in Iberian A. thaliana . Conclusions Overall, the interplay between the influence of precipitation seasonality on genetic diversity and the effect of restricted dispersal and genetic drift on genetic differentiation emerges as the major forces underlying the evolutionary trajectory of Iberian A. thaliana . 
    more » « less
  4. ; ; (, Ecology and Evolution)
    ABSTRACT Hybrid zones, where genetically distinct groups of organisms meet and interbreed, offer valuable insights into the nature of species and speciation. Here, we present a new R package,bgchm, for population genomic analyses of hybrid zones. This R package extends and updates the existingbgcsoftware and combines Bayesian analyses of hierarchical genomic clines with Bayesian methods for estimating hybrid indexes, interpopulation ancestry proportions, and geographic clines. Compared to existing software,bgchmoffers enhanced efficiency through Hamiltonian Monte Carlo sampling and the ability to work with genotype likelihoods combined with a hierarchical Bayesian approach, enabling inference for diverse types of genetic data sets. The package also facilitates the quantification of introgression patterns across genomes, which is crucial for understanding reproductive isolation and speciation genetics. We first describe the models underlyingbgchmand then provide an overview of the R package and illustrate its use through the analysis of simulated and empirical data sets. We show thatbgchmgenerates accurate estimates of model parameters under a variety of conditions, especially when the genetic loci analyzed are highly ancestry informative. This includes relatively robust estimates of genome‐wide variability in clines, which has not been the focus of previous models and methods. We also illustrate how both selection and genetic drift contribute to variability in introgression among loci and how additional information can be used to help distinguish these contributions. We conclude by describing the promises and limitations ofbgchm, comparingbgchmto other software for genomic cline analyses, and identifying areas for fruitful future development. 
    more » « less
  5. ; ; ; ; ; ; ; (, Evolutionary Applications)
    ABSTRACT Coastal organisms live in a dynamic environment where a myriad of environmental stressors, including climate change, ocean acidification, and human harvesting, act on variable spatio‐temporal scales. Each of these stressors may impose unique selective forces on a population, shaping a species' adaptive potential and its ability to persist under future climatic conditions. Genomic investigations of adaptive responses to environmental and anthropogenic disturbances remain rare, especially in marine systems. Here, we use whole genome sequencing data from the owl limpet,Lottia gigantea, and outlier detection methods to pinpoint signals of selection (1) across long‐standing environmental gradients spanning the species' distribution, (2) at the poleward edge of the species' range where it experienced a recent expansion, and (3) between sites vulnerable to or protected from human size‐selective harvesting within California. Loci associated with environmental gradients across the entire range show the strongest differentiation at the southern end of the species' range, potentially driven by adaptation to sea surface temperature and pH. Additionalad‐hocoutlier analyses revealed a distinct set of loci potentially under selection in the expanded range, with different functional roles than the range‐wide outliers. Despite demographic models suggesting that protection from harvesting has a positive impact on the abundance of large individuals, we did not find strong signals of selection or changes in genetic diversity between sites differing in harvesting vulnerability. Our findings suggest that range‐wide environmental selective signals established over longer time scales are distinct from those imposed by climatic anomalies at finer spatio‐temporal scales. We found that climatic variation has a stronger selective imprint than human harvesting, and thus conservation interventions should consider prioritizing the maintenance of climate‐related adaptive potential. Understanding how climatic trends and anomalies interact with anthropogenic pressures will allow us to make more informed decisions to sustain the evolutionary capacity ofL. giganteaand other key coastal species. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Text Encoded Conversion
  • XML
  • Save / Share this Record
  • Send to Email