Reference-free SNP calling: improved accuracy by preventing incorrect calls from repetitive genomic regions

Abstract

BackgroundSingle nucleotide polymorphisms (SNPs) are the most abundant type of genetic variation in eukaryotic genomes and have recently become the marker of choice in a wide variety of ecological and evolutionary studies. The advent of next-generation sequencing (NGS) technologies has made it possible to efficiently genotype a large number of SNPs in the non-model organisms with no or limited genomic resources. Most NGS-based genotyping methods require a reference genome to perform accurate SNP calling. Little effort, however, has yet been devoted to developing or improving algorithms for accurate SNP calling in the absence of a reference genome.ResultsHere we describe an improved maximum likelihood (ML) algorithm called iML, which can achieve high genotyping accuracy for SNP calling in the non-model organisms without a reference genome. The iML algorithm incorporates the mixed Poisson/normal model to detect composite read clusters and can efficiently prevent incorrect SNP calls resulting from repetitive genomic regions. Through analysis of simulation and real sequencing datasets, we demonstrate that in comparison with ML or a threshold approach, iML can remarkably improve the accuracy of de novo SNP genotyping and is especially powerful for the reference-free genotyping in diploid genomes with high repeat contents.ConclusionsThe iML algorithm can efficiently prevent incorrect SNP calls resulting from repetitive genomic regions, and thus outperforms the original ML algorithm by achieving much higher genotyping accuracy. Our algorithm is therefore very useful for accurate de novo SNP genotyping in the non-model organisms without a reference genome.ReviewersThis article was reviewed by Dr. Richard Durbin, Dr. Liliana Florea (nominated by Dr. Steven Salzberg) and Dr. Arcady Mushegian.

Reference-Free Population Genomics from Next-Generation Transcriptome Data and the Vertebrate–Invertebrate Gap

“…Dou et al [59] recently highlighted this problem, and proposed a method to overcome it, based on the idea that sequencing coverage is expected to be higher in repeated than in unique genomic regions. This approach does not apply to transcriptomic data, in which coverage primarily reflects the level of gene expression, which is not only determined by gene copy number.…”

Widespread intersex differentiation across the stickleback genome – The signature of sexually antagonistic selection?

“…This scenario leads to two testable predictions (see also Dou et al, 2012;McKinney, Waples, Seeb, & Seeb, 2017;Tsai, Evans, Noorai, Starr-Moss, & Clark, 2019). We here considered that both the reference genome and the new de novo genome are derived from a female (XX) individual.…”

Whole-Genome Restriction Mapping by “Subhaploid”-Based RAD Sequencing: An Efficient and Flexible Approach for Physical Mapping and Genome Scaffolding

“…For a 2bRAD data set, the read depth k for each unique restriction site theoretically follows the Poisson distribution, assuming that all sites across the genome are evenly sequenced (Dou et al 2012):$\text{Poisson}\left(k|C\right)\sim \frac{C^k{e}^{-C}}{k!}.$If the observed depth of one site is no more than 2, we regard the genotype of this site as 0 or “undetermined,” indicating that this marker information cannot be captured using the threshold method. Therefore, the FNR can be represented as:…”