Alignment-free population genomics: an efficient estimator of sequence diversity

Area of Science:

• Genomics
• Bioinformatics
• Computational Biology

Background:

• Comparative sequencing is crucial for genome functional annotation.
• Efficient computational tools are needed for analyzing vast amounts of genomic data.

Purpose of the Study:

• To present a novel strategy for comparing unaligned genomes.
• To develop an efficient and accurate estimator for the average number of mismatches per site (π).

Main Methods:

• Utilizing a coalescent approach combined with advanced sequence indexing algorithms.
• Implementing a maximum-likelihood estimator for π, improving upon previous alignment-free methods.
• Incorporating fluctuating coalescent times to enhance estimator accuracy.

Main Results:

• The new estimator demonstrates speed and accuracy, even with moderate recombination rates (ρ ≤ π).
• The method is effective for analyzing large genomes, with run times ideally growing linearly with sequence length.
• Sliding window analysis of Drosophila genomes identified a divergence minimum in the chromosome 3 pericentromeric region.

Conclusions:

• The developed computational strategy and estimator are valuable for comparative genomics.
• This approach facilitates the analysis of unaligned genomic sequences, particularly large genomes.
• The findings provide insights into genome divergence patterns in closely related species.