Inferring Gene Orders from Gene Maps using the Breakpoint Distance
Résumé
Preliminary to most comparative genomics studies is the annotation of chromosomes as ordered sequences of genes. Unfortunately, di fferent genetic mapping techniques usually give rise to di fferent maps with unequal gene content, and often containing sets of unordered neighboring genes. Only partial orders can thus be obtained from combining such maps. However, once a total order O is known for a given genome, it can be used as a reference to order genes of a closely related species characterized by a partial order P . In this paper, the problem is to fi nd a linearization of P that is as close as possible to O in term of the breakpoint distance. We first prove an NP-complete complexity result for this problem. We then give a dynamic programming algorithm whose running time is exponential for general partial orders, but polynomial when the partial order is derived from a bounded number of genetic maps. A time-effi cient greedy heuristic is then given for the general case, with a performance higher than 90% on simulated data. Applications to the analysis of grass genomes are presented.
Origine : Fichiers produits par l'(les) auteur(s)
Loading...