Uneven sequencing (coverage) depth can bias microbial intraspecies diversity estimates and how to account for it

An unbiased and accurate estimation of intraspecies diversity, i.e. the extent of genetic diversity within species (or microdiversity), is crucial for clinical and environmental microbiome studies. Although it is well appreciated that sequencing depth (or coverage depth) below 10X can provide biased estimates of microdiversity, typically underestimating diversity due to the random sampling of alleles, there is a widely accepted convention that microdiversity estimates tend to be relatively stable at sequencing depth exceeding 10X. Therefore, discarding species with <10X or rarefying to 10-20X sequencing depth are generally used to compare microdiversity among taxa and samples. Our findings showed that these biases may persist even at depth levels above 50-200X for all popular sequencing platforms, including Illumina, PacBio, and Oxford Nanopore. The biases mostly, but not always, represent an underestimation of diversity and were attributable to the incomplete recovery of Single Nucleotide Variants (SNVs) at lower sequencing depth levels. To address this issue, we recommend using rarefaction-based approaches to standardize data at least 50X, and ideally at 200X sequencing depth, which reduces differences between observed and expected microdiversity values to <0.5%. Furthermore, the Average Nucleotide Identity of reads (ANIr) metric is significantly less sensitive to sequencing depth variability than nucleotide diversity (π), making it a robust alternative for estimating microdiversity at sequencing depth close or exceeding 10X, without a need to rarefying data. Therefore, the sequencing depth thresholds proposed herein provide a more standardized framework for direct comparisons of microdiversity across samples and studies.