Abstract:Background and objective Short-read sequencing is becoming the standard of practice for the study of structural variants associated with disease. However, with the growth of sequence data largely surpassing reasonable storage capability, the biomedical community is challenged with the management, transfer, archiving, and storage of sequence data.Methods We developed Hierarchical mUlti-reference Genome cOmpression (HUGO), a novel compression algorithm for aligned reads in the sorted Sequence Alignment/Map (SAM) format. We first aligned short reads against a reference genome and stored exactly mapped reads for compression. For the inexact mapped or unmapped reads, we realigned them against different reference genomes using an adaptive scheme by gradually shortening the read length. Regarding the base quality value, we offer lossy and lossless compression mechanisms. The lossy compression mechanism for the base quality values uses k-means clustering, where a user can adjust the balance between decompression quality and compression rate. The lossless compression can be produced by setting k (the number of clusters) to the number of different quality values.Results The proposed method produced a compression ratio in the range 0.5-0.65, which corresponds to 35-50% storage savings based on experimental datasets. The proposed approach achieved 15% more storage savings over CRAM and comparable compression ratio with Samcomp (CRAM and Samcomp are two of the state-of-the-art genome compression algorithms). The software is freely available at https://sourceforge.net/projects/hierachicaldnac/with a General Public License (GPL) license.Limitation Our method requires having different reference genomes and prolongs the execution time for additional alignments.Conclusions The proposed multi-reference-based compression algorithm for aligned reads outperforms existing single-reference based algorithms.

Research on Multi-genome Index and Its Improved Sequence Alignment Algorithm

Gene Sequence Alignment on a Public Computing Platform

Indexing large genome collections on a PC

DNA sequences alignment method using sparse index on pan-genome graph

A novel fast multiple nucleotide sequence alignment method based on FM-index

Multiple-Searching Adaptive Genetic Algorithm for Multiple Sequence Alignment

Short Read Alignment Based On Maximal Approximate Match Seeds

Improvement in accuracy of multiple sequence alignment using novel group-to-group sequence alignment algorithm with piecewise linear gap cost

A Short Reads Alignment Algorithm Oriented To Massive Data

An algorithm for rapid noncoding RNA sequence-structure alignment

GPU Accelerated Biological Sequence Alignment

Investigating Memory Optimization of Hash-Index for Next Generation Sequencing on Multi-Core Architecture

Improvement of Performance of MegaBlast Algorithm for DNA Sequence Alignment

A heuristic algorithm for multiple sequence alignment base on progressive multiple alignment

Gene Position Index Mutation Detection Algorithm Based on Feedback Fast Learning Neural Network

AliBI: An Alignment-Based Index for Genomic Datasets

A local alignment of DNA based on parallelized MUMmer algorithm

Accelerating whole-genome alignment in the age of complete genome assemblies

High Performance Multiple Sequence Alignment Algorithms for Comparison of Microbial Genomes

HUGO: Hierarchical Multi-Reference Genome Compression for Aligned Reads.

GMAP and GSNAP for Genomic Sequence Alignment: Enhancements to Speed, Accuracy, and Functionality