PAPNC, a novel method to calculate nucleotide diversity from large scale next generation sequencing data

Area of Science:

• Virology
• Bioinformatics
• Genetics

Background:

• Estimating viral diversity is crucial for understanding pathogen evolution and drug resistance.
• Deep sequencing generates large datasets, necessitating efficient diversity estimation tools.
• Existing methods based on multiple alignments are impractical for large-scale deep sequencing data.

Purpose of the Study:

• To develop and validate a novel method for accurately calculating population diversity from next-generation sequencing data.
• To address the limitations of current diversity estimation methods for large datasets.

Main Methods:

• A novel method, Pairwise Alignment Positional Nucleotide Counting (PAPNC), was developed for diversity estimation.
• PAPNC was applied to 454 Titanium deep sequencing data from HIV patient plasma samples.
• Diversity measurements were compared with those obtained from single-genome sequencing (SGS) using average pairwise difference (APD) in MEGA5.

Main Results:

• PAPNC diversity estimates from deep sequencing (0.002-0.021) were comparable to APD from SGS (0.001-0.018).
• A high correlation (R²=0.96) was observed between PAPNC and SGS-derived APD, with an average difference of ~1% after PCR error correction.
• PCR error contributed a background diversity of 0.0016 in control samples.

Conclusions:

• PAPNC is a novel and effective method for calculating genetic diversity from large next-generation sequencing datasets.
• The method can be integrated into existing variation calling and haplotype reconstruction software.
• PAPNC facilitates the study of viral evolution and drug resistance through accurate diversity estimation.