A path integral approach for allele frequency dynamics under polygenic selection

Nathan W. Anderson,Lloyd Kirk,Joshua G. Schraiber,Aaron P. Ragsdale

DOI: https://doi.org/10.1101/2024.06.14.599114

2024-06-14

Abstract:Many phenotypic traits have a polygenic genetic basis, making it challenging to learn their genetic architectures and predict individual phenotypes. One promising avenue to resolve the genetic basis of complex traits is through evolve-and-resequence experiments, in which laboratory populations are exposed to some selective pressure and trait-contributing loci are identified by extreme frequency changes over the course of the experiment. However, small laboratory populations will experience substantial random genetic drift, and it is difficult to determine whether selection played a roll in a given allele frequency change. Predicting how much allele frequencies change under drift and selection had remained an open problem well into the 21st century, even those contributing to simple, monogenic traits. Recently, there have been efforts to apply the path integral, a method borrowed from physics, to solve this problem. So far, this approach has been limited to genic selection, and is therefore inadequate to capture the complexity of quantitative, highly polygenic traits that are commonly studied. Here we extend one of these path integral methods, the perturbation approximation, to selection scenarios that are of interest to quantitative genetics. In particular, we derive analytic expressions for the transition probability (i.e., the probability that an allele will change in frequency from , to in time ) of an allele contributing to a trait subject to stabilizing selection, as well as that of an allele contributing to a trait rapidly adapting to a new phenotypic optimum. We use these expressions to characterize the use of allele frequency change to test for selection, as well as explore optimal design choices for evolve-and-resequence experiments to uncover the genetic architecture of polygenic traits under selection.

Evolutionary Biology

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the dynamics of allele frequencies under polygenic selection. Specifically, the author focuses on how to identify gene loci that contribute to specific phenotypes through changes in allele frequencies in laboratory populations when these populations are exposed to certain selection pressures. However, because laboratory populations are small in size and will experience significant random genetic drift, it is difficult to determine whether a change in a certain allele frequency is caused by selection. Predicting changes in allele frequencies under the combined action of genetic drift and selection has been an unsolved problem for a long time, especially when it comes to complex, highly polygenic traits.

To this end, the author extended a path - integral - based method - the perturbation approximation method - and applied it to polygenic selection scenarios of interest in quantitative genetics. They derived analytical expressions for the transition probability of allele frequencies changing from \(x\) to \(y\) under the conditions of stabilizing selection and rapid adaptation to new phenotypic optima. These expressions are used to test whether allele - frequency changes are caused by selection and to explore the optimal design choices for evolution and re - sequencing experiments in order to reveal the genetic structure of polygenic traits under selection.

In short, this study aims to provide a theoretical framework to understand and predict the pattern of allele - frequency changes in the context of polygenic selection through analytical methods, thereby helping researchers more effectively identify gene loci affected by selection. This not only helps to understand the genetic basis of adaptive evolution, but also provides theoretical guidance for designing more effective experiments.