LDkit

LDkit performs linkage disequilibrium analysis to calculate LD decay, identify LD blocks, and assess linkage between specific sites and regions for applications in genetic association, population genetics, and evolutionary studies.

Key Features:

• Parallel Computing Support: Leverages parallel computing by distributing tasks across multiple threads; for example, analysis of chromosome 22 in the Human 1000 Genomes dataset (1,103,547 SNPs, 2,504 individuals) reduced runtime from ~77 minutes to under 6 minutes using 32 threads.
• Data Format Compatibility: Supports Variant Call Format (VCF) and PLINK 'ped + map' formats.
• Subset Analysis Capability: Enables analysis of specific subsets of individuals drawn from larger populations.
• Comprehensive LD Analysis Functions: Calculates and visualizes LD decay, LD blocks, and linkage between specific sites and regions.
• Implementation / Cross-Platform Support: Implemented in Java to allow execution across different operating systems.

Scientific Applications:

• Genetic Association Studies: Uses LD patterns to help identify genes linked to traits or diseases.
• Population Genetics: Analyzes patterns of genetic variation to infer population structure and demographic history.
• Evolutionary Biology: Examines changes in linkage disequilibrium to investigate evolutionary processes over time.

Methodology:

Parses VCF and PLINK 'ped + map' inputs, distributes computations across multiple threads for parallel processing, computes LD metrics including LD decay and LD blocks, analyzes linkage between sites and regions, performs subset analyses, and produces visualizations; implemented in Java.