Abstract
Background
Recovering the structure of ancestral genomes can be formalized in terms of properties of binary matrices such as the ConsecutiveOnes Property (C1P). The Linearization Problem asks to extract, from a given binary matrix, a maximum weight subset of rows that satisfies such a property. This problem is in general intractable, and in particular if the ancestral genome is expected to contain only linear chromosomes or a unique circular chromosome. In the present work, we consider a relaxation of this problem, which allows ancestral genomes that can contain several chromosomes, each either linear or circular.
Result
We show that, when restricted to binary matrices of degree two, which correspond to adjacencies, the genomic characters used in most ancestral genome reconstruction methods, this relaxed version of the Linearization Problem is polynomially solvable using a reduction to a matching problem. This result holds in the more general case where columns have bounded multiplicity, which models possibly duplicated ancestral genes. We also prove that for matrices with rows of degrees 2 and 3, without multiplicity and without weights on the rows, the problem is NPcomplete, thus tracing sharp tractability boundaries.
Conclusion
As it happened for the breakpoint median problem, also used in ancestral genome reconstruction, relaxing the definition of a genome turns an intractable problem into a tractable one. The relaxation is adapted to some biological contexts, such as bacterial genomes with several replicons, possibly partially assembled. Algorithms can also be used as heuristics for hard variants. More generally, this work opens a way to better understand linearization results for ancestral genome structure inference.
Introduction
Genomes, meant as the linear organization of genes along chromosomes, have been successively modelled by several mathematical objects. Sturtevant and Tan [1] first introduced permutations to study the evolution of genome structure. Starting in the 1980's [2], a large body of work focused on the mathematical and algorithmic properties of such models, including linear and circular genomes [3]. Multichromosomal linear genomes have been defined as generalizations of permutations: they are permutations cut in several pieces [4]. In this framework, hardness results of algorithmic complexity were ubiquitous as soon as three genomes were compared [5,6]. Even worse, if strings were used to model duplications and heterogeneous gene content, then even the basic problem of comparing two genomes proved to be hard [7].
In order to scale up and handle the dozens of available genomes, another model was needed. Bergeron, Mixtacki and Stoye [8] proposed to use a graph matching between gene extremities to define a genome. It simplified the presentation of the DoubleCut and Join (DCJ) theory [9] at the expense of relaxing the model of chromosomal structure as genomes could contain both linear and circular chromosomes. This can be seen as an unrealistic relaxation, as genomes are mostly either linear multichromosomal (eukaryotic nuclear genomes) or circular unichromosomal (bacterial or organelle genomes). But eukaryotes with organelles, or prokaryotes with several replicons, which have not yet been handled explicitly by a formal comparative genomics approach, arguably fit such a model. An unexpected consequence of this relaxation is that the comparison of three genomes with the breakpoint distance proved to be tractable, as an exact optimal median can be computed by solving a maximum weight perfect matching problem [10]. Moreover, the small parsimony problem, i.e., reconstructing the minimum number of evolutionary events along a given phylogeny, can be solved for any number of genomes for the SingleCut and Join (SCJ) distance by Fitch's parsimony algorithm on binary characters [11]. This opened the way to scalable methods at the level of large multispecies datasets.
An additional relaxation consists in allowing any graph, and not only a matching, to model genomes. Ancestral genome reconstruction methods often first compute sets of ancestral adjacencies (neighborhood relations between two genes) [1214], intervals (neighborhood relations between an arbitrary number of genes) [1518], which result in nonlinear structures. This, while also unrealistic at a first glance, allows computational breakthroughs, like incorporating duplications and heterogeneous gene content in the framework [19,20] with polynomial exact methods. Beyond the significant computational speedup, nonlinear genomes may help to understand the amount of error in the data [20].
Nevertheless, biological applications in general require linear genomes, which raises the question of linearizing a collection of adjacencies or intervals. The Linearization Problem is, given a set of weighted intervals (the weight indicates a confidence value based on phylogenetic conservation of intervals), to find a maximum weight subset which is compatible with a linear structure.
According to the definition of a linear structure, this can be described by some variant of the ConsecutiveOnes property (C1P) of binary matrices [10,15,17]: A binary matrix has the C1P if its columns can be ordered such that in each row, there is no 0 entry between two 1 entries. Here, each column is a gene or a gene extremity and each row is an interval. Adjacencies are a particular case of intervals of size two: In that case, the matrix, which has degree 2, can be identified with a graph (vertices are columns and edges are rows). In the case of adjacencies, the Linearization Problem translates to the Maximum Weight VertexDisjoint Path Cover Problem, so it is NPcomplete. A variant handles genomes with a single circular chromosome: A binary matrix has the circular C1P (Ci1P) if its columns can be ordered such that in each row, either there is no 0 entry between two 1 entries, or there is no 1 entry between two 0 entries. For adjacencies, the Linearization Problem contains the Maximum Weight Hamiltonian Cycle Problem, so it is also NPcomplete.
To the best of our knowledge, there is currently no tractability result known for the Linearization Problem. Currently all methods [1214,19] rely on heuristic or external Traveling Salesman Problem solvers, or branch and bound techniques [10,1518]. Moreover, none of the previously published methods is able to infer multichromosomal genomes, possibly with circular chromosomes, which is the natural model for bacterial genomes with plasmids.
In the present paper, we prove that the Linearization Problem for weighted adjacencies, when ancestral genomes can have several circular and linear chromosomes, is tractable. We prove this in a more general case, where multiple copies of columns are allowed. Here, instead of a permutation of the columns, one asks for a sequence on the alphabet of columns, containing at most m(c) occurrences of a column c. In the context of genome reconstruction, this allows to model genes with multiple copies in an ancestral genome [21] or to include telomere markers [22].
We show that this corresponds to finding a maximum weight fmatching, which, in turn, is reducible to finding a maximum weight matching. Also, following the complexity pattern already observed with the model of the C1P with multiplicities [21], we further show that the Linearization Problem for matrices with rows of degrees 2 and 3 is NPcomplete, even if all rows have the same weight and multiplicity one. We discuss the possibilities that our tractability result opens for ancestral genome reconstruction.
Results
A few definitions are needed to prove the two main results of this paper: (1) a polynomial algorithm for the linearization of degree 2 matrices with columns with multiplicity and weighted rows; and (2) an NPcompleteness proof for the linearization of matrices with rows of degrees 2 and 3, even if all multiplicities and row weights are equal to one.
The degree of a row of a binary matrix M (over {0, 1}) is the number of 1 entries in that row. The degree of M is the maximum degree over all its rows. In genomics, the columns of M are the genes, and its rows are the intervals of genes. If a row has degree 2, the interval is called an adjacency. A degree 2 matrix M can be identified with a graph, whose vertices are the columns of M, and edges are the adjacencies. We suppose that all rows are different (and in consequence the graph is simple: it has no multiedges).
A binary matrix (or submatrix) M has the Consecutive Ones Property (C1P) if its columns can be ordered such that in each row, the 1 entries are consecutive (there is no 0 entry between two 1 entries). If M has degree 2, it has the C1P if and only if the corresponding graph is a collection of vertexdisjoint paths. A matrix (or submatrix) M has the Circular Ones Property (Ci1P) if its columns can be ordered such that in each row, either the 1 entries are consecutive (there is no 0 entry between two 1 entries), or the 0 entries are consecutive (there is no 1 entry between two 0 entries); in other words the 1 entries are consecutive when the order of columns is viewed as a circle. If M has degree 2, it has the Ci1P if and only if the graph is a cycle or a collection of vertexdisjoint paths.
Given maximum copy numbers m(c) for each column c, M satisfies the Consecutive Ones Property with multiplicities (mC1P), if there is a sequence S of columns, containing at most m(c) occurrences of column c, and for each row r, the columns containing a 1 in r appear consecutively in S. The mCi1P is defined analogously.
The MAXROWC1P problem takes a binary matrix with weighted rows as input and asks for a subset of rows of maximum cumulative weight with the C1P. For graphs it is equivalent to the Maximum Weight VertexDisjoint Path Cover Problem and thus the MAXROWC1P is NPcomplete [23]. The MAXROWCi1P problem takes a matrix with weighted rows as input and asks for a subset of rows of maximum cumulative weight with the Ci1P. For graphs it can solve the Traveling Salesman Problem and thus the MAXROWCi1P is NPcomplete [23].
These two problems are classical and have been defined independently from comparative genomics, but model well the linearization of genomes with linear chromosomes, or a single circular chromosome, respectively. But the general case would better be modelled by the following. A matrix is componentmCi1P if there is a collection of cyclic sequences of columns that satisfy the following two conditions: (i) for each row r, the columns containing a 1 in r appear consecutively in at least one of the cyclic sequences; and (ii) the total number of occurrences of each column c in all cyclic sequences is at least one and at most m(c). In the particular case where m(c) = 1 for every column c, a matrix is componentmCi1P if its columns can be partitioned such that a row has 1s only in one part and each part is Ci1P. Here chromosomes are sequences, which mean possible ancestral gene orders. It is then the matter of solving the following problem.
MAXROWcomponentmCi1P
Input. A matrix with maximum copy numbers assigned to all columns and weighted rows;
Output. A subset of rows of maximum cumulative weight such that the obtained submatrix is componentmCi1P.
Note that it is equivalent if some sequences are not required to be circular, so it handles well the case where both circular and linear chromosomes are allowed. It is a relaxation of the previous problems, so the NPhardness does not follow from them. And in fact, the problem for degree 2 matrices (adjacencies) happens to be polynomial, as we now show in the next subsection.
A solution for matrices of degree two with weighted rows and multiplicites
For a degree 2 matrix M, let G_{M }be the corresponding graph with a node for each column and a weighted edge for each (weighted) row. Let m : V(G_{M}) → ℕ be the function specifying the maximum copy number of each column, i.e., the multiplicity limit for each vertex of G_{M}. We say that matrix M (resp., the corresponding graph G_{M}) is componentmCi1P for m if there exists a collection of cyclic walks (resp., corresponding cyclic sequences) that satisfies the following two conditions: (i) G_{M }is a subgraph of the union of cyclic walks; and (ii) the total number of occurrences of each vertex v in all cyclic walks is at most m(v).
A 2mmatching of a graph G is a spanning subgraph of G such that the degree of each vertex v ∈ V(G) is at most 2m(v). The following lemma shows the correspondence between spanning subgraphs of G that are componentmCi1P for m and 2mmatchings of G.
Lemma 1 A spanning subgraph of a graph G is componentmCi1P for m if and only if it is a 2mmatching of G.
We give a sketch of the proof. For more details, we refer the reader to [21], where a similar proof was given.
Proof. First, assume a spanning subgraph G' of G is componentmCi1P for m. Then there is a collection of cyclic walks satisfying conditions (i) and (ii). Since each vertex v appears at most m(v) times in these cyclic walks and each occurrence has only two neighbors, the degree of v in G' is at most 2m(v). Hence, G' is a 2mmatching of G.
Conversely, assume G' is a 2mmatching of G. If deg_{G' }(v) < 2m(v) for some v ∈ V(G'), then we add a new vertex v_{0 }and for each v such that deg_{G' }(v) < 2m(v), we add a new edge {v_{0}, v} with multiplicity 2m(v) − deg_{G' }(v) to G'. Since now every vertex of G' has even degree, each component C of G' is Eulerian, i.e., there is a cyclic walk which contains all edges of C, and each v ∈ V(C) appears exactly m(v) times in the walk. If C does not contain v_{0 }then this cyclic walk satisfies conditions (i) and (ii) for vertices in V(C). If C contains v_{0}, then after omitting all occurrences of v_{0 }we obtain a cyclic walk satisfying conditions (i) and (ii) for vertices in V(C). Hence, G' is componentmCi1P for m. QED
It follows that solutions to the MAXROWcomponentmCi1P for matrix M and m correspond to maximum weight 2mmatchings of G_{M}. Next, we give an algorithm for finding a maximum weight fmatching of G with running time (O((V(G) + E(G))^{3/2})), where f : V(G) → ℕ. We will use a more general form of Tutte's reduction for reducing the maximum weight fmatching problem to the maximum weight matching problem similar to the ones presented in [24,25].
Given an edge weighted graph G and function f, construct G' in the following way: For all x in V(G), let x_{1}, x_{2}, ..., x_{f(x) }be in V(G'); and for all e = {x, y}in E(G), let e_{x }and e_{y }be in V(G'). Now, for all e = {x, y} in E(G), let {x_{1}, e_{x}}, ..., {x_{f(x)}, e_{x}}, {e_{x}, e_{y}}, {y_{1}, e_{y}}, ...,{y_{f(y)}, e_{y}} be edges of G', and all these edges have weight w(e). This reduction is illustrated in Figure 1.
Figure 1. Reduction used to transform the maximum weight fmatching problem to the maximum weight matching problem. Edge weights are all one, unless otherwise indicated, and f is given by the white dots inside the nodes. The total edge weight in G is 8. The solid edges show a maximum weight fmatching in G (w = 6), and a corresponding maximum weight matching in G' (of weight 6 + 8 = 14).
Property 1 There is an fmatching in G with weight w if and only if there is a matching in G' with weight w + W, where .
An unweighted version of this property was shown in [25]. The weighted version can be shown in the same way, and hence, we omit the proof.
Since a maximum weight matching can be found in time [26], we have polynomial O((V(G) + E(G))^{3/2}) algorithms for the maximum weight fmatching problem and for the MAXROWcomponentmCi1P problem with multiplicities on matrices of degree 2.
Intractability for matrices of degree larger than two
The tractability does not generalize to matrices, that is, the MAXROWcomponentCi1P is already NPcomplete for unweighted matrices with rows of degrees 2 and 3. Note that the result for unweighted matrices implies NPcompleteness also for the cases when rows are weighted and/or columns have multiplicities.
We will first show that the following hypergraph covering problem is NPcomplete. Here we say that a hypergraph H = (V, E) is 2, 3uniform when all of its hyperedges are either 2edges or 3edges, that is, hyperedges that contain exactly two or three vertices. We will also denote 2,3uniform hypergraphs H = (V, E) as H = (V, E_{2}, E_{3}), where E_{2 }(resp., E_{3}) is its set of 2edges (resp., 3edges). We denote the power set of a set S with (also known as 2^{S}).
Definition 1 A graph covering of a 2, 3uniform hypergraph H = (V, E_{2}, E_{3}) is a graph G = (V, E') such that there exist a map , satisfying the following for every h ∈ E_{2 }∪ E_{3}:
(a) for every h ∈ E, and for every e ∈ c(h), e ⊆ h;
(b) c(h) = 1 if h ∈ E_{2 }and c(h) = 2 if h ∈ E_{3}; and
Here, we say that each set of edges c(h) covers the hyperedge h.
Informally, a graph covering of a 2,3uniform hypergraph is a graph constructed by picking an edge from each 2edge, and a pair of edges from each 3edge.
Problem 1 (The 2,3Uniform Hypergraph Covering by Cycles and Paths by Edge Removal Problem (23UCR Problem)) Given a 2, 3uniform hypergraph H = (V, E) and an integer k, is there a graph covering of H that consists of a collection of disjoint cycles and paths after removing at most k hyperedges from E?
Here we will show that Problem 1, the 23UCR Problem, is NPcomplete. Later in this section we will show that this implies that the MAXROWcomponentCi1P Problem is NPcomplete for matrices with rows of degrees 2 and 3. First, we must define the following NPcomplete version of 3SAT, which we will use to show NPcompleteness of Problem 1.
Problem 2 (The 3SAT(2,3) Problem) Given a CNF formula with the following three properties, is satisfiable?
(a) Formula has only 2clauses and 3clauses.
(b) Each variable x of has exactly two positive occurrences and one negative occurrence in the clauses.
(c) Exactly one positive occurrence of x appears in the 3clauses, while the other two occurrences appear in the 2clauses.
We show that this version of 3SAT is NPcomplete using a very similar proof to the one in [27], by reduction from 3SAT.
Theorem 1 The 3SAT(2,3) Problem is NPcomplete.
Proof. Clearly, the problem is in NP. We now show that it is NPhard by reduction from 3SAT by transforming a given formula that is an instance of 3SAT to a formula that is an instance of 3SAT(2,3) that is satisfiable if and only if is satisfiable. For each variable x of that has k occurrences, we first replace its k occurrences with x^{1}, x^{2},..., x^{k}, i.e., replace the ith occurrence of x (as literal x or ¬x) with x^{i}. We then add the following 2clauses: for i = 1, ..., k  1 and also , where for each i,
This "cycle" of implications (2clauses) on x^{1},..., x^{k}, ensures that for any truth assignment to the variables of , the values of are either all set to true or all set to false. In the first case, the x^{i}'s corresponding to the positive occurrences of x, are set to true and the x^{i}'s corresponding to the negated occurrence of x, are set to false. In the second case, the situation is reversed. Hence, any satisfying truth assignment to the variables of can be translated into a satisfying truth assignment to the variables of , and vice versa, i.e., is satisfiable if and only if is satisfiable. Since it is easy to verify that this transformation can be done in polynomial time, and that is indeed an instance of 3SAT(2,3), it follows that the 3SAT(2,3) Problem is NPcomplete. QED
We now show that the 23UCR Problem is NPcomplete by reduction from 3SAT(2,3).
Theorem 2 The 23UCR Problem is NPcomplete.
Proof. Clearly, the problem is in NP. We will show that it is also NPhard by reduction from 3SAT(2,3).
Given a 3SAT(2,3) formula with variables X = {x_{1}, ..., x_{n}} and set of 2clauses (resp., set of 3clauses), we construct a 2,3uniform hypergraph H_{ϕ }= (V, E). Hypergraph H_{ϕ }is composed of variable gadgets and clause gadgets which contains, among other vertices, a vertex for each literal of (what we will refer to as literal vertices: there are 3n = 2m_{2 }+ 3m_{3 }such vertices). The design of is such that there is a graph covering G of that consists of a collection of disjoint cycles and paths after removing at most m_{2 }+ n edges from E if and only if is satisfiable. For this proof, we call such a G a valid covering of H. Note that a valid covering does not contain any vertex of degree 3 or more.
Figure 2a depicts the variable gadget for variable x ∈ X with its two positive occurrences, labeled as x^{1 }and x^{2}, and its one negative occurrence ¬x in the clauses. We call the 3edge {x'', x''', x'''' } the auxiliary hyperedge, while the other two, {x^{1}, x^{2}, x'} and {¬x, x', x''}, are called the main hyperedges of the variable gadget.
Figure 2. (a) The variable gadget for variable x with literal vertices x^{1}, x^{2 }and ¬x, as well as the four auxiliary vertices x', x'', x''', x'''' that do not appear in any other hyperedge of . (b) 2clause gadget with literal vertices p and q, as well as the two auxiliary vertices c and c' that do not appear in any other hyperedge of H_{ϕ}. (c) 3clause gadget with literal vertices p, q and r.
Figure 2b (resp., 2c) depicts the clause gadget for the 2clause containing literals p, q (resp., and also r for a 3clause). For the 2clause gadget, we call the 2edge {c, c'} the auxiliary hyperedge. We will refer to literals of ϕ and the literal vertices of the gadgets of interchangeably when the context is clear. We have the following claim.
Lemma 2 Formula has a satisfying assignment if and only if has a valid covering.
Proof. "⇒" We first show that a satisfying assignment of can be used to construct a valid covering of .
For the variable gadget corresponding to x ∈ X, we first remove the main hyperedge that contains the literal(s) that is satisfied in the assignment, and then cover the remaining two edges as depicted in Figure 3: Figure 3a (resp., 3b) depicts how to cover the clause gadget when x is false (resp., true) in the assignment. In this figure (as in all remaining figures of this paper), hyperedges drawn with dashed lines are removed, while the straight lines are edges picked in the covering.
Figure 3. Two coverings of the variable gadget for x, when x is set to: (a) false, or (b) true in the assignment.
For a 2clause (resp., 3clause) c containing literals p, q (resp., and also r for a 3clause), without loss of generality let p be a literal that is satisfied in c (there has to be such a literal since it is a satisfying assignment). If c is a 2clause (resp., 3clause), we cover the corresponding gadget as depicted in Figure 4a (resp., 4b).
Figure 4. A covering of the (a) 2clause gadget; and (b) 3clause gadget, where literal p is satisfied.
In the above covering, since exactly m_{2 }+ n hyperedges were removed, and since it is easy to verify that each vertex has degree at most 2, it follows that it is a valid covering of .
"⇐" Now we show that if has a valid covering then is satisfiable.
For hypergraph , we say that a graph G = (V, E') selects a literal vertex for x ∈ X of if x is adjacent to two edges of G in some clause gadget of . Obviously, selected vertices of G correspond to a satisfying truth assignment of if and only if
(i) in every clause gadget, at least one literal vertex is selected, and
(ii) for every x ∈ X, at most one of x and ¬x is selected.
We call a graph G = (V, E') an expected behavior covering of when each variable (resp., clause) gadget of is covered in a way depicted in Figure 3 (resp., 4). It is easy to verify the following observation.
Observation 1 If a valid covering G = (V, E') of is also an expected behavior covering of , then G corresponds to a satisfying truth assignment of ϕ.
In the remainder of this lemma, we will give a set of transformations that converts a valid covering into an expected behavior covering while preserving the validity of the covering at each step. Assume that we have a valid covering of H_{ϕ}.
We say that a variable gadget is undecided in a valid covering of H_{ϕ }if neither of its two main hyperedges is removed. We first show that we can assume that there are no undecided variable gadgets.
Claim 1 We can transform a valid covering of H_{ϕ }into a valid covering that contains no undecided variable gadgets.
Proof. To prove this claim we do a case analysis on the possible configurations that an undecided variable gadget can have in a valid covering of H_{ϕ}, and show how we can locally transform each one to a decided configuration without affecting the validity of the covering of H_{ϕ}.
First, assume that the auxiliary hyperedge is removed in an undecided variable gadget. The set of possible configurations that the gadget can be in is depicted on the left in Figure 5. In this figure (as in all remaining figures) doubleheaded arrows pointing to two vertices in a 3edge represent the two coverings of this 3edge as explained in Figure 6.
Figure 5. The transformation in the case when the auxiliary hyperedge of the variable gadget is removed.
Figure 6. (a) a 3edge with a double arrow pointing to two vertices. (b)(c) the two configurations that are represented by (a).
We can transform any configuration of Figure 5 to the decided configuration on the right. It is easy to see that in the transformed configuration, the number of hyperedges removed is the same as in any initial configuration, and that vertices x', x'', x''' and x'''' (which do not intersect any vertex outside this variable gadget) have degree at most 2. Finally, since each initial configuration of Figure 5 is part of a valid covering of H_{ϕ}, and the degree of any literal vertex (x^{1}, x^{2 }and ¬x) affected by the transformation has only decreased or remained the same, it follows that the covering of H_{ϕ }remains valid after this local transformation.
Hence, we can assume that the auxiliary hyperedge is present in any undecided variable gadget. Without loss of generality we can then assume that any configuration of the undecided variable gadget must be in one of the two forms depicted on the left in Figure 7. In each case, we can perform the corresponding transformation shown in Figure 7. Again it is easy to see that the number of hyperedges removed has not increased, that vertices x', x'', x''', x'''', c and c' (which do not intersect any vertex outside of what is shown here) have degree at most 2, and that degree of any involved literal vertex has not increased. Hence, the covering remains valid. QED
Figure 7. Two sets of possible configurations of an undecided variable gadget and the corresponding transformation of the covering. (Note that if edge {c, c'} is also missing in either initial configuration on the left, that the corresponding configuration on the right still applies, since c and c' do not intersect any vertex outside this variable gadget, and the number of hyperedges removed has not increased.)
We have the following claim.
Claim 2 In any valid covering of H_{ϕ}, at least one hyperedge is removed from each 2clause gadget.
Proof. If no hyperedge is removed from the 2clause gadget (i.e., all hyperedges are covered) in a valid covering of H_{ϕ}, then vertex c (see Figure 2b) has degree 3, which contradicts the fact that this 2clause gadget is part of a valid covering of H_{ϕ}. QED
By Claims 1 and 2, at least n + m_{2 }hyperedges have been removed from the variable and 2clause gadgets, and since in any valid covering this is the maximum number of hyperedges which can be removed, we have the following corollary.
Corollary 1 We can transform a valid covering of H_{ϕ }into a valid covering where:
(a) exactly one hyperedge is removed from each variable gadget and each 2clause gadget of H_{ϕ}, and
(b) no hyperedge is removed from any 3clause gadget.
We have the following claim.
Claim 3 We can transform a valid covering of H_{ϕ }into an expected behavior valid covering.
Proof. Firstly, in the valid covering of H_{ϕ }we can assume, by Claim 1 and Corollary 1, that exactly one main hyperedge is removed and the auxiliary hyperedge is not removed from each variable gadget. However, this does not imply expected behavior. All possible configurations of a decided variable gadget without expected behavior and their corresponding transformations to the expected behavior covering are shown in Figure 8. Analogous to the proof of Claim 1, these local transformations do not affect validity of the covering.
Figure 8. Two sets of possible configurations of a decided variable gadget without expected behavior and the corresponding transformations.
Secondly, in the valid covering of H_{ϕ }we can assume, by Corollary 1, that exactly one hyperedge is removed from each 2clause gadget. Assume now that a 2clause gadget is not expected behavior covered. The only possible such configuration can be transformed to the expected behavior covering as shown in Figure 9. Again, these local transformations do not affect validity of the covering.
Figure 9. The only possible configuration of a 2clause gadget without expected behavior and the corresponding transformation.
Thirdly, in the valid covering of H_{ϕ }we can assume, again by Corollary 1, that the 3clause gadget is covered, and hence it is also expected behavior covered (see Figure 4b). Since all gadgets are expected behavior covered in this valid covering of H_{ϕ}, the claim holds. QED
It follows by Observation 1 and Claim 3, that if H_{ϕ }has a valid covering, then ϕ is satisfiable. This completes the proof of the lemma. QED
Finally, since the construction of H_{ϕ }is polynomial, then by Lemma 2 it follows that the 23UCR Problem is NPcomplete. QED
Let the componentCi1P by Row Removal Problem be the corresponding decision version of the MAXROWcomponentCi1P Problem as follows.
Problem 3 (The componentCi1P by Row Removal Problem) Given a binary matrix M and an integer k, can we obtain a submatrix that is componentmCi1P by removing at most k rows from M?
We now show that the componentCi1P by Row Removal Problem is NPcomplete for matrices with rows of degrees 2 and 3.
The following lemma shows the correspondence between the componentCi1P by Row Removal Problem for matrices with rows of degrees 2 and 3 and the 23UCR Problem. A 2,3uniform hypergraph H = (V, E) can be represented by a binary matrix B_{H }with V columns and E rows, where for each hyperedge h ∈ E, we add a row with 1's in the columns corresponding to the vertices in h and 0's everywhere else. Obviously, there is a onetoone correspondence between 2,3uniform hypergraphs and such matrices.
Lemma 3 A 2,3uniform hypergraph H = (V, E) can be covered by a collection of disjoint cycles and paths after removing at most k hyperedges from E if and only if matrix B_{H }has the componentCi1P after removing at most k rows.
Proof. Assume first that H has a covering G that consists of a collection of disjoint cycles and paths after removing at most k hyperedges from E. Remove the (at most k) rows from B_{H }that correspond to the hyperedges removed from E. Each path (resp., cycle) O of G defines a cyclic order on its set of vertices. Consider the cyclic ordering of the columns of each component of B_{H }corresponding to O. It is easy to see that each such cyclic ordering is a Ci1P ordering of its corresponding component, and hence B_{H }has the componentCi1P after removing at most k rows.
Conversely, assume that each component C = {v_{1},...,v_{C}} of the submatrix of B_{H }obtained by removing at most k rows is Ci1P with respect to cyclic order of its columns. Consider the following covering G of H, after removing the (at most k) hyperedges from E that correspond to the rows removed from B_{H}: for every hyperedge, pick the edge between two adjacent columns/vertices in π. Note that every picked edge is for some j, or . Hence, G consists of a collection of disjoint cycles and paths. QED
By Theorem 2 and Lemma 3 it follows that the componentCi1P by Row Removal Problem is NPcomplete for matrices with rows of degrees 2 and 3. Since this decision problem is NPcomplete, it follows that the MAXROWcomponentCi1P Problem is also NPcomplete for matrices with rows of degrees 2 and 3.
Theorem 3 The MAXROWcomponentCi1P Problem is NPcomplete.
Discussion/Conclusion
There are exact optimization [12,19,20] or empirical [1315] fast methods to construct ancestral adjacencies which do not necessarily form a linear signal. But to date, all linearization methods were heuristics or calls to Traveling Salesman Problem solvers [1214,19]. Moreover, no method is currently adapted to reconstruct bacterial ancestral genome with plasmid(s), while this situation is common in the living world.
We report here two results: (1) a polynomial variant of the Linearization Problem, when the output allows paths and cycles and a maximum number of copies per gene, in the case of degree 2 matrices with weighted rows; and (2) an NPcompleteness proof of the same problem for matrices with rows of degrees 2 and 3, even when multiplicities and weights are equal to one.
It is not the first time that a slight change in the formulation of a problem dramatically changes its computational status [10]. Even if such a relaxation is less realistic in certain contexts, solving the relaxation can also help to approach efficiently the constrained problem, like for DCJ and inversions/translocations for example [28,29]. More generally, 2factors (spanning subgraphs composed of collections of vertexdisjoint cycles) have been used to approximate Traveling Salesman solutions [24], so genomes composed of several circular chromosomes can be a way to approximate solutions for linear ones.
Moreover, considering genomes composed of linear and circular segments is appropriate for bacterial genomes where linear segments can be seen as segments of not totally recovered circular chromosomes. Currently no ancestral genome reconstruction method is able to handle bacterial genomes with plasmids, but rather they are restricted to eukaryotes or bacterial chromosomes with a single circular chromosome. For example, Darling et al. [30] reconstruct the ancestral genomes of Yersinia pestis strains but are limited to the main chromosome by their method, while there are 3 plasmids in most current species, and they are of capital importance since they are suspected to have provoked the pathogenicity of the plague agent. So it is crucial to include them in evolutionary studies, which justifies our model for future biological studies.
Furthermore, genes are often duplicated in genomes, and in the absence of a precise and efficient phylogenetic context, which is still absent for bacteria (no ancestral genome reconstruction method is able to handle horizontal transfers for example), a multicopy family translates into a multiplicity in the problem statements.
The ability to obtain such genomes in polynomial time from adjacencies also opens interesting perspectives for phylogenetic scaffolding of extant bacterial genomes [31] or more generally bacterial communities [32].
These applications are left as a future work.
Competing interests
The authors declare that they have no competing interests.
Authors' contributions
JM, MP, RW, CC and ET formalized and solved the linearization problems and wrote the paper.
Acknowledgements
We thank Jens Stoye for useful discussions. JM and CC are funded by NSERC Discovery Grants. MP is funded by a Marie Curie Fellowship from the Alain Bensoussan program of ERCIM. ET is funded by the Agence Nationale pour la Recherche, Ancestrome project ANR10BINF0101.
This article has been published as part of BMC Bioinformatics Volume 13 Supplement 19, 2012: Proceedings of the Tenth Annual Research in Computational Molecular Biology (RECOMB) Satellite Workshop on Comparative Genomics. The full contents of the supplement are available online at http://www.biomedcentral.com/bmcbioinformatics/supplements/13/S19
References

Sturtevant A, Tan C: The comparative genetics of Drosophila Pseudoobscura and Drosophila Melanogaster.
Journal of Genetics 1937, 34:415432. Publisher Full Text

Watterson GA, Ewens WJ, Hall TE, Morgan A: The chromosome inversion problem.
Journal of Theoretical Biology 1982, 99:17. Publisher Full Text

Fertin G, Labarre A, Rusu I, Tannier E, Vialette S: Combinatorics of genome rearrangements.

Hannenhalli S, Pevzner P: Transforming men into mice (polynomial algorithm for genomic distance problem).
36th Annual Symposium on Foundations of Computer Science, IEEE Comput. Soc. Press, Los Alamitos, CA 1995, 581592.

Bryant D: The complexity of the breakpoint median problem. In Tech Rep CRM2579. Centre de Recherches Mathématiques, Université de Montréal; 1998.

Caprara A: The reversal median problem.
INFORMS Journal on Computing 2003, 15:93113. Publisher Full Text

Blin G, Chauve C, Fertin G, Rizzi R, Vialette S: Comparing Genomes with Duplications: a Computational Complexity Point of View.
ACM/IEEE Transactions on Computational Biology and Bioinformatics 2007, 4:523534. PubMed Abstract  Publisher Full Text

Bergeron A, Mixtacki J, Stoye J: A Unifying View of Genome Rearrangements.
Algorithms in Bioinformatics, Proceedings of WABI'06, Volume 4175 of Lecture Notes in Computer Science 2006, 163173.

Yancopoulos S, Attie O, Friedberg R: Efficient sorting of genomic permutations by translocation, inversion and block interchange.
Bioinformatics 2005, 21:33403346. PubMed Abstract  Publisher Full Text

Tannier E: Yeast Ancestral Genome reconstruction: the possibilities of computational methods.
Comparative Genomics, Proceedings of RECOMBCG'09, Volume 5817 of Lecture Notes in Computer Science 2009, 112. PubMed Abstract

Feijao P, Meidanis J: SCJ: a breakpointlike distance that simplifies several rearrangement problems.
IEEE/ACM Transactions on Computational Biology and Bioinformatics 2011, 8:13181329. PubMed Abstract  Publisher Full Text

Ma J, Zhang L, Suh B, Raney B, Burhans R, Kent W, Blanchette M, Haussler D, Miller W: Reconstructing contiguous regions of an ancestral genome.
Genome Research 2006, 16:15571565. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Bertrand D, Gagnon Y, Blanchette M, ElMabrouk N: Reconstruction of Ancestral Genome subject to Whole Genome Duplication, Speciation, Rearrangement and Loss.
Algorithms in Bioinformatics, Proceedings of WABI'10, Volume 6293 of Lecture Notes in Bioinformatics 2010, 7889.

Muffato M, Louis A, Poisnel CE, Crollius HR: Genomicus: a database and a browser to study gene synteny in modern and ancestral genomes.
Bioinformatics 2010, 26:11191121. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Chauve C, Tannier E: A methodological framework for the reconstruction of contiguous regions of ancestral genomes and its application to mammalian genomes.
PLoS Computational Biology 2008, 4:e1000234. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Stoye J, Wittler R: A Unified Approach for Reconstructing Ancient Gene Clusters.
IEEE/ACM Trans Comput Biol Bioinf 2009, 6(3):387400. PubMed Abstract  Publisher Full Text

Chauve C, Gavranović H, Ouangraoua A, Tannier E: Yeast ancestral genome reconstructions: the possibilities of computational methods II.
Journal of Computational Biology 2010, 17:10971112. PubMed Abstract  Publisher Full Text

Jones BR, Rajaraman A, Tannier E, Chauve C: ANGES: Reconstructing ANcestral GEnomeS maps.
Bioinformatics 2012., 18 PubMed Abstract  Publisher Full Text

Ma J, Ratan A, Raney BJ, Suh BB, Zhang L, Miller W, Haussler D: DUPCAR: reconstructing contiguous ancestral regions with duplications.
Journal of Computational Biology 2008, 15:10071027. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Bérard S, Gallien C, Boussau B, Szollosi G, Daubin V, E T: Evolution of gene neighborhood within reconciled phylogenies.
Bioinformatics 2012. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Wittler R, Maňuch J, Patterson M, Stoye J: Consistency of sequencebased gene clusters.
Journal of Computational Biology 2011, 18(9):10231039. PubMed Abstract  Publisher Full Text

Chauve C, Maňuch J, Patterson M, Wittler R: Tractability results for the ConsecutiveOnes Property with multiplicity.
Combinatorial Pattern Matching, Proceedings of CPM'11, Volume 6661 of Lecture Notes in Computer Science 2011, 90103.

Garey M, Johnson D: Computers and Intractability: A Guide to the Theory of NPcompleteness. W. H. Freeman & Co; 1979.

Lovasz L, Plummer MD: Matching Theory, Volume 29 of Annals of Discrete Mathematics. North Holland; 1986.

Dessmark A, Lingas A, Garrido O: On Parallel Complexity of Maximum fmatching and the Degree Sequence Problem. In Proceedings of the 19th International Symposium on Mathematical Foundations of Computer Science 1994. MFCS '94, SpringerVerlag; 1994:316325.

Micali S, Vazirani VV: An algorithm for finding maximum matching in general graphs.

Papadimitriou C: Computational Complexity. Addison Wesley; 1994.

Miklós I, Tannier E: Bayesian sampling of genomic rearrangement scenarios via double cut and join.
Bioinformatics 2010, 26:30123019. PubMed Abstract  Publisher Full Text

Miklós I, Tannier E: Approximating the number of double cutanjoin scenarios.

Darling AE, Miklós I, Ragan MA: Dynamics of genome rearrangement in bacterial populations.
PLoS Genetics 2008, 4:e1000128. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Husemann P, Stoye J: Phylogenetic Comparative Assembly.
Algorithms for Molecular Biology 2010, 5:3. PubMed Abstract  BioMed Central Full Text  PubMed Central Full Text

Pop M: Genome assembly reborn: recent computational challenges.
Briefings in Bioinformatics 2009, 10:354366. PubMed Abstract  Publisher Full Text  PubMed Central Full Text