Abstract:1. Identifying genomic adaptation is key to understanding species' evolutionary responses to climate change. However, current methods to identify adaptive variation have two major limitations. First, when estimating genetic variation, most methods do not account for observational uncertainty in genetic data because of finite sampling and missing genotypes. Second, current methods use phenomenological models to partition genetic variation into adaptive and non-adaptive components. These phenomenological models are not mechanistic models of evolution and, therefore, do not faithfully capture the demographic history of the species. 2. We address these limitations by developing a hierarchical Bayesian model that explicitly accounts for both the observational uncertainty and underlying evolutionary processes. The first layer of the hierarchy is the data model that captures observational uncertainty by probabilistically linking RAD-sequence data to genetic variation. The second layer is a process model that provides a mechanistic explanation of how evolutionary forces, such as local adaptation, mutation, migration, and drift, maintain genetic variation. The third layer is the parameter model, which incorporates our knowledge about biological processes. For example, because most loci in the genome are expected to be neutral, the environmental sensitivity coefficients are assigned a regularized prior centered at zero. Together, the three models provide a rigorous probabilistic framework to identify local adaptation in wild organisms. 3. Analysis of simulated RAD-seq data shows that our statistical model can reliably infer adaptive genetic variation. To show the real-world applicability of our method, we re-analyzed RAD-seq data from Willow Flycatchers (Empidonax traillii) in the USA. We found 30 genes close to loci that showed a statistically significant association with temperature seasonality. Gene ontology suggests that several of these genes play a crucial role in egg mineralization, feather development, and the ability to withstand extreme temperatures. 4. Moreover, biogeographers can easily modify the data and process models to accommodate a wide range of genetic datasets (e.g., pool and low coverage genome sequencing) and demographic histories (e.g., range shifts), allowing them to construct statistical models specific to their study system.

On the plains and prairies of Minnesota: The role of mathematical statistics in biological explanation

A Unifying Model for the Analysis of Phenotypic, Genetic and Geographic Data

Mathematical modeling of evolution. Solved and open problems

A practical guide to selecting models for exploration, inference, and prediction in ecology

Random walk models in the life sciences: including births, deaths and local interactions

Local adaptation and genetic effects on fitness: Calculations for exponential family models with random effects

Unifying approaches from statistical genetics and phylogenetics for mapping phenotypes in structured populations

The Lost Art of Mathematical Modelling

Alternative study designs and nonparametric statistical methods for adaptive management studies of invasive plants

From Gregor Mendel to Eric Davidson: Mathematical Models and Basic Principles in Biology

Charles Darwin as a statistical thinker

Improved state-space models for inference about spatial and temporal variation in abundance from count data

Of mice and men: Sparse statistical modeling in cardiovascular genomics

A mathematical approach to territorial pattern formation

Social Network Analysis of a Grassland Rodent Community Using a Lotka-Volterra Modeling Approach

Statistical and Computational Challenges in Whole Genome Prediction and Genome-Wide Association Analyses for Plant and Animal Breeding

Open Problems in Mathematical Biology

A spatial approach to jointly estimate Wright's neighborhood size and long-term effective population size

Identifying genomic adaptation to local climate using a mechanistic evolutionary model

Algorithmic and Statistical Perspectives on Large-Scale Data Analysis