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.

How to handle high subgenome sequence similarity in allopolyploid Fragaria x ananassa: linkage disequilibrium based variant filtering

Defining Loci in Restriction-Based Reduced Representation Genomic Data from Nonmodel Species: Sources of Bias and Diagnostics for Optimal Clustering

Subgenome phasing for complex allopolyploidy: case-based benchmarking and recommendations

Blueprint for Phasing and Assembling the Genomes of Heterozygous Polyploids: Application to the Octoploid Genome of Strawberry

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

Ingap-Family:Accurate Detection of Meiotic Recombination Loci and Causal Mutations by Filtering out Artificial Variants Due to Genome Complexities

Large-scale identification of single-feature polymorphisms in complex genomes

Detection of Genomic Variations and DNA Polymorphisms and Impact on Analysis of Meiotic Recombination and Genetic Mapping.

Disambiguate: An open-source application for disambiguating two species in next generation sequencing data from grafted samples

Identifying Genome-Wide Sequence Variations And Candidate Genes Implicated In Self-Incompatibility By Resequencing Fragaria Viridis

Sensitive Long-Indel-Aware Alignment of Sequencing Reads

Variant calling in polyploids for population and quantitative genetics

Genotyping Polyploids from Messy Sequencing Data

Multiple Variant Calling Pipelines in Wheat Whole Exome Sequencing

A Novel Multi-Alignment Pipeline for High-Throughput Sequencing Data.

Positional‐based cloning ‘fail‐safe’ approach is overpowered by wheat chromosome structural variation

A general approach to single-nucleotide polymorphism discovery

polymapR—linkage analysis and genetic map construction from F1 populations of outcrossing polyploids

Variant filters using segregation information improve mapping of nectar-production genes in sunflower (Helianthus annuus L.)

HaploMerger: Reconstructing Allelic Relationships for Polymorphic Diploid Genome Assemblies.

Genome-Wide Association Mapping To Candidate Polymorphism Resolution In The Unsequenced Barley Genome