Long-Read Genome Assembly and Polishing Strategies

Introduction

The advent of third-generation sequencing technologies has enabled the generation of reads exceeding tens of kilobases in length, fundamentally transforming genome assembly [1, 2]. Unlike short-read platforms that produce highly accurate but short sequences, long-read platforms yield reads with substantially higher error rates (typically 5–15% per base) but provide the contiguity necessary to resolve repetitive regions, structural variants, and complex genomic architectures [3, 4]. In veterinary medicine, long-read assembly is critical for constructing complete genomes of livestock pathogens, zoonotic agents, and host species, facilitating accurate annotation of virulence factors, antimicrobial resistance genes, and mobile genetic elements [5, 1]. The assembly process typically involves two major phases: initial contig construction from raw long reads and subsequent polishing to correct residual errors [6, 7]. This article provides an exhaustive review of long-read genome assembly algorithms and polishing strategies, with emphasis on methodologies applicable to microbial and veterinary genomes.

Long-Read Assembly Algorithms

Overlap-Layout-Consensus (OLC) Approaches

Long-read assemblers predominantly employ the overlap-layout-consensus (OLC) paradigm, which is well suited for reads with high error rates and long lengths [8]. In OLC, all pairwise read overlaps are computed, a string graph is constructed, and consensus sequences are derived from the graph. The Miniasm assembler implements a minimal OLC approach that bypasses error correction, relying instead on raw read overlaps and subsequent consensus generation with tools such as Racon [8]. The FALCON assembler extends OLC with a hierarchical correction step and produces phased diploid assemblies through the FALCON-Unzip algorithm [9]. NextDenovo employs an efficient error correction module prior to OLC assembly, achieving high accuracy on noisy nanopore long reads [10]. Verkko2 integrates proximity-ligation data with long-read de Bruijn graphs to achieve telomere-to-telomere assembly, demonstrating substantial improvements in contiguity and phasing [11].

De Bruijn Graph and Hybrid Approaches

While de Bruijn graphs are standard for short-read assembly, their application to long reads is complicated by high error rates. However, hybrid assemblers combine short-read accuracy with long-read contiguity. The SPAdes assembler incorporates long reads via a hybrid mode that uses short-read de Bruijn graphs and long-read repeat resolution [3]. The Unicycler pipeline automates hybrid assembly for bacterial genomes, using short reads to correct long-read assemblies [1]. LRScaf and SSPACE-LongRead are scaffolding tools that leverage long reads to order and orient contigs produced by short-read assemblers [12, 13]. Gap filling is addressed by LR_Gapcloser, which uses raw long reads to close assembly gaps with high efficiency and low memory usage [14], and by gapFinisher, which processes SSPACE-LongRead output to fill gaps [15].

Multi-Sample and Haplotype-Aware Assembly

For population-scale studies, multi-sample assembly pipelines such as LORA enable simultaneous assembly of multiple genomes, leveraging shared haplotypes to improve contiguity [16]. StrainCascade provides an automated modular workflow for high-throughput bacterial genome reconstruction, integrating long-read assembly with characterization of plasmids and mobile elements [5]. Haplotype-aware error correction methods, such as those described by Barak et al., improve consensus accuracy in diploid or polyploid genomes by distinguishing allelic variants during correction [17].

Polishing Strategies

Polishing is the process of correcting base-level errors in a draft assembly. Errors in long-read assemblies are predominantly insertions and deletions (indels) in homopolymer regions and GC-rich sequences [7, 2]. Polishing strategies fall into two categories: self-polishing using only long reads, and hybrid polishing using short reads.

Self-Polishing with Long Reads

Self-polishing tools use the same long-read data to correct the assembly. Racon performs partial order alignment of reads against the assembly and generates a consensus [8]. Medaka uses a neural network model trained on specific sequencing platform error profiles to predict corrections [7]. NextPolish is a fast k-mer-based polisher that scores and counts k-mers from short reads but can also be used in a self-polishing mode with long reads [18]. Homopolish specifically targets homopolymer errors and has shown superior performance on nanopore assemblies [7, 2]. The PEPPER tool uses a deep learning model for margin polishing and is effective when combined with Medaka [7]. Evaluation of individual tools on microbial genomes indicates that Homopolish, PEPPER, and Medaka yield the best results, but no single tool is universally optimal [7].

Hybrid Polishing with Short Reads

Hybrid polishing leverages the high accuracy of short reads to correct long-read assemblies. Pilon is a widely used tool that aligns short reads to the assembly and identifies discrepancies [19]. Polypolish uses a conservative approach that only makes corrections supported by multiple short-read alignments, minimizing false positives [19]. Pypolca offers default and careful modes, with the careful mode avoiding false-positive errors at low coverage [19]. NextPolish outperforms Pilon in speed and correction accuracy when using short-read data [18]. The depth of short-read coverage significantly affects polishing performance: most benefits are achieved by 25x depth, and Polypolish-careful introduces no false-positive errors at any depth [19]. For very low depth (<5x), Polypolish-careful alone is recommended; for low depth (5–25x), a combination of Polypolish-careful and Pypolca-careful is optimal; for sufficient depth (>25x), Polypolish-default and Pypolca-careful are preferred [19].

Repeat-Aware and Iterative Polishing

Standard polishing tools can overcorrect in repetitive regions, introducing new errors. The telomere-to-telomere consortium developed a repeat-aware polishing strategy that uses a diverse panel of sequencing technologies to correct errors without overcorrection [6]. This approach improved the assembly quality value from 70.2 to 73.9 in the CHM13 human genome [6]. Iterative polishing, where multiple rounds of correction are applied, can further reduce error rates. For example, a second round of Homopolish after initial polishing improves accuracy [7]. However, iterative polishing must be carefully monitored to avoid introducing systematic biases [6].

Evaluation and Validation

Assessing assembly quality requires multiple metrics. BUSCO (Benchmarking Universal Single-Copy Orthologs) evaluates gene completeness [7]. Assembly quality value (QV) estimates per-base accuracy based on k-mer counts [6]. CloseRead is a tool that visualizes local assembly quality and diagnoses errors in structurally complex regions such as immunoglobulin loci, enabling targeted re-assembly [20]. Validation strategies include aligning short reads to the assembly and counting discordant k-mers, as well as comparing to reference genomes when available [6, 1]. For microbial genomes, gene prediction with Prokka and comparison to known sequences provides additional validation [7].

Workflow Diagram

The following Mermaid diagram illustrates a typical long-read genome assembly and polishing workflow.

flowchart TD
    A[Raw Long Reads], > B[Error Correction / Self-Correction]
    B, > C[Overlap-Layout-Consensus Assembly]
    C, > D[Draft Assembly]
    D, > E{Polishing Strategy}
    E, > F[Self-Polishing: Racon, Medaka, NextPolish]
    E, > G[Hybrid Polishing: Short Reads + Pilon, Polypolish, Pypolca]
    F, > H[Iterative Polishing]
    G, > H
    H, > I[Repeat-Aware Polishing]
    I, > J[Validation: BUSCO, QV, CloseRead]
    J, > K{Errors Detected?}
    K, >|Yes| L[Targeted Re-assembly / Gap Filling]
    L, > D
    K, >|No| M[Final Polished Assembly]

Frequently Asked Questions

What is the difference between self-polishing and hybrid polishing?

Self-polishing uses only long reads to correct the assembly, while hybrid polishing incorporates short reads to achieve higher accuracy, particularly in homopolymer regions [7, 19].

Which polishing tool is best for bacterial genomes?

The optimal tool depends on coverage and error profile. For nanopore-only assemblies, Homopolish, PEPPER, and Medaka perform well; for hybrid polishing, Pypolca-careful and Polypolish-careful are recommended [7, 19].

How much short-read coverage is needed for effective polishing?

Most polishing benefits are achieved at 25x short-read coverage; lower depths require conservative tools like Polypolish-careful to avoid false positives [19].

Can long-read assembly achieve telomere-to-telomere quality?

Yes, with advanced assemblers like Verkko2 and repeat-aware polishing, telomere-to-telomere assemblies are achievable for complex genomes [6, 11].

What are common errors in long-read assemblies?

The most common errors are indels in homopolymer runs and GC-rich regions, which can be addressed by specialized polishers like Homopolish [7, 2].

References

[1] Goldstein S, Beka L, Graf J, et al. Evaluation of strategies for the assembly of diverse bacterial genomes using MinION long-read sequencing. BMC Genomics. 2018. URL: https://www.semanticscholar.org/paper/edae24cff5f31d607a51012591831ab4703e4a0a

[2] Zhang P, Jiang D, Wang Y, et al. Comparison of De Novo Assembly Strategies for Bacterial Genomes. International Journal of Molecular Sciences. 2021. URL: https://www.semanticscholar.org/paper/7bb7227bf5c9b1828f20bfe44a4fb112c67b1370

[3] Liao YC, Lin SH, Lin HH. Completing bacterial genome assemblies: strategy and performance comparisons. Scientific Reports. 2015. URL: https://www.semanticscholar.org/paper/e681946228ec41ac44fb602f160c5382a3cc5ba1

[4] Lin HH, Liao YC. Evaluation and Validation of Assembling Corrected PacBio Long Reads for Microbial Genome Completion via Hybrid Approaches. PLoS ONE. 2015. URL: https://www.semanticscholar.org/paper/fa3bca72f71a74d6ff7b69292125baa42aca7df2

[5] Jordi SBU, Baertschi I, Li J, et al. StrainCascade: An automated, modular workflow for high-throughput long-read bacterial genome reconstruction and characterization. iScience. 2026. URL: https://pubmed.ncbi.nlm.nih.gov/42305602/

[6] Mc Cartney AM, Shafin K, Alonge M, et al. Chasing perfection: validation and polishing strategies for telomere-to-telomere genome assemblies. Nature Methods. 2021. URL: https://www.semanticscholar.org/paper/0c72c2614ccaaecab60e201312e605c002a69d2b

[7] Lee JY, Kong M, Oh JS, et al. Comparative evaluation of Nanopore polishing tools for microbial genome assembly and polishing strategies for downstream analysis. Scientific Reports. 2021. URL: https://www.semanticscholar.org/paper/bb772664b24c23d6ef4e5b9080fe678e797099a2

[8] Vaser R, Sovic I, Nagarajan N, et al. Fast and accurate de novo genome assembly from long uncorrected reads. bioRxiv. 2016. URL: https://www.semanticscholar.org/paper/9d62188a0c9895b4c8b6d9c8cbb76945ee1c2a7a

[9] Chin C, Peluso P, Sedlazeck F, et al. Phased Diploid Genome Assembly with Single Molecule Real-Time Sequencing. Nature Methods. 2016. URL: https://www.semanticscholar.org/paper/9af8db5ae66e5009b4deb8bcdc4beb93752e0e74

[10] Hu J, Wang Z, Sun Z, et al. NextDenovo: an efficient error correction and assembly tool for noisy long reads. Genome Biology. 2024. URL: https://www.semanticscholar.org/paper/493b019728735250b3d8c2cd6007276fabe71991

[11] Antipov D, Rautiainen M, Nurk S, et al. Verkko2 integrates proximity-ligation data with long-read De Bruijn graphs for efficient telomere-to-telomere genome assembly, phasing, and scaffolding. Genome Research. 2025. URL: https://www.semanticscholar.org/paper/793db66111b2edea42dbe3afff7edeeb206ac0c1

[12] Boetzer M, Pirovano W. SSPACE-LongRead: scaffolding bacterial draft genomes using long read sequence information. BMC Bioinformatics. 2014. URL: https://www.semanticscholar.org/paper/2399ca37840b4d6569d7f5de4f1a5ff499c89311

[13] Qin M, Wu S, Li A, et al. LRScaf: improving draft genomes using long noisy reads. BMC Genomics. 2018. URL: https://www.semanticscholar.org/paper/5a207c56224bd9d59a49936161d14e235cc6e141

[14] Xu GC, Xu TJ, Zhu R, et al. LR_Gapcloser: a tiling path-based gap closer that uses long reads to complete genome assembly. GigaScience. 2018. URL: https://www.semanticscholar.org/paper/4f5781c9031f14de77e4c13ca7fbaebef1a917d9 *** Disclaimer: This article is for educational and informational purposes only. It is not intended to substitute for professional veterinary advice, diagnosis, treatment, or regulatory guidance. Always consult a licensed veterinarian or qualified specialist regarding animal health, disease diagnosis, and therapeutic decisions.

[15] Kammonen JI, Smolander OP, Paulin L, et al. gapFinisher: A reliable gap filling pipeline for SSPACE-LongRead scaffolder output. PLoS ONE. 2019. URL: https://www.semanticscholar.org/paper/6ccd544c75307dc4b8d34306124a8ef54b707f2b

[16] Desvillechabrol D, Ouazahrou R, da Fonseca JP, et al. LORA: a polymorphic multi-sample long read assembly pipeline. NAR Genom Bioinform. 2026. URL: https://pubmed.ncbi.nlm.nih.gov/42344949/

[17] Barak P, Gibney D, Jain C. Haplotype-aware long-read error correction. Algorithms Mol Biol. 2026. URL: https://pubmed.ncbi.nlm.nih.gov/42374584/

[18] Hu J, Fan J, Sun Z, et al. NextPolish: a fast and efficient genome polishing tool for long-read assembly. Bioinformatics. 2019. URL: https://www.semanticscholar.org/paper/24c676f9011cf20b9ca590ffd5c55844612252d4

[19] Bouras G, Judd LM, Edwards RA, et al. How low can you go? Short-read polishing of Oxford Nanopore bacterial genome assemblies. bioRxiv. 2024. URL: https://www.semanticscholar.org/paper/d34c4c6073cb2b4c91ee3bede821b3e103c842ac

[20] Zhu Y, Watson CT, Safonova Y, et al. CloseRead: a tool for assessing assembly errors in immunoglobulin loci applied to vertebrate long-read genome assemblies. Genome Biology. 2025. URL: https://www.semanticscholar.org/paper/c03d06e2cbe5ac0badca663bf9d7067031e90991