MACSE: Multiple Alignment of Coding SEquences accounting for frameshifts and stop codons.

PLoS ONE. 2011;6(9):e22594 DOI 10.1371/journal.pone.0022594

 

Journal Homepage

Journal Title: PLoS ONE

ISSN: 1932-6203 (Online)

Publisher: Public Library of Science (PLoS)

LCC Subject Category: Medicine | Science

Country of publisher: United States

Language of fulltext: English

Full-text formats available: PDF, HTML, XML

 

AUTHORS

Vincent Ranwez
Sébastien Harispe
Frédéric Delsuc
Emmanuel J P Douzery

EDITORIAL INFORMATION

Peer review

Editorial Board

Instructions for authors

Time From Submission to Publication: 24 weeks

 

Abstract | Full Text

Until now the most efficient solution to align nucleotide sequences containing open reading frames was to use indirect procedures that align amino acid translation before reporting the inferred gap positions at the codon level. There are two important pitfalls with this approach. Firstly, any premature stop codon impedes using such a strategy. Secondly, each sequence is translated with the same reading frame from beginning to end, so that the presence of a single additional nucleotide leads to both aberrant translation and alignment.We present an algorithm that has the same space and time complexity as the classical Needleman-Wunsch algorithm while accommodating sequencing errors and other biological deviations from the coding frame. The resulting pairwise coding sequence alignment method was extended to a multiple sequence alignment (MSA) algorithm implemented in a program called MACSE (Multiple Alignment of Coding SEquences accounting for frameshifts and stop codons). MACSE is the first automatic solution to align protein-coding gene datasets containing non-functional sequences (pseudogenes) without disrupting the underlying codon structure. It has also proved useful in detecting undocumented frameshifts in public database sequences and in aligning next-generation sequencing reads/contigs against a reference coding sequence.MACSE is distributed as an open-source java file executable with freely available source code and can be used via a web interface at: http://mbb.univ-montp2.fr/macse.