Skip to main content
  • Methodology article
  • Open access
  • Published:

LTC: a novel algorithm to improve the efficiency of contig assembly for physical mapping in complex genomes

Abstract

Background

Physical maps are the substrate of genome sequencing and map-based cloning and their construction relies on the accurate assembly of BAC clones into large contigs that are then anchored to genetic maps with molecular markers. High Information Content Fingerprinting has become the method of choice for large and repetitive genomes such as those of maize, barley, and wheat. However, the high level of repeated DNA present in these genomes requires the application of very stringent criteria to ensure a reliable assembly with the FingerPrinted Contig (FPC) software, which often results in short contig lengths (of 3-5 clones before merging) as well as an unreliable assembly in some difficult regions. Difficulties can originate from a non-linear topological structure of clone overlaps, low power of clone ordering algorithms, and the absence of tools to identify sources of gaps in Minimal Tiling Paths (MTPs).

Results

To address these problems, we propose a novel approach that: (i) reduces the rate of false connections and Q-clones by using a new cutoff calculation method; (ii) obtains reliable clusters robust to the exclusion of single clone or clone overlap; (iii) explores the topological contig structure by considering contigs as networks of clones connected by significant overlaps; (iv) performs iterative clone clustering combined with ordering and order verification using re-sampling methods; and (v) uses global optimization methods for clone ordering and Band Map construction. The elements of this new analytical framework called Linear Topological Contig (LTC) were applied on datasets used previously for the construction of the physical map of wheat chromosome 3B with FPC. The performance of LTC vs. FPC was compared also on the simulated BAC libraries based on the known genome sequences for chromosome 1 of rice and chromosome 1 of maize.

Conclusions

The results show that compared to other methods, LTC enables the construction of highly reliable and longer contigs (5-12 clones before merging), the detection of "weak" connections in contigs and their "repair", and the elongation of contigs obtained by other assembly methods.

Background

Until very recently, genome sequencing projects such as the human ([1, 2]), mouse ([3], rice ([4]), or maize genome projects ([5, 6]) have relied on the construction of physical maps as a framework to support BAC-by-BAC or whole-genome shotgun sequencing [7]. Alternatively, genome sequencing could be conducted via whole genome shotgun approach ([8–11] etc.) or using the novel next generation sequencing technologies (e.g., [12–14]). Physical maps can be established via BAC clones fingerprinting using restriction enzyme profiling [15–22] or by digital fingerprinting [23–25], and subsequent assembly of the clones into contigs based on the systematic comparison of fingerprint profiles. BAC contigs are then ordered using molecular markers and genetic or radiation hybrid maps. In addition, to providing a framework for sequencing, such maps can be used for high-resolution gene mapping [26–28] and map-based gene cloning [29–34]. Despite significant progress in fingerprinting techniques such as High Information Content Fingerprinting (HICF) [22] and the development of efficient programs such as FPC (FingerPrinted Contigs) [35, 36], maximal likelihood-based reconstruction of physical map [37, 38], random cost algorithm minimizing triplewise linking distance [39], etc. to automatically cluster clones into contigs, physical mapping remains long, laborious, and expensive especially for large and complex genomes that contain a high amount of repeated sequences (e.g., maize or wheat genome). Therefore, the development of algorithms and methods making this process more cost effective is important in view of the increasing amount of non-model species that will be sequenced in the near future.

The basis of contig assembly is that the same DNA fragments in different clones are cut by a given restriction enzyme at the same sites. Hence, the presence of fragments with the same length in the fingerprints of two clones, c1 and c2, indicates a possible overlap between these clones. However, in large genomes, the abundance of repeated elements and the limited accuracy of scoring the band lengths may lead to the identification of shared bands for two clones that originate from different parts of a chromosome, thereby reducing the reliability of contig assembly. Thus, contig assembly relies on the identification of significant overlaps. This implies the calculation of p-values of clone overlap for any pair of clones c1 and c2. Namely, the p-value is the probability for two random clones to have the same or a higher number of shared bands by chance. The exact calculation of the p-value is usually problematic (reviewed in Wendl, 2005 [40]). The FPC software uses the Sulston approximation [17] that is based on the simplest model of taking into account tolerance (the accuracy of fragment length scoring), and the assumption that appearances of fragments (bands) of different lengths are independent and identically distributed (iid assumption). This approximation is valid in cases where a small number of clones are matching but can be very inaccurate in situations with an intermediate to large number of common bands [40, 41]. Moreover, it was observed that different bands can have very different abundances within a fingerprint database. To take this variation into account Nelson et al. (2005) [6] proposed to exclude the most abundant bands for reducing the proportion of false overlaps. Other ways of taking into account band frequencies based on Bayes Theorem were also proposed (e.g., [42, 43]), but such methods are cumbersome, especially with HICFs (e.g., [44]).

The FPC program package assembles clones into contigs based on fingerprints generated by either the end-labeled double digest method [15, 18] or the complete digest method [16, 45]. Because of technical difficulties related to the large amount of clones to order, FPC divides the clones into subsets of relatively small contigs in which clones are supposed to be highly significantly overlapping. Clones are ordered using local optimization and building band maps. To achieve a treatable size and ensure high accuracy of contigs, FPC users usually employ a very high initial threshold (cutoff) for the p-value of clone overlaps. However, for many clone pairs that do overlap physically, the p-value may not overcome such a cutoff. As a result, numerous short sub-contigs and singletons are produced calling for subsequent merging. In fact, the ordering of bands and highly overlapping clones in short contigs is questionable. The merging of short sub-contigs and singletons is also problematic.

Additional difficulties can be caused by the presence of "questionable" (chimerical or poorly fingerprinted) clones (referred to as Q-clones). The presence of Q-clones and false clone overlaps can result in a wrong clone order or even in the assembly of clones deriving from different parts of the genome into the same contig. The presence of bands with similar lengths also make a clone ordering complex thereby hampering the correct map assembly [46] and resulting in low-quality minimum tiling paths (MTP). In particular, unexpected gaps can arise when a MTP is checked by sequencing (e.g., via BAC end sequencing). Although many physical maps have been constructed with the standard FPC algorithm and successfully employed for genome sequencing, quite a lot of errors in contig assembling were also found [47]. The diversity of factors affecting the map quality in different situations calls for the development of new methods and tools complementing FPC and other existing packages. This problem is especially important for physical mapping of complex genomes with a high level of repeats, such as in wheat [48] and barley.

In this paper, we present a novel approach coordinating clone clustering and ordering. We also propose to use a new metric of clone overlap instead of the standard Sulston score. In contrast to FPC, the Linear Topology Contig (LTC, available upon request from the corresponding author) program starts clustering with a relatively relaxed cutoff and uses the topology of significant clone overlap to obtain longer contigs with a realistic (linear) structure. In each cluster, clones are ordered based on a global optimization procedure, and clones that disturb the order stability (assessed by re-sampling analysis) are excluded from the contig. Ordered contigs are then merged with a relaxed cutoff into longer contigs using the network representation of the significant clone overlaps as a control of the contig topology (similar ideas were used by Waterman et al., 1986 [49]; Cuticchia et al., 1992 [50]; Zhang et al., 1994 [51] and others). In addition to contig building, LTC can be used for verification, repairing, and elongation of contigs obtained by other methods (e.g., FPC).

The reliability of the proposed methodology was assessed with HICF data from the wheat 3B physical map [48]. The performance of LTC vs. FPC was compared also on the simulated BAC libraries based on the known genome sequences for chromosome 1 of rice and chromosome 1 of maize. The results demonstrate that contigs built by LTC are longer, better ordered, and more robust to errors caused by false and missing bands than those obtained by FPC and can make the MTP selection more effective, leading to more reliable physical maps and increased sequencing cost efficiency.

Methods

I. Contigs Construction

In LTC, the contig construction algorithm includes the following steps: (i) calculation of p-values for clone overlaps; (ii) temporary exclusion from the analysis of the clones and clone overlaps unproved by parallel paths; (iii) adaptively changing cutoff clustering; (iv) "non-linear" cluster splitting into sub-clusters with linear topological structure; (v) global optimization ordering; (vi) verification of the orders by re-sampling; (vii) sub-contig merging into contigs; and (viii) MTP construction.

(i) Calculation of p-values for clone overlaps

In the first step, LTC calculates all pair-wise p-values Pr(c1, c2) of clone overlaps and selects threshold Pr0 (cutoff) to declare clones c1 and c2 with Pr(c1, c2) < Pr0 as overlapping clones. A proper choice of the threshold Pr0 should provide a reasonable trade-off between two requirements: (a) providing a sufficient number of pairs of overlapping clones, and (b) reducing the proportion of false overlaps among selected clone pairs. Instead of the Sulston score Pr(Sulst)employed in the FPC package, LTC uses metrics Pr(Siid)and Pr(Sind)that estimate p-values more accurately (refer to Additional file 1, Section 1 for description of calculation) and the corresponding modifications Pr(SiidM)and Pr(SindM)that take into account the number of shared genetic markers. Metric Pr(Siid)is based on the Sulston model of "random" clones [17]. This model assumes that appearances of fragments (bands) of different lengths are independent and identically distributed (iid assumption). Unlike the Sulston score, Pr(Siid)provides a good approximation even in situations with relatively large clones and a high number of matching × bands. Metric Pr(Sind)is based on a model with similar assumptions, but bands may not be identically distributed (in contrast to iid assumption) (for details see Additional file 1, Section 1).

Let n1 and n2 be the numbers of bands in clones c1 and c2. LTC approximates the probability that the number of bands present in both clones c1 and c2 is equal or higher than k by exp{-(a0 + a1k + a2k2)}, for k > k0. Coefficients a0, a1, and a2 are estimated by the Monte-Carlo method, i.e., by simulating a large number of pairs of random clones (with n1 and n2 bands) and scoring the "observed" number of shared bands (refer to Additional file 1, Section 1 for details); value k0 is defined by inequality: exp{-(a0 + a1k + a2k2)} < 0.001 for all k > k0 (see also Additional file 1, Section 2). Shared genetic markers are taken into account by adding the term amarkm to a0 + a1k + a2k2, where m is the number of shared markers. We used amark = 100 ln 10 to make the overlap of clones with a shared marker more significant than 10-100. If bands are not assumed to be identically abundant, then the number of shared bands k is represented by -∑ lnf b , where the sum is taken over all shared bands and f b is the abundance of band b. Clearly, if all bands have equal abundances then k is proportional to the number of shared bands (i.e., both formulations will give the same result). Values f b can be estimated by the maximum likelihood method: the probability that a random clone with n bands contains band b is 1-(1-f b )n. [Logarithm of likelihood is then equal to ∑ c with b log(1-(1-f b )n(c)) + ∑ c without b log((1-f b )n(c)); for small f b (f b <<1/nmax, nmax= max c n(c)), the maximum likelihood estimation for f b is close to nmean-1π b /(1-π b ), where π b is a proportion of clones in the entire database having band b and nmean as the mean number of bands in the clone]. Taking band abundances into account satisfies the additive condition, when the sum of weights of two bands is equal to the weight of "pair of bands" that appears with an abundance equal to the product of abundances of the bands (because bands are supposed to be independent).

(ii) Temporal exclusion of clones and clone overlaps not proved by parallel paths (TENPP-procedure)

In contrast to FPC, LTC excludes putatively false significant overlaps and putatively problematic clones before clustering. The main idea behind the identification of problematic clones and clone overlaps is that each part of the chromosome is most probably covered by several clones (although in fact, some parts can be uncovered or poorly covered by clones of the BAC library). One can expect that chimerical clones and false clone overlaps are not proven by parallel clones [39, 52]. Thus, clustering should subdivide the clones into groups covering different parts of the chromosome. LTC clusters the clones in such a way that each position of a chromosomal region (without ends), covered by clones from the cluster, is covered by several (at least by three) significantly overlapping clones. Moreover, it requires that, even after excluding any single clone or clone overlap, for any pair of clones c i and c j from the cluster C0, a sequence of clones c(1),..,c( n ) from C0 should exist such that c(1) = c i , c( n ) = c j and overlap of clones c( k ) and c( k +1) is significant for all k = 1,..,n-1.

Let Pr0 be a liberal level of cutoff (we used 10-12≈10-3/N2, with N≈60,000 clones in our analyses). LTC calculates all pair-wise clone overlaps Pr and considers the results as a net of significant (relative to the selected cutoff Pr0) overlaps. For each edge e (significant overlap) of the net, parallel short paths (of 2-5 edges) connecting vertexes (clones) of edge e, but not going through edge e, are searched (see Additional file 2, Fig. AF2.1a). Edges having no such paths are temporarily excluded from the analysis. Then, the following procedure is applied to each vertex c0: for each pair of neighbor vertices c1 and c2 (connected with c0 by edge, i.e., Pr(c0,c1) < Pr0 and Pr(c0,c2) < Pr0), short "parallel" paths of 2-5 edges in this net, connecting vertices c1 and c2, but not going through vertex c0, are searched (Additional file 2, Fig. AF2.1b). Vertices having at least one pair of neighbors without such parallel paths are also temporarily excluded from the analysis. The excluded clones and overlaps can be used later to merge or elongate contigs (see below). For convenience, we refer to such a procedure of T emporal E xcluding from the analysis N ot proven by P arallel P aths as to the TENPP procedure with respective cutoff.

(iii) Clustering with "adaptively" varying cutoff

The ideal way to obtain a reliable contig assembly is to write down a full likelihood for the entire problem that allows a possibility for errors and repeats in data and find out all solutions with likelihood close to the maximum [37]. However, the expected length of band map for entire dataset is usually too long and finding maximal likelihood solution will be too hard even for modern computers. Moreover, estimating the probabilities included in the likelihood function is also not straightforward. Subdividing the data into clusters putatively covering non-overlapped parts of genome sequences by identification of putatively false clone overlaps and chimerical clones can simplify the problem [37]. In contrast to FPC that uses a stringent uniform cutoff LTC starts clustering with a liberal cutoff Pr0 by the single-linkage algorithm [53], and selects the resulting r easonably s ized clusters (rs-clusters). By "reasonably sized clusters", we arbitrarily (for certainty) define clusters with 6-500 clones. Actually, larger clusters also can be analyzed by the LTC program and subdivided into parts with linear topological structure (see below). Such a restriction was used to reduce the CPU time increasing cubically with number of clones in cluster. Additional reason for this restriction was the fact that the program Pajek [54] used for the visual control of the net structure of significant overlaps cannot represent large nets well. Using other programs for net visualization can help to overcome this problem. Artificial split of real contig can result in additional errors in ordering of obtained parts complicating their merging on the next steps.

For each cutoff level, LTC excludes putative false clone overlaps and putative chimerical clones in large clusters (with more than 500 clones) by the TENPP procedure described above. Then, the single linkage algorithm is run again. At the next step, LTC increases the stringency, but only after removing the selected reasonably sized clusters (i.e., protecting them from further "dissolving"). A schematic representation of the algorithm is provided in Figure 1. Refer to Additional file 1, Section 3 for simple example illustrating the efficiency of the proposed algorithm. Note that repeating the TENPP procedure will be effective only in cases of high proportions of false significant clone overlaps. If such proportion is not high, then it is reasonable to expect that most of the putatively chimerical clones and false clone overlaps will be excluded directly after the first TENPP procedure. Sizes of clusters can also be reduced by temporary exclusion of "buried" clones. These clones can be used later to obtain more reliable band maps.

Figure 1
figure 1

Scheme of clone clustering with adapting cutoff. Diamonds denote single-linkage clustering with corresponding cutoff. Circles denote procedure of excluding clones and clone overlaps not proved by parallel paths in the net of significant (relative to the corresponding cutoff) clone overlaps from the analysis.

(iv) Looking for linear topological structure

False clone overlaps and chimerical clones can lead to clusters with non-linear topological structures (Figure 2a) which is incompatible with the one-dimensional structure of eukaryotic chromosomes. To facilitate the detection and visualization of such clusters, LTC employs a representation of clusters as nets of significant overlaps. Obviously, ordering topologically non-linear clusters is problematic. To overcome this problem, we propose to split such clusters into sub-clusters with linear topological structure by excluding clones from the branching nodes from the analysis (Figure 2b). Non-linearity of the cluster structure can be detected by scoring ranks of vertices relative to the diametric path. By definition, a diametric path of the net is the longest, in terms of the number of edges, among the shortest paths between all possible pairs of vertices (e.g., [55]). Note that several diametric paths can exist for the same net. Any of diametric paths can be used for detection of non-linearity: The presence of vertices with rank 2 and higher indicates on non-short offshoots from the selected diametric path and, hence, a possibility of non-linear structure (Figure 2a).

Figure 2
figure 2

Problematic clusters with non-linear topological structure. (a) Identification of non-linearity of cluster topological structure: presence of vertices with not too small ranks (e.g., more than 2) points to non-short offshoots from the selected diametric path (marked in bold line). (b) Splitting clusters into sub-clusters having linear topological structure. See also Additional File 1, Section 4.

(v) Ordering using global optimization

LTC orders clusters with linear topological structures without constructing their band maps. Here, the ordering problems are formulated in terms of global optimization of some criterion (similar to Fickett and Cinkosky, 1992 [56]; Alizadah et al., 1993 [57]; Wang et al., 1994 [39]; Flibotte et al., 2004 [46]). For simplicity, let's consider only the situation where all clones are not buried. This can be achieved by the temporal exclusion of buried clones, although sometimes this can lead to loosing the contig connection. The criterion W(Ω) of clone order Ω = (cΩ(1),..., cΩ( n )) is calculated as:

W ( Ω ) = Σ k = 1 , .. , n − 1 W Ω ( k ) , Ω ( k + 1 ) − b ( Ω ) W 0 ,

where Wi,j= -log Pr(c i , c j ), b(Ω) is the number of adjacent (within ordering Ω) clones with Pr(cΩ( k ), cΩ( k +1)) > Pr0 and W0 is the penalty for non-significant overlap of adjacent clones.

The maximization of such criterion can be reformulated as the well-known and intensively studied Traveling Salesman Problem (TSP) without the need to return to the starting point. Let Wmax be the maximum of Wi,j= -log Pr(c i , c j ). LTC defines distance between two clones as

d ( c i , c j ) = W max − W i , j + W 0 1 { P r ( c i , c j ) > P r 0 } ,

where 1{Pr(c i , c j ) > Pr0} is indicator function equal to 1, if Pr(c i , c j ) > Pr0, and equal to zero otherwise. Global optimization is effective especially when marker information is also available [46]. The exact solution to the TSP is a computationally challenging problem. Nevertheless, good heuristics (e.g., based on evolution strategy optimization) for the solution of TSP were developed for situations where the number of vertices is up to 103 orders of magnitude [58]. Coordinates of clone ends within the contig based on the solution path Ω best can be calculated by methods proposed by Flibotte et al. (2004) [46]. Using such global optimization approach can result in a reduced number of Q-clones and of places in the contig where two neighbor clones have an unexpected non-significant overlap. More effective triplewise linking distance also can be used to reduce the effect of false positive and false negative bands [39].

(vi) Resampling verification of ordering

The quality of the clone order within a contig is characterized not only by the value of the chosen criterion, but also by its robustness to small uncertainty of band content of the clones, which can be referred to as "contig stability". To evaluate this stability, LTC uses jackknife re-sampling iterations (in contrast to bootstrap ones used in Cuticchia et al., 1993 [59] and Wang et al., 1994 [39], that can artificially increase the significance of clone overlaps). Namely, LTC first constructs the order using clone overlaps scored over all bands. In addition, it constructs orders using clone overlaps based on randomly selected subsets of bands (say, 95% of the total set of bands). Then, the identification of unstable regions is conducted based on the frequency distribution of the right-side and left-side neighbors for each clone in the contig order. The higher the deviation from 1 (i.e., from the "diagonal" pattern) is, the less certain the local order is ([58, 59]). One of the main reasons for the appearance of unstable orders is the high similarity of parallel clones that cannot be ordered properly due to missing and false bands. Excluding parallel clones allows the construction of a stable "skeleton" map, similarly to the approach suggested for building genetic maps (see [58]).

(vii) Merging of sub-contigs

After ordering, LTC tries to elongate contigs by merging those that display significant end-to-end overlaps (which may be also achievable by adding 1-2 supplementary intermediate clones), or by adding singletons. First, LTC re-analyzes all clones and clone overlaps temporally excluded at previous stages (see above). To elongate a concrete contig, LTC searches for all clones connected (by significant overlaps or via short paths of significant overlaps) with the clones from either of the contig ends. If adding all of the clones (for one of the two contig ends) does not lead to a violation of contig linearity, then such elongation is not problematic. If adding the clones does lead to branching (i.e., contradicts the linear structure of the chromosome), then each of the possibilities of linear elongations needs to be considered (e.g., see Additional File 2, Fig. AF2.2). The correct elongation can be detected by testing clone overlaps based on clone-end sequencing [23]. The same problem arises when clones from one contig significantly overlap with internal clones from another contig. The availability of DNA markers (in clones) with known chromosomal position helps to prevent the merging of contigs from different chromosomal zones. Contigs resulting from elongation should be reordered (see stage (v)).

(viii) MTP establishment

In the current version of LTC, MTP construction was based only on the topology of clone overlaps without building of the band map (in contrast to FPC). In our algorithms, we assume: (i) that MTP should include the terminal clones of the contig; (ii) the adjacent clones in MTP should significantly overlap (at the chosen threshold); and (iii) MTP should be of minimal length, i.e., include the minimal number of clones. To satisfy these conditions, we select clones for MTP from the diametric path of the net of the significant clone overlaps (see paragraph iv). In practice other criteria come into play, such as whether or not a particular BAC has a genetic marker on it and whether or not BAC end sequence is available. A further sophistication of the criterion may also the lengths of the clones (longer clones should be preferable) and overlaps with clones disturbing the contig linearity (presence of markers based on BAC-end-sequencing of such clones can clarify clone-overlaps).

II. Verification and improvement of FPC contigs

LTC can be used also to verify and improve contigs obtained with other methods and tools. To test the quality of a contig and understand the underlying reason(s) for the assembly problem, the following procedure can be applied: (i) calculate all pair-wise probabilities of clone overlaps and represent the contig as a net of significant overlaps; (ii) test connections within this net; and (iii) test for topological linearity of the contig. If poor overlap (at the sequence level) of adjacent clones in the MTP occurs together with low significance of the clones' overlap, but the contig is connected and has a linear structure, then it makes sense to attempt fixing the contig by reordering the clones and selecting an alternative MTP. In this case, one needs to take into account that if some part of the analysis was already done for the old MTP, then it can be cheaper to select a MTP with common parts with the old one, rather than to repeat the analysis for the optimal MTP. If the detected problems in the MTP can be explained by non-connectivity or non-linearity of the contig structure, then the contigs must be split into connected parts each with linear structures. LTC first temporally excludes clones and clone overlaps not proven by parallel paths from the analysis (see above). After contig splitting, reordering and verification, LTC attempts merging the contigs by decreasing cutoff stringency or via returning back the previously excluded clones and clone overlaps. LTC also checks (whenever possible) that the added clones are not specified as belonging to other parts of the chromosome.

Comparison of clone partitions

Different clone partitions obtained with different clustering schemes were compared for the number of clones covered by reasonable size clusters and by direct comparisons of clusters. The comparison was performed by scoring the Rand index R[60] and its modification R' (see Additional File 1, Secttion 5). Additionally to these characteristics, the mean number M of clusters of the one clone partition overlaps with each of the other clone partition clusters was scored; if the two partitions are very similar, M should be close to 1. Two modifications of the M value, M* and M**, that reduce the role played by small clusters and singletons were also scored (see Additional File 1, Section 5).

Graphical representation of the cluster topological structure

Topological structures of clusters are represented by the net of significant (relative to some specified cutoff) clone overlaps. The nets were drawn using the publicly available program Pajek [54]. Draft pictures were obtained with the Kamada-Kawai algorithm for drawing undirected graphs [61] performed in the Pajek program.

Simulation of BAC libraries based on the known genome sequences

To test the effectiveness of contig assembly algorithms we simulated clone libraries and clone fingerprints based on known genome sequences (analogously to Xu et al., 2004 [62] and Krzywinski et al., 2007 [63], instead of using artificial genome sequences employed by other authors, e.g. Cuticchia et al., 1992 [50]; Soderlund et al., 1997 [35]). By Renzyme we denote the sequentially ordered set of N restriction sites r i corresponding to the chosen enzyme sequence (we used HindIII). For convenience, we supplement this set by r0 and r N +1 corresponding to the start and end points of the sequence. Let L be the total length of the sequence. For each generating clone the program selects a start and end points from Renzyme. Index istart is selected randomly from 0 to N + 1. Index iend calculated by istart + s h, where s is equal to 1 or -1 with probabilities 0.5 (defines clone direction) and h is the integer part of a normally distributed random value with mean a and variance σ2 (defines clone sequence length distribution). Values a and σ2 are selected such that mean and standard deviation of clone length be about 120 kbp and 30 kbp, respectively (a = 120λ, σ = 30 λ, where λ = (N + 1)/L). Chimerical clones were generated as union of regular clones, occurring with a probability pchimer (= 0.05 in our simulations).

Fingerprint for generated clones were defined by the set of bands of K types (corresponding to K used protruding-end restriction enzymes) presented in the clone; we used enzymes BamHI, EcoRI, XbaI and XhoI (K = 4) analogously to Ding et al. (1999) [21]. Presence of a band b with length L b and type k b is defined by the presence of sequence part of length L b bp bounded by the restriction site of enzyme k b from the one side and restriction site of one of the K enzymes mentioned above, or restriction site of blunt-end restriction enzyme (we used HaeIII) or clone end from the another side. In fact, sequence parts bounded by two different protruding-end restriction enzymes are twice included into fingerprint. Bands with length <Lmin = 50 bp or >Lmin = 500 were filtered out. Clones with number of bands n < nmin = 50 or n > nmax = 250 were also filtered out.

In the preparation of input data for the FPC and LTC programs we introduced noise into sequence length scoring to simulate errors in wet fingerprinting process. "Observed" band length was calculated by L b + err, where err is a random value uniformly distributed within interval (-ε, ε) or within (-3ε, 3ε) with probabilities p ε = 0.90 and p3 ε = 1-p ε = 0.1 respectively. This imples that difference between two observations of the same band is less than constant tolerance value t with probability τ = p ε 2(1-(1-t/2ε)2) + p3 ε 2(1-(1-t/6ε)2) + 2p ε p3 ε (t/3ε). We used ε = t/ 2 = 0.2 implying τ = 0.94. To satisfy FPC format, band lengths were multiplied by 30 and rounded (hence tolerance t = 12). Type of band was taken into account by adding 0, 5000, 10000 or 15000 for the resulted band length for bands of type k = 1, 2, 3 and 4, respectively. Some bands of clones were excluded (with probability pmissing = 0.05) to simulate false negatives caused by problems in PCR reactions. Automatic FPC assembly was conducted using cutoffs 10-75 - 10-45 with step × 105 and DQer for contigs having more than 10% Q-clones (with step × 109) and ReBuild with corresponding cutoff if needed (similar to Paux et al., 2008 [48]).

Materials

To illustrate the advantages of the LTC analytical framework we employed the results obtained in the physical mapping of wheat chromosome 3B [48]. The corresponding database included High Information Content Fingerprinting results of 56,952 of BAC clones obtained from the chromosome 3B specific library. The initial input data for LTC included band lengths classified according to four dyes. For each dye, up to 4500 distinct points (band sizes) were obtained. Two bands were considered of the same if their sizes were within a tolerance = t (e.g., constant t = 4).

The elements of the LTC analytical framework were tested with fingerprinting data from two regions of chromosome 3B corresponding to clones located in the 3BL7-0.63-1.00 and 3BS1-0.33-0.55 deletion bins. Sixty-nine contigs comprising 3,606 clones that were found to belong to bin 3BL7-0.63-1.00 and 180 (8,167 clones) to bin 3BS1-0.33-0 [48] were used in the analysis. Using LTC, these 249 contigs were verified/corrected/extended; alternative contigs were also constructed and compared to the results obtained with FPC.

Simulated BAC libraries

Simulated BAC libraries LibRiceChr1 and LibMaizeChr1 were constructed based on the known genome sequences SeqRiceChr1 and SeqMaizeChr1 for chromosome 1 of rice (45 Mbp, available at http://rgp.dna.affrc.go.jp/whoga/download.html.en, file chr01.fa.gz) and chromosome 1 of maize (300 Mbp, available at http://ftp.maizesequence.org/current/assembly, we used version that was downloaded at June 10, 2010), respectively. Basic characteristics of the libraries are summarized in table Table 1.

Table 1 Basic characteristics of simulated BAC-libraries

Results and Discussion

Distribution of bands and clone lengths

In total, 30 to 275 bands were scored (out of 18,000 possible) for each of the N = 56,952 clones of the 3B BAC library (for more details of this data see [48]). The frequency of a band b was calculated as π b = N b /N, where N b is the number of clones containing bands with lengths different from b not more than the tolerance value. The bands showed quite variable abundances (Additional File 2, Fig. AF2.3a). In particular, one band was observed with π b = 0.93, one with π b = 0.78, six with π b from 0.2 up to 0.32, 220 with π b from 0.1 up to 0.2, 654 with π b from 0.05 up to 0.1, and others were observed in less than 5% of the clones. The minimal band abundance was 1.2%. We found more than 30 of the five-band combinations ("band haplotypes") present in multiple clones in a range of 1000-2271 (3.6% of total) clones, where π b of bands was in the range of 0.120 to 0.156, pointing to highly significant "linkage disequilibria" of the bands. This can result from physical overlap of considered clones and from repeats [6, 64]. The distribution of clone length (scored as number of bands) was bimodal (Additional File 2, Fig. AF2.3b). The origin of this bimodality is in the construction of the BAC library itself. Three sub-libraries were constructed including large, medium, and small fragments (see [65]). Fingerprinting results revealed that the medium length sub-library was actually a mix of large and small inserts, and not medium inserts. This resulted in two distinct populations of clone length.

LTC, a new algorithm for building contigs from fingerprinted clones

In this work we have developed a new analytical framework, LTC, for contig assembly of fingerprinted BAC clones that can be used as an alternative or as a complement to FPC. In contrast to FPC, the LTC program starts clustering with a relatively relaxed cutoff and uses the topology of significant clone overlaps to obtain longer and more realistic contig structures. Instead of a uniform cutoff, LTC uses a procedure that adaptively increases the cutoff stringency. Using stringent cutoffs only for large clusters generates fewer short contigs or singletons. After ordering, LTC merges the contigs by relaxing the cutoff (analogously to FPC), hence, this approach can be referred to as "up-down-up", in contrast to the FPC approach that can be referred to as "down-up" (stringent cutoffs applied from the beginning). LTC differs from FPC in a number of important features that are: (a) the metrics of clone overlap, (b) consideration of the band abundances, (c) the algorithm used for clone clustering, (d) the use of the topological structure of clusters for contig construction, (e) the application of global optimization methods for clone ordering, and (f) the assessment of the reliability of the mapping results by re-sampling. These parameters should enable the construction of more robust and longer contigs in particular when dealing with large and repetitive genomes.

To demonstrate the efficiency of LTC in contig assembly, we compared contigs obtained by LTC with contigs obtained using the standard FPC package for BAC clones from the wheat chromosome 3B [48]. The results are presented for the three main analytical stages shared between the LTC and FPC packages, i.e., (i) the identification of significant clone overlaps; (ii) the subdivision of clones into clusters of significantly overlapped clones; and (iii) the ordering of clones within clusters. To assess the capacity of LTC to verify, repair, and elongate contigs obtained by FPC, we also ran LTC on the 3B dataset previously obtained by FPC for the same contigs.

(i) Identification of significant clone overlaps by LTC

Both FPC and LTC identify significant clone overlaps in a similar way: for each pair of clones, c1 and c2, the p-value Pr(c1, c2) of clone overlap is calculated and compared with a predefined threshold Pr0 (cutoff). Clones c1 and c2 with Pr(c1, c2) < Pr0 are declared as overlapping clones. In contrast to FPC that uses Sulston score Pr(Sulst)[17], LTC uses new alternative metrics Pr(Siid)and Pr(Sind)(see paragraph (i) in Methods) for more accurate approximation of p-value. To compare metrics used in FPC and LTC, values Pr(Siid), Pr(Sulst), and Pr(Snid)were calculated for each pair of clones from the 3B dataset (56,952 clones). The results indicate that Pr(Siid)generally overestimates the significance of clone overlap compared to Pr(Snid). For highly significant clone overlaps, Pr(Siid)also provides a substantially higher p-value than the standard Sulston score Pr(Sulst)(see Figure 3). However, logarithms of Pr(Siid), Pr(Sulst), and Pr(Snid)were not absolutely correlated, especially for strongly overlapping clones. For example, liberal cutoffs P0(Snid)= 10-12, P0(Sulst)= 10-17.2, and P0(Siid)= 10-19.4 gave comparable numbers of significant clone overlaps (~ 12.7 per clone, see Figure 3a). We found that 2.1% of clone-overlaps, which was significant with Pr(Snid), was not significant with Pr(Sulst)(i.e., were identified as false positives) and vice versa (see Figure 3b). Similarly, 2.5% and 2.7% of false positives were found for pairs based on the Pr(Siid)& Pr(Snid), and Pr(Siid)& Pr(Sulst)criteria for clone overlaps, respectively. Such inconsistencies of the criteria can be explained by non-identical band abundances and by insufficient accuracy of the Sulston score to estimate clone overlap probability (see [40]).

Figure 3
figure 3

Comparison of metrics: (a) Number of significant clone overlaps K per clone for different cutoffs. Using liberal cutoffs P0(Snid)= 10-12, P0(Sulst)= 10-17.2, and P0(Siid)= 10-19.4 give about the same number of significant clone overlaps (~ 12.7 per clone). (b) Pr(Sulst)vs. Pr(Snid). Dots in rectangle I correspond to clone overlaps significantly relative to Pr(Sulst)= 10-17.2 but not significantly relative to Pr(Snid)= 10-12 (false positives). Dots in rectangle I correspond to clone overlaps significantly relative to Pr(Sulst)= 10-17.2 but not significantly relative to Pr(Snid)= 10-12 (false positives). Dots in rectangle IV correspond to clone overlaps significantly relative to Pr(Snid)= 10-12 but not significant relatively to Pr(Sulst)= 10-17.2 (false negatives).

A more accurate estimation of p-values for clone overlaps used in LTC reduces the proportion of falsely significant clone overlaps and increases the proportion of real clone overlaps that have a significant number of common bands. This leads to a reduction of clustering errors. Such errors can result in two undesired outcomes: (a) a wrong contig ordering with poor overlap of some adjacent clones from MTP, referred to as a gap that will call for (non-natural) splitting the contig into shorter contigs; and (b) a wrong partition of contigs into "independent" contigs that could not be merged. Both outcomes will yield shorter contigs than one could obtain using more correct clone clustering. Hence, using more accurate metrics for clone overlaps can result in longer contigs even with FPC algorithms for contig assembly.

Analogously to FPC, choosing cutoff stringency for clone clustering should be based on a tradeoff between the advantages of stringent and liberal cutoff values. Indeed, on the one hand, clustering with a liberal cutoff results in large but unreliable clusters, where ordering and detecting problematic clones and clone overlaps is challenging. On the other hand, clustering with a stringent cutoff results in many small clusters from which the ordering and identification of problematical clones is presumably "easy", while merging is difficult. Therefore, using liberal cutoffs may be a reasonable strategy if and only if powerful tools for multipoint ordering are available, such as with LTC.

(ii) Clustering the 3B fingerprinting dataset: comparing different procedures

Partition of clones using the LTC adaptive clustering procedure

We applied the LTC main adaptive clustering procedure to the 56,952 fingerprints obtained for chromosome 3B [48]. Out of N(N-1)/2 = 1.6·109 possible clone overlaps, 361,571 (i.e., 0.02%) were significant at Pr(SnidM)< Pr0 = 10-12 (Figure 4a). From these, 2,155 clone overlaps (0.6% of the significant ones) and 4,097 clones (7.2% of total) have not been confirmed by parallel paths (Figure 4a). After temporal exclusion of these clones and clone overlaps, the entire database was subdivided into 828 clusters with a minimal number of clones per cluster of 6. In total, 6,550 clones (11.5% of total) were found in clusters with smaller sizes or as singletons. Among the 828 clusters, one was large (9,386 clones), three were of intermediate size (912, 617, and 525 clones) whereas the remaining ones included less than 336 clones per cluster (Figure 4a).

Figure 4
figure 4

Adaptive clustering using different scenarios: (a) full scheme (see Fig. 1); (b) reduced scheme (TENPP procedure used only once); (c) clustering with adaptively changing cutoffs (for simplicity, only cutoffs 10-12, 10-25, 10-35, 10-45, and 10-50 were used in the example); (d) clustering with uniform stringent cutoff 10-50 (non-adaptive clustering analogous to FPC). Clustering (a) and (b) were based on P(SnidM); clustering (c) and (d) were based on P(Snid).K - number of significant clone overlaps; n - number of clones in clusters of reasonable size (from 6 to 500); ∑n - total number of clones in clusters of reasonable size; N ¯ - number of clones in clusters with less than 6 clones, singletons, or clones excluded by TENPP procedures; K ¯ - number of significant clone overlaps excluded by TENPP procedures; k - number of clusters of reasonable size; M = ∑n/k - mean number of clones in clusters of reasonable size; N = 56,952 is the total number of clones in the database; N i -number of clones in the large or intermediate size cluster i.

The TENPP procedure was used again on the entire database with cutoffs Pr1 = 10-3Pr0 = 10-15 which found additional significant clone overlaps and 243 clones not proven by short parallel paths. After temporal exclusion of these clones and clone overlaps, four large and intermediate clusters were subdivided into 108 clusters with at least 6 clones per cluster. In this stage, additional 339 clones were found either in smaller size clusters, or as singletons, or were excluded by the last TENPP procedure. Among the remaining 108 clusters, one was large (1,494 clones), one was intermediate (622 clones), and the others contained up to 402 clones per cluster (Figure 4a).

Another round of the TENPP procedure with cutoffs Pr2 = 10-3Pr1 = 10-18 identified additional significant clone overlaps and 225 clones not proven by short parallel paths. Large and intermediate clusters were subdivided into 13 clusters with a minimal number of clones per cluster of at least 6 after temporal exclusion. This resulted in 150 clones that were found in clusters with smaller sizes or as singletons, or that were excluded by the last application of TENPP procedure. Thirteen clusters were obtained with one intermediate containing 794 clones while the rest had up to 395 clones per cluster (Figure 4a).

A final TENPP procedure with cutoffs Pr3 = 10-3Pr2 = 10-21 found additional significant clone overlaps and 219 clones not proven by short parallel paths. After temporal exclusion, the remaining clusters were subdivided into 9 clusters with 6 to 216 clones; 78 clones were found in clusters with smaller sizes and singletons or were excluded by the last TENPP procedure (Figure 4a).

Thus, in total, 56,952 fingerprints were subdivided into 951 clusters comprised of 6 to 402 clones; 11,214 clones were found in clusters with less than 6 clones, were singletons, or were excluded by the last TENPP procedures. The average number of clones per cluster was 48.3. The partition of clones obtained with this first approach will be referred hereafter as CLTC.

Skipping some components of the LTC clustering algorithm

To demonstrate the power of the LTC adaptive clustering procedure, we skipped some of the clustering scheme components described above (Figure 4a). Three versions were considered: (i) using TENPP only once; (ii) clustering with adaptively changing cutoff stringency but without the TENPP procedure; and (iii) clustering with a uniformly stringent cutoff (analogous to the FPC algorithm).

  1. (i)

    Figure 4b shows that if the TENPP is used only once (with cutoff Pr 0 = 10-12) then in addition to reasonable size clusters, clustering with cutoff Pr 1 = 10-15 results also in two large (with 3,428 and 1,399 clones) and two intermediate size (with 880 and 615 clones) clusters. Repeated clustering with the more stringent cutoff Pr 2 = 10-18 gives additional reasonable size clusters plus one large cluster with 2,878 clones and two intermediate size clusters (with 640 and 615 clones) (Figure 4b). Further increase of the cutoff stringency up to Pr 3 = 10-21 splits these clusters into one large cluster with 1,900 clones and two intermediate size clusters (with 625 and 613 clones, Figure 4b). The clone partition obtained with this scenario was called C 1.

  2. (ii)

    Without using the TENPP procedure (Figure 4c), one very large cluster is obtained even for very stringent cutoffs (e.g., 51,780 for Pr 0 = 10-12, 48,577 for Pr 1 = 10-25, 41,527 for Pr 2 = 10-35, and 19,590 for Pr 3 = 10-45). This clone partition was called C 2.

  3. (iii)

    Clustering with a stringent uniform cutoff (10-50) (Figure 4d) resulted in more singletons and smaller clusters compared to clustering with adaptively changing cutoff (Figure 4c). The clone partition obtained with the uniform cutoff will be referred to as C 3, whereas the partition obtained with FPC (after ordering and merging) will be called C FPC. In partition C FPC, 41,295 of the 56,952 BAC clones were automatically subdivided into 1,995 contigs (with 6 and more clones) that were manually merged into 1,036 contigs with 6 to 290 clones per contig; the remaining 15,675 clones were found in clusters with less than 6 clones or were singletons [48].

Among the 951 clusters of CLTC (Figure 4a) only 47 have non-linear topological structures, and 31 of them have only one branching point. The branching 47 non-linear clusters were split into about 171 "topologically linear" clusters. Note that among the clusters obtained with the other scenarios (Figure 4c-d) the proportion of clusters with non-linear topological structure was much higher (up to 70% in C3, not shown). Thus, these results demonstrate that the TENPP procedure enables the use of much more liberal cutoffs to obtain clusters of reasonable size. Repeating the procedure with more stringent cutoffs helps to split large clusters into reasonable size clusters. Clustering with adaptively changing cutoffs helps to protect from "dissolving" the reasonably sized clusters obtained at liberal cutoff stringency under more stringent cutoffs (see Additional File 2, Fig. AF2.4).

Comparing LTC vs. FPC clustering

The CLTC clustering contains more clones in the reasonable size clusters than CFPC: 80% vs. 73% of the 56,952 clones. As described in the Methods section (paragraph iii), "reasonable size of clusters" were arbitrary (for certainty) defined as clusters with 6 to 500 clones; obviously, other ranges can also be successfully used. Although the number of clusters in CFPC was less than in CLTC after the split into clusters with linear topological structure (1,036 vs. about 1,200), CLTC clusters are proven by parallel clones and hence should be more reliable. Moreover, many of the FPC clusters were obtained by manual merging smaller clusters, whereas 1200 clusters were obtained by LTC before any merging was undertaken. The Rand index was high, R(CLTC, CFPC) = 0.997 because most of the clone pairs appeared in different clusters. The modification of R: R'(CLTC, CFPC) = 0.56 was more informative, reflecting that large clusters in CLTC and CFPC were rather different (one of the FPC clusters overlapped with eight LTC clusters). Nevertheless, the values of M1(CLTC, CFPC) = 3.7, M1*(CLTC, CFPC) = 2.15, and M1**(CLTC, CFPC) = 1.25 for the mean number of clusters from CLTC overlapped with CFPC (see Comparison of clone partitions in the Methods section) indicating that most of the differences between CFPC and CLTC originate from the difference between Q-clones excluded by FPC and those excluded by the TENPP procedure in LTC. If one excludes the clones and clone overlaps from the CFPC clusters that were excluded by TENPP in LTC and a clustering is applied with a liberal cutoff Pr0 = 10-12, then a new clustering CFPC(TENPP) is obtained with M1**(CLTC, CFPC(TENPP)) = 1.3 and M1**(CFPC(TENPP), CLTC) = 1.0. This demonstrates that each of CFPC(TENPP) clusters is actually a sub-cluster of CLTC. Thus, all of LTC contigs are longer than the corresponding FPC sub-contigs proven by parallel clones.

(iii) Verification of FPC contigs

The LTC adaptive clustering procedure was also used to test 249 FPC contigs (each containing 6 to 290 clones; the corresponding clone partition was mentioned in the previous section and denoted by CFPC) that were assigned to the 3BL7-0.63-1.00 and 3BS1-0.33-0.55 deletion bins of wheat chromosome 3B [48]. The LTC analysis indicated that 201 of the contigs (80%) had regions where the overlap between two adjacent clones was not significant with the Pr0(Siid)= 10-12 criterion and therefore were considered as gaps by LTC. The maximum number of gaps per contig was 44 (found in a contig of 283 clones). The gaps detected in the FPC-based contigs were of five main origins: (i) Contig consisting of two (or more) non-connected parts; (ii) Buried clones causing conflicts between requests of linear ordering and contig connectivity; (iii) Weak power of ordering tools in FPC based on local optimization; (iv) Non-synchronous utilization of information on shared markers and common bands; (v) Topologically non-linear structure of the contig, mostly caused by the presence of clones and clone overlaps not proven by short parallel paths.

Contig consisting of two non-connected parts

Only one contig located in the 3BL7-0.63-1.00 bin corresponded to this category. The net of significant clone overlaps for this contig consisted of two non-connected parts (even when p-value of clone overlaps was calculated using standard and the very liberal cutoff Pr0(Sulst)= 10-10). These parts can be connected by adding only one clone (see Additional File 2, Fig. AF2.5). Such a situation may result from the exclusion of one or two connecting clones (3B_043_I24 and 3B_073_B21) presumably considered as Q-clones by FPC at the stage of cluster ordering.

Buried clones causing conflicts between ordering and contig connectivity

The significance of clone overlap is determined by the number of common bands and clone lengths. Hence, if clone c i is buried in clone c j and significantly overlaps with clone c k ,, then it is not necessary that the overlap of clones c i and c k is also significant (see Figure 5a). Thus, after excluding buried clones, the net of significant clone overlaps can lose its connectivity (see Figure 5b). On the other hand, it is likely that each path along the net of significant connections, that visits all vertices, will visit the vertex corresponding to clone c i , buried in clone c j , before and after visiting the vertex corresponding to clone c j (see Figure 5c). In other words, vertex corresponding to clone c j is visited twice, and therefore, such paths cannot be found directly by tools constructed for solving the classical TSP. Such a situation was found in 12 contigs and could be solved by temporal exclusion of only a subset of buried clones from the cluster before ordering.

Figure 5
figure 5

Examples of situations where buried clones cause problems in clone ordering by LTC. (a) Clone c i is buried in clone c j and significantly overlaps with clone c k , but overlap of clones c i and c k is not significant. (b) The net of significant clone overlaps for the situation (a). Exclusion of buried clone c i leads to a loss of connection between the left and right parts of the net. (c) Clones cannot be ordered to fit the requirement that adjacent clones overlap significantly (because the buried clone significantly overlapped with one clone only).

Weak power of ordering tools based on local optimization

Gaps that can be repaired by clone reordering within clusters were found in 198 of the 249 considered contigs (data not shown). Such reordering of clones usually changes clone positions in the contig only locally (see Additional File 2, Fig. AF2.6). Band maps based on LTC clone ordering were compared with FPC band maps. It was found that adjacent clones in LTC ordering cover on average more common bands in a corresponding band map than adjacent clones in the FPC band map (data not shown). This leads to better correspondence between number of common bands and number of common positions in the band map for significantly overlapped clones. Similar results were obtained by Flibotte et al. (2004) [46]. Increased coverage of the band map leads to higher robustness of LTC contigs to technical errors in fingerprinting. It also results in a lower number of band repeats. Hence, band maps based on LTC ordering are somewhat simpler than band maps obtained by FPC.

Topologically non-linear structure of the contig

A non-linear contig structure was found in 26 of the considered contigs (see examples in Figure 6). This was mostly due to the presence of clones and clone overlaps not proven by short parallel paths. Excluding such clones and clone overlaps (by TENPP procedure) and clustering by the single linkage algorithm subdivided all of the 26 contigs with non-linear structure into linear parts, but also splits some other (linear) contigs into smaller parts. Hence, some of LTC contigs are shorter than corresponding FPC contigs, but are much more robust to the presence of chimerical clones and false significant clone overlaps. It should be noted that end-to-end merging of LTC contigs seems to be a less problematic procedure than merging of FPC contigs because LTC does not place most problematic clones to the contig ends.

Figure 6
figure 6

Net representation of problematic FPC contigs with non-linear topological structure. Presented four contigs from the 3BL7-0.63-1.00 and 3BS1-0.33-0.55 deletion bins display non-linear topological structures. Vertices represent clones; edges represent significant (relative to arbitrary selected cutoff Pr0(SnidM)= 10-12) clone overlaps. Branching regions are marked by solid black circle.

Non-synchronous utilization of information on common markers and common bands

In certain cases, the anchoring of markers to physical contigs leads to situations where a marker is found in two different contigs or even in two non-overlapping clones of the same contig. Such a situation can arise in cases when a single marker is situated on the overlap of two clones having only a non-significant number of common bands. It also can be a result of marker duplication or errors in marker amplification or contig assembly. Using FPC tools, such a situation cannot be recognized (because any pair of clones usually has several falsely common bands that can be erroneously considered as real overlap) and generally results in the artificial fusion of the two contigs. In contrast, LTC is able to detect such a situation because marker information is used synchronously with shared-band information in the analysis. In this case, the presence of duplicated markers and/or errors in marker amplification will usually lead to non-linearity in the topological structure of the cluster that will be detected by LTC (see Additional File 2, Fig. AF2.7).

(iv) Elongation of FPC contigs by using LTC

FPC contigs from the 3BL7-0.63-1.00 and 3BS1-0.33-0.55 deletion bins were compared with the contigs that were constructed from the entire fingerprinting dataset using LTC (without initial clustering/ordering by FPC). It appeared that short FPC contigs (with 6-30 clones) were usually parts of LTC contigs. In such situations, LTC allowed the elongation of FPC contigs in a natural way. Several examples of such elongation of FPC contigs are presented in Figure 7. Long FPC contigs (with 50 and more clones) were usually different from those obtained with LTC. We assume that the main reason for such a difference is that clones and clone overlaps excluded by the TENPP procedure do not fully coincide with Q-clones and false significant clone overlaps identified by FPC. This leads to differences in clone partitions into clusters (see above) and can result in different contigs. In such cases, PCR amplification of BAC-end-sequences should be performed to assess the robustness of the contigs and determine which assembly is correct. For example, in the 3B dataset, some overlaps of adjacent clones from the MTP were not confirmed by sequencing and no satisfactory alternative MTPs could be found using FPC for these contigs (see Additional File 2, Fig. AF2.8a and b). For such situations alternative LTC contig assemblies were proven useful (e.g., to construct alternative elongation of the verified MTP (Additional File 2, Fig. AF2.8b). Note, sometimes LTC contigs included parts of different FPC contigs that were considered as belonging to different chromosomal zones based on marker assignment information (Additional File 2, Fig. AF2.9). We explain such observations by the fact that the location of the contig in the chromosomes is determined by position markers that were found in the clones of the contig. Hence, if contig assembly was made incorrectly, then the position of clones on the chromosomes is also determined with errors. Verification of the contig using BAC-end-sequencing can clarify the contig assembly and help to determine the position of clones on the chromosome.

Figure 7
figure 7

Examples of de novo assembled contigs that proved to contain contigs constructed by FPC. Black vertices reflect clones from the corresponding FPC contigs. White vertices correspond to clones which were not specified to belong to the zones 3BS1 or 3BL7.

Comparison of LTC vs. FPC on simulated BAC libraries based on known genome sequences

Simulated BAC libraries LibRiceChr1 and LibMaizeChr1 were automatically assembled using FPC program and our LTC package (see sections Methods and Materials). Resulted contigs were tested using methods from LTC program and compared with real position of simulated clones in the genome sequence. Results of this analysis are briefly summarized in the Table 2.

Table 2 Contig assembly with FPC and LTC programs

The results presented in Table 2 indicate that LTC contigs are longer and more reliable than those obtained by FPC: (a) average contig length 1,731 (LTC) vs. 1,045 kbp (FPC) for rice, and 1,532 vs. 1,128 kbp for maize; (b) number of contigs with non-linear topological structure was 0 out of 25 (LTC) vs. 8 out of 42 (FPC) observed in rice analysis, and 0 out of 190 (LTC) vs. 23 out of 257 (FPC) in maize. Low efficiency of chimerical clones identification seems to be the main source of the errors in the FPC contig assembly: only 9% (rice) or 7% (maize) of chimerical clones were not included into FPC contigs. Excluding clones and clone overlaps not proven by parallel clones (see also [52]) used in LTC algorithm enabled us to identify and exclude 89% (rice) and 81% (maize) of chimerical clones even using much more liberal cutoff: 10-25- 10-30 (LTC) vs. 10-45-10-75 (FPC). Diagnosis of contig linearity implemented in LTC can also assist in detecting problematic contigs assembled using FPC: 8 out of 22 chimerical rice FPC contigs and 23 out of 108 chimerical maize contigs were detected as having non-linear structure. LTC was also able to identify problems in clone positioning within FPC contigs leading to non-significant overlaps of adjacent clones. Reordering of clones within contigs without chimerical clones was usually able to fix these problems.

Using these simulated BAC-libraries we found that LTC can encounter some difficulties with ordering contigs that contain long repeats causing false significant overlap of clones that belong local neighborhoods. Likewise, we found one pair regions separated by ~27 Mbp in the sequence of maize chromosome 1 (positioned at 52 Mbp and 79.5 Mbp from the sequence start) with several significant (up to 10-32) overlaps between corresponding clones. In fact, these overlaps caused non-linear cluster topology (after TENPP procedure), but overlaps of clones within each of these two regions were by far more significant (10-70-10-140) and increasing the cutoff stringency resulted in split of this cluster into two with linear topological structure. LTC can also encounter difficulties in situations where chimerical clone consists of two parts that in fact belong to non-overlapping regions of a small neighborhood in the sequence. In addition, we found that most of chimerical clones that have not been detected by LTC actually consisted of one long part (100-150 kbp) belonging to the region of the contig and shorter part(s) (30-60 kbp) outside the region. We think that these chimerical clones can be identified by constructing more accurate maximal likelihood band map for LTC contig [37, 38].

General discussion

Although many physical maps constructed with the standard FPC algorithm have been successfully employed for genome sequencing, a substantial amount of errors were found in contig assembling [47] and subsequently corrected using complementary methods. Various factors, such as the genome composition (e.g., abundance of repeats), employed wet strategies and technologies of DNA cloning and clone fingerprinting, insufficient genome coverage of the BAC libraries, and a low accuracy of band scoring can affect the efficiency and accuracy of BAC assembly, hence the quality of the physical map. The power and flexibility of bioinformatics tools and human factors may also play a role: (i) from the beginning, reading the band sizes is somewhat subjective because it is very sensitive to selection of the threshold parameters for automatic peak detection; (ii) clustering (contig assembly) can be highly affected by false significant clone overlaps and chimerical clones; (iii) cluster ordering depends on the selection of an initial clone and the presence of equal-sizes but not common bands in different clones; (iv) decision-making in contig merging is difficult to formalize and the result depends on the choices made among a high number of possibilities; (v) anchoring (linking genetic and physical maps) is also mostly made by hand and verification of mapping results (at each stage) is needed.

Our results show that many of the previously mentioned difficulties can become less problematic with the use of new physical mapping algorithms and empowering of standard algorithms by additional tools. The quality of physical maps can be improved by a more accurate identification of chimerical clones and false clone overlaps. In particular this can be achieved by a more reliable scoring of the clone overlap p-value, by the utilization of "information content" of the bands, and by investigating the topology of nets of significant clone overlaps. Clustering can be improved by the utilization of clone overlaps with more liberal p-values: such overlaps result in larger clusters and hence simplify contig elongation and merging. Using more effective tools of contig ordering and ordering verification also improves contig quality and leads to shorter MTPs. More reliable multilocus genetic maps allow aligning the physical contigs to the correct chromosomal position more accurately thereby improving the quality of anchoring.

Conclusions

The proposed LTC methodology helps in obtaining more realistic (relevant) clusters of clones that are expected to have linear topological structure corroborating the linear structure of the eukaryotic chromosome. Analytical tools for clone ordering based on global optimization methods allow achieving coordinated clone orders and, presumably, shorter band map. Tools of LTC can also be used to "cure" and elongate contigs obtained with FPC or other methods. In particular, the proposed framework proved effective for fixing gaps in MTPs detected at the sequence level as well as for choosing alternative MTPs to increase the efficiency of sequencing.

We have shown that the LTC program has a few advantages over standard FPC in contig assembly. The disadvantages of LTC are in its less friendly interface, which lead to the possibility to work only with a Windows operating system and to a rather primitive output of the results. Some steps, such as, the building of band maps and automatic merging of contigs, the possibility of compound alignment and optimization of physical and genetic maps, are still not implemented in LTC. At this stage, it would make sense to utilize the advantages of these two packages in the following way: (i) construct the contigs with LTC at the beginning of a physical mapping project; (ii) for each resulting contig separately construct band maps using FPC; (iii) test FPC orders using LTC; (iv) resolve the detected problems by excluding clones disturbing the LTC-order in FPC-contig; (v) select a MTP for the verified contigs; (vi) if some gaps are found in the resulting MTP, try to cure them by changing or adding clones to the MTP using LTC; (vii) align contigs to a consensus genetic map [66, 67] using genetic markers using FPC. In addition, LTC can be used to check the results of physical maps already constructed by FPC and provide additional information or corrections by curing, elongating, and merging the FPC contigs. This can be useful in the final steps of a FPC-based project in which LTC can be used for selecting or curing the MTP.

Although our approximation of p-value seems to be more accurate than the standard Sulston score, better approximations are needed to estimate p-values for highly significant overlaps. This can be important for very large databases of clones (e.g., when mapping is based on genome-wide rather than chromosome-wide fingerprinting) especially for organisms with highly repeated genomes. To achieve stability of clone ordering we need to construct effective algorithms of excluding extra parallel clones for highly covered regions and buried clones without endangering contig connectivity. Ordering of clones is associated with the bands' relative positions in the chromosome. Hence, LTC procedure used for clone ordering can be improved by the parallel construction of a globally optimal map of bands. The process of merging contigs can be further optimized by better coordinating the "wet" and "dry" tools. We also think that criteria for the MTP selection can be improved by taking into account the number of bands and band abundances in clones and clone overlaps.

References

  1. Lander ES, Linton LM, Birren B, Nusbaum C, Zody MC, Baldwin J, Devon K, Dewar K, Doyle M, FitzHugh W, Funke R, Gage D, Harris K, Heaford A, Howland J, Kann L, Lehoczky J, LeVine R, McEwan P, McKernan K, Meldrim J, Mesirov JP, Miranda C, Morris W, Naylor J, Raymond C, Rosetti M, Santos R, Sheridan A, Sougnez C, et al.: Initial sequencing and analysis of the human genome. Nature 2001, 409(6822):860–921. 10.1038/35057062

    Article  CAS  PubMed  Google Scholar 

  2. McPherson JD, Marra M, Hillier L, Waterston RH, Chinwalla A, Wallis J, Sekhon M, Wylie K, Mardis ER, Wilson RK, Fulton R, Kucaba TA, Wagner-McPherson C, Barbazuk WB, Gregory SG, Humphray SJ, French L, Evans RS, Bethel G, Whittaker A, Holden JL, McCann OT, Dunham A, Soderlund C, Scott CE, Bentley DR, Schuler G, Chen HC, Jang W, Green ED, et al.: A physical map of the human genome. Nature 2001, 409(6822):934–941. 10.1038/35057157

    Article  CAS  PubMed  Google Scholar 

  3. Gregory SG, Sekhon M, Schein J, Zhao S, Osoegawa K, Scott CE, Evans RS, Burridge PW, Cox TV, Fox CA, Hutton RD, Mullenger IR, Phillips KJ, Smith J, Stalker J, Threadgold GJ, Birney E, Wylie K, Chinwalla A, Wallis J, Hillier L, Carter J, Gaige T, Jaeger S, Kremitzki C, Layman D, Maas J, McGrane R, Mead K, Walker R, et al.: A physical map of the mouse genome. Nature 2002, 418(6899):743–750. 10.1038/nature00957

    Article  CAS  PubMed  Google Scholar 

  4. Chen M, Presting G, Barbazuk WB, Goicoechea JL, Blackmon B, Fang G, Kim H, Frisch D, Yu Y, Sun S, Higingbottom S, Phimphilai J, Phimphilai D, Thurmond S, Gaudette B, Li P, Liu J, Hatfield J, Main D, Farrar K, Henderson C, Barnett L, Costa R, Williams B, Walser S, Atkins M, Hall C, Budiman MA, Tomkins JP, Luo M, et al.: An integrated physical and genetic map of the rice genome. Plant Cell 2002, 14(3):537–545. 10.1105/tpc.010485

    Article  PubMed  PubMed Central  Google Scholar 

  5. Coe E, Cone K, McMullen M, Chen SS, Davis G, Gardiner J, Liscum E, Polacco M, Paterson A, Sanchez-Villeda H, Soderlund C, Wing R: Access to the maize genome: an integrated physical and genetic map. Plant Physiol 2002, 128(1):9–12. 10.1104/pp.010953

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Nelson WM, Bharti AK, Butler E, Wei F, Fuks G, Kim H, Wing RA, Messing J, Soderlund C: Whole-genome validation of high-information-content fingerprinting. Plant Physiol 2005, 139(1):27–38. 10.1104/pp.105.061978

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Lewin HA, Larkin DM, Pontius J, O'Brien SJ: Every genome sequence needs a good map. Genome Res 2009, 19(11):1925–1928. 10.1101/gr.094557.109

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Staden R: A strategy of DNA sequencing employing computer programs. Nucleic Acids Res 1979, 6(7):2601–2610. 10.1093/nar/6.7.2601

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Roach JC, Boysen C, Wang K, Hood L: Pairwise end sequencing: a unified approach to genomic mapping and sequencing. Genomics 1995, 26(2):345–353. 10.1016/0888-7543(95)80219-C

    Article  CAS  PubMed  Google Scholar 

  10. Fleischmann RD, Adams MD, White O, Clayton RA, Kirkness EF, Kerlavage AR, Bult CJ, Tomb JF, Dougherty BA, Merrick JM: Whole-genome random sequencing and assembly of Haemophilus influenzae Rd. Science 1995, 269(5223):496–512. 10.1126/science.7542800

    Article  CAS  PubMed  Google Scholar 

  11. Adams MD, Celniker SE, Holt RA, Evans CA, Gocayne JD, Amanatides PG, Scherer SE, Li PW, Hoskins RA, Galle RF, George RA, Lewis SE, Richards S, Ashburner M, Henderson SN, Sutton GG, Wortman JR, Yandell MD, Zhang Q, Chen LX, Brandon RC, Rogers YH, Blazej RG, Champe M, Pfeiffer BD, Wan KH, Doyle C, Baxter EG, Helt G, Nelson CR, et al.: The genome sequence of Drosophila melanogaster. Science 2000, 287(5461):2185–2195. 10.1126/science.287.5461.2185

    Article  PubMed  Google Scholar 

  12. Voelkerding KV, Dames SA, Durtschi JD: Next-generation sequencing: from basic research to diagnostics. Clin Chem 2009, 55(4):641–658. 10.1373/clinchem.2008.112789

    Article  CAS  PubMed  Google Scholar 

  13. Li R, Fan W, Tian G, Zhu H, He L, Cai J, Huang Q, Cai Q, Li B, Bai Y, Zhang Z, Zhang Y, Wang W, Li J, Wei F, Li H, Jian M, Li J, Zhang Z, Nielsen R, Li D, Gu W, Yang Z, Xuan Z, Ryder OA, Leung FC, Zhou Y, Cao J, Sun X, Fu Y, et al.: The sequence and de novo assembly of the giant panda genome. Nature 2010, 463(7279):311–317. 10.1038/nature08696

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Metzker ML: Sequencing technologies - the next generation. Nat Rev Genet 2010, 11(1):31–46. 10.1038/nrg2626

    Article  CAS  PubMed  Google Scholar 

  15. Coulson A, Sulston J, Brenner S, Karn J: Toward a physical map of the genome of the nematode Caenorhabditis elegans. Proc Natl Acad Sci USA 1986, 83(20):7821–7825. 10.1073/pnas.83.20.7821

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Olson MV, Dutchik JE, Graham MY, Brodeur GM, Helms C, Frank M, MacCollin M, Scheinman R, Frank T: Random-clone strategy for genomic restriction mapping in yeast. Proc Natl Acad Sci USA 1986, 83(20):7826–7830. 10.1073/pnas.83.20.7826

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Sulston J, Mallett F, Staden R, Durbin R, Horsnell T, Coulson A: Software for genome mapping by fingerprinting techniques. Comput Appl Biosci 1988, 4(1):125–132.

    CAS  PubMed  Google Scholar 

  18. Gregory SG, Howell GR, Bentley DR: Genome mapping by fluorescent fingerprinting. Genome Res 1997, 7(12):1162–1168.

    CAS  PubMed  PubMed Central  Google Scholar 

  19. Marra MA, Kucaba TA, Dietrich NL, Green ED, Brownstein B, Wilson RK, McDonald KM, Hillier LW, McPherson JD, Waterston RH: High throughput fingerprint analysis of large-insert clones. Genome Res 1997, 7(11):1072–1084.

    CAS  PubMed  PubMed Central  Google Scholar 

  20. Zhang HB, Wing RA: Physical mapping of the rice genome with BACs. Plant Mol Biol 1997, 35(1–2):115–127. 10.1023/A:1005733909099

    Article  CAS  PubMed  Google Scholar 

  21. Ding Y, Johnson MD, Colayco R, Chen YJ, Melnyk J, Schmitt H, Shizuya H: Contig assembly of bacterial artificial chromosome clones through multiplexed fluorescence-labeled fingerprinting. Genomics 1999, 56(3):237–246. 10.1006/geno.1998.5734

    Article  CAS  PubMed  Google Scholar 

  22. Luo MC, Thomas C, You FM, Hsiao J, Ouyang S, Buell CR, Malandro M, McGuire PE, Anderson OD, Dvorak J: High-throughput fingerprinting of bacterial artificial chromosomes using the snapshot labeling kit and sizing of restriction fragments by capillary electrophoresis. Genomics 2003, 82(3):378–389. 10.1016/S0888-7543(03)00128-9

    Article  CAS  PubMed  Google Scholar 

  23. Venter JC, Smith HO, Hood L: A new strategy for genome sequencing. Nature 1996, 381(6581):364–366. 10.1038/381364a0

    Article  CAS  PubMed  Google Scholar 

  24. Prade RA, Griffith J, Kochut K, Arnold J, Timberlake WE: In vitro reconstruction of the Aspergillus (= Emericella) nidulans genome. Proc Natl Acad Sci USA 1997, 94(26):14564–14569. 10.1073/pnas.94.26.14564

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Mahairas GG, Wallace JC, Smith K, Swartzell S, Holzman T, Keller A, Shaker R, Furlong J, Young J, Zhao S, Adams MD, Hood L: Sequence-tagged connectors: a sequence approach to mapping and scanning the human genome. Proc Natl Acad Sci USA 1999, 96(17):9739–9744. 10.1073/pnas.96.17.9739

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Faris JD, Gill BS: Genomic targeting and high-resolution mapping of the domestication gene Q in wheat. Genome 2002, 45(4):706–718. 10.1139/g02-036

    Article  CAS  PubMed  Google Scholar 

  27. Weikard R, Goldammer T, Laurent P, Womack JE, Kuehn C: A gene-based high-resolution comparative radiation hybrid map as a framework for genome sequence assembly of a bovine chromosome 6 region associated with QTL for growth, body composition, and milk performance traits. BMC Genomics 2006, 7: 53. 10.1186/1471-2164-7-53

    Article  PubMed  PubMed Central  Google Scholar 

  28. Le Cunff L, Garsmeur O, Raboin LM, Pauquet J, Telismart H, Selvi A, Grivet L, Philippe R, Begum D, Deu M, Costet L, Wing R, Glaszmann JC, D'Hont A: Diploid/polyploid syntenic shuttle mapping and haplotype-specific chromosome walking toward a rust resistance gene (Bru1) in highly polyploid sugarcane (2n approximately 12x approximately 115). Genetics 2008, 180(1):649–660. 10.1534/genetics.108.091355

    Article  PubMed  PubMed Central  Google Scholar 

  29. Wang GL, Holsten TE, Song WY, Wang HP, Ronald PC: Construction of a rice bacterial artificial chromosome library and identification of clones linked to the Xa-21 disease resistance locus. Plant J 1995, 7(3):525–533. 10.1046/j.1365-313X.1995.7030525.x

    Article  CAS  PubMed  Google Scholar 

  30. Tanksley SD, Ganal MW, Martin GB: Chromosome landing: a paradigm for map-based gene cloning in plants with large genomes. Trends Genet 1995, 11(2):63–68. 10.1016/S0168-9525(00)88999-4

    Article  CAS  PubMed  Google Scholar 

  31. Huang L, Brooks SA, Li W, Fellers JP, Trick HN, Gill BS: Map-based cloning of leaf rust resistance gene Lr21 from the large and polyploid genome of bread wheat. Genetics 2003, 164(2):655–664.

    CAS  PubMed  PubMed Central  Google Scholar 

  32. Ma J, SanMiguel P, Lai J, Messing J, Bennetzen JL: DNA rearrangement in orthologous orp regions of the maize, rice and sorghum genomes. Genetics 2005, 170(3):1209–1220. 10.1534/genetics.105.040915

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Tameling WI, Baulcombe DC: Physical association of the NB-LRR resistance protein Rx with a Ran GTPase-activating protein is required for extreme resistance to Potato virus X. Plant Cell 2007, 19(5):1682–1694. 10.1105/tpc.107.050880

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Monden Y, Naito K, Okumoto Y, Saito H, Oki N, Tsukiyama T, Ideta O, Nakazaki T, Wessler SR, Tanisaka T: High potential of a transposon mPing as a marker system in japonica × japonica cross in rice. DNA Res 2009, 16(2):131–140. 10.1093/dnares/dsp004

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Soderlund C, Longden I, Mott R: FPC: a system for building contigs from restriction fingerprinted clones. Comput Appl Biosci 1997, 13(5):523–535.

    CAS  PubMed  Google Scholar 

  36. Soderlund C, Humphray S, Dunham A, French L: Contigs built with fingerprints, markers, and FPC V4.7. Genome Res 2000, 10(11):1772–1787. 10.1101/gr.GR-1375R

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Bhandarkar SM, Machaka SA, Shete SS, Kota RN: Parallel computation of a maximum-likelihood estimator of a physical map. Genetics 2001, 157(3):1021–1043.

    CAS  PubMed  PubMed Central  Google Scholar 

  38. Bhandarkar S, Huang J, Arnold J: Information theoretic approach to genome reconstruction. In Handbook of computational molecular biology. Edited by: Aluru S. Chapman & Hall/CRC Computer and Information Science Series; 2006:1–26.

    Google Scholar 

  39. Wang Y, Prade RA, Griffith J, Timberlake WE, Arnold J: A fast random cost algorithm for physical mapping. Proc Natl Acad Sci USA 1994, 91(23):11094–11098. 10.1073/pnas.91.23.11094

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Wendl MC: Probabilistic assessment of clone overlaps in DNA fingerprint mapping via a priori models. J Comput Biol 2005, 12(3):283–297. 10.1089/cmb.2005.12.283

    Article  CAS  PubMed  Google Scholar 

  41. Wendl MC: Algebraic correction methods for computational assessment of clone overlaps in DNA fingerprint mapping. BMC Bioinformatics 2007, 8: 127. 10.1186/1471-2105-8-127

    Article  PubMed  PubMed Central  Google Scholar 

  42. Balding DJ, Torney DC: Statistical analysis of DNA fingerprint data for ordered clone physical mapping of human chromosomes. Bull Math Biol 1991, 53(6):853–879.

    Article  CAS  PubMed  Google Scholar 

  43. Mayraz G, Shamir R: Construction of physical maps from oligonucleotide fingerprints data. J Comput Biol 1999, 6(2):237–252. 10.1089/cmb.1999.6.237

    Article  CAS  PubMed  Google Scholar 

  44. Meyers BC, Scalabrin S, Morgante M: Mapping and sequencing complex genomes: let's get physical! Nat Rev Genet 2004, 5(8):578–588. 10.1038/nrg1404

    Article  CAS  PubMed  Google Scholar 

  45. Marra M, Kucaba T, Sekhon M, Hillier L, Martienssen R, Chinwalla A, Crockett J, Fedele J, Grover H, Gund C, McCombie WR, McDonald K, McPherson J, Mudd N, Parnell L, Schein J, Seim R, Shelby P, Waterston R, Wilson R: A map for sequence analysis of the Arabidopsis thaliana genome. Nat Genet 1999, 22(3):265–270. 10.1038/10327

    Article  CAS  PubMed  Google Scholar 

  46. Flibotte S, Chiu R, Fjell C, Krzywinski M, Schein JE, Shin H, Marra MA: Automated ordering of fingerprinted clones. Bioinformatics 2004, 20(8):1264–1271. 10.1093/bioinformatics/bth079

    Article  CAS  PubMed  Google Scholar 

  47. Salzberg SL, Yorke JA: Beware of mis-assembled genomes. Bioinformatics 2005, 21(24):4320–4321. 10.1093/bioinformatics/bti769

    Article  CAS  PubMed  Google Scholar 

  48. Paux E, Sourdille P, Salse J, Saintenac C, Choulet F, Leroy P, Korol A, Michalak M, Kianian S, Spielmeyer W, Lagudah E, Somers D, Kilian A, Alaux M, Vautrin S, Berges H, Eversole K, Appels R, Safar J, Simkova H, Dolezel J, Bernard M, Feuillet C: A physical map of the 1-gigabase bread wheat chromosome 3B. Science 2008, 322(5898):101–104. 10.1126/science.1161847

    Article  CAS  PubMed  Google Scholar 

  49. Waterman MS, Griggs JR: Interval graphs and maps of DNA. Bull Math Biol 1986, 48(2):189–195.

    Article  CAS  PubMed  Google Scholar 

  50. Cuticchia AJ, Arnold J, Timberlake WE: The use of simulated annealing in chromosome reconstruction experiments based on binary scoring. Genetics 1992, 132(2):591–601.

    CAS  PubMed  PubMed Central  Google Scholar 

  51. Zhang P, Schon EA, Fischer SG, Cayanis E, Weiss J, Kistler S, Bourne PE: An algorithm based on graph theory for the assembly of contigs in physical mapping of DNA. Comput Appl Biosci 1994, 10(3):309–317.

    CAS  PubMed  Google Scholar 

  52. Arratia R, Lander ES, Tavare S, Waterman MS: Genomic mapping by anchoring random clones: a mathematical analysis. Genomics 1991, 11(4):806–827. 10.1016/0888-7543(91)90004-X

    Article  CAS  PubMed  Google Scholar 

  53. Jain AK, Dubes RC: Algorithms for clustering data. N.J.: Prentice Hall: Englewood Cliffs; 1988.

    Google Scholar 

  54. Batagelj V, Mrvar A: Pajek - analysis and visualization of large networks. In Graph drawing software. Edited by: Jünger M, Mutzel P. Berlin: Springer; 2003:77–103.

    Google Scholar 

  55. Bollobas B: Modern graph theory. 1st edition. Springer; 2002.

    Google Scholar 

  56. Fickett J, Cinkosky M: A genetic algorithm for assembling chromosome physical maps. In The Second International Conference on Bioinformatics. Supercomputing and Complex Genomic Analysis. Edited by: Lim H, Fickett J, Cantor C, Robbins R. New Jersey: World Scientific; 1992:273–285.

    Google Scholar 

  57. Alizadah F, Karp R, Newberg LA, Weisser D: Physical mapping of chromosomes: A combinatorial problem in molecular biology. Algorithmica 1995, 13(Numbers 1–2):52–76. 10.1007/BF01188581

    Article  Google Scholar 

  58. Mester D, Ronin Y, Minkov D, Nevo E, Korol A: Constructing large-scale genetic maps using an evolutionary strategy algorithm. Genetics 2003, 165(4):2269–2282.

    CAS  PubMed  PubMed Central  Google Scholar 

  59. Cuticchia AJ, Arnold J, Timberlake WE: ODS: ordering DNA sequences--a physical mapping algorithm based on simulated annealing. Comput Appl Biosci 1993, 9(2):215–219.

    CAS  PubMed  Google Scholar 

  60. Rand WM: Objective criteria for the evaluation of clustering methods. Journal of the American Statistical Association 1971, 66: 846–850. 10.2307/2284239

    Article  Google Scholar 

  61. Kamada T, Kawai S: An algorithm for drawing general undirected graphs. Information Processing Letters 1989, 31: 7–15. 10.1016/0020-0190(89)90102-6

    Article  Google Scholar 

  62. Xu Z, Sun S, Covaleda L, Ding K, Zhang A, Wu C, Scheuring C, Zhang HB: Genome physical mapping with large-insert bacterial clones by fingerprint analysis: methodologies, source clone genome coverage, and contig map quality. Genomics 2004, 84(6):941–951. 10.1016/j.ygeno.2004.08.014

    Article  CAS  PubMed  Google Scholar 

  63. Krzywinski M, Bosdet I, Mathewson C, Wye N, Brebner J, Chiu R, Corbett R, Field M, Lee D, Pugh T, Volik S, Siddiqui A, Jones S, Schein J, Collins C, Marra M: A BAC clone fingerprinting approach to the detection of human genome rearrangements. Genome Biol 2007, 8(10):R224. 10.1186/gb-2007-8-10-r224

    Article  PubMed  PubMed Central  Google Scholar 

  64. Nelson WM, Dvorak J, Luo MC, Messing J, Wing RA, Soderlund C: Efficacy of clone fingerprinting methodologies. Genomics 2007, 89(1):160–165. 10.1016/j.ygeno.2006.08.008

    Article  CAS  PubMed  Google Scholar 

  65. Safar J, Bartos J, Janda J, Bellec A, Kubalakova M, Valarik M, Pateyron S, Weiserova J, Tuskova R, Cihalikova J, Vrana J, Simkova H, Faivre-Rampant P, Sourdille P, Caboche M, Bernard M, Dolezel J, Chalhoub B: Dissecting large and complex genomes: flow sorting and BAC cloning of individual chromosomes from bread wheat. Plant J 2004, 39(6):960–968. 10.1111/j.1365-313X.2004.02179.x

    Article  CAS  PubMed  Google Scholar 

  66. Mester DI, Ronin YI, Korostishevsky MA, Pikus VL, Glazman AE, Korol AB: Multilocus consensus genetic maps (MCGM): formulation, algorithms, and results. Comput Biol Chem 2006, 30(1):12–20. 10.1016/j.compbiolchem.2005.09.007

    Article  CAS  PubMed  Google Scholar 

  67. Korol AB, Mester D, Frenkel Z, Ronin YI: Methods for genetic analysis in the Triticeae . In Genetics and Genomics of the Triticeae. Edited by: Feuillet C, Muehlbauer GJ. Springer + Business Media LLC; 2009:163–199. full_text

    Chapter  Google Scholar 

Download references

Acknowledgements

The research leading to these results has received funding from the European Community's Seventh Framework Programme (FP7/2007-2013) under the grant agreement FP7-212019 (TriticeaeGenome project, http://www.triticeaegenome.eu). It was also partially supported by Binational Agricultural Research and Development Fund (BARD research project US-3873-06) and by the Israeli Ministry of Absorption.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Zeev Frenkel.

Additional information

Authors' contributions

ZF designed algorithms, wrote the code for LTC package, conducted the analysis using LTC (for all presented datasets) and FPC packages (for simulated data only) and wrote the manuscript. EP provided the HICF data on 3B, conducted the FPC analysis for 3B, suggested the algorithm for simulations, proposed ideas for comparison of LTC and FPC methodologies and participated in preparing the manuscript. DM developed algorithms and wrote code for fast global optimization of clones ordering within the contigs. CF contributed to the discussion and participated in preparing the manuscript. AK supervised the project, proposed ideas for LTC methodologies, contributed to developing of the algorithms, helped to conduct the analysis and wrote the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

12859_2010_4167_MOESM1_ESM.DOC

Additional file 1: More detailed description of LTC approach. We have included the following items in an additional file named add_file_1.doc: (i) Description of p-value approximation for clone overlapping; (ii) Paragraph about comparison of LTC-approximation of p-value for clone overlapping, Sulston score, and mutual overlap statistic; (iii) Example illustrating clustering with "adaptively" varying cutoff; (iv) Examples of complications that one can meet in dealing with clusters having non-linear topological structure; and (v) A paragraph about comparison of clustering results obtained by different methods. (DOC 128 KB)

12859_2010_4167_MOESM2_ESM.DOC

Additional file 2: Supplemental figs. In this file we present figures that were not included into the main text of the paper. These figures clarify ideas of LTC approach and illustrate some of our results. (DOC 242 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution 2.0 International License (https://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Frenkel, Z., Paux, E., Mester, D. et al. LTC: a novel algorithm to improve the efficiency of contig assembly for physical mapping in complex genomes. BMC Bioinformatics 11, 584 (2010). https://doi.org/10.1186/1471-2105-11-584

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2105-11-584

Keywords