Computational population genetics, clearly explained

Recently, in April 2026, new AADR versions were released on Harvard Dataverse. Among the more important additions are the new compatibility datasets introduced for reducing platform-specific bias when co-analyzing ancient DNA generated with different experimental setups. In practice, this matters when combining data produced with different capture reagents such as Agilent (AG), Twist (TW), and shotgun (SG), because these can introduce systematic differences that may affect downstream population genetic analyses. The compatibility panels were added to minimize that problem and make mixed-platform datasets more directly comparable. ...

I recently published dt, a modern data transformation tool designed to make the awk workflows commonly used on this blog more intuitive, expressive, and fast. Dt is written in Rust because it compiles to a single binary that runs anywhere, and it uses Polars for the actual data processing, giving you columnar operations that handle large files efficiently. The syntax uses explicit functions (filter(), select(), mutate()) that you chain together with pipes, so transformations read like a recipe instead of a regex puzzle. There’s also an interactive REPL that shows you the result after each operation, letting you build complex pipelines step-by-step, catch mistakes early, and undo errors with .undo. ...

In this post, I’ll cover how to build a transparent, deterministic admixture modeler from scratch using R. The main advantage of this approach is automated combinatorial testing: instead of manually specifying and running each potential ancestry model one-by-one in qpAdm, this script systematically tests hundreds or even thousands of source combinations in seconds and ranks them by fit quality. While qpAdm is the “industry standard” and provides sophisticated covariance weighting through block jackknife procedures, it requires manual intervention for each model specification and can often feel like a “black box.” This script offers a transparent, complementary approach that uses the exact same mathematical foundation (f4-statistic regression with quadratic programming constraints) to rapidly explore the models. If you have 10 potential sources (or even more) and want to test all 2-way, 3-way, and 4-way combinations, that’s 375 separate qpAdm runs against a single automated run with this script. ...

mergeit is part of the EIGENSOFT package and can be used to merge exactly two EIGENSTRAT datasets without converting to PACKEDPED format first. In this post, I’ll show how to merge the sample we created in Pseudohaploid Genotyping for Ancient DNA: BAM to EIGENSTRAT with the AADR dataset. Setting up EIGENSOFT mergeit is part of the EIGENSOFT package. You can install it via conda: conda install -c bioconda eigensoft If you haven’t installed conda yet, see the Miniconda setup in the pseudohaploid genotyping post. ...

This is a follow-up to my previous post Processing Ancient DNA: From FASTQ to Aligned BAM, where I aligned an ancient DNA sample against the hs37d5 reference genome, producing a filtered BAM compatible with the AADR dataset. In this post, I’ll cover pseudohaploid genotype calling using pileupCaller and converting the output to EIGENSTRAT format for use with ADMIXTOOLS. Since we just created this BAM ourselves in the previous post, we already know it’s aligned to hs37d5. However, if you’re starting with a BAM file, you’ll need to verify the reference genome first. I’ll start by showing how to check BAM headers to identify the reference genome. ...