European LD proxies lose tagging for 43% of GWAS leads in the African panel: a chromosome 22 analysis with SwiftLD, an 826-fold faster, byte-identical vcftools LD
Abstract
SwiftLD is a drop-in replacement for the pairwise linkage disequilibrium (LD) modes of vcftools 0.1.17, --hap-r2 and --geno-r2, which population-genetic pipelines use to compute r2, D and D' between nearby variants. On three regions of 1000 Genomes chromosome 22 with 503 to 2,504 samples it wrote the same output as vcftools, byte for byte, and was 826-fold faster (geometric mean at two threads on an AVX-512 CPU). It was also 7.6-fold faster than PLINK 2, the fastest alternative, from VCF input to LD output on all five regions tested; PLINK 2 writes a different file, although its r2 values agree with vcftools' within 10^-6. One megabase with 503 samples took vcftools 29 s and SwiftLD 0.049 s, and on an ordinary AVX2 laptop SwiftLD was 102-fold faster than vcftools at 12 threads. vcftools was slow because it repeatedly re-reads genotypes from a temporary file, counts haplotypes individually and flushes every output row. I built SwiftLD with correctness-gated coding agents (Claude Opus 5.5): every change had to keep the output byte-identical to vcftools. No output byte differed in 167,749 differential tests on three machines. On chromosome 22, the best European LD proxy of a GWAS lead kept r2 >= 0.8 for 57.4% of leads in the African panel and 85.6% to 92.7% in other panels (98.1% in a European split-half control). Code and a two-command reproduction: https://github.com/joe-b-christensen/swiftld.
02 /Specification & Classification
Click any tag to browse matching papers- Research Focus
- Approach
Reproducibility ledger
Nobody has reported trying to reproduce this paper yet. Did you replicate an experiment or run the code? Your report saves the next lab a week.
Sign in to file a reproducibility report.
Related papers
- emeraLD: Rapid Linkage Disequilibrium Estimation with Massive Data Sets
Corbin Quick, Christian Fuchsberger, … Hyun Min Kang·BioinformaticsPeer-reviewed·Jun 28, 2018
- ARCHER-LD: Rapid Long-Range Linkage Disequilibrium Calculations at Biobank Scale using GPU Acceleration
Rachit Kumar, Pankhuri Singhal, … Anurag Verma·bioRxiv·Sep 21, 2026
- vcf2gwas - Python API for comprehensive GWAS analysis using GEMMA
Frank Vogt, Gautam S. Shirsekar, Detlef Weigel·BioinformaticsPeer-reviewed·Oct 8, 2021
- LDBlockShow: a fast and convenient tool for visualizing linkage disequilibrium and haplotype blocks based on variant call format files
Shan‐Shan Dong, Weiming He, … Tie‐Lin Yang·Briefings in BioinformaticsPeer-reviewed·Aug 24, 2020
- SAFE-LD: A novel method for the estimation of linkage disequilibrium from summary statistics
Giulia Elizabeth de Sanctis, Sodbo Sharapov, … Nicola Pirastu·bioRxiv·Sep 29, 2025
- Approximately independent linkage disequilibrium blocks in human populations
Tomaz Berisa, Joseph K. Pickrell·BioinformaticsPeer-reviewed·Sep 22, 2015