View source for Tree alignment

In [[computational phylogenetics]], '''tree alignment''' is a [[computational problem]] concerned with producing [[multiple sequence alignment]]s, or alignments of three or more sequences of [[DNA]], [[RNA]], or [[protein]]. Sequences are arranged into a [[phylogenetic tree]], modeling the evolutionary relationships between [[species]] or [[taxa]]. The [[edit distance]]s between sequences are calculated for each of the tree's internal vertices, such that the sum of all edit distances within the tree is minimized. Tree alignment can be accomplished using one of several algorithms with various trade-offs between manageable tree size and computational effort.

==Definition==
'''Input''': A set <math>S</math> of sequences, a [[phylogenetic tree]] <math>T</math> leaf-labeled by <math>S</math> and an [[edit distance]] function <math>d</math> between sequences.

'''Output''': A labeling of the internal vertices of <math>T</math> such that <math>\Sigma_{e \in T} d(e)</math> is minimized, where <math>d(e)</math> is the [[edit distance]] between the endpoints of <math>e</math>.

The task is [[NP-hard]].<ref>{{Citation | last1=Elias | first1=Isaac | title=Settling the intractability of multiple alignment | doi=10.1089/cmb.2006.13.1323 | year=2006 | journal=J Comput Biol | volume=13 | issue=7 | pages=1323–1339 | pmid=17037961| citeseerx=10.1.1.6.256 }}</ref>

==Background==

===Sequence alignment===
{{Main|Sequence alignment}}
[[File:Sequence Alignment Otakuftf&Kingftf.png|thumb|550px|upright|This is a simple Sequence Alignment of Insulin gene between rat, human and chicken. The labeled nucleotides are the different nucleotides with rats I and --- means the missing nucleotides]]
In [[bioinformatics]], the basic method of information processing is to contrast the sequence data. [[Biologist]]s use it to discover the function, structure, and evolutionary information in biological sequences. The following analyses are based on the [[sequence assembly]]: the [[phylogenetic]] analysis, the [[haplotype]] comparison, and the prediction of [[RNA]] structure. Therefore, the efficiency of sequence alignment will directly affect the efficacy of solving these problems. In order to design a rational and efficient sequence alignment, the algorithm derivation becomes an important branch of research in the field of bioinformatics.

Generally, sequence alignment means constructing a string from two or more given strings with the greatest similarity by adding letters, deleting letters, or adding a space for each [[Connection string|string]]. The multiple sequence alignment problem is generally based on pairwise sequence alignment and currently, for a pairwise sequence alignment problem, biologists can use a [[dynamic programming]] approach to obtain its optimal solution. However, the multiple sequence alignment problem is still one of the more challenging problems in bioinformatics. This is because finding the optimal solution for multiple sequence alignment has been proven as an [[NP-complete]] problem and only an approximate optimal solution can be obtained.<ref>L Wang, T Jiang. On the complexity of multiple sequence alignment[J]. Journal of Computational Biology, 194,1(4):337— 34.</ref>

===Distance matrix method===
[[Distance method]] measures the minimum operation number of character [[insertion (genetics)|insertion]]s, [[deletion (genetics)|deletion]]s, and [[Base-pair substitution|substitution]]s that are required to transform one sequence ''u'' to the other sequence ''v'' when being operated on a pair of strings. The calculation of edit distance can be based on [[dynamic programming]], and the equation is in O(|u|×|v|) time, where |u| and |v| are the lengths of u and v.<ref>Yen Hung Chen, On the bottleneck tree alignment problems, INFORMATION SCIENCES; JUN 1, 2010; 180; 11; p2134-p2141</ref> The efficient estimation of edit distance is essential as [[Distance method]] is a basic principle in [[computational biology]]<ref>{{cite journal | last1=Ostrovsky | first1=Rafail | last2=Rabani | first2=Yuval | title=Low distortion embeddings for edit distance | journal=Journal of the ACM | publisher=Association for Computing Machinery (ACM) | volume=54 | issue=5 | date=2007-10-01 | issn=0004-5411 | doi=10.1145/1284320.1284322 | pages=23–es | s2cid=15341956 }}</ref> For functions of hereditary properties "symmetrization" can be used. Due to a series of functions being used to calculate edit distance, different functions may result in different results. Finding an optimal edit distance function is essential for the tree alignment problem.

==The problem of tree alignment==
{{Unreferenced section|date=January 2019}}
[[File:Exponential.svg|thumbnail|right|This figure indicates the growth rate about the exponential time, the polynomial time and the linear time]]
Tree alignment results in a [[NP-hard]] problem, where scoring modes and alphabet sizes are restricted. It can be found as an algorithm, which is used to find the optimized solution. However, there is an exponential relationship between its efficiency and the number sequences, which means that when the length of the sequence is very large, the computation time required to get results is enormously long. Using star alignment to get the approximate optimized solution is faster than using tree alignment. However, whatever the degree of multiple-sequence similarity is, the time complexity of star alignment has a proportional relationship with the square of the sequence number and the square of the sequence average length. As usual, the sequence in MSA is so long that it is also inefficient or even unacceptable. Therefore, the challenge of reducing the time complexity to linear is one of the core issues in tree alignment.

==Combinatorial optimization strategy==
[[Combinatorial optimization]] is a good strategy to solve MSA problems. The idea of combinatorial optimization strategy is to transform the multiple sequence alignment into pair sequence alignment to solve this problem. Depending on its transformation strategy, the combinatorial optimization strategy can be divided into the tree alignment algorithm and the star alignment algorithm. For a given multi-sequence set <math>S</math> ={<math>s_1</math>,...,<math>s_n</math> }, find an [[evolutionary tree]] which has n leaf nodes and establishing one to one relationship between this evolutionary tree and the set <math>S</math>. By assigning the sequence to the internal nodes of the evolutionary tree, we calculate the total score of each edge, and the sum of all edges' scores is the score of the evolutionary tree. The aim of tree alignment is to find an assigned sequence, which can obtain a maximum score, and get the final matching result from the evolutionary tree and its nodes' assigned sequence. Star alignment can be seen as a special case of the tree alignment. When we use star alignment, the evolutionary tree has only one internal node and n leaf nodes. The sequence, which is assigned to the internal node, is called the core sequence.<ref>Serafim Batzoglou. The many faces of sequence alignment[J]. Briefings in Bioinformatics. [http://bib.oxfordjournals.org/content/6/1/6.full.pdf+html 2005,6(1):6—22]</ref>

===The keyword tree theory and the Aho-Corasick search algorithm===
When the [[Combinatorial game theory|combinatorial]] optimization strategy is used to transform the multiple sequence alignment into pair sequence alignment, the main problem is changed from "How to improve the efficiency of multiple sequence alignment" to "How to improve the efficiency of pairwise sequence alignment." The Keyword Tree Theory and the Aho-Corasick search algorithm is an efficient approach to solve the pairwise sequence alignment problem. The aim of combining the keyword tree theory and the Aho-Corasick search algorithm is to solve this kind of problem: for a given long string <math>T</math> and a set of short strings <math>P</math>={<math>p_1</math>,<math>p_2</math>,... ,<math>p_z</math>} (z∈N,z>1), find the location of all <math>P_i</math> in <math>T</math>. The keyword tree produced by set <math>P</math> is used, and then searched for in <math>T</math> with this keyword tree through the Aho-Corasick search algorithm.<ref>Aho A V, Corasick M J. Efficient string matching: an aid to bibliographic search[J]. Communications of ACM, [ftp://163.13.200.222/assistant/bearhero/prog/%A8%E4%A5%A6/ac_bm.pdf 1975,18(6): 333—340]{{Dead link|date=July 2018 |bot=InternetArchiveBot |fix-attempted=no }}.</ref> The total time complexity of using this method to find all <math>P_i</math>'s location in the T is O(<math>m</math>+<math>n</math>+<math>k</math>), where <math>m</math>=|<math>T</math>| (the length of <math>T</math>), <math>n</math>=Σ|<math>P_i</math>| (the sum of all <math>P_i</math>'s lengths) and <math>k</math> means the sum of occurrence for all <math>P_i</math> in <math>T</math>.

====Keyword tree theory====
The keyword tree of the set <math>P</math>={<math>p_1</math>,<math>p_2</math>,... ,<math>p_z</math> }  (z∈N,z>1) is a rooted tree, whose root denoted by K, and this keyword tree satisfies:

(1): Each edge clearly demarcates one letter.

(2): Any two edges separated from the same node are to correspond to different letters.

(3) Each pattern <math>P_i</math> (i=1,2,...,z) corresponds to a node <math>v</math>, and the path from the root K to the node <math>v</math> can exactly correctly spell the string <math>P_i</math>.

For each leaf node of this K tree, it corresponds to one of the certain patterns of set <math>P</math>.

<math>L(v)</math> is used to represent the STRING which is connected from the root node to the node <math>v</math>. <math>Lp(v)</math> will then be used to represent the length of the longest suffix (also, this suffix is the prefix of one of patterns in the set <math>P</math>). Searching this prefix from the root node in the keyword tree, and the last node denoted by <math>n_v</math> when the search is over.{{clarify|date=January 2019|reason=unclear phrasing, "Searching this prefix from the root node in the keyword tree, and the last node denoted by n_v when the search is over. When Lp(v)=0, n_v =K. The ordered pair (v, n_v) is called a failure link." This is unclear, and I don't know what it means.}}<ref>D Gusfield. Algorithms on strings, trees and sequences: computer science and computational biology[M]. [http://www.cambridge.org/us/academic/subjects/computer-science/algorithmics-complexity-computer-algebra-and-computational-g/algorithms-strings-trees-and-sequences-computer-science-and-computational-biology Cambridge: Cambridge University Press.1997.]</ref>

For example, the set <math>P</math>={potato, tattoo, theater, other}, and the keyword tree is shown on the right. In that example, if <math>L(v)</math>=potat, then <math>Lp(v)</math>=|tat|=3, and the failure link of the node <math>v</math> is shown in that figure.

Establishing a failure link is the key to improve the time complexity of the Aho-Corasick algorithm. It can be used to reduce the original polynomial time to the linear time for searching. Therefore, the core of keyword tree theory is to find all failure links (which also means finding all <math>n_v</math>s) of a keyword tree in the linear time. It is assumed that every <math>n_v</math> of all nodes <math>v</math>, whose distance from the root node is less than or equal <math>k</math>, can be found. The <math>n_v</math> of the node <math>v</math> whose distance from the root node is <math>k</math> + 1 can then be sought for. Its parent node is <math>v'</math>, and the letter represented by the node <math>v</math> and <math>v'</math>, is <math>x</math>.

(1): If the next letter of the node <math>n_v' </math> is <math>x</math>, the other node of this edge can be set as <math>w</math>, and <math>n_v</math>=<math>w</math>.

(2): If all letters are not <math>x</math> by searching all edges between <math>n_v' </math> and its child nodes, <math>L(n_v)</math> is a suffix of <math>L(n_v' )</math> plus <math>x</math>. Because this suffix matches the STRING beginning with the root node (similar to prefix), the <math>x</math> after <math>n_v' </math> can be detected or not. If not, this process can be continued until <math>x</math> or the root node is found.

====Aho-Corasick search algorithm====
{{Main|Aho–Corasick string matching algorithm}}
After establishing all failure links in the keyword tree, the Aho-Corasick search algorithm is used to find the locations of all <math>P_i</math> (i=1,2,...,z) in the linear time. In this step, the time complexity is O(m+k).

==Other strategies==
{{Main|Multiple sequence alignment}}
In MSA, DNA, RNA, and proteins, sequences are usually generated and they are assumed to have an evolutionary relationship. By comparing generated maps of RNA, DNA, and sequences from evolutionary families, people can assess conservation of proteins and find functional gene domains by comparing differences between evolutionary sequences. Generally, heuristic algorithms and tree alignment graphs are also adopted to solve multiple sequence alignment problems.

===Heuristic algorithm===
{{Main|Heuristic (computer science)}}
Generally, [[heuristic algorithm]]s rely on the [[iterative]] strategy, which is to say that based on a comparison method, optimizing the results of multiple sequence alignment by the iterative process. Davie M. proposed using the [[particle swarm optimization]] algorithm to solve the multiple sequence alignment problem; Ikeda Takahiro proposed a heuristic algorithm which is based on the [[A* search algorithm]]; E. Birney first proposed using the [[hidden Markov model]] to solve the multiple sequence alignment problem; and many other biologists use the [[genetic algorithm]] to solve it.<ref>RobertC Edgar, Serafim Batzoglou. Multiple sequence alignment[J]. Current Opinion in Structural Biology. [http://robotics.stanford.edu/~serafim/Publications/2006_MSA_COSB.pdf 2006,16(3):368— 373] {{Webarchive|url=https://web.archive.org/web/20131023062245/http://robotics.stanford.edu/~serafim/Publications/2006_MSA_COSB.pdf |date=2013-10-23 }}.</ref><ref>Notredame C, Higgins D.G. SAGA:sequence alignment by genetic algorithm [J]. Nucleic Acids Research. [http://www.tcoffee.org/Publications/Ps_pdf/saga_paper.pdf 1996,24(8):1515-1524].</ref> All these algorithms generally are robust and insensitive to the number of sequences, but they also have shortcomings. For example, the results from the particle swarm optimization algorithm are unstable and its merits depend on the selection of random numbers, the runtime of the A* search algorithm is too long, and the genetic algorithm is easy to fall into local excellent.{{clarify|date=January 2019|reason=what does "local excellent" mean?}}

===Tree alignment graph===
Roughly, tree alignment graph aims to align trees into a graph and finally synthesize them to develop statistics. In biology, tree alignment graphs (TAGs) are used to remove the evolutionary conflicts or overlapping taxa from sets of trees and can then be queried to explore uncertainty and conflict. By integrating methods of aligning, synthesizing and analyzing, the TAG aims to solve the conflicting relationships and partial overlapping [[taxon]] sets obtained from a wide range of sequences. Also, the tree alignment graph serves as a fundamental approach for [[supertree]] and [[grafting]] exercises, which have been successfully tested to construct supertrees by Berry.<ref>Wilkinson M, Pisani D, Measuring support and finding unsupported relationships in supertrees, Systematic Biology 54:823-831.</ref> Because the transformation from trees to a graph contain similar [[node (graph theory)|node]]s and [[Glossary of graph theory#Basics|edges]] from their source trees, TAGs can also provide extraction of original source trees for further analysis. TAG is a combination of a set of aligning trees. It can store conflicting hypotheses evolutionary relationship and synthesize the source trees to develop evolutionary hypotheses. Therefore, it is a basic method to solve other alignment problems.<ref>Stephen A. Smith, Joseph W. Brown, analyzing and synthesizing phylogenies using tree alignment graphs, PLoS Computational Biology 9(9).</ref>

==See also==
* [[Generalized tree alignment]]

==References==
{{reflist|20em}}

{{DEFAULTSORT:Tree Alignment}}
[[Category:Computational phylogenetics]]
[[Category:NP-complete problems]]