Abstract:The coalescent is a powerful statistical framework that allows us to infer past population dynamics leveraging the ancestral relationships reconstructed from sampled molecular sequence data. In many biomedical applications, such as in the study of infectious diseases, cell development, and tumorgenesis, several distinct populations share evolutionary history and therefore become dependent. The inference of such dependence is a highly important, yet a challenging problem. With advances in sequencing technologies, we are well positioned to exploit the wealth of high-resolution biological data for tackling this problem. Here, we present adaPop , a probabilistic model to estimate past population dynamics of dependent populations and to quantify their degree of dependence. An essential feature of our approach is the ability to track the time-varying association between the populations while making minimal assumptions on their functional shapes via Markov random field priors. We provide nonparametric estimators, extensions of our base model that integrate multiple data sources, and fast scalable inference algorithms. We test our method using simulated data under various dependent population histories and demonstrate the utility of our model in shedding light on evolutionary histories of different variants of SARS-CoV-2. Genomic data provide information about evolutionary dynamics—such as an evolving epidemic and tumorgenesis—that is difficult to infer from other source of data. One of the main computational challenges in inferring past population histories is to jointly model dependent subpopulations and correctly quantifying their dependencies over time. When distinct subpopulations have common ancestry in the past and evolve under shared environmental pressure, their population dynamics become dependent. In this work, we propose an efficient inference method for studying dependent population dynamics from genetic data in the coalescent framework: an approach that considers the stochastic process of the "coalescence" of genealogical lineages traveling back in time to explain the statistical properties of a sample's genetic variation. We also extend our framework to jointly model the ancestral and sampling processes incorporating sampling times as an additional source of information. We validate our methods via extensive simulations and demonstrate that our methods provide new insights into the evolutionary dynamics of SARS-CoV-2 novel variants.

Minimal-assumption inference from population-genomic data

Decoding coalescent hidden Markov models in linear time

MAGIC: A diffusion-based imputation method reveals gene-gene interactions in single-cell RNA-sequencing data

Estimation of demography and mutation rates from one million haploid genomes

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

Limits and convergence properties of the sequentially Markovian coalescent

Exact Limits of Inference in Coalescent Models

Improved inference of population histories by integrating genomic and epigenomic data

A Principled Approach to Deriving Approximate Conditional Sampling Distributions in Population Genetics Models with Recombination

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

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

adaPop: Bayesian inference of dependent population dynamics in coalescent models

Coalescence 2.0: a multiple branching of recent theoretical developments and their applications

An empirical approach to demographic inference with genomic data

Inference of Locus-Specific Population Mixtures From Linked Genome-Wide Allele Frequencies

Accurate inference of population history in the presence of background selection

Fast and accurate haplotype inference with hidden markov model

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

Accelerated Bayesian inference of population size history from recombining sequence data

Ancestral inference from haplotypes and mutations

An Efficient Bayesian Inference Framework for Coalescent-Based Nonparametric Phylodynamics