AnchoRNA: Full virus genome alignments through conserved anchor regions

Abstract

Background: Multiple sequence alignment of full viral genomes can be challenging due to factors such as long sequences, large insertions/deletions (spanning several 100 nucleotides), large number of sequences, sequence divergence, and high computational complexity in particular when computing alignments based on RNA secondary structures. Standard alignment methods often face these issues, in particular when processing highly variable sequences or when specific phylogenetic analysis is required on selected subsequences. Results: We present an algorithm to determine high quality anchors that define partitions of sequences and guide the alignment of viral genomes to respect well conserved, and therefore functionally significant, regions. This new approach is implemented in the Python-based command line tool AnchoRNA, which is designed to identify conserved regions, or anchors, within coding sequences. By default, anchors are searched in translated coding sequences accounting for high mutation rates in viral genomes. AnchoRNA enhances the accuracy and efficiency of full-genome alignment by focusing on these crucial conserved regions. AnchoRNA anchors are compared with those produced by three other anchor-finding tools, and anchor-guided alignments are evaluated using three alignment programs. Utilizing a dataset of 55 representative Pestivirus genomes, AnchoRNA identified 55 anchors that are used for guiding the alignment process. The incorporation of these anchors led to improvements across tested alignment tools, highlighting the effectiveness of AnchoRNA in enhancing alignment quality, especially in viral genomes. Conclusions: We demonstrate that AnchoRNA effectively improves the alignment of complete viral genomes by identifying conserved anchor regions that guide multiple sequence alignments.