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.

Accurate detection of identity-by-descent segments in human ancient DNA

Estimating effective population size trajectories from time-series Identity-by-Descent (IBD) segments

Characterizing identity by descent segments in Chinese interpopulation unrelated individual pairs

Personalized genealogical history of UK individuals inferred from biobank-scale IBD segments

Haplotype-based inference of recent effective population size in modern and ancient DNA samples

Inferring demographic history from a spectrum of shared haplotype lengths

Placing Ancient DNA Sequences into Reference Phylogenies

Detecting population continuity and ghost admixture among ancient genomes

Parental relatedness through time revealed by runs of homozygosity in ancient DNA

Reconstruct recent multi-population migration history by using identical-by-descent sharing

Accurate genome-wide phasing from IBD data

Conflation of short identity-by-descent segments bias their inferred length distribution

KIN: a method to infer relatedness from low-coverage ancient DNA

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

A Probabilistic Method for Estimating the Sharing of Identity by Descent for Populations with Migration

Identifying and Interpreting Apparent Neanderthal Ancestry in African Individuals

A rapid, accurate approach to inferring pedigrees in endogamous populations

Rapid detection of identity-by-descent tracts for mega-scale datasets

New Insights Into Mitochondrial DNA Reconstruction and Variant Detection in Ancient Samples

Inferring Admixture Histories of Human Populations Using Linkage Disequilibrium

The length of haplotype blocks and signals of structural variation in reconstructed genealogies