DipSkmer: Reference-free population genomics with diploid genome skims
Discuss this preprint
Start a discussion What are Sciety discussions?Listed in
This article is not in any list yet, why not save it to one of your lists.Abstract
Ecologists and conservation biologists rely on genetic diversity as a key essential biodiversity variable (EBV) used to track population health and dynamics, and utilize the population parameter θ (estimated by the average pairwise genomic distance) as a key metric of diversity. While whole-genome-sequencing (wgs) is increasingly affordable, it will be considerable time before the full diversity of life is represented by high-quality assembled genomes; even then, constant monitoring will still require repeated sampling of populations. In contrast, genome skimming (low-coverage, short-read wgs) is highly cost-effective but challenging to analyze because the coverage is too low for assembly and reliable error correction. Mature methods, such as Mash, exist for estimating pairwise genomic distances based on the Jaccard similarity of k -mer sets computed using sketching techniques. Some, such as Skmer, additionally model the impacts of low coverage. These methods have been successfully applied to assembly-free species identification and phylogenetics; however, their use in population genetics has been limited. This is because these methods implicitly treat genomes as haploid and heterozygosity confounds true estimates of genomic distance for diploid organisms. In this paper, we address this problem through a number of technical advances. First, we use coalescent theory to mathematically derive how the Jaccard index between two diploid samples changes with the scaled population size parameter ( θ ). Next, we derive an estimator that computes θ from the Jaccard index, in addition to several auxiliary variables, which we also estimate from the genome skims. The resulting method, DipSkmer, enables more accurate estimates of coverage, sequencing error, and pairwise nucleotide distance for diploid samples. Analyses of both simulated and empirical datasets show that for diploids and low distances (e.g., < 2%), Dip-Skmer produces the most accurate pairwise distance estimates, outperforming existing alignment-free methods such as Mash and Skmer, and closely approximates ANGSD, a reference and alignment-based tool.
Availability
The code for DipSkmer is available at https://github.com/echarvel3/ReSkmer/tree/DipSkmer-REFACTOR . Simulation scripts and environments are available at https://github.com/echarvel3/dipskmer_scripts .
Author Summary
The process of obtaining full-genome population genomic measurements for biodiversity monitoring remains expensive due to the need for high-coverage sequencing and reference assemblies. Genome skimming has been shown to be a viable, low-coverage alternative for obtaining genomic distances, and alignment- and assembly-free methods exist for analyzing nuclear data from skimming data to estimate the distance between samples. However, existing methods fail to model within-sample heterozygosity, expected for diploid organisms. Given the dominance of diploidy among species of interest to ecologists, the implications of these simplifying assumptions warrant further study. Here, we present a mathematical model of the k -mer sets sampled from two diploid genomes from a Wright-Fisher population. We use the model to develop DipSkmer, a k -mer-based, reference-free method for estimating nucleotide diversity and population divergence that, unlike its predecessors, models within-sample heterozygosity. Benchmarking shows more accurate genomic diversity estimates compared to existing reference-free, genome-skimming methods and comparable performance to the popular high-coverage, reference-based method, ANGSD. Thus, DipSkmer enables accessible, less expensive population monitoring through genetic diversity estimates.
