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.

Multi-CSAR: a web server for scaffolding contigs using multiple reference genomes

CSA: A High-Throughput Chromosome-Scale Assembly Pipeline for Vertebrate Genomes.

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

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

RegScaf: a regression approach to scaffolding

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

Reference-Assisted Chromosome Assembly

MeDuSa: a multi-draft based scaffolder

MOMS: A pipeline for scaffolding using multi‐optical maps

A Web-Based Multi-Genome Synteny Viewer for Customized Data

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

A Scalable and Accurate Targeted Gene Assembly Tool (SAT-Assembler) for Next-Generation Sequencing Data

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

Multi-genome synteny detection using minimizer graph mappings

CSA: a Web Service for the Complete Process of ChIP-Seq Analysis

Cbrowse: A Sam/Bam-Based Contig Browser For Transcriptome Assembly Visualization And Analysis

CAFU: a Galaxy Framework for Exploring Unmapped RNA-Seq Data

PEP_scaffolder: Using (homologous) Proteins to Scaffold Genomes

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

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

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