PyBSASeq

PyBSASeq identifies likely trait-associated single nucleotide polymorphisms (ltaSNPs) from Bulked Segregant Analysis coupled with Next Generation Sequencing (BSA-Seq) data to map genomic regions conditioning traits while reducing required sequencing coverage.

Key Features:

• Algorithmic Efficiency: Implemented in Python, PyBSASeq uses a novel algorithm applying Fisher's exact test to identify likely trait-associated SNPs (ltaSNPs) and localize associated genomic regions.
• Cost-Effectiveness: Demonstrates over five-fold higher sensitivity compared to SNP index and G-statistic approaches, enabling approximately an 80% reduction in sequencing coverage and costs.
• Reduced Coverage Requirement: Requires significantly lower sequencing coverage than SNP index and G-statistic methods, facilitating BSA-Seq analysis in species with large genomes.

Scientific Applications:

• DNA polymorphism–trait association mapping: Rapid identification of DNA polymorphism–trait associations from BSA-Seq data.
• QTL mapping: Mapping of qualitative and quantitative trait loci (QTL) using SNP frequency differences between bulks.
• Large-genome and budget-constrained studies: Application to species with large genomes or projects constrained by sequencing budget where lower coverage is advantageous.

Methodology:

PyBSASeq identifies ltaSNPs using Fisher's exact test and computes the ratio of ltaSNPs to total SNPs within chromosomal intervals to detect genomic regions associated with traits.