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.

Sequencing accuracy and systematic errors of nanopore direct RNA sequencing

Sequencing DNA with nanopores: Troubles and biases

Direct RNA sequencing (RNA004) allows for improved transcriptome assessment and near real-time tracking of methylation for medical applications

Advances in nanopore direct RNA sequencing

Error correction enables use of Oxford Nanopore technology for reference-free transcriptome analysis

Oxford Nanopore MinION Direct RNA-Seq for Systems Biology

From squiggle to basepair: computational approaches for improving nanopore sequencing read accuracy

Comparison of the Two Up-to-date Sequencing Technologies for Genome Assembly: HiFi Reads of Pacific Biosciences Sequel II System and Ultralong Reads of Oxford Nanopore

Comparative analysis of 43 distinct RNA modifications by nanopore tRNA sequencing

Applications and Potentials of Nanopore Sequencing in the (epi)genome and (epi)transcriptome Era

Analysis of bacterial transcriptome and epitranscriptome using nanopore direct RNA sequencing

Nanopore sequencing of RNA and cDNA molecules in Escherichia coli

Systematic comparison of tools used for m6A mapping from nanopore direct RNA sequencing

An Error Correction Method of Nanopore Sequencing Data Using Deep Learning

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

Decoding the epitranscriptional landscape from native RNA sequences

Correcting modification-mediated errors in nanopore sequencing by nucleotide demodification and reference-based correction

Benchmarking of the Oxford Nanopore MinION sequencing for quantitative and qualitative assessment of cDNA populations

On Nanopore DNA Sequencing by Signal and Noise Analysis of Ionic Current.

The newest Oxford Nanopore R10.4.1 full-length 16S rRNA sequencing enables the accurate resolution of species-level microbial community profiling

Quantitative profiling of pseudouridylation dynamics in native RNAs with nanopore sequencing