# A novel hierarchical clustering algorithm for gene sequences

- Dan Wei
^{1, 2}, - Qingshan Jiang
^{2}Email author, - Yanjie Wei
^{2}Email author and - Shengrui Wang
^{3}

**13**:174

https://doi.org/10.1186/1471-2105-13-174

© Wei et al.; licensee BioMed Central Ltd. 2012

**Received: **5 November 2011

**Accepted: **30 June 2012

**Published: **23 July 2012

## Abstract

### Background

Clustering DNA sequences into functional groups is an important problem in bioinformatics. We propose a new alignment-free algorithm, mBKM, based on a new distance measure, DMk, for clustering gene sequences. This method transforms DNA sequences into the feature vectors which contain the occurrence, location and order relation of *k*-tuples in DNA sequence. Afterwards, a hierarchical procedure is applied to clustering DNA sequences based on the feature vectors.

### Results

The proposed distance measure and clustering method are evaluated by clustering functionally related genes and by phylogenetic analysis. This method is also compared with BlastClust, CD-HIT-EST and some others. The experimental results show our method is effective in classifying DNA sequences with similar biological characteristics and in discovering the underlying relationship among the sequences.

### Conclusions

We introduced a novel clustering algorithm which is based on a new sequence similarity measure. It is effective in classifying DNA sequences with similar biological characteristics and in discovering the relationship among the sequences.

## Keywords

## Background

With the development of advanced biotechnology, more and more biological sequence information has been generated. The amount of genetic data is growing faster than the rate at which it can be analyzed. Clustering techniques provide a viable solution for handling and analyzing such rapidly growing genetic data. Clustering algorithms partition sequences into different biologically meaningful groups, facilitating therefore the prediction of functions of genes [1]. When a new gene is assigned to a cluster, the biological function of this cluster can be attributed to this gene with high confidence. On the other hand, clustering gene sequences into groups may also help with analyzing evolutionary relationships among the sequences in a cluster [2].

Clustering of gene sequences requires calculation of similarity between sequences. There are two clustering approaches according to the similarity measure used in a clustering method. One is based on sequence alignment. The similarity between two gene sequences is measured by the scores obtained from an alignment algorithm such as BLAST [3] or FASTA [4]. Although sequence alignment gives good solutions, it is relatively difficult to cluster a large number of sequences because of its computational complexity. Moreover, if the sequences in the set vary in length, a satisfactory alignment is hard to achieve, resulting in a low accuracy of clustering.

The other approach for similarity measure is to use alignment-free methods [5–10]. In recent years, several alignment-free measures have been proposed. The word-based measure is one of the most widely used methods [11–14]. This method chooses a short word length *k*, maps each sequence onto an *n*-dimensional vector according to its *k*-length tuple (also called *k*-tuple or *k*-word) properties, and then assesses the similarity of any two vectors by measures such as Euclidean distance [15], Mahalanobis distance [16], Kullback–Leibler discrepancy [17], cosine distance [18] or Pearson’s correlation coefficient [19]. In recent years, several novel alignment-free measures [20, 21] have been designed for DNA sequences analysis. Yang et al. [22] extended the *k*-tuple distance, which is based on the difference in tuple frequencies, to clustering gene sequences. Their tuple-based method determines the similarity of sequences by considering only tuple frequencies and ignoring the positional information within a sequence.

Major algorithms used in gene sequence clustering can be divided into two categories according to the result format: hierarchical clustering algorithms and partitional clustering algorithms [23]. Hierarchical clustering is widely used for detecting clusters in genomic data. It generates a set of partitions forming a cluster hierarchy. According to linkage criteria, there are three hierarchical clustering methods including single-linkage clustering (SL), complete-linkage clustering (CL) and average-linkage clustering (AL) [24]. With SL, clusters may be merged together due to single sequences being close to each other, even though many of the sequences in each cluster may be very distant to each other [25]. CL tends to find compact clusters of approximately equal diameters [25]. With CL, all objects in a cluster are similar to each other. AL can be seen as an intermediate between single and complete linkage clustering, resulting in more homogeneous clusters than those obtained by the single-linkage method [26]. For instance, BlastClust [27] and GeneRage [28] employ single linkage clustering approach; SWORDS [29] is based on word frequencies as profiles to merge clusters hierarchically; and Uchiyama [30] use average linkage clustering algorithm to classify genes. Hierarchical approaches may yield fairly good results, but they require the similarity of all pairs of sequences and quickly arrive at a bottleneck in terms of computational time and memory usage for large-scale data sets [31].

Partitioning algorithms have also been used. Partitional clustering obtains a partition of data objects by optimizing some clustering criterion. Partitional clustering algorithms are simple and well-suited for clustering large datasets [32]. K*-*means (KM) [33, 34] is a commonly used method of partitional clustering methods. KM has a lower order of computational complexity and demands less physical memory than the hierarchical method. It is suitable for clustering large gene data. Some KM-based algorithms, such as those introduced by Wan et al. [33], Kelarev et al. [34], Tseng et al. [35] and Ashlock et al. [36], have been developed to group DNA sequences. The major drawback of KM compared to hierarchical clustering algorithms is the lack of hierarchical relationships in its results. To remedy the problem, bisecting K-means (BKM), a hierarchical variation of KM, was proposed to build a tree of clusters in a top-down fashion by splitting the least homogeneous cluster into two more homogeneous ones. BKM can produce either a flat clustering or a hierarchical clustering by recursively applying KM. It has a linear complexity and is relatively efficient and scalable. Recent study [37] concluded that BKM outperforms KM and performs equally well or better than hierarchical methods when it partitions the dataset based on a homogeneity criteria. The bisecting approach is very attractive for genomic studies [38].

Hierarchical clustering produces a nested series of partitions, where the results are usually depicted as a dendrogram while partitional clustering produces a flat partition. BlastClust [27] is a hierarchical clustering method based on BLAST scores as the measure of sequence similarity. BlastClust computes pairwise similarity of all sequences by BLAST alignment and then clusters sequences by the single linkage clustering method which produces clusters of linear topology. The performance of BlastClust is limited by the size of the input data. CD-HIT-EST [39], a partitional approach, is also widely used to cluster DNA sequences. CD-HIT-EST uses an incremental clustering process and avoids the unnecessary alignments by a short word filtering mechanism, which detects similar sequences by counting the number of identical short words between them. The purpose of filters is to decide whether the identity between two sequences is above or below a threshold without aligning them, therefore speeding up the clustering process. Though CD-HIT-EST is based on alignment, it can avoid too many pairwise alignments by using a filter, thus it is faster than BlastClust, and can handle larger datasets.

Recent studies reveal also that BlastClust is less effective for clustering divergent sequences [40], and its performance strongly depends on the choice of optimal BLAST parameters including similarity threshold, percent identity, and alignment length [41]. CD-HIT-EST, on the other hand, does not provide hierarchical relationships between clusters of sequences. In many situations both CD-HIT-EST and BlastClust yield clusters with only one sequence [41]. All the traditional clustering methods based on sequence alignment encounter computational difficulties in dealing with large biological databases.

The approach presented in this paper involves a new alignment-free distance measure based on *k*-tuples, DMk (Distance Measure based on *k*-tuples) [42], and a modified bisecting K-means clustering algorithm, mBKM (modified Bisecting K-Means algorithm). mBKM aims to speed up the clustering process by using the alignment-free similarity measure, and is able to produce either a hierarchical clustering or a partition clustering result. We have applied mBKM with DMk in clustering gene sequences and performing phylogenetic analysis. DMk shows better performance than the *k*-tuple distance in our experiments, and mBKM outperforms SL, CL, AL, BKM and KM when tested on public gene sequence datasets. Furthermore, the proposed method also outperforms alignment-based methods such as BlastClust and CD-HIT-EST.

## Methods

A gene is a stretch of DNA that codes for a single polypeptide chain [43]. A gene sequence is a succession of four symbols {A, C, G, T}. Because the similarity between the genes of two species indicates their evolutionary relationship, it is used in many clustering algorithms. The goal of sequence clustering is to partition biological sequences into meaningful/functional groups according to the similarity information, which is calculated using either an alignment-based method or an alignment-free method.

*S*

_{1}= AGCACACA and

*S*

_{2}= ACACAGTA,

*S*

_{1}

^{ P }and

*S*

_{2}

^{ P }are used to represent the

*p*

^{th}characters in

*S*

_{1}and

*S*

_{2}, respectively. The alignment score [45] for (

*S*

_{1}

*S*

_{2}) is given by

where *E* is the cost of an alignment operation: deletion, substitution, or insertion. However this distance measure relies on sequence alignment. Since sequence alignment suffers in computational aspect with regard to large biological databases, clustering methods relying on sequence alignment have difficulties in dealing with the large gene data. An alignment-free similarity measure helps avoid the computational complexity of multiple sequence alignment for similarity computation. In this paper we propose a new alignment-free similarity measure, DMk, based on which we developed mBKM to cluster gene sequences.

In the follows, we will present DMk first, and then describe mBKM algorithms.

### A new similarity measure: DMk

In this section, we introduce a new similarity measure which takes into account the occurrence, location and order relation of *k*-tuple in a DNA sequence.

Sequences are numerically transformed to feature vectors that can be processed by data mining algorithms. Let Σ be the alphabet set of nucleotides (Σ = {A, C, G, T}). A sequence of length *s*, *S*, is defined as a linear succession of *s* symbols from Σ. A segment of *k* consecutive symbols in sequence *S* (*k ≤ s*) is designated as a *k*-tuple. There is a set of 4^{
k
} possible *k*-tuples, *W*_{
k
}*.* The number of occurrences of a *k*-tuple *w*, *N*_{
w
}, is counted by moving a sliding window of length *k* over the sequence with *k* - 1 bp overlapping step size.

*k*-tuple occur in the sequence to explore the sequence structure. For a DNA sequence

*S*

*p*

_{ r }is the location of the

*r*

^{ th }occurrence of

*k*-tuple

*w*, where

*p*

_{0}= 0. And

*α*

_{ r }is given as,

in which *m* stands for the number of occurrences of *w*. *α*_{
r
} reflects the density of *w* and is closely related to the location where *w* occurs in the sequence. Each *w* begins at the 1/*α*_{1} position, and {*α*_{
1,
}*α*_{
2
}*,…,α*_{
m
}} for repetition of *w* forms an array whose *r*^{
th
} element indicates the relative position of two neighboring *w* in the sequence. This array allows us to find all subsequent repeats of *w*.

*α*

_{ r }, we define

*β*

_{ j }as a partial sum of {

*α*

_{ r }}.

*β*

_{ j }is calculated by the following formula:

{*α*_{
r
}} is a list of non-negative real numbers, and *β*_{
j
} is totally ordered by ≤, so *β*_{1}, *β*_{2}, …, *β*_{
m
} is also an ordered set. {*α*_{1}, *α*_{2},…, *α*_{
m
}} and {*β*_{1}, *β*_{2},…, *β*_{
m
}} determine each other uniquely. *β*_{
j
} is only dependent of the number and positions of *w* and independent on other *k*-tuples. Given the set of {*β*_{1}, *β*_{2},…, *β*_{
m
}}, one can obtain where *w* occurs and how many times *w* occurs in the sequence.

*Shannon’s entropy*[48], which illuminates the total information measure of source on the average, is a measure of order/disorder. According to [49], when using the totally ordered set {

*β*

_{1}

*β*

_{2},…,

*β*

_{ m }} to calculate the probabilities, the Shannon entropy reflects the degree of importance of position in a sequence. We construct a discrete probability distribution $Q=({q}_{1},{q}_{2},\dots ,{q}_{m})$, where ${q}_{i}={\beta}_{i}/\sum _{i=1}^{m}{\beta}_{i}$, and $\sum _{i=1}^{m}{q}_{i}=1$. The

*Shannon entropy*of the discrete probability distribution is calculated by

For each *k*-tuple *w* in the sequence, not only the information of tuple numbers but also the information of tuple positions is involved in the definition of *H*. We take *H* as the feature of *w* in the sequence, and then construct a vector consisted of *H* of all possible *k*-tuples in the given sequence.

For a fixed *k*, there are 4^{
k
} distinct *k*-tuples to be considered. These *k*-tuples in a fixed 4^{
k
}-dimension feature vector are denoted by $({H}_{1},{H}_{2},\dots ,{H}_{{4}^{k}})$, where *H*_{
i
} means the feature representation of the *i*th *k*-tuple. This feature vector based on *H* can be regarded as an index for its corresponding sequence.

*k*-tuple distance is the sum of the differences in frequency over all possible

*k*-tuples; on the other hand, we use Euclidean distance between Shannon entropy of

*k*-tuples in sequences to measure the similarity. This distance measure method is referred as DMk. For any two sequences

*X*and

*Y*, DMk can be calculated as:

where ${h}_{{w}_{i}}^{X}$ and ${h}_{{w}_{i}}^{Y}$ represent the Shannon entropy values of the *i*^{
th
}*k*-tuple in sequences *X* and *Y*, respectively. DMk can be calculated from following algorithm:

Algorithm Name: DMk for similarity measure

Input: sequences {*S*_{1}, *S*_{2},…, *S*_{
N
}}.

Output: similarity matrix, (*d*(*X*,*Y*))_{
N*N
}.

- 1.
For each sequence, search and locate each

*k*-tuple;1.1 For each

*k*-tuple, use Equation (1) to calculate ${\alpha}_{r}(1\le r\le m)$1.2 For each

*k*-tuple, use Equation (2) to calculate ${\beta}_{j}(1\le j\le m)$;1.3 For each

*k*-tuple, use Equation (3) to calculate*H*; - 2.
For each sequence, construct 4

^{ k }-component vector by*H*of all*k*-tuples. - 3.
For any two sequences, use Equation (4) to calculate the distance between the two sequences.

- 4.
Return {

*d*}.

### A new clustering algorithm: mBKM

KM can be used to obtain a hierarchical clustering solution using a repeated bisecting approach [50, 51]. BKM is such an algorithm and it can produce either a partitional or a hierarchical clustering.

BKM has a linear time complexity in each bisecting step. Recent study [51] concludes BKM outperforms KM as well as the agglomerative approach in terms of accuracy and efficiency. Consequently, the bisecting approach is very attractive in many applications for clustering and genomic data analysis.

*K*clusters are obtained. The outcome is structured as a binary tree. There are two key steps in a typical BKM. The first one is the selection of initial centroids. Generally the initial centroids are chosen randomly in BKM. The second key step is the rule,

*ζ*, for selection of a existing cluster to be split in each bisecting step.

*ζ*is typically given by the following three approaches [50]:

- 1)
Choosing the cluster with largest size;

- 2)Selecting the cluster with the overall similarity$\frac{1}{{\left|C\right|}^{2}}\sum _{\begin{array}{l}s\in C\\ s\text{'}\in C\end{array}}d(s,s\text{'})$(6)

*d*(

*s*,

*s*

^{’}).

*C*is a cluster;

- 3)
Using a criterion based on both size and overall similarity.

Because the differences between these methods are small in terms of the final clustering result, the way of splitting the largest remaining cluster is recommended [50].

- 1.
Randomly choosing the initial centroids in BKM may result in too adjacent elements selected. If the initial centroids are too close, the algorithm will reach a local optimization. Moreover, different sets of initial cluster centroids can lead to different final clustering results.

- 2.
The algorithm for choosing one existing cluster to split in each bisecting step usually selects the cluster with the largest size. Although this leads to reasonably good and balanced clustering solution, it cannot gracefully work for datasets where the natural clusters are of different sizes, as it will tend to partition larger clusters first. In real biological data, the number of elements in every cluster may not always be similar.

- 1)
Selecting Initial Cluster Centroids

*s*

_{1},

*s*

_{2}, …,

*s*

_{ N }}, let

*d*(

*s*

_{ i },

*s*

_{ j })(

*i*,

*j*= 1, 2, …

*N*) be the distance between any two sequences in the dataset. We choose the sequence ${s}_{{c}_{1}}$ and ${s}_{{c}_{2}}$ as the cluster centroid according the following rule:

- 2)
Selecting the Cluster to Split

BKM algorithm usually partitions the largest size cluster into two smaller ones and yields clusters with similar size. However, a cluster with large number is not always the loose one. If one existing cluster is a loose one, in which its members are not closely related to each other, the cluster will be selected to be split.

*C*

_{ j }is defined as following:

where *μ*_{
j
} is the centroid of sequences in *C*_{
j
}, *d* (*s*_{
i
}, *μ*_{
j
}) is the distance between *s*_{
i
} and *μ*_{
j
}, and *n*_{
j
} is the number of sequences in the cluster.

A small variance of a cluster indicates that the members in the cluster tend to be closely related to the mean. In other words, the smaller the variance is, the more compact the cluster is, and vice versa.

Based on the above idea, we outline mBKM algorithm as follows.

Algorithm Name: mBKM for clustering sequences

Input: sequences {*s*_{1}, *s*_{2}, …, *s*_{
N
}}, a distance function *d* between sequences, the number of clusters *K*.

Output: Set of *K* clusters.

- 1.
Initialization: Regard the whole dataset {

*s*_{1},*s*_{2}, …,*s*_{ N }} as a single cluster. - 2.
Pick a cluster to split.

- 3.
Find two sub-clusters:

3.1 Select two initial centroids using Equation (6);

3.2 Assign the sequences to the closest centroid;

3.3 Recalculate two centroids based on the sequences assigned to the cluster;

3.4 Repeat steps 3.2 and 3.3 until no change in cluster centroid calculation.

- 4.
Calculate the variance of each cluster according Equation (7) and take the split that produces the clustering result with the highest variance.

- 5.
Repeat steps 2, 3 and 4 until the desired number

*K*is reached.

This algorithm outputs a binary tree of sequences, where each leaf represents a sequences and each node represents a sequence collection.

## Results and discussion

The proposed method is evaluated by clustering functionally related gene sequences and by phylogenetic analysis. We present our evaluation results in two parts. The first one aims at testing the efficiency of our similarity measure, DMk. The second one is to illustrate the efficiency of the proposed clustering method, mBKM.

*j*and class

*i*

*F*(

*i*

*j*) is defined as:

where *i* =1, 2, …, *e* *j* = 1, 2, …, *f* *precision*(*i*, *j*) = *n*_{
ij
}/*n*_{
j
}*recall*(*i*, *j*) = *n*_{
ij
}/*n*_{
i
}*e* is the number of classes, and *f* is the number of clusters. *n*_{
ij
} is the number of the sequences of class *i* in cluster *j* *n*_{
i
} is the number of the sequences of class *i*, and *n*_{
j
} is the number of the sequences of cluster *j*.

where *N* is the total number of sequences in the data set. Clearly, an F-measure has a value between 0 and 1. The larger the F-measure is, the better the clustering result is.

### Evaluation of similarity measure

To evaluate the proposed similarity measure, we test DMk on gene sequence data sets and compare it with the *k*-tuple distance. We also verify the effectiveness of DMk by assessing how well it performs on phylogenetic analysis.

#### Gene sequences clustering

Genes of the same family usually share similar sequences, functional domains, and even interacting partners. When a new gene is assigned to a cluster, the biological function of this cluster can be attributed to this gene with high confidence.

**Description for the Data Sets**

Data | Name | Number | Average length (bp) | Description |
---|---|---|---|---|

DS1 | beta-globin | 176 | 1531 | Cytochrome P450 |

beta-Hemoglobin | 89 | 448 | Hemoglobin subunit | |

integrin_alpha | 142 | 3360 | Integrin, alpha | |

ketoacyl-synt1 | 43 | 754 | Estradiol 17-beta-dehydrogenase 8 | |

myoglobin | 55 | 478 | Cytoglobin Myoglobin | |

RWD | 93 | 825 | RWD domain-containing protein | |

VCL | 92 | 2746 | Vinculin | |

Histone | 81 | 668 | Histone | |

DS2 | HBG106679 | 22 | 446 | Copper uptake protein 2 |

HBG108349 | 49 | 718 | Prolactin | |

HBG079775 | 26 | 3152 | Transcription elongation factor SPT5 | |

HBG058842 | 34 | 1351 | TNFR superfamily member 1A | |

HBG002834 | 92 | 951 | Calumenin/Reticulocalbin | |

HBG050441 | 58 | 1899 | ATP-binding cassette sub-family G member | |

DS3 | HBG093787 | 32 | 1769 | Hypothetical membrane proteins |

HBG099893 | 34 | 430 | Putative membrane protein precursor | |

HBG415481 | 65 | 557 | Phasin like/family protein | |

HBG423057 | 32 | 236 | Hypothetical proteins | |

HBG050644 | 99 | 3129 | Beta galactosidase, beta glucuronidase, Evolved beta-D-galactosidase alpha subunit | |

HBG364776 | 48 | 1069 | Formate dehydrogenase gamma subunit precursor | |

DS4 | HBG000080 | 29 | 674 | BWK-1,CG6617-PA , Zgc:73100 C20orf11 homolog , RH01588p |

HBG060165 | 28 | 163 | ATP synthase, H + transporting mitochondrial F1 complex/epsilon subunit | |

HBG010471 | 48 | 1802 | Hypothetical Glycosyl transferase, family 25/Endoplasmic reticulum targeting sequence containing protein | |

HBG000013 | 70 | 318 | 60 S ribosomal protein L36a-like, 60 S ribosomal protein L42, L44, IP15820p, RPL | |

HBG000026 | 18 | 3157 | Eukaryotic translation initiation factor 2-alpha kinase 3 precursor, Eukaryotic translation initiati | |

HBG065748 | 48 | 1238 | AT20832p,AT27361p, CG10513-PA, CG10514-PA, CG10550-PA, isoform A, CG10553-PA,CG10559-PA,CG10560-P |

Four widely used clustering algorithms, including KM, single-linkage clustering (SL), complete-linkage clustering (CL) and average-linkage clustering (AL), have been chosen in the experiments. For comparison, we perform the clustering tests on all data sets using the *k*-tuple distance and DMk distance. In this paper, we set *k* value to 3. For protein coding genes, a tuple size of 3 is a good choice according to reference [22]. We also tested the clustering performance on different *k* values, and the result confirms that a small *k* value is preferred, see Additional file 1: Table S1. For larger *k* values, there are more tuples with zero frequencies and less information is captured by the algorithm.

KM algorithm would yield different results during multiple executions due to its stochastic feature for initialization. We examine KM in ten runs and report the average performance. The AL, CL and SL hierarchical algorithms generate one solution for each of them. We obtain the result of hierarchical clustering algorithms by analyzing the hierarchical tree using the expected number of cluster as input parameters.

*k*-tuple distance. In our experiments, on average, the value of the F-measure given by DMk is 18% better than by the

*k*-tuple distance (p = 0.0165, one-sided paired

*t*-test) in KM, 49.7% better in SL (p = 0.0028), 24.9% better in CL (p = 0.016), and 35.8% better in AL (p = 0.01885). Clearly, DMk provides a significant improvement in clustering sequences. On the four data sets, the F-measure of DMk is improved more than 20% compared with that of the

*k*-tuple distance during the same clustering process in most cases. DMk outperforms the

*k*-tuple distance in the experiments. This is because DMk considers the occurrence, location and order relation of tuples in sequence and can capture more information in the sequence, while the

*k*-tuple distance considers frequency alone and ignore the position of tuples in a sequence. In addition, we have tested DMk and

*k*-tuple measures on protein sequences with a

*k*value of 2, and the results indicate that DMk performs better than

*k*-tuple distance (data not shown). Thus in practical DMk measure can also be applied in clustering protein sequences after tuning current algorithm.

**The F-measures of the Data Sets**

Method | DS1 | DS2 | DS3 | DS4 |
---|---|---|---|---|

KM with | 0.5738 | 0.7828 | 0.5543 | 0.6532 |

SL with | 0.3544 | 0.4148 | 0.3307 | 0.3244 |

CL with | 0.5153 | 0.7253 | 0.5588 | 0.516 |

AL with | 0.5113 | 0.6956 | 0.5578 | 0.3185 |

BKM with | 0.5725 | 0.7876 | 0.5498 | 0.6551 |

mBKM with | 0.5882 | 0.7913 | 0.5691 | 0.6722 |

KM with DMk | 0.7 | 0.8261 | 0.7716 | 0.8284 |

SL with DMk | 0.601 | 0.7948 | 0.8188 | 0.6535 |

CL with DMk | 0.7172 | 0.9295 | 0.6868 | 0.7468 |

AL with DMk | 0.7898 | 0.9365 | 0.6963 | 0.8498 |

BKM with DMk | 0.7346 | 0.8511 | 0.8044 | 0.8813 |

mBKM with DMk | 0.808 | 0.9645 | 0.9143 | 0.9587 |

#### Phylogenetic analysis

In this experiment, the proposed similarity measure DMk is further tested by phylogenetic analysis. In order to evaluate the similarity measures, we use UPGMA in the PHYLIP package, a widely used clustering algorithm in phylogenetic analysis. The tree is drawn by TREEVIEW program [53].

**The similarity/dissimilarity matrix for the 10 full β-globin gene sequences based on DMk**

Species | Human | Goat | Opossum | Gallus | Lemur | Mouse | Rat | Gorilla | Bovine | Chimpanzee |
---|---|---|---|---|---|---|---|---|---|---|

Human | 0 | 22.95 | 37.65 | 111.47 | 14.02 | 35.21 | 20.68 | 3.42 | 25.07 | 3.54 |

Goat | 0 | 41.22 | 65.70 | 18.80 | 35.05 | 33.93 | 32.36 | 6.04 | 33.05 | |

Opossum | 0 | 42.54 | 33.29 | 64.03 | 51.64 | 46.35 | 40.41 | 49.73 | ||

Gallus | 0 | 90.93 | 80.07 | 95.26 | 121.09 | 61.69 | 122.65 | |||

Lemur | 0 | 21.39 | 18.50 | 17.19 | 18.12 | 18.74 | ||||

Mouse | 0 | 16.04 | 33.64 | 27.60 | 37.59 | |||||

Rat | 0 | 17.69 | 30.53 | 20.58 | ||||||

Gorilla | 0 | 33.66 | 0.80 | |||||||

Bovine | 0 | 35.46 | ||||||||

Chimpanzee | 0 |

In Table 3, the most similar species pairs are human-gorilla, human-chimpanzee and gorilla-chimpanzee, which are expected from their evolutionary relationship. A slightly less similar species pair is goat-bovine. On the other hand, gallus is separated from the rest, this coincides with the fact that gallus is the only nonmammalian species among these 10 species. We can also find that opossum is far away from the remaining mammals. These results are consistent with biological morphology.

*k*-tuple distance is shown in Figure 1(a).

The tree in Figure 1 (a) has some consistencies with biological morphology. Although it supports the separation of gallus relative to other species, its obvious drawback is that it fails to separate (mouse, rat) and (goat, bovine) from opossum. From Figure 1 (b), gallus is separated from the rest and opossum is far away from the other species. This topology is in good agreement with that presented by Feng et al. [54] and Cao et al. [55] except for the relative position of rodents.

DMk measures the similarity between DNA sequences more effective than the *k*-tuple distance. This is because DMk measures the distance between DNA sequences based on sequence structure and composition. Through evaluation on gene families and constructing phylogenetic trees of full gene sequences of 10 species, we find that DMk gives more competitive results compared to the *k*-tuple distance.

### Evaluation of clustering methods

To evaluate the effectiveness of the proposed clustering algorithm, mBKM, we apply mBKM in clustering gene sequences and compare it with several clustering algorithms. Moreover, we use our method, mBKM with similarity measure DMk, in phylogenetic analysis to show how well the genes are grouped together and how well the resulting trees agree with existing phylogenies.

#### Performance comparison of clustering methods

In order to illustrate the efficiency of mBKM in gene sequence clustering, we ran mBKM with the *k*-tuple distance and DMk on real data sets listed in Table 1. The clustering results are compared with those of KM, SL, CL, AL and BKM algorithms. For BKM, the number of iterations for each bisecting step is set to 5. We ran BKM 10 times to obtain the average F-measure. By combing the six clustering algorithms with two similarity measures, we have 12 combinations of clustering algorithm for performance assessment. The combinations are KM with *k*-tuple, SL with *k*-tuple, CL with *k*-tuple, AL with *k*-tuple, BKM with *k*-tuple, mBKM with *k*-tuple, KM with DMk, SL with DMk, CL with DMk, AL with DMk, BKM with DMk and mBKM with DMk.

The clustering performance of different clustering methods is the result of a combination of factors, including the types of sequence distances used for clustering and the choice of clustering algorithms. Table 2 shows the clustering performance on the data sets for all 12 clustering methods. For each data set, we set the number of cluster as the real number of class during the clustering run. For example, the real number of cluster is 8 in DS1 and 6 in DS2.

From Table 2, we observe that mBKM using DMk achieves best result and clearly outperforms other methods for the four data sets. The average F-measure of mBKM with *k*-tuple is about 2.2% higher than KM with *k*-tuple (p = 0.036), 45% higher than SL with *k*-tuple (p = 0.00195), 11.4% higher than CL with *k*-tuple (p = 0.0424), 19% higher than AL with *k*-tuple (p = 0.08615) and 2.3% higher than BKM (p = 0.0141). For mBKM with DMk, F-measures for DS1, DS2, DS3, and DS4 are 0.808, 0.9645, 0.9143, and 0.9587 respectively. On average, the value of F-measure given by mBKM is 14.2% better than KM (p = 0.00025), 21.3% better than SL (p = 0.0105), 15.4% better than CL (p = 0.02835), 10.1% better in AL (p = 0.0686), and 2.3% higher than BKM (p = 0.0015) respectively. These results show that our method, combining mBKM with DMk, is able to achieve high quality results on all the data sets.

*k*-tuple distance and DMk, respectively. The data used for generating these figures are included in Additional file 1: Tables S3-S10.

Figure 2 illustrates the results of the six clustering algorithms with the *k*-tuple distance. From Figure 2 and Additional file 1: Tables S3-S6, mBKM achieves better F-measures than other five clustering algorithms for the real number of clusters on all the data sets. Although the other clustering algorithms give slightly better results in terms of F-measure in some cases, mBKM performs better than the other clustering algorithms in terms of the average of the F-measures values (average values are shown in Additional file 1: Tables S3-S6). This result shows that on average, mBKM performs better than other clustering algorithms for a range of cluster numbers, in the vicinity of real number of clusters. It also implies that varying the number of clusters as input for these clustering algorithms does not affect the performance.

Figure 3 shows the results of clustering algorithms with DMk. mBKM obtains the highest F-measure values among the six clustering algorithms at the real number of clusters. On average, mBKM achieves better results than the other clustering algorithms for DS2, DS3, and DS4. For DS1, the average value of mBKM is very close to that of AL and higher than those of the other clustering algorithms. Overall mBKM produces consistently high quality clusters in the neighborhood of the real number of cluster (data shown in Additional file 1: Tables S7-S10). The F-measures given by mBKM are higher than those of other clustering methods at the corresponding number of clusters in most cases.

From Figures 2 and 3, we can see that DMk achieves better cluster quantity than the *k*-tuple distance in terms of F-measure. Using same clustering algorithm on the same data set, DMk achieves higher average of the F-measure values than the *k*-tuple distance, and DMk also obtains higher F-measures at corresponding number of clusters (data shown in Additional file 1: Tables S3-S10). From both Figures, we find that F-measure changes as the number of cluster changes. As it is known, F-measure is a balanced measure of precision and recall. It is an ideal condition when the number of cluster is equal to the real number. When the number of cluster is greater than or less than the real number, the F-measure will be affected.

With regard to clustering algorithms, SL performs poorly in many cases, and this may be because that SL uses the nearest pair of sequences and may lead to bad splits of one cluster if two or more clusters show different pattern densities. For KM and BKM, the results of many runs are lower than those of mBKM. On the whole, mBKM achieves better results than other clustering algorithms, and mBKM combining with DMk achieves best results among these clustering methods in our experiments.

The task of sequence clustering is to group given sequences into clusters. The similarity measure, DMk, measures the similarity between DNA sequences based solely on the *k*-tuple. It is more effective than the *k*-tuple distance, which is one of the most widely used methods. The clustering algorithm, mBKM, can obtain better clustering results and can reveal the relationships among clusters in hierarchical manner. In the next experiments, we combine mBKM with DMk to clustering DNA sequences.

In order to further illustrate the efficiency of our method, combining mBKM and DMk, we compare mBKM with DMk to two other clustering programs: BlastClust [27] and CD-HIT-EST [39]. BlastClust is an alignment-dependent clustering algorithm. BlastClust is from NCBI Blast package. BlastClust accepts a number of parameters that can be used to control the clustering stringency including thresholds for score density (−*S* parameter), and alignment length (−*L* parameter). CD-HIT-EST is a popular DNA clustering program based on greedy incremental clustering method. CD-HIT-EST groups DNA sequences into clusters that meet a user-defined similarity threshold (−*c* parameter) and uses short-word filters to rapidly determine that if two sequences are similar, which reduces the number of full alignments necessary.

*p*as

*F*(input type is nucleotide sequence) and vary the input parameters, -

*S*and –

*L*, to evaluate the results. The score density, –

*S*parameter, varies between 10 and 90 with step size 10, and the alignment length, –

*L*parameter, varies between 0.1 and 0.9 with step size 0.1. Other parameters are kept default. For CD-HIT-EST, because the sequence identity threshold, -

*c*parameter, should be greater than or equal to 0.8 in the program, we vary -

*c*parameter between 0.8 and 1 with step size 0.02, and set the word length as default value. The best results from different parameter combination are recorded. For mBKM with DMk, we set the size of

*k*-tuple as 3 and use the real number of clusters as input. As BlastClust and CD-HIT-EST do not use the number of clusters as input, we choose the resulting class

*i*, which has the max

*F(i,j)*for cluster

*j*, to calculate the F-measures. The results, which contain the corresponding F-measures and the execution time, are summarized in Table 4.

**Clustering results on the data sets listed in Table** 1

mBKM with DMk | BlastClust | CD-HIT-EST | ||||
---|---|---|---|---|---|---|

Data | F-measure | Time(s) | F-measure | Time(s) | F-measure | Time(s) |

DS1 | 0.8080 | 6.875 | 0.4525 | 48 | 0.2713 | 39.8 |

DS2 | 0.9645 | 1.844 | 0.7515 | 13.6 | 0.5924 | 6.4 |

DS3 | 0.9143 | 2.375 | 0.3693 | 12.7 | 0.3157 | 17.1 |

DS4 | 0.9587 | 1.328 | 0.5224 | 9.3 | 0.4007 | 6.8 |

Table 4 demonstrates that mBKM with DMk produces good results relative to each original cluster set in terms of F-measure. Every F-measure of mBKM with DMk is higher than 0.8 and the highest is 0.9645. It is also seen in the table that mBKM with DMk outperforms BlastClust and CD-HIT-EST on all the data sets. BlastClust and CD-HIT-EST tend to give more clusters than the real numbers of classes, therefore, BlastClust and CD-HIT-EST give high precision and low recall value. But neither of these two performs well in terms of F-measure. The execution times reported in Table 4 for algorithm comparison show mBKM with DMk is faster than BlastClust and CD-HIT-EST.

For the cases that the real number of clusters is unknown, the performance of our algorithm will be affected. In order to compare with BlastClust and CD-HIT-EST on a relatively fair ground, we can vary the number of clusters and take the average of the F-measure values over the different numbers of clusters. For instance, we run mBKM with DMk with the range of 3–20 numbers and the average values of F-measure are 0.7065, 0.8533, 0.8205 and 0.8429 for DS1, DS2, DS3 and DS4, respectively. As shown in Additional file 1: Tables S7-S10, these values are also higher than the corresponding F-measure of BlastClust and CD-HIT-EST.

#### Phylogenetic analysis

- 1)
The clustering result of 10 species

*rodents*relative to

*ferungulates*and

*primates*. Moreover, the tree in Figure 4(a) is identical to the tree in Figure 4(b)4(c) and the tree built MUSCLE with UPGMA. In experiment, the branch (bovine, goat) is not classified well by CLUSTALW with ML. Furthermore, it took about 0.1 second for our method. However, UPGMA with CLUSTALW and MUSCLE for the same data set took 5.1 and 1.2 seconds to build the tree, respectively, and ML with CLUSTALW and MUSCLE took 8 and 4.1 seconds to build the tree, respectively.

- 2)
The Clustering result of 60 H1N1 viruses

As is seen from Figure 5(a), 60 H1N1 viruses are distinctly divided into four main groups using our method. The four groups, include European swine older than 2009 (G1), the avian older than 2009 (G2), American swine older than 2009 (G3) and the new 2009 viruses from human, swine and avian (G4). The result shows that the new 2009 human H1N1 viruses have closer relationship with old American swine than old avian and European swine. This grouping result is generally consistent with the topology given by CLUSTALW with UPGMA, which is shown in Figure 5(b), and the one presented by MUSCLE with UPGMA, which is provided in the Additional file 1, as well as the result suggested by zhao et al. [58]. Figure 5(c), built by MUSCLE using ML method, also shows the new 2009 human H1N1 viruses have close relationship with old American swine except the position of the group (old avian swine, European swine) is different from the positions in Figure 5(a) and 5(b). CLUSTALW with ML (in Additional file 1) also classifies the 60 H1N1 viruses into four groups except that swine/Wisconsin/1961 and swine/Wisconsin/1961 are not classified well.

Our method analyzed the 60 H1N1 viruses within 1 second, while UPGMA with CLUSTALW and MUSCLE of the same data set took 460 and 60.1 seconds to build the tree, and ML with CLUSTALW and MUSCLE took 571 and 188.1 seconds to build the tree, respectively.

Our method, mBKM with DMk, performs well when clustering 10 species and 60 H1N1 viruses. It obtains similar results to the alignment-based method. Furthermore, our method is much faster than the alignment-based methods.

#### Scalability test

For DMk, the time complexity of transforming the gene sequence *s*_{1}⋯*s*_{
l
} to a vector is *O* (*l*4^{
K
}), thus the time complexity of generating the vectors for the whole sequence database is $O(N\overline{l}{4}^{k})$, where $\overline{l}$ is the average length of the sequences and *N* is the number of sequences. The value of *k* set to 3 yields good results in our experiments, and we fix *k* to 3 as the size of *k*-tuple. DMk have linear time complexity with respect to both $\overline{l}$ and *N*.

The time consumed for mBKM calculation is primarily determined by choosing the initial cluster centroids. For *N* sequences, this step has a time complexity of *O* (*N*^{2}). The time complexity of clustering step in mBKM is *O* (*N* log *K*). The following scalability test on our method, mBKM with DMk, confirms that our method has linear time complexity with respect to the average length of the sequences. The scalability test uses theoretical model sequences composed of the four symbols ‘A’, ‘C’, ’G’ and ‘T’. The method is implemented in Java and on a computer with 3.00 GHz CPU and 2 GB RAM.

*O*(

*N*

^{2}). The scalability with respect to the length of sequences was tested on five datasets with five different sequence lengths: 10000, 20000, 30000, 40000, 50000 and each set consists of 4 clusters and 100 sequences. The sensitivity with respect to the length of the sequence is illustrated in Figure 7(b), from which we can see that the time of our method increases linearly when the length of sequences increases.

## Conclusions

In this paper, we presented a novel approach for DNA sequence clustering, mBKM, based on a new sequence similarity measure, DMk, which is extracted from DNA sequences based on the position and composition of oligonucleotide pattern. The experimental results show the method of combining mBKM with DMk is effective in classifying DNA sequences with similar biological characteristics and in discovering the underlying relationship among the sequences. In addition, DMk can achieve comparable or better accuracy than the frequency-based distance measure. Our proposed method can be applied to study gene families and it can also help with the prediction of novel genes. Furthermore, mBKM with DMk can generate cluster trees that are useful to understand the processes governing the gene evolution. In addition, our method may be extended for protein sequence analysis and metagenomics of identifying source organisms of metagenmic data. Our method has limitations too. For example, the method did not consider edge length, and has not address problems with long repeated sequences or long insertions. In future we will try to address these problems.

## Declarations

### Acknowledgements

This work is supported by the National Natural Science Foundation of China under Grant No.61175123 and No.10771176, and the Shenzhen New Industry Development Fund under grant No.CXB201005250021A.

## Authors’ Affiliations

## References

- Demuth JP, De Bie T, Stajich JE, Cristianini N, Hahn MW: The evolution of mammalian gene families. PLoS One. 2006, 1: 1-10. 10.1371/journal.pone.0000001.View ArticleGoogle Scholar
- Zhao B, Duan V, Yau SS: A novel clustering method via nucleotide-based Fourier power spectrum analysis. JTheor Biol. 2011, 279: 83-89. 10.1016/j.jtbi.2011.03.029.View ArticleGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: A basic local alignment search tool. JMol Biol. 1990, 215: 403-410.View ArticleGoogle Scholar
- Pearson WR, Lipman DJ: Improved tools for biological sequence comparison. ProcNatlAcad Sci USA. 1988, 85 (8): 2444-2488. 10.1073/pnas.85.8.2444.View ArticleGoogle Scholar
- Vinga S, Almeida J: Alignment-free sequence comparison-a review. Bioinformatics. 2003, 19 (4): 513-523. 10.1093/bioinformatics/btg005.View ArticlePubMedGoogle Scholar
- Haubold B, Reed FA, Pfaffelhuber P: Alignment-free estimation of nucleotide diversity. Bioinformatics. 2011, 27 (4): 449-455. 10.1093/bioinformatics/btq689.View ArticlePubMedGoogle Scholar
- Liu Z, Meng J, Sun X: A novel feature-based method for whole genome phylogenetic analysis without alignment: application to HEV genotyping and subtyping. Biochem Biophys Res Commun. 2008, 368 (2): 223-230. 10.1016/j.bbrc.2008.01.070.View ArticlePubMedGoogle Scholar
- Domazet-Loso M, Haubold B: Efficient estimation of pairwise distances between genomes. Bioinformatics. 2009, 25 (24): 3221-3227. 10.1093/bioinformatics/btp590.View ArticlePubMedGoogle Scholar
- Domazet-Loso M, Haubold B: Alignment-free detection of local similarity among viral and bacterial genomes. Bioinformatics. 2011, 27 (11): 1466-1472. 10.1093/bioinformatics/btr176.View ArticlePubMedGoogle Scholar
- Kelil A, Wang S, Brzezinski R, Fleury A: CLUSS: Clustering of protein sequences based on a new similarity measure. BMC Bioinformatics. 2007, 8: 286-10.1186/1471-2105-8-286.PubMed CentralView ArticlePubMedGoogle Scholar
- Reinert G, Chew D, Sun FZ, Waterman MS: Alignment-free sequence comparison (I): statistics and power. JComput Biol. 2009, 16 (12): 1615-1634.View ArticleGoogle Scholar
- Dai Q, Liu X, Yao Y, Zhao F: Numerical characteristics of word frequencies and their application to dissimilarity measure for sequence comparison. JTheor Biol. 2011, 276 (1): 174-180. 10.1016/j.jtbi.2011.02.005.View ArticleGoogle Scholar
- Lu G, Zhang S, Fang X: An improved string composition method for sequence comparison. BMC Bioinformatics. 2008, 9 (Suppl 6): S15-10.1186/1471-2105-9-S6-S15.PubMed CentralView ArticlePubMedGoogle Scholar
- Aita T, Husimi Y, Nishigaki K: A mathematical consideration of the word-composition vector method in comparison of biological sequences. BioSystems. 2011, 106: 67-75. 10.1016/j.biosystems.2011.06.009.View ArticlePubMedGoogle Scholar
- Blaisdell BE: A measure of the similarity of sets of sequences not requiring sequence alignment. ProcNatlAcad Sci USA. 1986, 83: 5155-5159. 10.1073/pnas.83.14.5155.View ArticleGoogle Scholar
- Wu TJ, Burke JP, Davison DB: A measure of DNA sequence dissimilarity based on Mahalanobis distance between frequencies of words. Biometrics. 1997, 53 (4): 1431-1439. 10.2307/2533509.View ArticlePubMedGoogle Scholar
- Wu TJ, Hsieh YC, Li LA: Statistical measures of DNA dissimilarity under Markov chain models of base composition. Biometrics. 2001, 57 (2): 441-448. 10.1111/j.0006-341X.2001.00441.x.View ArticlePubMedGoogle Scholar
- Stuart GW, Moffett K, Baker S: Integrated gene and species phylogenies from unaligned whole genome protein sequences. Bioinformatics. 2002, 18 (1): 100-108. 10.1093/bioinformatics/18.1.100.View ArticlePubMedGoogle Scholar
- Fichant G, Gautier C: Statistical method for predicting protein coding regions in nucleic acid sequences. ComputAppl Biosci. 1987, 3 (4): 287-295.Google Scholar
- Wang J, Zheng X: WSE, a new sequence distance measure based on word frequencies. Math Biosci. 2008, 215 (1): 78-83. 10.1016/j.mbs.2008.06.001.View ArticlePubMedGoogle Scholar
- Zheng X, Qin Y, Wang J: A Poisson model of sequence comparison and its application to coronavirus phylogeny. Math Biosci. 2009, 217 (2): 159-166. 10.1016/j.mbs.2008.11.006.View ArticlePubMedGoogle Scholar
- Yang K, Zhang L: Performance comparison of gene family clustering methods with expect curated gene family data set in Arabidposis thaliana. Planta. 2008, 228: 439-447. 10.1007/s00425-008-0748-7.View ArticlePubMedGoogle Scholar
- Dong G, Pei J: Classification, clustering, features and distances of sequence Data. Sequence Data Mining. 2007, 33: 47-65. 10.1007/978-0-387-69937-0_3.View ArticleGoogle Scholar
- Sokal RR, Rohlf FJ: Biometry: The Principles and Practice of Statistics in Biological Research. 1995, W. H. Freeman and Company, New York, 3Google Scholar
- Everitt BS, Landau S, Leese M: Cluster Analysis. 2001, Oxford University Press, LondonGoogle Scholar
- Loewenstein Y, Portugaly E, Fromer M, Linial M: Efficient algorithms for accurate hierarchical clustering of huge datasets: tackling the entire protein space. Bioinformatics. 2008, 24 (13): i41-i49. 10.1093/bioinformatics/btn174.PubMed CentralView ArticlePubMedGoogle Scholar
- National Center for Biotechnology Information (NCBI): Documentation of the BLASTCLUST-algorithm. ftp://ftp.ncbi.nih.gov/blast/documents/blastclust.html,
- Enright AJ, Ouzounis CA: GeneRAGE: a robust algorithm for sequence clustering and domain detection. Bioinformatics. 2000, 16 (5): 451-457. 10.1093/bioinformatics/16.5.451.View ArticlePubMedGoogle Scholar
- Chaudhuri P, Das S: SWORDS: A statistical tool for analyzing large DNA sequences. J Biosci. 2002, 27 (1): 1-6. 10.1007/BF02703678.View ArticlePubMedGoogle Scholar
- Uchiyama I: Hierarchical clustering algorithm for comprehensive orthologous-domain classification in multiple genomes. Nucleic Acids Res. 2006, 34 (2): 647-658. 10.1093/nar/gkj448.PubMed CentralView ArticlePubMedGoogle Scholar
- Hao X, Jiang R, Chen T: Clustering 16 S rRNA for OTU prediction: a method of unsupervised Bayesian clustering. Bioinformatics. 2011, 27 (5): 611-618. 10.1093/bioinformatics/btq725.PubMed CentralView ArticlePubMedGoogle Scholar
- Sperisen P, Pagni M: JACOP: a simple and robust method for the automated classification of protein sequences with modular architecture. BMC Bioinformatics. 2005, 6: 216-10.1186/1471-2105-6-216.PubMed CentralView ArticlePubMedGoogle Scholar
- Wan X, Bridges SM, Boyle JA, Boyle AP: Interactive Clustering for Exploration of Genomic Data. SmartEng Design. 2002, 12: 753-758.Google Scholar
- Kelarev A, Kang B, Steane D: Clustering Algorithms for ITS Sequence Data with Alignment Metrics. Lect Notes ComputSci. 2006, 4304: 1027-1031. 10.1007/11941439_116.View ArticleGoogle Scholar
- Tseng GC: Penalized and weighted K-means for clustering with scattered objects and prior information in high-throughput biological data. Bioinformatics. 2007, 23 (17): 2247-2255. 10.1093/bioinformatics/btm320.View ArticlePubMedGoogle Scholar
- Ashlock D, Warner E: Classifying Synthetic and Biological DNA Sequences with Side Effect Machines. IEEE Symposium on Computational Intelligence in Bioinformatics and Computational Biology: 15–17 Sept. 2008; Sun Valley, ID. 2008, 22-29.View ArticleGoogle Scholar
- Zhao Y, Karypis G: Criterion functions for document clustering: experiments and analysis. 2001, Technical ReportGoogle Scholar
- Kashef R, Kamel MS: Enhanced bisecting k-means clustering using intermediate cooperation. Pattern Recognit. 2009, 42 (11): 2557-2569. 10.1016/j.patcog.2009.03.011.View ArticleGoogle Scholar
- Li W, Godzik A: Cd-hit: a Fast Program for Clustering and Comparing Large Sets of Protein or Nucleotide Sequences. Bioinformatics. 2006, 22 (13): 1658-1659. 10.1093/bioinformatics/btl158.View ArticlePubMedGoogle Scholar
- Alam I, Cornell M, Soanes DM, Hedeler C, Wong HM, Rattray M, Hubbard SJ, Talbot NJ, Oliver SG, Paton NW: A Methodology for Comparative Functional Genomics. JIntegr Bioinform. 2007, 4 (3): 69-Google Scholar
- Picardi E, Mignone F, Pesole G: EasyCluster: a fast and efficient gene-oriented clustering tool for large-scale transcriptome data. BMC Bioinformatics. 2009, 10 (Suppl 6): S10-10.1186/1471-2105-10-S6-S10.PubMed CentralView ArticlePubMedGoogle Scholar
- Wei D, Jiang Q: A DNA Sequence Distance Measure Approach for Phylogenetic Tree Construction. 5th IEEE International Conference on Bio-Inspired Computing: Theories and Applications: 23–26 Sept. 2010 Changsha. Edited by: Li K, Tang Z, Li R, Nagar AK, Thamburaj R. 2010, 204-212.Google Scholar
- Neumann-Held EM: The gene is dead-Long live the gene: Conceptualizing genes the constructionist way. Sociobiology and Bioeconomics: the Theory of Evolution in Biological and Economic Theory. Edited by: Koslowski P. 1999, Springer, Berlin, 105-137.View ArticleGoogle Scholar
- White JR, Navlakha S, Nagarajan N, Ghodsi M, Kingsford C, Pop M: Alignment and clustering of phylogenetic markers- implications for microbial diversity studies. BMC bioinformatics. 2010, 11: 152-10.1186/1471-2105-11-152.PubMed CentralView ArticlePubMedGoogle Scholar
- Waterman MS: Introduction to Computational Biology: Maps, Sequences, and Genomes. 1995, Chapman and Hall, LodonView ArticleGoogle Scholar
- Durbin R, Eddy SR, Krogh A, Mitchison G: Biological Sequence Analysis: probabilistic models of proteins and nucleic acids. 1998, Cambridge University Press, CambridgeView ArticleGoogle Scholar
- Nair ASS, Mahalakshmi T: Visualization of genomic data using inter-nucleotide distance signals. Proceedings of IEEE Genomic Signal Processing: 11–13 July 2005; Bucharest. 2005Google Scholar
- Schmitt AO, Herzel H: Estimating the Entropy of DNA sequences. JTheor Biol. 1997, 188 (3): 369-377. 10.1006/jtbi.1997.0493.View ArticleGoogle Scholar
- Li C, Wang J: Relative entropy of DNA and its application. Physica A: Stat Mech Appl. 2005, 347: 465-471.View ArticleGoogle Scholar
- Steinbach M, Karypis G, Kumar V: A comparison of document clustering techniques. KDD Workshop on Text Mining: 20–23. 2000, August ; BostonGoogle Scholar
- Zhao Y, Karypis G: Hierarchical Clustering Algorithms for Document Datasets. Data Mining Knowl Discov. 2005, 10: 141-168. 10.1007/s10618-005-0361-3.View ArticleGoogle Scholar
- Larsen B, Aone C: Fast and effective text mining using linear-time document clustering. Proceedings of the fifth ACM SIGKDD international conference on Knowledge discovery and data mining: 15–18 August 1999; San Diego. Edited by: Fayyad U, Chaudhuri S, Madigan D. 1999, 16-22.View ArticleGoogle Scholar
- Page RD: TreeView: an application to display phylogenetic trees on personal computers. Bioinformatics. 1996, 12: 357-358. 10.1093/bioinformatics/12.4.357.View ArticleGoogle Scholar
- Feng J, Hu Y, Wan P, Zhang A, Zhao W: New method for comparing DNA primary sequences based on a discrimination measure. JTheor Biol. 2010, 266 (4): 703-707. 10.1016/j.jtbi.2010.07.040.View ArticleGoogle Scholar
- Cao Y, Janke A, Waddell PJ, Westerman M, Takenaka O, Murata S, Okada N, Paabo S, Hasegawa M: Conflict among individual mitochondrial proteins in resolving the phylogeny of eutherian orders. JMol Evol. 1998, 47 (3): 307-322.View ArticleGoogle Scholar
- Chenna R, Sugawara H, Koike T, Lopez R, Gibson TJ, Higgins DG, Thompson JD: Multiple sequence alignment with the Clustal series of programs. Nucleic Acids Res. 2003, 31 (13): 3497-3500. 10.1093/nar/gkg500.PubMed CentralView ArticlePubMedGoogle Scholar
- Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32 (5): 1792-1797. 10.1093/nar/gkh340.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhao B, He RL, Yau SS: A new distribution vector and its application in genome clustering. MolPhylogenet Evol. 2011, 59 (2): 438-443.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.