Comparing segmentations by applying randomization techniques
- Niina Haiminen^{1}Email author,
- Heikki Mannila^{1, 2} and
- Evimaria Terzi^{1}
DOI: 10.1186/1471-2105-8-171
© Haiminen et al; licensee BioMed Central Ltd. 2007
Received: 17 January 2007
Accepted: 23 May 2007
Published: 23 May 2007
Abstract
Background
There exist many segmentation techniques for genomic sequences, and the segmentations can also be based on many different biological features. We show how to evaluate and compare the quality of segmentations obtained by different techniques and alternative biological features.
Results
We apply randomization techniques for evaluating the quality of a given segmentation. Our example applications include isochore detection and the discovery of coding-noncoding structure. We obtain segmentations of relevant sequences by applying different techniques, and use alternative features to segment on. We show that some of the obtained segmentations are very similar to the underlying true segmentations, and this similarity is statistically significant. For some other segmentations, we show that equally good results are likely to appear by chance.
Conclusion
We introduce a framework for evaluating segmentation quality, and demonstrate its use on two examples of segmental genomic structures. We transform the process of quality evaluation from simply viewing the segmentations, to obtaining p-values denoting significance of segmentation similarity.
Background
Segmental structure of various scales exists in genomic sequences. Many evolutionary and genetic mechanisms leading to variation in DNA operate on segments of the genome (duplications and inversions of segments, recombinations). Furthermore, eukaryotic chromosomes consist of alternating regions of gene-rich and gene-poor regions. A gene-rich region can be further decomposed into non-coding segments, segments that contain regulatory information, and genes, which in turn consist of introns and exons. Also, remnants of viral or microbial inserts in a genome form a type of segmental structure.
There are many types of features with which one can segment the sequences. For any given technique, there may exist alternative biological features to segment on. For example, if the goal is to identify coding and noncoding segments in a sequence, one may study the distribution of three-letter words (codons) along the sequence to determine where to set the segment boundaries. An alternative would be, for example, to segment on the frequency of nucleotides in each third position, as done in [1].
There are also many possible segmentation techniques for discovering the segmental structure in DNA sequences. The techniques include recursive segmentation methods [1–7], Bayesian methods [8, 9], hidden Markov models [10–13], and wavelet analysis [14], among others. For a review of various approaches, see [15].
If a reliable annotation of the underlying segmental structure exists, it is of interest to find out which feature set and segmentation method give a result closest to the true segmental structure. This can give insight into the biological process that is responsible for creating and maintaining the segmental structure. Furthermore, the segmentation technique and feature set that is the best on an annotated sequence may also yield results close to true segmental structure on unannotated data.
In order to find differences between segmentations one has to define a notion of similarity or distance between them. In this paper we describe the distance between two segmentations P and Q by two numbers – conditional entropy of P given Q, and vice versa. For a similar approach to comparing clusterings see Meilă [26]. Conditional entropy is an information theoretic measure that quantifies the amount of information that one segmentation gives about the other. The sum of these two conditional entropies also defines the entropy distance between the two segmentations. Ideally, we want both terms in the sum to be small, instead of only requiring their sum to be small. We give an example of how this makes a difference when comparing segmentations. By using this measure we can easily rank a set of segmentations with respect to their distance from the underlying true segmentation, if one is known.
The sum of these two conditional entropies also defines the entropy distance between the two segmentations. Ideally, we want both terms in the sum to be small, instead of only requiring their sum to be small. We give an example of how this makes a difference when comparing segmentations. By using this measure we can easily rank a set of segmentations with respect to their distance from the underlying true segmentation, if one is known.
Knowing the best segmentation technique and feature set for a given sequence is still not enough. Namely, one candidate segmentation could be better than the other candidates, but all could still be quite far from the true one. That is, we want to find out if the best result is in some sense significant. The problem of deciding segmentation significance has been addressed before in the case where the sequence data itself is known [9]. Our approach does not rely on the sequence data, since it takes as input only the set of segmentations we want to evaluate. Our technique is therefore more general, and also applies to cases where the segmentations are obtained by using alternative biological features. To our knowledge, the issue of significance of segmentation similarity has not been considered before.
We test the significance of segmentation similarity by generating random segmentations and computing the distances between the underlying true segmentation and the randomized segmentations. If a random segmentation is about as close to the true one as our candidate segmentation is, then the agreement between the true segmentation and our candidate one is not very interesting. Note that in our randomization approach we do not make any assumptions regarding the properties of a good segmentation, but only consider the values of conditional entropies summarizing the similarity between two segmentations. To define a randomization procedure, we have to specify the class of segmentations from which we sample random elements. In this paper we use two classes of segmentation: (a) segmentations that have a given k number of segments, and (b) segmentations that have the same number of segments and the same segment-length distribution as the true segmentation.
We apply this randomization technique to examples on coding-noncoding structure and to isochore detection. The results show that the small distances obtained by some segmentation techniques and biological features are indeed significant, while for others, the obtained segmentations are only as similar to the ground truth as a majority of the randomly generated segmentations.
Results and discussion
In this section we show the results of our randomization techniques for evaluating the discovered segmentations in two examples of genomic sequences for which a segmental structure is already known.
Example 1: Coding-noncoding structure
Discovering the locations of genes in a DNA sequence is an important task to which computational methods give different predictions. In this example, we evaluate the closeness of different segmentation results to known gene boundaries.
We used a dataset consisting of a 25 kb region of bacterium Rickettsia prowazekii [16] (positions 535 – 560 kb), containing 13 coding segments. The correct number of segments in the data is 25 (13 genes separated by 12 noncoding segments [GenBank:AJ235269]). We denote the correct underlying segmentation into coding and noncoding regions by T.
We applied the encoding scheme described in [1] to transform the DNA sequence into a 12-dimensional signal. This encoding of the sequence captures the codon usage, see [1] and references therein. The segmentation techniques we applied are the entropic segmentation (E) by Bernaola-Galván et al. [1] and the least squares segmentation method (L). In the entropic segmentation the input is split recursively until no more significant splits can be found. The optimal least squares segmentation with k segments is computed with a dynamic-programming algorithm [17]. Note that the entropic segmentation decides on the number of segments, while the least-squares method takes the number of segments as input. The entropic segmentation method found 19 segments (at 99% significance level) for this dataset, outputting segmentation E_{ cd }. The least squares method for k = 25 output segmentation L_{ cd }.
We additionally applied the least squares method on six other features: frequencies of 2-letter words (AA,AC,...,TT), 1-letter words (A,C,G,T), and each nucleotide separately (A;C;G;T). We denote these features by {2,1, A, C, G, T} respectively. The sequences we constructed in such a way correspond to a 16-dimensional, a 4-dimensional and four 1-dimensional signals. Each point in these sequences corresponds to the frequency of a feature in an 10 bp window. We divided the sequence into 2500 non-overlapping windows.
Conditional entropies and their significances for segmentations on coding-noncoding data
P | H (T|P) | p _{ ℓ } | p _{ k } | H (P|T) | p _{ ℓ } | p _{ k } |
---|---|---|---|---|---|---|
L _{ cd } | 0.219 | 0 | 0 | 0.284 | 0 | 0 |
E _{ cd } | 0.286 | 0 | 0 | 0.143 | 0 | 0 |
L _{2} | 0.473 | 0 | 0 | 0.353 | 0 | 0 |
L _{1} | 0.566 | 0.0018 | 0.0005 | 0.364 | 0 | 0 |
L _{ A } | 0.879 | 0.1642 | 0.4219 | 0.466 | 0 | 0 |
L _{ C } | 1.212 | 1 | 0.9770 | 0.512 | 0 | 0.0001 |
L _{ G } | 0.738 | 0.0522 | 0.0732 | 0.849 | 0.1195 | 0.0661 |
L _{ T } | 0.455 | 0 | 0 | 0.289 | 0 | 0 |
Example 2: Isochore structure
Isochores are large (>> 300 kb) DNA segments fairly homogeneous in their guanine and cytosine (G+C) content. Exactly defining isochore borders in the human genome remains an open problem, but different computational approaches exist. Isochores are discussed in biological literature already since [18]. In our experiments we used a dataset consisting of the 100 Mb short arm telomeric region of human chromosome 1 ([19]). Unlike in the previous example, there is no well defined biological annotation of the segmental (isochore) structure of this genomic region. We consider as the ground truth segmentation T results by Costantini et al. [20] and alternatively results by IsoFinder [7]. These results differ considerably in their segment number (k = 114 vs k = 1305), reflecting genomic structures at different granularities. We also study segmentations on the major histocompatibility (MHC) region in chromosome 6, for which some biologically validated isochore boundaries are known. To generate candidate segmentations, we use the least squares dynamic programming algorithm (L) described in the previous example, with the same number of segments as the in respective T.
Conditional entropies and their significances for isochore segmentations with 114 segments
P | H (T|P) | p _{ ℓ } | p _{ k } | H (P|T) | p _{ ℓ } | p _{ k } |
---|---|---|---|---|---|---|
L _{ gc } | 0.740 | 0 | 0 | 0.683 | 0 | 0 |
L _{2} | 0.688 | 0 | 0 | 0.771 | 0 | 0 |
L _{1} | 0.750 | 0 | 0 | 0.778 | 0 | 0 |
L _{ A } | 0.937 | 0.0352 | 0.0389 | 0.876 | 0.0010 | 0.0090 |
L _{ C } | 0.797 | 0 | 0 | 0.684 | 0 | 0 |
L _{ G } | 0.836 | 0.0002 | 0.0004 | 0.821 | 0.0001 | 0.0005 |
L _{ T } | 0.984 | 0.1758 | 0.1748 | 0.840 | 0.0002 | 0.0010 |
Next we consider IsoFinder [7] results as the ground truth segmentation T, with k = 1305. We generated least squares candidate segmentations L_{ f }with this k and the same features f as in the previous case, the only difference being that here the sequence was aggregated into 10 kb windows (obtaining 10,000 data points). The segmentations are shown in Additional File 1, along with the chosen reference segmentation T. It is not clear from the figure which of the segmentations, if any, are close to T. The distributions of conditional entropies H (R|T) and H (T|R) for R ∈ ${C}_{N,k}$ for and R ∈ ${C}_{N,k,\ell}$ in 10,000 randomizations are also shown in the Additional File 1. In this case, the significance results depend on the choice of k-versus ℓ-randomization: if the former is chosen, all the candidate segmentations are significantly close to T. This means that the least squares segmentations are indeed significantly closer to T than segmentations with randomly assigned k boundaries. In the case of ℓ-randomization, we find that H (R|T) < H (L_{ f }|T) for all our candidate segmentations L_{ f }. Thus all segmentations L_{ f }are far from T. The reason is that the segment length distribution of T is far from the distribution for any L_{ f }, and this prevents the candidate segmentations from being significantly close to T.
We also study the major histocompatibility (MHC) region in chromosome 6 ([21]). The 3.67 Mb region contains two experimentally validated isochore structures [22]. The known MHC isochore boundaries are around 1.8 Mb, 2.5 Mb and 3.4 Mb. We studied predictions from wavelet multiresolution analysis by Wen and Zhang [14] (k = 7), results by Costantini et al. [20] (k = 6), and predictions by the tool IsoFinder [23] (k = 8), choosing in turn each one of these as the ground truth and thus performing three randomization experiments. We used least squares segmentations (with k = 7) as alternative candidate segmentations (as discussed earlier, with window size 100 kb). The segmentations and results from the randomization tests are shown in Additional File 2. The randomization tests show that all segmentations, except in some cases those for features A and T, are significantly close to the chosen ground truth segmentation. In particular, all the ground truth segmentations are found to be significantly close to each other.
Conclusion
In biological sequence analysis, there exist situations where many alternative segmentations for the underlying biological structure are proposed. We give a framework for evaluating the quality of results produced by different segmentation methods. Our approach also applies to cases where the segmentations are obtained by using alternative biological features, as we base our analysis only on the segment boundaries. Applicable segmentation distance measures and randomization tests are discussed, and results are shown for two applications of segmenting genomic data.
Methods
Comparing segmentations
We first give some intuition for our segmentation-comparison method and provide some basic definitions. Consider a sequence S with length |S| = N and a segmentation P of S. The segmentation partitions the sequence into k non-overlapping and contiguous intervals that span the whole sequence and they are called segments. A segmentation P with k segments can be fully defined using (k + 1) segment boundaries p_{0},...,p_{ k }, where p_{ i }∈ S, p_{ i }<p_{i+1}for every i, and p_{0} = 0 and p_{ k }= N. The i-th segment of P, denoted by ${\overline{p}}_{i}$, is defined to be the interval ${\overline{p}}_{i}$ = (p_{i-1}, p_{ i }]. Each segment ${\overline{p}}_{i}$ consists of |${\overline{p}}_{i}$| points that correspond to the length of the segment.
The maximum value that the entropy of a segmentation can have is log N.
That is, the conditional entropy of segmentation P given segmentation Q is the expected amount of information we need to identify the segment of P a point belongs to, given that we know the segment of this point in Q.
The following lemma gives an efficient algorithm for computing the conditional entropies between two segmentations. The algorithm runs in time O (k_{ p }+ k_{ q }).
Lemma 1. Let P and Q be two segmentations. Denote by U their union, i.e., the segmentation defined by the segment boundaries that appear in P or in Q. The conditional entropy of P given Q, H (P|Q), can be computed using the following closed formula
H (P|Q) = H (U) - H (Q).
Intuitively, the entropy of P given Q tells us how much information we obtain about P, if we know that that we are in a specific segment of Q. The more information Q reveals about the structure P, the more similar segmentations P and Q are, and the smaller the value of H (P|Q).
The single value H (P|Q) does not give the whole picture of segmentation similarity, however. For example, consider the case where segmentation Q consists of a single segment. Then, using Lemma 1 we can verify our intuition that knowledge about Q gives us no information about P, i.e., H (P|Q) = H (P). However, notice that H (Q|P) = 0, for any P.
Consider also the case where Q consists of N segments where each segment has length 1. In this case, H (P|Q) = 0, that is, Q gives lots of information about P, irrespective of the structure of P. As before, observe that H (Q|P) = log N - H (P). Thus, if P has low entropy, the value of H (Q|P) is large.
The above examples show that the similarity of two segmentations P and Q cannot be judged just by the single value H (P|Q) or H (Q|P). Rather, we can conclude that segmentations P and Q are similar only if both H (P|Q) and H (Q|P) are small. Even when using the entropy distance between two segmentations (see [25]), defined as D_{ H }(P, Q) = H (P|Q) + H (Q|P), we can get small values to segmentations that are quite different. We show in the experimental section that considering the two conditional entropies separately gives more accurate results than using their sum.
Randomization techniques
Consider a segmentation algorithm that given as input sequence S outputs a segmentation P. The plethora of segmentation algorithms and segmentation criteria naturally raises the question of how good and how informative segmentation P is. Assume that we a priori know a ground-truth segmentation T of S. Then, we can say that segmentation P is good if P is similar to T. Thus, using the definitions in the Methods section, P is a good segmentation if H (P|T) and H (T|P) are small. However, a natural question is how small is small enough? Or, is there a threshold in the values of the conditional entropies below which we can characterize segmentation P as being correct or interesting? Finally, can we set this threshold universally for all segmentations? In this section we describe a set of randomization techniques that we devise in order to provide an answer to these questions.
Our generic methodology is the following. Given a segmentation P and a ground-truth segmentation T of the same sequence, we first compute H (P|T) and H (T|P). We compare the values of these conditional entropies with the values of the conditional entropies H (R|T) and H (T|R) for a random segmentation R. We conclude that P is similar to T, and thus interesting, if the values of H (P|T) (and H (T|P)) are small compared to the values of H (R|T) (and H (T|R)) for a large majority of random segmentations R.
Consider a class ${C}_{N}$ of segmentations for sequences of length N. Then, the randomization test is conducted as follows. Pick random segmentations R ∈ ${C}_{N}$. For each such R compute H (T|R), and compare H (T|P) against the distribution of the values H (T|R). Similarly, compute the values of H (R|T) for a large number of segmentations R ∈ ${C}_{N}$ and compare these values with the value of H (P|T). In general, if segmentations T and P have a very different number of of segments, one of H (T|R) and H(R|T) will be large for any R from ${C}_{N}$. The randomization method we describe is best suited for the case when T and P have about the same number of segments.
We still need to specify the class ${C}_{N}$ of segmentations from which the random segmentations are picked. We define two classes of segmentations. Intuitively, the first class is used for checking if the candidate segmentation P is significantly closer to T than random segmentations R with the same number of segments as in T. Imagine that the segmentation procedure that generated P has knowledge of the segment number in T. By using this class we find out if it is enough to guess a segmentation as close to T as P is, by just randomly assigning a correct number of segment boundaries. The other class is used similarly, for checking if the knowledge of T's segment length distribution is enough to generate segmentations as close to T as P. An analog is found in classification problems, where the true class labels in T are permuted to check if the candidate classification P offers more insight into T than we would expect from guessing a random classification R.
In the first case, if T has k segments, then we restrict the random segmentations to those that have k segments as well. We denote by ${C}_{N,k}$ the class of all segmentations with k segments that partition sequences of length N. We have $\left|{C}_{N,k}\right|=\left(\begin{array}{c}N\\ (k-1)\end{array}\right)$, since there are $\left(\begin{array}{c}n\\ (k-1)\end{array}\right)$ ways to choose k - 1 segment boundaries from the N points of the sequence (the first and the last boundary are always fixed). We call the randomization test in which the random segmentations R ∈ ${C}_{N,k}$ a k-randomization test.
We also introduce the ℓ-randomization test, specified as follows. Consider a segmentation T with segments $\{{\overline{t}}_{1},\mathrm{...},{\overline{t}}_{k}\}$. Each segment ${\overline{t}}_{i}$ has length |${\overline{t}}_{i}$| and these lengths define a distribution of the segment lengths of the segmentation. There are a total of at most k! segmentations that have the same segment lengths as T does (maximized when all segments in T have different lengths). For a given distribution of segment lengths ℓ we denote by ${C}_{N,k,\ell}$ the class of segmentations with k segments and lengths ℓ. Obviously, and since the segmentations in . differ only in the order in which the segments with different lengths appear. Note that for a random segmentation R ∈ ${C}_{N,k,\ell}$, the conditional entropies w.r.t. segmentation T are equal, i.e., H (T|R) = H (R|T). This follows from the fact that H (T) = H (R), since both segmentations contain exactly the same segments. The ℓ-randomization restricts our attention to random segmentations with segment length distribution being the same as the segment-length distribution of the ground-truth segmentation T. That is, the significance of H (T|P) and H (P|T) for a candidate segmentation P are evaluated under the assumption that the segment-length distribution is known. In the special case where all segments in T have the same length, the segment length distribution uniquely characterizes the segmentations R ∈ ${C}_{N,k,\ell}$. In this case H (T|R) = H (R|T) = 0. Moreover, any segmentation P that does not have equal-length segments has H (P|T) > 0 and H (T|P) > 0 and thus is considered far from the ground truth w.r.t. the ℓ-randomization test.
Declarations
Authors’ Affiliations
References
- Bernaola-Galván P, Grosse I, Carpena P, Oliver JL, Román-Roldán R, Stanley HE: Finding borders between coding and noncoding DNA regions by an entropic segmentation method. Phys Rev Letters. 2000, 85 (6): 1342-1345.View ArticleGoogle Scholar
- Bernaola-Galván P, Román-Roldán R, Oliver JL: Compositional segmentation and long-range fractal correlations in DNA sequences. Phys Rev E. 1996, 53 (5): 5181-5189.View ArticleGoogle Scholar
- Oliver JL, Román-Roldán R, Pérez J, Bernaola-Galván P: SEGMENT: identifying compositional domains in DNA sequences. Bioinformatics. 1999, 15 (12): 974-979.View ArticlePubMedGoogle Scholar
- Li W: Delineating relative homogeneous G+C domains in DNA sequences. Gene. 2001, 276 (1–2): 57-72.View ArticlePubMedGoogle Scholar
- Li W: DNA segmentation as a model selection process. Proc RECOMB '01, ACM Press. 2001, 204-210.View ArticleGoogle Scholar
- Li W, Bernaola-Galván P, Haghighi F, Grosse I: Applications of recursive segmentation to the analysis of DNA sequences. Computers and Chemistry. 2002, 26 (5): 491-510.View ArticlePubMedGoogle Scholar
- Oliver JL, Carpena P, Hackenberg M, Bernaola-Galván P: IsoFinder: computational prediction of isochores in genome sequence. Nucl Acids Res. 2004, 32: W287-W292.PubMed CentralView ArticlePubMedGoogle Scholar
- Liu JS, Lawrence CE: Bayesian inference on biopolymer models. Bioinformatics. 1999, 15 (1): 38-52.View ArticlePubMedGoogle Scholar
- Husmeier D, Wright F: A Bayesian approach to discriminate between alternative DNA sequence segmentations. Bioinformatics. 2002, 18 (2): 226-234.View ArticlePubMedGoogle Scholar
- Churchill GA: Stochastic models for heterogeneous DNA sequences. Bull Math Biol. 1989, 51 (1): 79-94.View ArticlePubMedGoogle Scholar
- Peshkin L, Gelfand MS: Segmentation of yeast DNA using hidden Markov models. Bioinformatics. 1999, 15 (12): 980-986.View ArticlePubMedGoogle Scholar
- Melodelima C, Guéguen L, Piau D, Gautier C: Prediction of human isochores using a hidden Markov model. Proc JOBIM. 2005, 427-434.Google Scholar
- Guéguen L: Sarment: Python modules for HMM analysis and partitioning of sequences. Bioinformatics. 2005, 21 (16): 3427-3428.View ArticlePubMedGoogle Scholar
- Wen SY, Zhang CT: Identification of isochore boundaries in the human genome using the technique of wavelet multiresolution analysis. Biochem and Biophys Res Comm. 2003, 311 (1): 215-222.View ArticleGoogle Scholar
- Braun JV, Müller HG: Statistical methods for DNA sequence segmentation. Statistical Science. 1998, 13 (2): 142-162.View ArticleGoogle Scholar
- Andersson SG, Zomorodipour A, Andersson JO, Sicheritz-Ponten T, Alsmark UC, Podowski RM, Naslund AK, Eriksson AS, Winkler HH, Kurland CG: The genome sequence of Rickettsia prowazekii and the origin of mitochondria. Nature. 1998, 396 (6707): 133-140.View ArticlePubMedGoogle Scholar
- Bellman R: On the approximation of curves by line segments using dynamic programming. Commun ACM. 1961, 4 (6): 284-View ArticleGoogle Scholar
- Macaya G, Thiery JP, Bernardi G: An approach to the organization of eukaryotic genomes at a macromolecular level. J Mol Biol. 1976, 108 (1): 237-254.View ArticlePubMedGoogle Scholar
- Lander ES, Linton LM, Birren B: Initial sequencing and analysis of the human genome. Nature. 2001, 409 (6822): 860-921.View ArticlePubMedGoogle Scholar
- Costantini M, Clay O, Auletta F, Bernardi G: An isochore map of human chromosomes. Genome Research. 2006, 16 (4): 536-541.PubMed CentralView ArticlePubMedGoogle Scholar
- The MHC sequencing consortium: Complete sequence and gene map of a human major histocompatibility complex. Nature. 1999, 401 (6756): 921-923.Google Scholar
- Fukagawa T, Sugaya K, Matsumoto KI, Okomura K, Ando A, Inoko H, Ikemura T: A boundary of long-range G+C% mosaic domains in the human MHC locus: pseudoautosomal boundary-like sequence exists near the boundary. Genomics. 1995, 25 (1): 184-191.View ArticlePubMedGoogle Scholar
- Oliver JL, Bernaola-Galvan P, Carpena P, Roman-Roldan R: Isochore chromosome maps of eukaryotic genomes. Gene. 2001, 276 (1–2): 47-56.View ArticlePubMedGoogle Scholar
- Cover TM, Thomas JA: Elements of Information Theory. 1991, Wiley Series in TelecommunicationsView ArticleGoogle Scholar
- Mielikäinen T, Terzi E, Tsaparas P: Aggregating time partitions. Proc KDD '06, ACM Press. 2006, 347-356.View ArticleGoogle Scholar
- Meilă M: Comparing clusterings: an axiomatic view. Proc. ICML '05. 2005, ACM Press, 577-584.View ArticleGoogle Scholar
Copyright
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.