Identifying genomic adaptation to local climate using a mechanistic evolutionary model
Nikunj Goel,Christen M. Bossu,Justin J. Van Ee,Erika Zavaleta,Kristen C. Ruegg,Mevin B. Hooten
DOI: https://doi.org/10.1101/2024.11.21.624744
2024-11-22
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.
Biology