Abstract:Reconstructing the demographic history of populations and species is one of the greatest challenges facing population geneticists. Mazet et al (2016) introduced, for a sample of size k=2 haploid genomes, a time- and sample-dependent parameter which they called the IICR (inverse instantaneous coalescence rate). Here we extend their work to larger sample sizes and focus on Tk, the time to the first coalescence event in a haploid sample of size k where k>=2. We define the IICRk as the Inverse Instantaneous Coalescence Rate among k lineages. We show that (i) under a panmictic population k(k-1)/2 * IICRk is equivalent to Ne, (ii) the IICRk can be obtained by either simulating Tk values or by using the Q-matrix approach of Rodriguez et al (2018) and we provide the corresponding Python and R scripts. We then study the properties of the k(k-1)/2 * IICRk under a limited set of n-island and stepping-stone models. We show that (iii) in structured models the k(k-1)/2 * IICRk is dependent on the sample size and on the sampling scheme, even when the genomes are sampled in the same deme. For instance, we find that k(k-1)/2 * IICRk plots for individuals sampled in the same deme will be shifted towards recent times with a lower plateau as k increases. We thus show that (iv) the k(k-1)/2 * IICRk cannot be used to represent ``the demographic history'' in a general sense, (v) the k(k-1)/2 * IICRk can be estimated from real or simulated genomic data using the PSMC/MSMC methods (Li and Durbin, 2011; Schiffels and Durbin, 2013) (vi) the MSMC2 method produces smoother curves that infer something that is not the k(k-1)/2 * IICRk, but are close to the IICR in the recent past when all samples are obtained from the same deme. Altogether we argue that the PSMC, MSMC and MSMC2 plots are not expected to be identical even when the genomes are sampled from the same deme, that none can be said to represent the ``demographic history of populations'' and that they should be interpreted with care. We suggest that the PSMC, MSMC and MSMC2 could be used together with the k(k-1)/2 * IICRk to identify the signature of population structure, and to develop new strategies for model choice.

Inferring the demographic history from DNA sequences: An importance sampling approach based on non-homogeneous processes

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

Comparison of Bayesian Coalescent Skyline Plot Models for Inferring Demographic Histories

Accurate inference of population history in the presence of background selection

Geometry of the sample frequency spectrum and the perils of demographic inference

Improved inference of population histories by integrating genomic and epigenomic data

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

Sequential importance sampling for multi-resolution Kingman-Tajima coalescent counting

Sampling through time and phylodynamic inference with coalescent and birth-death models

Limits and convergence properties of the sequentially Markovian coalescent

Accurate and flexible estimation of effective population size history

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

Accelerated Bayesian inference of population size history from recombining sequence data

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

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

The Confounding Effects of Population Structure, Genetic Diversity and the Sampling Scheme on the Detection and Quantification of Population Size Changes

Ancestral inference from haplotypes and mutations

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

Stairway Plot 2: demographic history inference with folded SNP frequency spectra

Integrating Pool-seq uncertainties into demographic inference

The Expected Sample Allele Frequencies from Populations of Changing Size via Orthogonal Polynomials