Generalized gap model for bacterial artificial chromosome clone fingerprint mapping and shotgun sequencing

Area of Science:

• Genomics
• Bioinformatics
• Computational Biology

Background:

• The Lander-Waterman theory is foundational for analyzing genome mapping projects, particularly bacterial artificial chromosome (BAC) fingerprinting and shotgun sequencing.
• Characterizing gaps and coverage is crucial for assembling complete genomes, but existing models have limitations at higher coverages.

Purpose of the Study:

• To develop an extended Lander-Waterman theory applicable to a wider range of bacterial artificial chromosome (BAC) project parameters.
• To provide enhanced descriptive statistics for characterizing gaps and coverage in genome mapping projects.
• To address and resolve inconsistencies and paradoxes previously reported with the Lander-Waterman theory at higher coverages.

Main Methods:

• Developed a theoretical extension to the Lander-Waterman model.
• Analyzed the model's behavior across a broader spectrum of project parameters and coverage levels.
• Investigated the relationship between the extended model and other existing genome mapping models.
• Applied the extended theory to calculate closure probabilities and determine optimal stopping points in BAC shotgun sequencing.

Main Results:

• The extended Lander-Waterman theory supports a larger set of descriptive statistics and is applicable to more project parameters.
• Previous assertions of Lander-Waterman theory inconsistency at higher coverages were found to be incorrect.
• A well-known, seemingly different genome mapping model was identified as mathematically equivalent to the extended Lander-Waterman theory.
• The paradox of infinite island lengths in genome mapping was resolved.

Conclusions:

• The enhanced Lander-Waterman theory provides a more robust framework for analyzing bacterial artificial chromosome (BAC) fingerprinting and shotgun sequencing projects.
• The corrected understanding of the theory's behavior at higher coverages improves gap characterization and genome assembly strategies.
• Applications include improved probability density function evolution, closure probability calculations, and probabilistic stopping point determination for BAC sequencing.