Estimating Temporally Variable Selection Intensity from Ancient DNA Data.

Novel technologies for recovering DNA information from archaeological and historical specimens have made available an ever-increasing amount of temporally spaced genetic samples from natural populations. These genetic time series permit the direct assessment of patterns of temporal changes in allele frequencies and hold the promise of improving power for the inference of selection. Increased time resolution can further facilitate testing hypotheses regarding the drivers of past selection events such as the incidence of plant and animal domestication. However, studying past selection processes through ancient DNA (aDNA) still involves considerable obstacles such as postmortem damage, high fragmentation, low coverage, and small samples. To circumvent these challenges, we introduce a novel Bayesian framework for the inference of temporally variable selection based on genotype likelihoods instead of allele frequencies, thereby enabling us to model sample uncertainties resulting from the damage and fragmentation of aDNA molecules. Also, our approach permits the reconstruction of the underlying allele frequency trajectories of the population through time, which allows for a better understanding of the drivers of selection. We evaluate its performance through extensive simulations and demonstrate its utility with an application to the ancient horse samples genotyped at the loci for coat coloration. Our results reveal that incorporating sample uncertainties can further improve the inference of selection.

Estimating temporally variable selection intensity from ancient DNA data with the flexibility of modelling linkage and epistasis.

Innovations in ancient DNA (aDNA) preparation and sequencing technologies have exponentially increased the quality and quantity of aDNA data extracted from ancient biological materials. The additional temporal component from the incoming aDNA data can provide improved power to address fundamental evolutionary questions like characterizing selection processes that shape the phenotypes and genotypes of contemporary populations or species. However, utilizing aDNA to study past selection processes still involves considerable hurdles like how to eliminate the confounding factor of genetic interactions in the inference of selection. To address this issue, we extend the approach of He et al., 2023 to infer temporally variable selection from the aDNA data in the form of genotype likelihoods with the flexibility of modelling linkage and epistasis in this work. Our posterior computation is carried out by a robust adaptive version of the particle marginal Metropolis-Hastings algorithm with a coerced acceptance rate. Our extension inherits the desirable features of He et al., 2023 such as modelling sample uncertainty resulting from the damage and fragmentation of aDNA molecules and reconstructing underlying gamete frequency trajectories of the population. We evaluate its performance through extensive simulations and show its utility with an application to the aDNA data from pigmentation loci in horses.