Abstract:Background Many medicinal plants are known for their complex genomes with high ploidy, heterozygosity, and repetitive content which pose severe challenges for genome sequencing of those species. Long reads from Oxford nanopore sequencing technology (ONT) or Pacific Biosciences Single Molecule, Real-Time (SMRT) sequencing offer great advantages in de novo genome assembly, especially for complex genomes with high heterozygosity and repetitive content. Currently, multiple allotetraploid species have sequenced their genomes by long-read sequencing. However, we found that a considerable proportion of these genomes (7.9% on average, maximum 23.7%) could not be covered by NGS (Next Generation Sequencing) reads (uncovered region by NGS reads, UCR) suggesting the questionable and low-quality of those area or genomic areas that can’t be sequenced by NGS due to sequencing bias. The underlying causes of those UCR in the genome assembly and solutions to this problem have never been studied. Methods In the study, we sequenced the tetraploid genome of Veratrum dahuricum (Turcz.) O. Loes (VDL), a Chinese medicinal plant, with ONT platform and assembled the genome with three strategies in parallel. We compared the qualities, coverage, and heterozygosity of the three ONT assemblies with another released assembly of the same individual using reads from PacBio circular consensus sequencing (CCS) technology, to explore the cause of the UCR. Results By mapping the NGS reads against the three ONT assemblies and the CCS assembly, we found that the coverage of those ONT assemblies by NGS reads ranged from 49.15 to 76.31%, much smaller than that of the CCS assembly (99.53%). And alignment between ONT assemblies and CCS assembly showed that most UCR can be aligned with CCS assembly. So, we conclude that the UCRs in ONT assembly are low-quality sequences with a high error rate that can’t be aligned with short reads, rather than genomic regions that can’t be sequenced by NGS. Further comparison among the intermediate versions of ONT assemblies showed that the most probable origin of those errors is a combination of artificial errors introduced by “self-correction” and initial sequencing error in long reads. We also found that polishing the ONT assembly with CCS reads can correct those errors efficiently. Conclusions Through analyzing genome features and reads alignment, we have found the causes for the high proportion of UCR in ONT assembly of VDL are sequencing errors and additional errors introduced by self-correction. The high error rates of ONT-raw reads make them not suitable for self-correction prior to allotetraploid genome assembly, as the self-correction will introduce artificial errors to > 5% of the UCR sequences. We suggest high-precision CCS reads be used to polish the assembly to correct those errors effectively for polyploid genomes.

Integration of Hybrid and Self-Correction Method Improves the Quality of Long-Read Sequencing Data.

A Crowdsourcing Method For Correcting Sequencing Errors For The Third-Generation Sequencing Data

Pseudo-Sanger Sequencing: Massively Parallel Production of Long and Near Error-Free Reads Using NGS Technology

Hybrid-hybrid correction of errors in long reads with HERO

Comprehensive assessment of error correction methods for high-throughput sequencing data

HALC: High throughput algorithm for long read error correction

Bi-Level Error Correction for PacBio Long Reads

Turn ‘noise’ to signal: accurately rectify millions of erroneous short reads through graph learning on edit distances

Instance-based Error Correction for Short Reads of Disease-Associated Genes.

Quality-Score Guided Error Correction for Short-Read Sequencing Data Using Cuda

An Approach to Correcting DNA Sequencing Error

How Error Correction Affects PCR Deduplication: A Survey Based on UMI Datasets of Short Reads

Comprehensive Assessment of Isoform Detection Methods for Third-Generation Sequencing Data

MapReduce for Accurate Error Correction of Next-Generation Sequencing Data

A Sequence-Based Novel Approach for Quality Evaluation of Third-Generation Sequencing Reads.

DBG2OLC: Efficient Assembly of Large Genomes Using Long Erroneous Reads of the Third Generation Sequencing Technologies

Comparison of ONT and CCS sequencing technologies on the polyploid genome of a medicinal plant showed that high error rate of ONT reads are not suitable for self-correction

Opportunities and challenges in long-read sequencing data analysis

Error filtering, pair assembly and error correction for next-generation sequencing reads

NmTHC: a hybrid error correction method based on a generative neural machine translation model with transfer learning

ResSeq: Enhancing Short-Read Sequencing Alignment by Rescuing Error-Containing Reads