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.

Optimal assembly strategies of transcriptome related to ploidies of eukaryotic organisms

TRITEX: Chromosome-Scale Sequence Assembly of Triticeae Genomes with Open-Source Tools

Optimizing de novo transcriptome assembly from short-read RNA-Seq data: a comparative study

Comparison of De Novo Transcriptome Assemblers and k-mer Strategies Using the Killifish, Fundulus heteroclitus

Combining independent de novo assemblies optimizes the coding transcriptome for nonconventional model eukaryotic organisms

Optimized Sequencing Depth and De Novo Assembler for Deeply Reconstructing the Transcriptome of the Tea Plant, an Economically Important Plant Species

Performance of Gene Expression Analyses Using De Novo Assembled Transcripts in Polyploid Species

Augmenting transcriptome assembly combinatorially

De Novo Transcriptome Assembly of RNA-Seq Reads with Different Strategies

Comparative study of de novo assembly and genome-guided assembly strategies for transcriptome reconstruction based on RNA-Seq

Construction of a de novo assembly pipeline using multiple transcriptome data sets from Cypripedium macranthos (Orchidaceae)

Evaluation of Assembly Strategies Using RNA-Seq Data Associated with Grain Development of Wheat (triticum Aestivum L.)

Single-pollen-cell Sequencing for Gamete-Based Phased Diploid Genome Assembly in Plants

A high-throughput SNP discovery strategy for RNA-seq data

The Impacts of Read Length and Transcriptome Complexity for De Novo Assembly: A Simulation Study

A Transcriptome Post-Scaffolding Method for Assembling High Quality Contigs

Computational Approaches for Transcriptome Assembly Based on Sequencing Technologies

Gcapda: a Haplotype-Resolved Diploid Assembler

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

BinPacker: Packing-Based De Novo Transcriptome Assembly from RNA-seq Data

Enhancing transcriptome analysis in medicinal plants: multiple unigene sets in Astragalus membranaceus