Inference of population genetic structure from temporal samples of DNA
Francois, O.; Demaille, B.; Severine, L.; Jay, F.
Show abstract
The recent years have seen a growing number of studies investigating evolutionary questions using ancient DNA techniques and temporal samples of DNA. To address these questions, one of the most frequently-used algorithm is based on principal component analysis (PCA). When PCA is applied to temporal samples, the sample dates are, however, ignored during analysis, which could lead to some misinterpretations of the results. Here we introduce a new factor analysis (FA) method for which individual scores are corrected for the effect of allele frequency drift through time. Based on a diffusion approximation, our approach approximates allele frequency drift in a random mating population by a Brownian process. Exact solutions for estimates of corrected factors are obtained, and a fast estimation algorithm is presented. We compared data representations obtained from the FA method with PCA and with PC projections in simulations of divergence and admixture scenarios. Then we applied FA with correction for temporal drift to study the evolution of hepatitis C virus in a patient infected by multiple strains, and to describe the population structure of ancient European samples.
Matching journals
The top 5 journals account for 50% of the predicted probability mass.
Similar papers in this journal
Similar papers in this journal
- Mutation rate, selection, and epistasis inferred from RNA virus haplotypes via neural posterior estimation 94%
- Impact and mitigation of sampling bias to determine viral spread: evaluating discrete phylogeography through CTMC modeling and structured coalescent model approximations 93%
- Global estimates of the fitness advantage of SARS-CoV-2 variant Omicron 93%
Similar papers in this journal
Similar papers in this journal
"Similar papers" are the closest papers from that journal in the model's embedding space. They show what the match is built on, but the ranking comes mostly from a classifier over the whole training set, not from these examples alone.