Abstract:BackgroundRecent advances in sequencing technologies have allowed for the retrieval of ancient DNA data (aDNA) from skeletal remains, providing direct genetic snapshots from diverse periods of human prehistory. Comparing samples taken in the same region but at different times, hereafter called “serial samples”, may indicate whether there is continuity in the peopling history of that area or whether an immigration of a genetically different population has occurred between the two sampling times. However, the exploration of genetic relationships between serial samples generally ignores their geographical locations and the spatiotemporal dynamics of populations. Here, we present a new coalescent-based, spatially explicit modelling approach to investigate population continuity using aDNA, which includes two fundamental elements neglected in previous methods: population structure and migration. The approach also considers the extensive temporal and geographical variance that is commonly found in aDNA population samples.ResultsWe first showed that our spatially explicit approach is more conservative than the previous (panmictic) approach and should be preferred to test for population continuity, especially when small and isolated populations are considered. We then applied our method to two mitochondrial datasets from Germany and France, both including modern and ancient lineages dating from the early Neolithic. The results clearly reject population continuity for the maternal line over the last 7500 years for the German dataset but not for the French dataset, suggesting regional heterogeneity in post-Neolithic migratory processes.ConclusionsHere, we demonstrate the benefits of using a spatially explicit method when investigating population continuity with aDNA. It constitutes an improvement over panmictic methods by considering the spatiotemporal dynamics of genetic lineages and the precise location of ancient samples. The method can be used to investigate population continuity between any pair of serial samples (ancient-ancient or ancient-modern) and to investigate more complex evolutionary scenarios. Although we based our study on mitochondrial DNA sequences, diploid molecular markers of different types (DNA, SNP, STR) can also be simulated with our approach. It thus constitutes a promising tool for the analysis of the numerous aDNA datasets being produced, including genome wide data, in humans but also in many other species.

Improved inference of population histories by integrating genomic and epigenomic data

Limits and convergence properties of the sequentially Markovian coalescent

Accurate inference of population history in the presence of background selection

Inferring demographic and selective histories from population genomic data using a two-step approach in species with coding-sparse genomes: an application to human data

Esmc: a Statistical Model to Infer Admixture Events from Individual Genomics Data

Inferring the joint demographic history of multiple populations from multidimensional SNP frequency data

Inferring demographic history from a spectrum of shared haplotype lengths

Investigating population continuity with ancient DNA under a spatially explicit simulation framework

Extending the IICR to multiple genomes and identification of limitations of some demographic inferential methods

Simultaneous Inference of Past Demography and Selection from the Ancestral Recombination Graph under the Beta Coalescent

Demographic inference using genetic data from a single individual: Separating population size variation from population structure

The interplay of recombination landscape, a transposable element and population history in European populations of

Efficient moment-based inference of admixture parameters and sources of gene flow

Decoding coalescent hidden Markov models in linear time

Methylomes reveal recent evolutionary changes in populations of two plant species

Faster inference of complex demographic models from large allele frequency spectra

Statistical inference of transmission fidelity of DNA methylation patterns over somatic cell divisions in mammals

Evaluation of ancient DNA imputation: a simulation study

Integrating Pool-seq uncertainties into demographic inference

A Fast and Flexible Statistical Model for Large-Scale Population Genotype Data: Applications to Inferring Missing Genotypes and Haplotypic Phase

Ignoring population structure in hominin evolutionary models can lead to the inference of spurious admixture events