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.

Hybrid-hybrid correction of errors in long reads with HERO

Telomere-to-Telomere Phased Genome Assembly Using HERRO-Corrected Simplex Nanopore Reads

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

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

Bi-Level Error Correction for PacBio Long Reads

Pre-Assembly NGS Correction of ONT Reads Achieves HiFi-Level Assembly Quality

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

Correcting Illumina sequencing errors for human data

Hybrid Assembly of Ultra-Long Nanopore Reads Augmented with 10X-Genomics Contigs: Demonstrated with a Human Genome.

HALC: High throughput algorithm for long read error correction

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

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

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

Efficient assembly of nanopore reads via highly accurate and intact error correction

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

Repeat and haplotype aware error correction in nanopore sequencing reads with DeChat

MapReduce for Accurate Error Correction of Next-Generation Sequencing Data

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

Joint Haplotype Phasing and Genotype Calling of Multiple Individuals Using Haplotype Informative Reads

An Approach to Correcting DNA Sequencing Error

An Optimization Approach to the Reconstruction of Positional DNA Sequencing by Hybridization with Errors