Ant Colony Optimization based Founder Sequence Reconstruction

7 61 0
Ant Colony Optimization based Founder Sequence Reconstruction

Đang tải... (xem toàn văn)

Thông tin tài liệu

In this paper we propose an ant colony optimization algorithm (ACO) based method, equipped with some important improvements, for the founder DNA sequence reconstruction problem. The proposed method yields excellent performance when validating on 108 test sets from three benchmark datasets. Comparing with the best by far corresponding method, our proposed method performs better in 45 test sets, equally well in 44 and worse only in 19 sets. These experimental results demonstrate the efficacy and perspective of our proposed method.

VNU Journal of Science: Comp Science & Com Eng., Vol 33, No (2017) 59-65 Ant Colony Optimization based Founder Sequence Reconstruction Anh Vu Thi Ngoc1, Dinh Phuc Thai2, Hoang Duc Nguyen2, Thanh Hai Dang2,∗, Dong Do Duc2 1The 2Faculty Hanoi college of Industrial Economics of Information Technology, VNU University of Engineering and Technology Abstract Reconstruction of a set of genetic sequences (founders) that can combine together to form given genetic sequences (e.g DNA) of individuals of a population is an important problem in evolutionary biology Such reconstruction can be modeled as a combinatorial optimization problem, in which we have to find a set of founders upon that genetic sequences of the population can be generated using a smallest number of recombinations In this paper we propose an ant colony optimization algorithm (ACO) based method, equipped with some important improvements, for the founder DNA sequence reconstruction problem The proposed method yields excellent performance when validating on 108 test sets from three benchmark datasets Comparing with the best by far corresponding method, our proposed method performs better in 45 test sets, equally well in 44 and worse only in 19 sets These experimental results demonstrate the efficacy and perspective of our proposed method Received 11 Sep 2017; Revised 31 Dec 2017; Accepted 31 Dec 2017 Keywords: Founder sequence reconstruction (FSR), Ancestor genes, Ant colony optimization (ACO) * To this end, the main challenge is at the problem of determining the plausible number of founder (ancestor) sequences and of finding themselves for a given finite offspring sequences It is well known as the founder sequence reconstruction problem Various methods have been recently proposed for reconstructing founder sequences, such as those based on dynamic programming [2], tree search [3], neighboring search [4] and metaheuristics [5] In this paper we propose a ant colony optimization (ACO) based method for the founder sequence reconstruction problem The manuscript is structured as follows: • Section first formulates the problem of founder sequence reconstruction and Section then presents related works that have been Introduction Today we have been observing a huge amount of biological sequences (e.g DNA/genes, proteins) steadily being generated thanks to the unprecedentedly fast development of bio-technologies Having genetic sequences of a population, researchers are often interested in the evolution history of the population, which can be traced back by re-constructing such given sequences from a small number of not-yet identified ancestors (namely founder sequences) using some genetic operators Many biological studies have demonstrated the efficacy of this approach [1] * Corresponding author E-mail.: hai.dang@vnu.edu.vn https://doi.org/10.25073/2588-1086/vnucsce.170 59 A.V.T Ngoc et al / VNU Journal of Science: Comp Science & Com Eng., Vol 33, No (2017) 59-65 60 successfully applied to the problem with good results reported • Our proposed algorithm, experimental results and comparisons with previously proposed state-of-the-art related methods are described in Section • Section gives some conclusions for the proposed method It also suggests some potential follow-ups to improve the method further Problem statement length m defined over a finite set S , i.e., Ci = Ci1 , Ci , with Cij  S (which can be A, C, G, T if recombinants of interest are DNA sequences), we need to find a set of k founder sequences F = ( F1 , F2 , , each of length m defined over the set S A set F is considered valid if the set of recombinants C can be reconstructed from F This means that, each recombinant Ci can be decomposed into pi components (  pi  m ) Fr , Fr ,  , Fr Founder Sequences Reconstruction Problem (FSRP) is defined as follows: Given a set of n recombinants C = (C1 , C2 ,  , Cn ) , each Ci is a sequence of i1 i2 ip so that each piece Fr ( j = 1,2,  , pi ) appears at ij least once at the same position as in Ci K L Figure Haloptye sequences as recombinants, which are supposed to be originated from a set of predefined founder sequences using a decomposition with breakpoints A valid decomposition is considered reducible if two consecutive pieces not appear in the same founder sequence Among such reducible ones the FSRP aims to find out the optimal decompositions with a minimum number of required breakpoints The number of breakpoints for a solution F can be calculated using the formula: n i =1 i  p m In this paper we consider a common biological application in that each recombinant is a haplotype sequence, i.e S = {0,1} , where and are the two possible common alleles On the left side of Figure is an example of a set C of haplotype sequences, which is presented in form of a matrix In the middle part is a valid founder sequences ( a , b and c ) assuming that the number of founder sequences is set to The optimal decomposition with breakpoints on the recombinants into sections, which are part of the founder sequences, is shown on the right-hand side Breakpoints are marked with vertical bars The FSRP was first introduced by Ukkonen [2] and has been proven NP-Hard [6] with k > Related work This section introduces two state-of-the-art algorithms proposed for the FSR problem, namely Recblock [3] and LNS [4], which have achieved excellent results on benchmark datasets 3.1 RecBlock algorithm RecBlock [3] is a FSR algorithm based on tree search Given k founder sequences each of length m , the algorithm encodes them as a matrix with k rows and m columns RecBlock A.V.T Ngoc et al / VNU Journal of Science: Comp Science & Com Eng., Vol 33, No (2017) 59-65 reviews the columns of the matrix from left to right Vertex Vl at the depth l of the search tree is part of a solution for the prefix part of the founders till the column l Each vertex Vl is labeled with a number of breakpoints BP( Vl ) in the process of reconstructing recombinants by far Recblock uses some strategies to speed up the reconstruction: • Only consider the founder sequences in the alphabet order to avoid revisiting permutations • A vertex is not extended further if its breakpoint number greater than that of the best solution so far Given two vertices Vl and Vl at the depth 61 solution found in the current episode is used to learn (tune  ) and go for the next turn Our proposed method for FSR has input and output as follows: Input: binary matrix C of size n * m representing a recombinant set and k is the number of the founder sequences to be found Output: binary matrix F of size k * m string representing the founder sequences so that BP (C , F ) is minimal Here, BP (C , F ) is the number of breakpoints required to obtain C from F In general, our ACO based method for FSR works as depicted in Algorithm 1: and respectively, if l1 l2 , BP(Vl )  BP(Vl )  n (where n is the of number of recombinants), we may ignore Vl for downstream analysis 3.2 Large neighborhood search algorithm LNS-1c is empirically considered the best algorithm proposed by far for solving the FSR problem [4] This algorithm uses the nearestneighbor search strategy over a large neighborhood of constructed solutions During searching the neighborhood, the algorithm picks out a set F free  F beforehand, then uses the algorithm Recblock to search for alternative founder sequences in FF free Whenever a better solution is found out, LNS1c performs local search over neighborhood from scratch Proposed method 4.1 Ant colony optimization based FSR Ant colony optimization [7] (ACO) is a metaheuristic method simulating how ants in nature find paths from their nest to food sources, which turn out to be a reinforcement learning method ACO solves optimization problems throughout many episodes, in each of which every ant travels to find solutions based on heuristic information and pheromone matrix  containing information learned The best 4.2 Structure graph for the FSR problem For the sake of visualization, we simulate the FSR problem as the problem of finding paths on a corresponding structure graph (see Figure 2) This structure graph includes a start, an end node and m columns Each column has k vertices, of which each corresponds to a state of the corresponding column in the matrix F of founder sequences In particularly, each state is a binary string of length k Each vertex has edges connecting to all ones in the next column We can see all paths starting from the start to the end node has to go through every column once, at which one state is chosen Each journey of ants travelling from the start to the end node therefore corresponds to a complete matrix of founder sequences 4.3 How ants travel on the structure graph When travelling on the structure graph, ants chose a next vertex to visit at random The 62 A.V.T Ngoc et al / VNU Journal of Science: Comp Science & Com Eng., Vol 33, No (2017) 59-65 algorithm is described in pseudo code in Algorithm ?? The probability at which a vertex is chosen is proportional to its level of compatibility to the matrix constructed by ants so far This level is calculated through heuristic and pheromone information  Particularly, the j vertex in the column i will be visited by an ant with a probability Pi , j = [ i , j ] [ a , j ] [ i ,l ] [ a ,l ] l Where: •  a, j is the heuristic value (see 4.3.1) a, j = BP(Ci , Fa  j ) where: • Ci is the matrix of the first i columns of matrix C • Fa is the solution that ant a has built (with i  columns) • Fa  j is the matrix resulted when ant a intends to visit vertex j To give an example, when i = we have the structure graph as in Figure •  i, j is the pheromone information (see 4.3.2) •  ,  are two parameters of an ACO determining the correlation between the heuristic value and the pheromone information Figure Structure graph when i = 4.3.2 Pheromone information In the FSR problem, we denote  ij as the 4.3.1 Heuristic information While constructing the optimal solution, heuristic information is calculated according to the level of compatibility to the matrix that is yielded with the next moves of ants In more details, when an ant is going to the j vertex in the column i the heuristic information is calculated as follows pheromone information of the j th vertex in the column i in the graph Vertices being visited in the optimal solutions found in every searching phase by ants so far will be learnt such that they are of high priority to be visited in next phases There are various pheromone updating methods that have been proposed for ACO We select the Smoothed Max-Min Ant system [8] because it yields the best results in our experiments In this regard, the pheromone information is updated after each loop as follows:  ij = (1   ) ij   ij where:   if (i, j )  T  ij =    maxif (i, j )  T and T is the optimal solution that ants found after the loop and (i, j ) is the vertex j in the column i of the structure graph 4.4 Improved ACO for FSRP Figure Structure graph for the ACO-based founder sequence reconstruction 4.4.1 Ants find solutions synchronously Note that the problem solution space is extremely large, if working independently with A.V.T Ngoc et al / VNU Journal of Science: Comp Science & Com Eng., Vol 33, No (2017) 59-65 each other ants could hardly to concentrate on potential regions of the searching space We therefore propose a search strategy for ants as follows: We let ants (in the set Ants) find solutions in parallel When moving to the next column, instead of letting each ant choose the next vertex to go, we create a new ant set (called NewAnts) to prolong paths created by ants in the set Ants In particular, if an ant a prolongs the path for an ant a , it means that ant a will go over the similar journey as ant a before moving to the next vertex in the next column When having NewAnts with the same size as Ants, we move to the next column and repeat such a new ant set building procedure from NewAnts until having a complete solution set This procedure is depicted in pseudo code in Algorithm For more details, when going from the column i  to the column i , each ant a  NewAnts will randomly choose an ant a  Ants to prolong its path and a vertex j in the column i to move forward The ant a is chosen with a probability also based on the heuristic and pheromone information, as follows: Pa , j = [ i , j ] [ a , j ] [ ax l i ,l ] [ a ,l ] x 4.4.2 Other improvements Neighborhood search: To lower the probability of missing good solutions while searching, we recommend using the reduced version of the algorithm RecBlock (3.2) to find other better solutions within the vicinity of the best by far solution found by ants Instead of browsing the whole founder sequences, for each founder in 63 the optimal solution found by far we use RecBlock to find another alternative better one Searching along two dimensions: With the newly proposed search strategy, ants will quickly converge onto some solution regions, leading to a low diversity of found solutions To improve this problem, apart from searching forward from the start to the end vertex, we also let ants search backward along the opposite direction (i.e from the end back to start vertex) The search direction is periodically changed When searching backward, the complete different heuristic information is used, leading to the potential of finding new solutions Experimental results We compare our proposed FSR algorithm called ACOFSRP with the best corresponding one by far, i.e LNS-1c [4] on benchmark data sets, namely rnd (random), evo and ms (each contains test set) All sequences in the first data set is randomly generated while those in the two latter ones are generated according to evolutionary models All three are used in the study of LNS-1c We experiments with the founder sequence length k  5,6,7,8,9,10 for each of such test sets, leading to a total of 108 tests We also experiments with different variants of ACOFSRP by not using either one of two improvements or both on the same three benchmark sets Experimental results show that ACOFSRP outperforms its two variants, demonstrating the power of two proposed improvements in ACOFSRP (data not shown) Due to the random nature of ACOFSRP, we perform each test 20 times and the run time of each is limited to 10 hours These numbers are and 72, respectively, in the study of LNS-1c [4] The program is run on a CPU with 12GB RAM and 4GHz processor Table ?? shows the detailed performance, in terms of the solution quality (number of required breakpoints) and the running time, of ACOFSRP and LNS-1c on three benchmark data sets Note that the values for ACOFSRP are the averages of those from 20 running times A.V.T Ngoc et al / VNU Journal of Science: Comp Science & Com Eng., Vol 33, No (2017) 59-65 64 Table Detailed performance of our ACOFSRP and LNS-1c on three benchmark sets # founders 10 10 10 10 10 10 ACOFSRP LNS-1c ACOFSRP LNS-1c Value Time(s) Value Time(s) Value Time(s) Value Time(s) rnd-30_60 evo-30_60 48427 145 3996 145 372 4501 372 44255 5394 324 5695 324 94 94 53 293 7644 289 8136 906 65 65 86 96096 12502 263 12361 268 45 45 353 240 22388 246 175659 36 27293 36 51 90559 36041 221 34456 229 28 28 rnd-30_90 evo-30_90 72903 203 6222 585 6753 585 203 60 516 79754 118 7491 514 8501 118 52 55418 12225 461 12506 472 69 69 19 07173 20652 417 19270 426 43 43 382 31562 399 12679 35 35383 35 69 36056 353 36055 370 244167 31 31 28 rnd-30_150 evo-30_150 976 11244 976 134777 381 10419 381 893 858 14045 865 216875 230 13178 230 72 766 20532 778 140918 131 21422 131 72 30531 698 31618 710 250463 63 63 59 639 36054 666 87405 39 36071 39 38 36120 591 36094 619 21046 35 12 rnd-50_100 evo-50_100 8644 1211 9290 1213 65968 368 368 145 1084 12766 1097 60881 250 12072 250 113 985 20193 1009 8769 174 21207 174 14706 44145 123 34994 124 910 31773 928 149 845 36063 875 113792 99 36061 99 2507 84 36128 794 36098 830 221118 83 3696 rnd-50_150 evo-50_150 1797 14459 1800 195873 522 12464 522 132 1606 19572 1622 144474 319 19894 319 109 1466 31384 1484 221180 205 33503 205 1354 36044 1385 85140 135 36059 135 169 1262 36130 1320 222181 101 36116 101 108 83 36174 1194 36122 1240 244166 82 291 rnd-50_250 evo-50_250 3031 26742 3043 101246 1126 21491 1126 3060 2698 34085 2725 172785 726 29774 726 1060 2461 36056 2508 251951 450 36042 450 259 2276 36090 2330 176486 258 36072 258 603 2133 36137 2204 244380 141 36186 141 12100 85 36269 2012 36256 2097 257557 83 275 ACOFSRP LNS-1c Value Time(s) Value Time(s) ms-30_60 4520 124 124 209 100 98859 99 5871 17273 81 7194 81 70 54798 69 11135 60 59 17377 2002 38579 50 33364 50 ms-30_90 8933 167 167 747 136 10240 136 768 114 12369 114 30934 97 126402 96 16197 85 83 32062 216 36057 74 73 1648 ms-30_150 252 11476 251 4986 189 16279 189 1421 154 24401 153 25361 125 32750 125 7590 103 36050 103 106022 36118 88 88 22794 ms-50_100 310 12258 310 2192 251 16089 251 18039 210 25576 212 442 177 34846 178 51495 156 36056 155 38758 138 36137 137 30080 ms-50_150 430 18911 429 48449 346 25681 346 26957 287 30661 286 1958 240 36047 241 130741 201 36072 203 170493 175 36120 174 8253 ms-50_250 615 23672 613 2171 482 33887 479 48013 396 36050 396 16430 338 36076 336 23916 288 36121 283 243608 257 36228 248 7413 A.V.T Ngoc et al / VNU Journal of Science: Comp Science & Com Eng., Vol 33, No (2017) 59-65 On the random data set ( rnd ), ACOFSRP could procedure solutions better than LNS-1c for 32 among total 36 cases On-par solutions are observed in the remaining cases Regarding the running time, ACOFSRP requires shorter time than LNS-1c for 32 cases while longer only for remaining cases On the data set evo , ACOFSRP is beated by LNS-1c in terms of excution time for all cases Nevertheless, solutions yielded by ACOFSRP are on-par with those of LNS-1c for 32 out of 36 cases For the remaining cases, the solution goodness scores by ACOFSRP are worse than those by LNS-1c (The small differences are observed, i.e up to breakpoints) On the data set ms , ACOFSRP produced solutions are better than and equal to those yielded by LNS-1c for 12 and 10 cases, respectively Interestingly, among such 22, ACOFSRP requires remarkably shorter runing time than LNS-1c for 12 cases For the remaining 14 cases, ACOFSRP produce solutions worse than LNS1c ./table_combine_all.tex Conclusion Founder gene sequence reconstruction (FSR) for a given population can be modeled as a combinatorial optimization problem, which has been proven NP-hard In this paper we propose a novel method based on ant colony optimization algorithms (ACO) coupled with two other important improvements (i.e local search and back forward search) to solve the founder gene sequence reconstruction problem Experiments on the benchmark data sets show better or equal results for almost sets when comparing to the best corresponding method, demonstrating the efficacy and future perspectives of our proposed method G g 65 Acknowledgments This work has been supported by Vietnam National University, Hanoi (VNU), under Project No QG.15.21 References [1] G Tyson, J Chapman, H Philip, E Allen, R Ram, P M Richardson, V Solovyev, E M Rubin, D Rokhsar, J F Banfield, Community structure and metabolism through reconstruction of microbial genomes from the environment, Nature 428 (2004) 37–43 [2] E Ukkonen, Finding Founder Sequences from a Set of Recombinants, Springer Berlin Heidelberg, Berlin, Heidelberg, 2002, pp 277–286 [3] A Roli, C Blum, Tabu Search for the Founder Sequence Reconstruction Problem: A Preliminary Study, Springer Berlin Heidelberg, Berlin, Heidelberg, 2009, pp 1035–1042 [4] A Roli, S Benedettini, T Stützle, C Blum, Large neighbourhood search algorithms for the founder sequence reconstruction problem, Computers Operations Research 39 (2) (2012) pp 213–224 [5] C Blum, A Roli, Metaheuristics in combinatorial optimization: Overview and conceptual comparison, ACM Comput Surv 35 (3) (2003) 268–308 [6] P Rastas, E Ukkonen, Haplotype inference via hierarchical genotype parsing, in: Proceedings of the 7th International Conference on Algorithms in Bioinformatics, WABI’07, Springer-Verlag, Berlin, Heidelberg, 2007, pp 85–97 [7] M Dorigo, T Stützle, Ant Colony Optimization, Bradford Company, Scituate, MA, USA, 2004 [8] D Do Duc, H Hoang Xuan, Smooth and threelevels ant systems: Novel aco algorithms for solving traveling salesman problem, in: Ad Cont to the International Conference: IEEERIVF 2010, pp 33–37 ... alternative founder sequences in FF free Whenever a better solution is found out, LNS1c performs local search over neighborhood from scratch Proposed method 4.1 Ant colony optimization based FSR Ant colony. .. novel method based on ant colony optimization algorithms (ACO) coupled with two other important improvements (i.e local search and back forward search) to solve the founder gene sequence reconstruction. .. letting each ant choose the next vertex to go, we create a new ant set (called NewAnts) to prolong paths created by ants in the set Ants In particular, if an ant a prolongs the path for an ant a ,

Ngày đăng: 29/01/2020, 23:53

Tài liệu cùng người dùng

  • Đang cập nhật ...

Tài liệu liên quan