{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,9,2]],"date-time":"2026-09-02T06:50:33Z","timestamp":1788331833861,"version":"build-2803163510"},"reference-count":39,"publisher":"Public Library of Science (PLoS)","issue":"10","license":[{"start":{"date-parts":[[2025,10,16]],"date-time":"2025-10-16T00:00:00Z","timestamp":1760572800000},"content-version":"vor","delay-in-days":0,"URL":"http:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"content-domain":{"domain":["www.plosone.org"],"crossmark-restriction":false},"short-container-title":["PLoS One"],"abstract":"<jats:p>Protocols and pipelines for SARS-CoV-2 genome sequencing were rapidly established when the COVID-19 outbreak was declared a pandemic. The most widely used approach for sequencing SARS-CoV-2 includes targeted enrichment by PCR, followed by shotgun sequencing and reference-based genome assembly. As the continued surveillance of SARS-CoV-2 worldwide is transitioning towards a lower level of intensity, it is timely to re-visit the sequencing protocols and pipelines established during the acute phase of the pandemic. In the current study, we have investigated the impact of primer scheme and reference genome choice by sequencing samples with multiple primer schemes (Artic V3, V4.1 and V5.3.2) and re-processing reads with multiple reference genomes. We have also analysed the temporal development in ambiguous base calls during the emergence of the BA.2.86.x variant. We found that the primers used for targeted enrichment can result in recurrent ambiguous base calls, which can accumulate rapidly in response to the emergence of a new variant. We also found examples of consistent base calling errors, associated with PCR artifacts and amplicon drop-out. Similarly, misalignments and partially mapped reads on the reference genome resulted in ambiguous base calls, as well as defining mutations being omitted from the assembly. These findings highlight some key limitations of using targeted enrichment by PCR and reference-based genome assembly for sequencing SARS-CoV-2, and the importance of continuously monitoring and updating primer schemes and bioinformatic pipelines.<\/jats:p>","DOI":"10.1371\/journal.pone.0334009","type":"journal-article","created":{"date-parts":[[2025,10,16]],"date-time":"2025-10-16T17:23:01Z","timestamp":1760635381000},"page":"e0334009","update-policy":"https:\/\/doi.org\/10.1371\/journal.pone.corrections_policy","source":"Crossref","is-referenced-by-count":3,"title":["SARS-CoV-2 sequencing artifacts associated with targeted PCR enrichment and read mapping"],"prefix":"10.1371","volume":"20","author":[{"ORCID":"https:\/\/orcid.org\/0000-0003-1093-1645","authenticated-orcid":true,"given":"Kirsten Maren","family":"Ellegaard","sequence":"first","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-2020-9739","authenticated-orcid":true,"given":"Vithiagaran","family":"Gunalan","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-0372-128X","authenticated-orcid":true,"given":"Raphael","family":"Sieber","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0001-7929-2998","authenticated-orcid":true,"given":"Sharmin Jamshid","family":"Baig","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0001-5105-7056","authenticated-orcid":true,"given":"Nicolai Balle","family":"Larsen","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-2753-5991","authenticated-orcid":true,"given":"Marc","family":"Bennedb\u00e6k","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0003-1705-4008","authenticated-orcid":true,"given":"Jonas","family":"Bybjerg-Grauholm","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-5671-5609","authenticated-orcid":true,"given":"Leandro Andr\u00e9s","family":"Escobar-Herrera","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0009-0005-5281-211X","authenticated-orcid":true,"given":"Tobias","family":"Gress","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0009-0006-6920-8231","authenticated-orcid":true,"given":"Theis Hass","family":"Thorsen","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Anders","family":"Krusager","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0009-0004-5515-9964","authenticated-orcid":true,"given":"Gitte Nygaard","family":"Aasbjerg","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Nour Saad","family":"Al-Tamimi","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0009-0001-5040-2513","authenticated-orcid":true,"given":"Casper","family":"Westergaard","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-5755-5520","authenticated-orcid":true,"given":"Christina Wiid","family":"Svarrer","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0001-7491-4137","authenticated-orcid":true,"given":"Morten","family":"Rasmussen","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0003-0321-1180","authenticated-orcid":true,"given":"Marc","family":"Stegger","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"340","published-online":{"date-parts":[[2025,10,16]]},"reference":[{"key":"pone.0334009.ref001","unstructured":"Ghebreyesus, Tedros Adhanom. WHO media briefing [Internet]. https:\/\/www.who.int\/director-general\/speeches\/detail\/who-director-general-s-opening-remarks-at-the-media-briefing-on-covid-19---11-march-2020"},{"key":"pone.0334009.ref002","doi-asserted-by":"crossref","unstructured":"Carabelli AM, Peacock TP, Thorne LG, Harvey WT, Hughes J, COVID-19 Genomics UK Consortium, et al. SARS-CoV-2 variant biology: immune escape, transmission and fitness. Nat Rev Microbiol [Internet]. 2023 Jan 18 [cited 2025 Jun 17]. https:\/\/www.nature.com\/articles\/s41579-022-00841-7","DOI":"10.1038\/s41579-022-00841-7"},{"issue":"6","key":"pone.0334009.ref003","doi-asserted-by":"crossref","first-page":"361","DOI":"10.1038\/s41579-023-00878-2","article-title":"The evolution of SARS-CoV-2","volume":"21","author":"PV Markov","year":"2023","journal-title":"Nat Rev Microbiol"},{"issue":"3","key":"pone.0334009.ref004","doi-asserted-by":"crossref","first-page":"183","DOI":"10.1038\/nrmicro.2016.182","article-title":"Clinical and biological insights from viral genome sequencing","volume":"15","author":"CJ Houldcroft","year":"2017","journal-title":"Nat Rev Microbiol"},{"issue":"6","key":"pone.0334009.ref005","doi-asserted-by":"crossref","first-page":"1261","DOI":"10.1038\/nprot.2017.066","article-title":"Multiplex PCR method for MinION and Illumina sequencing of Zika and other virus genomes directly from clinical samples","volume":"12","author":"J Quick","year":"2017","journal-title":"Nat Protoc"},{"key":"pone.0334009.ref006","unstructured":"Artic Network [Internet]. https:\/\/artic.network"},{"key":"pone.0334009.ref007","unstructured":"SARS-CoV-2 version 5.3.2 scheme release [Internet]. https:\/\/community.artic.network\/t\/sars-cov-2-version-5-3-2-scheme-release\/462"},{"issue":"1","key":"pone.0334009.ref008","doi-asserted-by":"crossref","first-page":"8","DOI":"10.1186\/s13059-018-1618-7","article-title":"An amplicon-based sequencing framework for accurately measuring intrahost virus diversity using PrimalSeq and iVar","volume":"20","author":"ND Grubaugh","year":"2019","journal-title":"Genome Biol"},{"issue":"21","key":"pone.0334009.ref009","article-title":"Sources of PCR-induced distortions in high-throughput sequencing data sets","volume":"43","author":"JM Kebschull","year":"2015","journal-title":"Nucleic Acids Res"},{"issue":"7","key":"pone.0334009.ref010","doi-asserted-by":"crossref","first-page":"1687","DOI":"10.1093\/nar\/18.7.1687","article-title":"DNA recombination during PCR","volume":"18","author":"A Meyerhans","year":"1990","journal-title":"Nucleic Acids Res"},{"key":"pone.0334009.ref011","doi-asserted-by":"crossref","first-page":"106","DOI":"10.1016\/j.virusres.2016.12.009","article-title":"Error rates, PCR recombination, and sampling depth in HIV-1 whole genome deep sequencing","volume":"239","author":"F Zanini","year":"2017","journal-title":"Virus Research"},{"issue":"2","key":"pone.0334009.ref012","doi-asserted-by":"crossref","first-page":"189","DOI":"10.1111\/tbed.14339","article-title":"SARS-CoV-2 RNA stability in dry swabs for longer storage and transport at different temperatures","volume":"69","author":"A Alfaro-N\u00fa\u00f1ez","year":"2022","journal-title":"Transbound Emerg Dis"},{"issue":"3","key":"pone.0334009.ref013","doi-asserted-by":"crossref","first-page":"2000045","DOI":"10.2807\/1560-7917.ES.2020.25.3.2000045","article-title":"Detection of 2019 novel coronavirus (2019-nCoV) by real-time RT-PCR","volume":"25","author":"VM Corman","year":"2020","journal-title":"Euro Surveill"},{"issue":"1","key":"pone.0334009.ref014","doi-asserted-by":"crossref","DOI":"10.1128\/spectrum.03591-22","article-title":"Rapid and Flexible RT-qPCR Surveillance Platforms To Detect SARS-CoV-2 Mutations","volume":"11","author":"K Spiess","year":"2023","journal-title":"Microbiol Spectr"},{"issue":"10","key":"pone.0334009.ref015","doi-asserted-by":"crossref","DOI":"10.1371\/journal.pone.0274889","article-title":"A RT-qPCR system using a degenerate probe for specific identification and differentiation of SARS-CoV-2 Omicron (B.1.1.529) variants of concern","volume":"17","author":"R Jessen","year":"2022","journal-title":"PLOS ONE"},{"key":"pone.0334009.ref016","unstructured":"Krueger F. Trim Galore [Internet]. Babraham Institute; https:\/\/github.com\/FelixKrueger\/TrimGalore"},{"key":"pone.0334009.ref017","unstructured":"Li H. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. 2013 [cited 2024 May 2]. https:\/\/arxiv.org\/abs\/1303.3997"},{"issue":"16","key":"pone.0334009.ref018","doi-asserted-by":"crossref","first-page":"2078","DOI":"10.1093\/bioinformatics\/btp352","article-title":"The Sequence Alignment\/Map format and SAMtools","volume":"25","author":"H Li","year":"2009","journal-title":"Bioinforma Oxf Engl"},{"issue":"21","key":"pone.0334009.ref019","doi-asserted-by":"crossref","first-page":"2987","DOI":"10.1093\/bioinformatics\/btr509","article-title":"A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data","volume":"27","author":"H Li","year":"2011","journal-title":"Bioinformatics"},{"issue":"67","key":"pone.0334009.ref020","doi-asserted-by":"crossref","first-page":"3773","DOI":"10.21105\/joss.03773","article-title":"Nextclade: clade assignment, mutation calling and quality control for viral genomes","volume":"6","author":"I Aksamentov","year":"2021","journal-title":"JOSS"},{"key":"pone.0334009.ref021","unstructured":"Posit team. RStudio: Integrated Development Environment for R [Internet]. Posit Software; http:\/\/www.posit.co\/"},{"issue":"43","key":"pone.0334009.ref022","doi-asserted-by":"crossref","first-page":"1686","DOI":"10.21105\/joss.01686","article-title":"Welcome to the Tidyverse","volume":"4","author":"H Wickham","year":"2019","journal-title":"JOSS"},{"key":"pone.0334009.ref023","unstructured":"Baptiste A. gridExtra: Miscellaneous Functions for \u201cGrid\u201d Graphics [Internet]. https:\/\/CRAN.R-project.org\/package=gridExtra"},{"key":"pone.0334009.ref024","unstructured":"Wickham H. stringr: Simple, Consistent Wrappers for Common String Operations [Internet]. https:\/\/CRAN.R-project.org\/package=stringr"},{"key":"pone.0334009.ref025","unstructured":"Wickham H. ggplot2: Elegant Graphics for Data Analysis [Internet]. https:\/\/ggplot2.tidyverse.org"},{"key":"pone.0334009.ref026","unstructured":"Neuwirth E. RColorBrewer: ColorBrewer Palettes [Internet]. https:\/\/CRAN.R-project.org\/package=RColorBrewer"},{"issue":"2","key":"pone.0334009.ref027","doi-asserted-by":"crossref","first-page":"178","DOI":"10.1093\/bib\/bbs017","article-title":"Integrative Genomics Viewer (IGV): high-performance genomics data visualization and exploration","volume":"14","author":"H Thorvaldsd\u00f3ttir","year":"2013","journal-title":"Brief Bioinform"},{"issue":"14","key":"pone.0334009.ref028","doi-asserted-by":"crossref","first-page":"3059","DOI":"10.1093\/nar\/gkf436","article-title":"MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform","volume":"30","author":"K Katoh","year":"2002","journal-title":"Nucleic Acids Res"},{"issue":"9","key":"pone.0334009.ref029","doi-asserted-by":"crossref","first-page":"1189","DOI":"10.1093\/bioinformatics\/btp033","article-title":"Jalview Version 2\u2014a multiple sequence alignment editor and analysis workbench","volume":"25","author":"AM Waterhouse","year":"2009","journal-title":"Bioinformatics"},{"issue":"12","key":"pone.0334009.ref030","doi-asserted-by":"crossref","first-page":"3291","DOI":"10.1093\/bioinformatics\/btac311","article-title":"ReadItAndKeep: rapid decontamination of SARS-CoV-2 sequencing reads","volume":"38","author":"M Hunt","year":"2022","journal-title":"Bioinformatics"},{"issue":"12","key":"pone.0334009.ref031","doi-asserted-by":"crossref","DOI":"10.1016\/S1473-3099(23)00682-5","article-title":"Virus isolation and neutralisation of SARS-CoV-2 variants BA.2.86 and EG.5.1","volume":"23","author":"R Lassauni\u00e8re","year":"2023","journal-title":"Lancet Infect Dis"},{"issue":"36","key":"pone.0334009.ref032","doi-asserted-by":"crossref","first-page":"2300460","DOI":"10.2807\/1560-7917.ES.2023.28.36.2300460","article-title":"First cases of SARS-CoV-2 BA.2.86 in Denmark, 2023","volume":"28","author":"M Rasmussen","year":"2023","journal-title":"Euro Surveill"},{"key":"pone.0334009.ref033","unstructured":"Hodcroft, Emma. Covariants [Internet]. Variant: 23I (Omicron). https:\/\/covariants.org\/variants\/23I.Omicron"},{"key":"pone.0334009.ref034","article-title":"Addressing pandemic-wide systematic errors in the SARS-CoV-2 phylogeny","author":"M Hunt","year":"2024","journal-title":"bioRxiv"},{"key":"pone.0334009.ref035","doi-asserted-by":"crossref","first-page":"45","DOI":"10.1016\/j.gene.2010.08.009","article-title":"Reducing chimera formation during PCR amplification to ensure accurate genotyping","volume":"469","author":"RP Smyth","year":"2010","journal-title":"Gene"},{"issue":"9","key":"pone.0334009.ref036","doi-asserted-by":"crossref","first-page":"2083","DOI":"10.1093\/nar\/30.9.2083","article-title":"Heteroduplexes in mixed-template amplifications: formation, consequence and elimination by \u201creconditioning PCR\u201d","volume":"30","author":"JR Thompson","year":"2002","journal-title":"Nucleic Acids Res"},{"key":"pone.0334009.ref037","author":"N De Maio","year":"2024","journal-title":"Rate variation and recurrent sequence errors in pandemic-scale phylogenetics"},{"key":"pone.0334009.ref038","unstructured":"Russell Corbett-Detig. Issues with SARS-CoV-2 sequencing data [Internet]. 2021. https:\/\/virological.org\/t\/issues-with-sars-cov-2-sequencing-data\/473\/16"},{"issue":"1","key":"pone.0334009.ref039","doi-asserted-by":"crossref","first-page":"708","DOI":"10.1186\/s12864-016-3030-6","article-title":"Viral deep sequencing needs an adaptive approach: IRMA, the iterative refinement meta-assembler","volume":"17","author":"SS Shepard","year":"2016","journal-title":"BMC Genomics"}],"container-title":["PLOS One"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/dx.plos.org\/10.1371\/journal.pone.0334009","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2025,10,16]],"date-time":"2025-10-16T17:23:11Z","timestamp":1760635391000},"score":1,"resource":{"primary":{"URL":"https:\/\/dx.plos.org\/10.1371\/journal.pone.0334009"}},"subtitle":[],"editor":[{"given":"Nihad A.M","family":"Al-Rashedi","sequence":"first","affiliation":[],"role":[{"vocabulary":"crossref","role":"editor"}]}],"short-title":[],"issued":{"date-parts":[[2025,10,16]]},"references-count":39,"journal-issue":{"issue":"10","published-online":{"date-parts":[[2025,10,16]]}},"URL":"https:\/\/doi.org\/10.1371\/journal.pone.0334009","relation":{},"ISSN":["1932-6203"],"issn-type":[{"value":"1932-6203","type":"electronic"}],"subject":[],"published":{"date-parts":[[2025,10,16]]}}}