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.

LDscaff: LD-based scaffolding of de novo genome assemblies

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

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

PEP_scaffolder: using (homologous) proteins to scaffold genomes.

PEP_scaffolder: Using (homologous) Proteins to Scaffold Genomes

L_RNA_scaffolder: scaffolding genomes with transcripts

RegScaf: a regression approach to scaffolding

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

Identifying the Causes and Consequences of Assembly Gaps Using a Multiplatform Genome Assembly of a Bird-of-paradise.

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

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

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

EPGA: de novo assembly using the distributions of reads and insert size

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

EPGA-SC : A Framework for <italic>de novo</italic> Assembly of Single-Cell Sequencing Reads

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

Paired-End Sequencing of Long-Range DNA Fragments for De Novo Assembly of Large, Complex Mammalian Genomes by Direct Intra-Molecule Ligation

Genome sequence assembly evaluation using long-range sequencing data

Scalable De Novo Genome Assembly Using Pregel

Distributed-Memory Parallel Contig Generation for De Novo Long-Read Genome Assembly

Genome analysis EPGA 2 : memory-efficient de novo assembler