ARCHER-LD: Rapid Long-Range Linkage Disequilibrium Calculations at Biobank Scale using GPU Acceleration

Clicks: 1
ID: 332846
2026
Free PDF
Article Quality & Performance Metrics
Overall Quality
Not rated
Combines reader engagement with the AI quality analysis. This article has not been analysed, so there is no overall score — reader engagement is measured and shown alongside.
AI Quality Assessment
Not analyzed
Readership in this journal

Ranked #955 of 957 articles by views in BMC Bioinformatics

Most read Least read

Bar heights use a square-root scale. Only the 120 most-read articles are drawn; the journal has 957 in total.

Mint this article as an NFT
Not yet minted

Create a permanent, verifiable on-chain record of this article on the Scimatic Network. The NFT is held in your Journament account, and you can withdraw it to your own wallet at any time.

5 SUSD one-off · no wallet required
Sign in to mint
Abstract
Abstract Summary Linkage disequilibrium (LD) information from one’s own dataset is considered optimal for downstream analyses such as statistical fine-mapping. However, the computational complexity (≈N22 computations for N variants) lead most studies to use external reference panels, such as from 1000 Genomes. To capture LD across all variant pairs in biobank-scale whole-genome sequencing (WGS) datasets with hundreds of millions of variants, new computational strategies are essential. We present a novel approach that uses multi-GPU distributed computing to compute R2 for every variant pair in a dataset. On chromosome 22 (1.8 million variants) of the 30x WGS 1000 Genomes dataset, our method using eight consumer-level 12GB GPUs (NVIDIA RTX 2080TIs) takes <10 minutes, while the same calculation with PLINK using a high-end 64-threaded CPU (Intel Xeon Gold 6338) takes > 80 minutes, corresponding to a ≈8x speedup. We further demonstrate true biobank-scale performance in the Penn Medicine Biobank (PMBB; 57,170 samples), computing chromosome 1 LD (≈1.38 million variants) in 49 minutes versus 22.7 hours for PLINK, a ≈28x speedup. With this method, we successfully computed on the aforementioned 30x WGS 1000 Genomes dataset (≈120 million variants and ≈2500 samples) the entire LD matrix ($>$1e16, or 10 quadrillion elements) in under 6 hours using 512 NVIDIA 40GB A100 GPUs on the Department of Energy Argonne Leadership Computing Facility Polaris Supercomputer. We make this tool, coded in Python using CuPy, publicly available. Using this tool, researchers can leverage the full extent of their genomic data without relying on external LD reference panels and acquire more accurate, population-specific findings, particularly for groups underrepresented in existing databases. Availability and Implementation ARCHER-LD is an open-source tool available on GitHub under the Apache License 2.0 at https://github.com/rachitk/ARCHER-LD. The version of record used in this manuscript is archived via Zenodo at https://doi.org/10.5281/zenodo.23399887. Supplementary Information Supplementary material is available at Bioinformatics online.
Reference Key
openalex_W7222915871 Use this key to autocite in the manuscript while using SciMatic Manuscript Manager or Thesis Manager
Authors Rachit Kumar, Pankhuri Singhal, David Zhang, Christopher C. Carson, Mitchell Conery, Álex Rodríguez, Tarak Nandi, Marylyn DeRiggi Ritchie, Mathialakan Thavappiragasam, Benjamin F. Voight, K Madduri Ravi, Anurag Verma
Journal BMC Bioinformatics
Year 2026
DOI
10.1093/bioinformatics/btag783
URL
https://doi.org/10.1093/bioinformatics/btag783
Keywords Keywords not found

Citations

No citations found. To add a citation, contact the admin at info@scimatic.org

Comments

No comments yet. Be the first to comment on this article.