AI Insight
Researchers developed ARCHER-LD, a GPU-accelerated computational tool that dramatically speeds up the calculation of linkage disequilibrium (LD) patterns across genetic variants in large-scale genomic datasets. The method achieves 8-28x faster performance compared to existing CPU-based tools like PLINK, completing analysis of 120 million variants from whole-genome sequencing data in under 6 hours using supercomputing resources. This enables researchers to calculate LD directly from their own biobank-scale datasets rather than relying on external reference panels.
Why it matters
This tool allows researchers to perform more accurate, population-specific genetic analyses by using LD patterns from their actual study populations rather than potentially mismatched reference databases. This is particularly important for improving genetic research in populations currently underrepresented in existing genomic databases, potentially reducing health disparities in precision medicine applications.
Understand the Science
⚠️ Preprint – Noch nicht peer-reviewed
Dieser Artikel wurde noch nicht von unabhängigen Experten begutachtet. Die Ergebnisse sind vorläufig und sollten mit Vorsicht interpretiert werden.
Linkage disequilibrium (LD) information from one’s own dataset is considered optimal for downstream analyses such as statistical fine-mapping. However, the computational complexity ({approx}N2 / 2 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 {approx}8x speedup. We further demonstrate true biobank-scale performance in the Penn Medicine Biobank (PMBB; 57,170 samples), computing chromosome 1 LD ({approx}1.38 million variants) in 49 minutes versus 22.7 hours for PLINK, a {approx}28x speedup. With this method, we successfully computed on the aforementioned 30x WGS 1000 Genomes dataset ({approx}120 million variants and {approx}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.