Abstract:Selecting the values of parameters used by de novo genomic assembly programs, or choosing an optimal de novo assembly from several runs obtained with different parameters or programs, are tasks that can require complex decision-making. A key parameter that must be supplied to typical next generation sequencing (NGS) assemblers is the k-mer length, i.e., the word size that determines which de Bruijn graph the program should map out and use. The topic of assembly selection criteria was recently revisited in the Assemblathon 2 study (Bradnam et al., 2013). Although no clear message was delivered with regard to optimal k-mer lengths, it was shown with examples that it is sometimes important to decide if one is most interested in optimizing the sequences of protein-coding genes (the gene space) or in optimizing the whole genome sequence including the intergenic DNA, as what is best for one criterion may not be best for the other. In the present study, our aim was to better understand how the assembly of unicellular fungi (which are typically intermediate in size and complexity between prokaryotes and metazoan eukaryotes) can change as one varies the k-mer values over a wide range. We used two different de novo assembly programs (SOAPdenovo2 and ABySS), and simple assembly metrics that also focused on success in assembling the gene space and repetitive elements. A recent increase in Illumina read length to around 150 bp allowed us to attempt de novo assemblies with a larger range of k-mers, up to 127 bp. We applied these methods to Illumina paired-end sequencing read sets of fungal strains of Paracoccidioides brasiliensis and other species. By visualizing the results in simple plots, we were able to track the effect of changing k-mer size and assembly program, and to demonstrate how such plots can readily reveal discontinuities or other unexpected characteristics that assembly programs can present in practice, especially when they are used in a traditional molecular microbiology laboratory with a 'genomics corner'. Here we propose and apply a component of a first pass validation methodology for benchmarking and understanding fungal genome de novo assembly processes.

Methodology for Assessing the Quality of Genomic Assembly Based on the Analysis of K-Mers Frequency in a Parallel Sequencing Sequencer

METHODOLOGY FOR ASSESSING THE QUALITY OF GENOMIC ASSEMBLY BASED ON THE ANALYSIS OF THE FREQUENCY OF K-MERS IN A PARALLEL SEQUENCING SEQUENCER

An efficient parallel algorithm for multiple sequence similarities calculation using a low complexity method.

Nanopore Sequencing Technology and Tools for Genome Assembly: Computational Analysis of the Current State, Bottlenecks and Future Directions

Estimation of genomic characteristics by analyzing k-mer frequency in de novo genome projects

Merqury: reference-free quality, completeness, and phasing assessment for genome assemblies

Efficient De Novo Assembly and Recovery of Microbial Genomes from Complex Metagenomes Using a Reduced Set of k-mers

KAT: a K-mer analysis toolkit to quality control NGS datasets and genome assemblies

Benchmarking genome assembly methods on metagenomic sequencing data

Estimation of genome size using k-mer frequencies from corrected long reads

Assessment of Metagenomic Assemblers Based on Hybrid Reads of Real and Simulated Metagenomic Sequences

The complex task of choosing a de novo assembly: lessons from fungal genomes

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

Accelerating K-mer Frequency Counting with GPU and Non-Volatile Memory

These are not the k-mers you are looking for: efficient online k-mer counting using a probabilistic data structure

Assembly algorithms for next-generation sequencing data

Pseudo-Sanger Sequencing: Massively Parallel Production of Long and Near Error-Free Reads Using NGS Technology

Gene sequence analysis model construction based on k-mer statistics

Simple, reference-independent analyses help optimize hybrid assembly of microbial community metagenomes

Exploring Genome Characteristics and Sequence Quality Without a Reference

A Novel Approach for Accurate Sequence Assembly Using de Bruijn graphs