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.

L_RNA_scaffolder: scaffolding genomes with transcripts

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

PEP_scaffolder: Using (homologous) Proteins to Scaffold Genomes

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

PEP_scaffolder: using (homologous) proteins to scaffold genomes.

L_RNA_scaffolder: Scaffolding Genomes with Transcripts (vol 14, 604, 2013)

Correction To: P_RNA_scaffolder: a Fast and Accurate Genome Scaffolder Using Paired-End RNAsequencing Reads

LDscaff: LD-based scaffolding of de novo genome assemblies

Bridger: A New Framework for De Novo Transcriptome Assembly Using Rna-Seq Data

Rascaf: Improving Genome Assembly with RNA Sequencing Data.

Scallop2 Enables Accurate Assembly of Multiple-End RNA-seq Data

A Transcriptome Post-Scaffolding Method for Assembling High Quality Contigs

Whole-Genome Restriction Mapping by "Subhaploid"-Based RAD Sequencing: An Efficient and Flexible Approach for Physical Mapping and Genome Scaffolding.

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

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

RegScaf: a regression approach to scaffolding

IsoTree: A New Framework for De Novo Transcriptome Assembly from RNA-seq Reads

Deep annotation of long noncoding RNAs by assembling RNA-seq and small RNA-seq data

IsoTree: De Novo Transcriptome Assembly from RNA-Seq Reads

Agtar: A Novel Approach For Transcriptome Assembly And Abundance Estimation Using An Adapted Genetic Algorithm From Rna-Seq Data

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