Attack of the PCR clones: Rates of clonality have little effect on RAD-seq genotype calls

Molecular Ecology Resources
By: , and 

Links

Abstract

Interpretation of high-throughput sequence data requires an understanding of how decisions made during bioinformatic data processing can influence results. One source of bias that is often cited is PCR clones (or PCR duplicates). PCR clones are common in restriction site-associated sequencing (RAD-seq) data sets, which are increasingly being used for molecular ecology. To determine the influence PCR clones and the bioinformatic handling of clones have on genotyping, we evaluate four RAD-seq data sets. Data sets were compared before and after clones were removed to estimate the number of clones present in RAD-seq data, quantify how often the presence of clones in a data set causes genotype calls to change compared to when clones were removed, investigate the mechanisms that lead to genotype call changes and test whether clones bias heterozygosity estimates. Our RAD-seq data sets contained 30%–60% PCR clones, but 95% of RAD-tags had five or fewer clones. Relatively few genotypes changed once clones were removed (5%–10%), and the vast majority of these changes (98%) were associated with genotypes switching from a called to no-call state or vice versa. PCR clones had a larger influence on genotype calls in individuals with low read depth but appeared to influence genotype calls at all loci similarly. Removal of PCR clones reduced the number of called genotypes by 2% but had almost no influence on estimates of heterozygosity. As such, while steps should be taken to limit PCR clones during library preparation, PCR clones are likely not a substantial source of bias for most RAD-seq studies.

Publication type Article
Publication Subtype Journal Article
Title Attack of the PCR clones: Rates of clonality have little effect on RAD-seq genotype calls
Series title Molecular Ecology Resources
DOI 10.1111/1755-0998.13087
Volume 20
Issue 1
Year Published 2020
Language English
Publisher Wiley
Contributing office(s) Coop Res Unit Leetown
Description 13 p.
First page 66
Last page 78
Google Analytic Metrics Metrics page
Additional publication details