Abstract:Admixture graphs are mathematical structures that describe the ancestry of populations in terms of divergence and merging (admixing) of ancestral populations as a graph. An admixture graph consists of a graph topology, branch lengths, and admixture proportions. The branch lengths and admixture proportions can be estimated using numerous numerical optimization methods, but inferring the topology involves a combinatorial search for which no polynomial algorithm is known. In this paper, we present a reversible jump MCMC algorithm for sampling high-probability admixture graphs and show that this approach works well both as a heuristic search for a single best-fitting graph and for summarizing shared features extracted from posterior samples of graphs. We apply the method to 11 Native American and Siberian populations and exploit the shared structure of high-probability graphs to characterize the relationship between Saqqaq, Inuit, Koryaks, and Athabascans. Our analyses show that the Saqqaq is not a good proxy for the previously identified gene flow from Arctic people into the Na-Dene speaking Athabascans. One way of summarizing historical relationships between genetic samples is by constructing an admixture graph. An admixture graph describes the demographic history of a set of populations as a directed acyclic graph representing population splits and mergers. The greedy search algorithms that are typically used to infer admixture graphs may fail to find the globally optimal graph. We here improve on these approaches by developing a novel MCMC sampling method, AdmixtureBayes , that can sample from the posterior distribution of admixture graphs. This enables an effective search of the entire state space as well as the ability to report a level of confidence in the sampled graphs. We apply AdmixtureBayes to a set of Native American and Arctic genomes to reconstruct the demographic history of these populations and report posterior probabilities of specific admixture events. While some previous studies have identified the ancient Saqqaq culture as a source of introgression into Athabascans, we instead find that it is the Siberian Koryak population, not the Saqqaq, that serves as the best proxy for gene flow into Athabascans.

A central limit theorem concerning uncertainty in estimates of individual admixture

Revealing the range of maximum likelihood estimates in the admixture model.

Understanding admixture fractions: theory and estimation of gene-flow

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

Interpreting forensic DNA mixtures: allowing for uncertainty in population substructure and dependence

Association Studies with Imputed Variants Using Expectation-Maximization Likelihood-Ratio Tests

An unconditional exact test for the Hardy-Weinberg equilibrium law: sample-space ordering using the Bayes factor.

An Empirical Bayes Approach for Constructing the Confidence Intervals of Clonality and Entropy

Sensitivity of inferences in forensic genetics to assumptions about founding genes

MaCH-admix: Genotype Imputation for Admixed Populations.

A Comparison of Approaches to Account for Uncertainty in Analysis of Imputed Genotypes

Modeling population structure under hierarchical Dirichlet processes

Population Admixture: Detection by Hardy-Weinberg Test and Its Quantitative Effects on Linkage-Disequilibrium Methods for Localizing Genes Underlying Complex Traits.

Faster model-based estimation of ancestry proportions

Posterior predictive checks to quantify lack-of-fit in admixture models of latent population structure

Confidence intervals for nonhomogeneous branching processes and polymerase chain reactions

Assessing the limits of local ancestry inference from small reference panels

Bayesian inference of admixture graphs on Native American and Arctic populations

Needles in the Haystack: Identifying Individuals Present in Pooled Genomic Data

Approximations of the allelic frequency spectrum in general supercritical branching populations

Quantifying uncertainty of molecular mismatch introduced by mislabeled ancestry using haplotype-based HLA genotype imputation