Accurate sequencing of DNA motifs able to form alternative (non-B) structures
- Matthias H. Weissensteiner1,11,12,
- Marzia A. Cremona2,3,4,11,
- Wilfried M. Guiblet1,5,
- Nicholas Stoler6,
- Robert S. Harris1,
- Monika Cechova1,7,13,
- Kristin A. Eckert4,8,
- Francesca Chiaromonte4,9,10,
- Yi-Fei Huang1,4 and
- Kateryna D. Makova1,4
- 1Department of Biology, The Pennsylvania State University, University Park, Pennsylvania 16802, USA;
- 2Department of Operations and Decision Systems, Université Laval, Quebec, Quebec G1V0A6, Canada;
- 3Population Health and Optimal Health Practices, CHU de Québec–Université Laval Research Center, Québec, Quebec G1V4G2, Canada;
- 4Center for Medical Genomics, The Pennsylvania State University, University Park, Pennsylvania 16802, USA;
- 5Laboratory of Cell Biology, NCI-CCR, National Institutes of Health, Bethesda, Maryland 20892, USA;
- 6Department of Biochemistry and Molecular Biology, The Pennsylvania State University, University Park, Pennsylvania 16802, USA;
- 7Faculty of Informatics, Masaryk University, 60200 Brno, Czech Republic;
- 8Department of Pathology, The Pennsylvania State University, College of Medicine, Hershey, Pennsylvania 17033, USA;
- 9Department of Statistics, The Pennsylvania State University, University Park, Pennsylvania 16802, USA;
- 10Institute of Economics and L'EMbeDS, Sant'Anna School of Advanced Studies, Pisa 56127, Italy
-
↵11 These authors contributed equally to this work.
Abstract
Approximately 13% of the human genome at certain motifs have the potential to form noncanonical (non-B) DNA structures (e.g., G-quadruplexes, cruciforms, and Z-DNA), which regulate many cellular processes but also affect the activity of polymerases and helicases. Because sequencing technologies use these enzymes, they might possess increased errors at non-B structures. To evaluate this, we analyzed error rates, read depth, and base quality of Illumina, Pacific Biosciences (PacBio) HiFi, and Oxford Nanopore Technologies (ONT) sequencing at non-B motifs. All technologies showed altered sequencing success for most non-B motif types, although this could be owing to several factors, including structure formation, biased GC content, and the presence of homopolymers. Single-nucleotide mismatch errors had low biases in HiFi and ONT for all non-B motif types but were increased for G-quadruplexes and Z-DNA in all three technologies. Deletion errors were increased for all non-B types but Z-DNA in Illumina and HiFi, as well as only for G-quadruplexes in ONT. Insertion errors for non-B motifs were highly, moderately, and slightly elevated in Illumina, HiFi, and ONT, respectively. Additionally, we developed a probabilistic approach to determine the number of false positives at non-B motifs depending on sample size and variant frequency, and applied it to publicly available data sets (1000 Genomes, Simons Genome Diversity Project, and gnomAD). We conclude that elevated sequencing errors at non-B DNA motifs should be considered in low-read-depth studies (single-cell, ancient DNA, and pooled-sample population sequencing) and in scoring rare variants. Combining technologies should maximize sequencing accuracy in future studies of non-B DNA.
Footnotes
-
↵Present addresses: 12Institute of Avian Research, D-26386 Wilhelmshaven, Germany;
-
↵13Department of Biomolecular Engineering, University of California Santa Cruz, Santa Cruz, CA 95064, USA
-
[Supplemental material is available for this article.]
-
Article published online before print. Article, supplemental material, and publication date are at https://www.genome.org/cgi/doi/10.1101/gr.277490.122.
-
Freely available online through the Genome Research Open Access option.
- Received November 9, 2022.
- Accepted May 4, 2023.
This article, published in Genome Research, is available under a Creative Commons License (Attribution-NonCommercial 4.0 International), as described at http://creativecommons.org/licenses/by-nc/4.0/.