AluMine

AluMine identifies and genotypes polymorphic Alu element insertions in sequencing data using an alignment-free 32-mer k-mer frequency approach.

Key Features:

• Alignment-free 32-mer k-mer frequency analysis: Counts frequencies of 32-mer sequences (k-mers) from raw sequencing reads rather than performing sequence alignment.
• REF-plus Discovery: Identifies k-mers present in the sample but absent from the reference genome to detect candidate polymorphic insertions.
• REF-minus Discovery: Detects k-mers present in the reference but missing from the sample to indicate possible deletions or absence variants.
• Merging and filtering workflow: Integrates REF-plus and REF-minus results to produce a consolidated list of candidate polymorphic insertions.
• Genotyping with 32-mer pairs: Uses insertion/deletion-specific 32-mer pairs to call genotypes directly from sequencing reads without alignment.
• Genotyping concordance: Reports high concordance with experimental genotypes (98.7%).
• Performance: Discovery pipeline runtime reported as approximately 2 hours per individual and genotyping runtime between 0.4 and 4 hours per individual depending on hardware.
• Implementation languages: Analysis scripts implemented in PERL and bash.

Scientific Applications:

• Discovery of polymorphic Alu insertions: Detection and genotyping of novel and known polymorphic Alu element insertions from sequencing data for studies of genomic variation and genetic diversity.
• Disease association studies: Generation of Alu insertion genotypes to support investigations of potential links between Alu variation and disease phenotypes.
• Population genetics and evolutionary analysis: Analysis of Alu insertion patterns across samples and populations to study population structure and evolutionary dynamics.

Methodology:

Counts k-mer (32-mer) frequencies from raw sequencing reads; performs REF-plus (sample-present/reference-absent) and REF-minus (reference-present/sample-absent) comparisons; merges and filters discovery results; genotypes loci using insertion/deletion-specific 32-mer pairs to call genotypes directly from sequencing reads without alignment.