Abstract:Abstract Motivation Crucial to the correctness of a genome assembly is the accuracy of the underlying scaffolds that specify the orders and orientations of contigs together with the gap distances between contigs. The current methods construct scaffolds based on the alignments of ‘linking’ reads against contigs. We found that some ‘optimal’ alignments are mistaken due to factors such as the contig boundary effect, particularly in the presence of repeats. Occasionally, the incorrect alignments can even overwhelm the correct ones. The detection of the incorrect linking information is challenging in any existing methods. Results In this study, we present a novel scaffolding method RegScaf. It first examines the distribution of distances between contigs from read alignment by the kernel density. When multiple modes are shown in a density, orientation-supported links are grouped into clusters, each of which defines a linking distance corresponding to a mode. The linear model parameterizes contigs by their positions on the genome; then each linking distance between a pair of contigs is taken as an observation on the difference of their positions. The parameters are estimated by minimizing a global loss function, which is a version of trimmed sum of squares. The least trimmed squares estimate has such a high breakdown value that it can automatically remove the mistaken linking distances. The results on both synthetic and real datasets demonstrate that RegScaf outperforms some popular scaffolders, especially in the accuracy of gap estimates by substantially reducing extremely abnormal errors. Its strength in resolving repeat regions is exemplified by a real case. Its adaptability to large genomes and TGS long reads is validated as well. Availability and implementation RegScaf is publicly available at https://github.com/lemontealala/RegScaf.git. Supplementary information Supplementary data are available at Bioinformatics online.

dScaff – an automatic bioinformatics framework for scaffolding draft de novo assemblies based on reference genome data

P_RNA_scaffolder: a Fast and Accurate Genome Scaffolder Using Paired-End RNA-sequencing Reads

Efficient Iterative Hi-C Scaffolder Based on N-best Neighbors

LDscaff: LD-based scaffolding of de novo genome assemblies

ntLink: a toolkit for de novo genome assembly scaffolding and mapping using long reads

PEP_scaffolder: Using (homologous) Proteins to Scaffold Genomes

RegScaf: a regression approach to scaffolding

SCOP: a Novel Scaffolding Algorithm Based on Contig Classification and Optimization.

L_RNA_scaffolder: scaffolding genomes with transcripts

MeDuSa: a multi-draft based scaffolder

SLR-superscaffolder: a de novo scaffolding tool for synthetic long reads using a top-to-bottom scheme

Comparison of Hi-C-Based Scaffolding Tools on Plant Genomes

BOSS: a Novel Scaffolding Algorithm Based on an Optimized Scaffold Graph

PEP_scaffolder: using (homologous) proteins to scaffold genomes.

Benchmarking of Hi-C tools for scaffolding de novo genome assemblies

New algorithms for accurate and efficient de novo genome assembly from long DNA sequencing reads

De Novo Assembly of Human Genomes with Massively Parallel Short Read Sequencing.

Exploring Genome Characteristics and Sequence Quality Without a Reference

Dime: A Novel Framework for De Novo Metagenomic Sequence Assembly

Benchmarking of next and third generation sequencing technologies and their associated algorithms for de novo genome assembly

De Novo Assembly Methods for Next Generation Sequencing Data