MOSAL: software tools for multiobjective sequence alignment
© Paquete et al.; licensee BioMed Central Ltd. 2014
Received: 3 September 2013
Accepted: 5 January 2014
Published: 8 January 2014
Multiobjective sequence alignment brings the advantage of providing a set of alignments that represent the trade-off between performing insertion/deletions and matching symbols from both sequences. Each of these alignments provide a potential explanation of the relationship between the sequences. We introduce MOSAL, a software tool that provides an open-source implementation and an on-line application for multiobjective pairwise sequence alignment.
KeywordsPairwise sequence alignment Multiobjective optimization Dynamic programming
Sequence alignment is in the core of many bioinformatics applications. It aims to identify regions of similarity in sequences of biological data, such as nucleotide and amino acid residues. The procedure consists of inserting gaps between the residues so that similar symbols from several sequences become aligned. For two sequences, dynamic programming algorithms can compute the optimal alignment in an efficient manner . However, for very large DNA or protein databases, heuristic approaches like FASTA and BLAST have been used [2, 3]. See  for an extensive review from a computational point of view.
Any of these approaches rely on the a priori definition of coefficients that are assigned to the components of the score function. These weights are usually defined by default in most of the software packages for sequence alignment and are usually not modified by the practitioner. However, there is a considerable disagreement about how to weight each coefficient. A small change in the weights can lead to a completely different alignment.
One way of overcoming the problem of setting weights is to consider a multiobjective formulation, where the practitioner is provided a set of optimal alignments representing the trade-off between components of the score function, for instance, substitution score given by a substitution matrix and the number of gaps; in this case, an alignment is optimal if there is no other alignment with better substitution score value and lesser number of gaps. Usually, there is not only one optimal alignment but several for which this notion of optimality holds; such set of all optimal alignments is called the Pareto optimal alignment set.
Under a multiobjective formulation, no weights are needed to be set up. Moreover, according to a classical result in the multiobjective optimization field , this optimal set contains not only all of the optima of a weighted sum formulation, but also many other alignments that are not possible to find at all by the weighted sum approach. Each of these alignments can be seen as a potential explanation of the relationship between the sequences and may be of interest for the practitioner for a more in-depth analysis. In fact, several other problems in bioinformatics have been already reformulated from a multiobjective point of view .
A multiobjective approach to pairwise sequence alignment has been explored by several researchers, both from a problem formulation and algorithmic point of view [7–11]. Recently, it has been applied to the construction of phylogenetic trees, which has shown to provide complementary information to that obtained by common methods .
MOSAL is a software tool that results from the problem formulation given in  with the aim of providing an open-source implementation and an on-line application where this implementation can be tested. The web-server is available at http://mosal.dei.uc.pt and physically located at the Department of Informatics Engineering, University of Coimbra, and is one of the outcomes of a national funded research project on multiobjective sequence alignment.
The code is written in C and provided under a GNU General Public License. A makefile is available for compilation under GNU/Linux. The implementation can be setup for several multiobjective score functions as described in : maximization of the number of matches or substitution score and minimization of gaps or indels.
Speed-up techniques described in  are also implemented and can be parameterized, in particular, the maximum size of the lower bound set for the pruning technique. This parameter should be defined with some care; if too small, the pruning has a reduced effect, and if too large, a excessive number of comparisons may reduce the advantage of pruning in terms of CPU-time. For most of the benchmarks tested, a value of 10 seems to be the most appropriate .
Command line options
Path to the 1st sequence file (FASTA)
Path to the 2nd sequence file (FASTA)
Use indels or gaps
dp | dpp -b=N
Do not use or use pruning technique.
If yes, specify the size of lower bound (N)
Use substitution score instead of matches
(F is the path to the matrix file)
Output only the scores without the alignments
- Step 1:
Insertion of each sequence in FASTA format in a text box. The user can choose either Protein or DNA type of sequence in a switch button.
- Step 2:
Choice of the score function with switch buttons. The user can choose either matches or substitution score for the first score function component and either indels or gaps for the second score function component. If substitution score is chosen, the user can choose a substitution score matrix (PAM 100, 250 and BLOSUM 62, 75, 80, 85 if Protein option is chosen in the previous step) or can even provide one in a predefined text format.
- Step 3:
Choice of the sequence alignment options: with or without the alignments and with or without pruning technique. If pruning is chosen, the number of bounds must be provided (10 is given by default). The option without alignment provides only the score function values of the alignments.
- Step 4:
Submit to the server, with the option of sending an e-mail to the user with the output files.
A visualization tool in the on-line application allows to visualize all the alignments and the corresponding score function values produced by the implementation or by the on-line application. The coloring scheme used in the Sequence Manipulation Suite (see http://www.bioinformatics.org/sms2/) is also applied here to allow the identification of potential regions of interest in the several alignments.
MOSAL provides a set of tools for the practitioner to perform a more in-depth analysis on the relation between a pair of biological sequences. The multiobjective formulation that is explored by the framework provides further insight into the confidence of the alignments obtained by common methods; for instance, a large number of optimal scores suggests that a single alignment may be insufficient to understand the relation between the sequences and that further investigation is required. Moreover, the output can be used to construct phylogenetic trees as suggested in .
This work was support by the Fundação para a Ciência e Tecnologia, project MOSAL - Multiobjective sequence alignment (PTDC/EIA-CCO/098674/2008) and by FEDER, Programa Operacional Factores de Competitividade do QREN, ref. COMPETE: FCOMP-01-0124-FEDER-010024.
- Needleman SB, Wunsch CD: A general method applicable to the search for similarities in the amino acid sequence of two proteins. J Mol Biol. 1970, 48 (3): 443-453. 10.1016/0022-2836(70)90057-4.View ArticlePubMedGoogle Scholar
- Altschul S, Gish W, Miller W, Myers E, Lipman D: A basic local alignment search tool. J Mol Biol. 1990, 215: 403-410.View ArticlePubMedGoogle Scholar
- Lipman D, Pearson W: Rapid and sensitive protein similarity searches. Science. 1985, 227: 1435-1441. 10.1126/science.2983426.View ArticlePubMedGoogle Scholar
- Gusfield D: Algorithms on Strings, Trees, and Sequences. Computer Science and Computational Biology. 1997, New York: Cambridge University PressView ArticleGoogle Scholar
- Ehrgott M: Multicriteria optimization. 2005, Berlin: SpringerGoogle Scholar
- Handl J, Kell DB, Knowles JD: Multiobjective Optimization in Bioinformatics and Computational Biology. IEEE/ACM Trans Comput Biol Bioinform. 2007, 4 (2): 279-292.View ArticlePubMedGoogle Scholar
- Roytberg M, Semionenkov M, Tabolina O: Pareto-optimal alignment of biological sequences. Biophysics. 1999, 44 (4): 565-577.Google Scholar
- Taneda A: Multi-objective pairwise RNA sequence alignment. Bioinformatics. 2010, 26 (19): 2383-2390. 10.1093/bioinformatics/btq439.View ArticlePubMedGoogle Scholar
- Abbasi M, Paquete L, Liefooghe A, Pinheiro M, Matias P: Improvements on bicriteria pairwise sequence alignment: algorithms and applications. Bioinformatics. 2013, 29 (8): 996-1003. 10.1093/bioinformatics/btt098.View ArticlePubMedGoogle Scholar
- DeRonne K, Karypis G: Pareto optimal pairwise sequence alignment. IEEE/ACM Trans Comput Biol Bioinform. 2013, 10 (2): 481-493.View ArticlePubMedGoogle Scholar
- Schnattinger T, Schöning U, Kestler H: Structural RNA alignment by multi-objective optimization. Bioinformatics. 2013, 29 (13): 1607-1613. 10.1093/bioinformatics/btt188.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.