- Open Access
Protein sequence classification using feature hashing
© Caragea et al; licensee BioMed Central Ltd. 2012
- Published: 21 June 2012
Recent advances in next-generation sequencing technologies have resulted in an exponential increase in the rate at which protein sequence data are being acquired. The k-gram feature representation, commonly used for protein sequence classification, usually results in prohibitively high dimensional input spaces, for large values of k. Applying data mining algorithms to these input spaces may be intractable due to the large number of dimensions. Hence, using dimensionality reduction techniques can be crucial for the performance and the complexity of the learning algorithms. In this paper, we study the applicability of feature hashing to protein sequence classification, where the original high-dimensional space is "reduced" by hashing the features into a low-dimensional space, using a hash function, i.e., by mapping features into hash keys, where multiple features can be mapped (at random) to the same hash key, and "aggregating" their counts. We compare feature hashing with the "bag of k-grams" approach. Our results show that feature hashing is an effective approach to reducing dimensionality on protein sequence classification tasks.
- Hash Function
- Support Vector Machine Classifier
- Latent Dirichlet Allocation
- Average Mutual Information
- Latent Semantic Indexing
Many problems in computational biology, e.g., protein function prediction, subcellular localization prediction, etc., can be formulated as sequence classification tasks , where the amino acid sequence of a protein is used to classify the protein in functional and localization classes.
Protein sequence data contain intrinsic dependencies between their constituent elements. Given a protein sequence x = x 0, ⋯, x n-1over the amino acid alphabet, the dependencies between neighboring elements can be modeled by generating all the contiguous (potentially overlapping) sub-sequences of a certain length k, x i-k , ⋯, x i-1, i = k, ⋯, n, called k-grams, or sequence motifs. Because the protein sequence motifs may have variable lengths, generating the k-grams can be done by sliding a window of length k over the sequence x, for various values of k. Exploiting dependencies in the data increases the richness of the representation. However, the fixed or variable length k-gram representations, used for protein sequence classification, usually result in prohibitively high-dimensional input spaces, for large values of k. Applying data mining algorithms to these input spaces may be intractable due to the large number of dimensions. Hence, using dimensionality reduction techniques can be crucial for the performance and the complexity of the learning algorithms.
In this paper, we study the applicability of feature hashing to protein sequence classification and address the following main questions: (i) How effective is feature hashing on prohibitively high dimensional k-gram representations?; (ii) What is the influence of the hash size (i.e., the reduced dimension) on the performance of protein sequence classifiers that use hash features, and what is the hash size at which the performance starts degrading, due to hash collisions?; and (iii) How does the performance of feature hashing compare to that of the "bag of k-grams" approach? The results of our experiments on three protein subcellular localization data sets show that feature hashing is effective at reducing dimensionality on protein sequence classification tasks.
The paper is organized as follows. In Section 2, we discuss the related work. We provide background on feature hashing in Section 3. Section 4 presents experiments and results, and Section 5 concludes the paper.
Feature selection [5, 7, 12] is a dimensionality reduction technique, which attempts to remove redundant or irrelevant features in order to improve classification performance of learning algorithms. Feature selection methods have been widely used in Bioinformatics for tasks such as protein function prediction and gene prediction, where the features could be k-grams; microarray analysis; mass spectra analysis; single nucleotide polymorphisms (SNPs) analysis, among others (see  for a review).
Topic models, such as Latent Dirichlet Allocation (LDA) , Probabilistic Latent Semantic Analysis (PLSA) , and Latent Semantic Indexing (LSI)  are dimensionality reduction models, designed to uncover hidden topics, i.e., clusters of semantically related words that co-occur in text documents. LSI uses singular value decomposition to identify topics, which are then used to represent documents in a low dimensional "topic" space. LDA models each document as a mixture of topics (drawn from a conjugate Dirichlet prior), and each topic as a distribution over the words in the vocabulary. LDA has recently emerged as an important tool for modeling protein data. For example, Airoldi et al.  proposed the mixed membership stochastic block models to learn hidden protein interaction patterns. Pan et al.  used LDA to discover latent topic features, which correspond to hidden structures in the protein data, and input these features to random forest classifiers to predict protein interactions. However, topic models are computationally expensive, for example, LDA requires inference at runtime to estimate the topic distribution.
Feature abstraction methods  are designed to reduce a model input size by grouping "similar" features into clusters of features. Specifically, feature abstraction learns an abstraction hierarchy over the set of features using hierarchical agglomerative clustering, based on the Jensen-Shannon divergence. A cut or level of abstraction through the resulting abstraction hierarchy specifies a compressed model, where the nodes (or abstractions) on the cut are used as "features" in a classification model. Silvescu et al.  used feature abstraction to simplify the data representation provided to a learner on biological sequence classification tasks.
Shi et al.  and Weinberger et al.  presented hash kernels to map the high-dimensional input spaces into low-dimensional spaces for large scale classification and large scale multitask learning (i.e., personalized spam filtering for hundreds of thousands of users), respectively. Ganchev and Dredze  empirically showed that hash features can produce accurate results on various NLP applications. Forman and Kirshenbaum  proposed a fast feature extraction approach by combining parsing and hashing for text classification and indexing. Hashing techniques have been also used in Bioinformatics. For example, Wesselink et al.  applied hashing to find the shortest contiguous subsequence that uniquely identifies a DNA sequence from a collection of DNA sequences. Buhler and Tompa  applied Locality-Sensitive Hashing (LSH) , a random hashing/projection technique, to discover transcriptional regulatory motifs in eukaryotes and ribosome binding sites in prokaryotes. Furthermore, Buhler  applied LSH to find short ungapped local alignments on a genome-wide scale. Shi et al.  used hashing to compare all subgraph pairs on biological graphs.
In the context of protein sequence classification, it is worth mentioning the fixed and variable-order Markov models (MMs), which capture dependencies between neighboring sequence elements. MMs are among the most widely used generative models of sequence data . In a kth order MM, the sequence elements satisfy the Markov property: each element is independent of the rest given the k preceding elements. One main disadvantage of MMs in practice is that the number of parameters increases exponentially with the range k of direct dependencies, thereby increasing the risk of overfitting. Begleiter et al.  (and papers cited therein) have examined methods for prediction using variable order MMs, including probabilistic suffix trees, which can be viewed as variants of abstraction wherein the abstractions are constrained to share suffixes.
In contrast to the approaches above, we used feature hashing, a very inexpensive approach, to reduce dimensionality on protein sequence classification tasks, and compared it with the "bag of k-grams" approach.
The traditional k-gram approaches construct a vocabulary of size d, which contains all k-grams in a protein data set. A protein sequence is represented as a vector x with as many entries as the number of k-grams in the vocabulary. For a protein sequence, an entry i in x can record the frequency of k-gram i in the sequence, denoted by x i . Because only a small number of k-grams (compared to the vocabulary size) occur in a particular sequence, the representation of x is very sparse, i.e., only a small number of entries of x are non-zero. However, storing the parameter vectors in the original input space requires O(d) numbers, which
Algorithm 1 Feature Hashing
Input: Protein sequence x; hash functions h and ξ, .
Output: Hashed feature vector x h .
x h : = [0, ⋯, 0];
for all k-gram ∈ x do
i = h (k-gram) % b; //Places k-grams into hash bins, from 0 to b-1.
; //Updates the i th hash feature value.
return x h //Records values of hash features.
may become difficult given today's very large collections of protein and DNA sequence data. Feature hashing eliminates the need for such a requirement by implicitly encoding the mapping into a hash function. Next, we briefly overview feature hashing.
Feature hashing [8–11] is a dimensionality reduction technique, in which high-dimensional input vectors x of size d are hashed into low-dimensional feature vectors x h of size b. The procedure for hashing a protein sequence x is shown in Algorithm 1 and is briefly described next (see also Figure 1). Let denote the set of all possible strings (or k-grams) and h and ξ be two hash functions, such that and , respectively. For a protein sequence x, each k-gram in x is directly mapped, using h, into a hash key, which represents the index of the k-gram in the feature vector x h , such that the hash key is a number between 0 and b - 1. Note that h can be any hash function, e.g. hashCode() of the Java String class, or murmurHash function available online at http://sites.google.com/site/murmurhash/. Each index in x h stores the value ("frequency counts") of the corresponding hash feature. The hash function ξ indicates whether to increment or decrement the hash dimension of the k-gram, which renders the hash feature vector x h to be unbiased (see  for more details).
However, for many practical applications, the value of b can be smaller than the theoretical lower bound. This may be problematic as the smaller the size of the hash vector x h becomes, the more collisions occur in the data. Even a single collision of very high frequency words with different class distributions, can result in significant loss of information. Next, we empirically study the applicability of feature hashing on a protein subcellular localization prediction task.
Experiments and results
We used three protein subcellular localization data sets in our study: psortNeg introduced in  and available online at http://www.psort.org/dataset/datasetv2.html, and plant, and non-plant introduced in  and available online at http://www.cbs.dtu.dk/services/TargetP/datasets/datasets.php. The psortNeg data set is extracted from PSORTdb v.2.0 Gram-negative sequences, which contains experimentally verified localization sites. Our data set consists of all proteins that belong to exactly one of the following five classes: cytoplasm (278), cytoplasmic membrane (309), periplasm (276), outer membrane (391) and extracellular (190). The total number of proteins in this data set is 1444. The plant data set contains 940 proteins belonging to one of the following four classes: chloroplast (141), mitochondrial (368), secretory pathway/signal peptide (269) and other (consisting of 54 proteins with label nuclear and 108 examples with label cytosolic). The non-plant data set contains 2738 proteins, each in one of the following three classes: mitochondrial (361), secretory pathway/signal peptide (715) and other (consisting of 1224 proteins labeled nuclear and 438 proteins labeled cytosolic).
Our experiments are designed to explore the following questions: (i) How effective is feature hashing on prohibitively high-dimensional k-gram representations?; (ii) What is the influence of the hash size on the performance of biological sequence classifiers that use hash features, and what is the hash size at which the performance starts degrading, due to hash collisions?; and (iii) How does the performance of feature hashing compare to that of the "bag of k-grams" approach?
To answer these questions, we proceeded with the following steps. We first preprocessed the data by generating all the k-grams from each collection of sequences, i.e., generating all the contiguous (potentially overlapping) sub-sequences of length k, for various values of k. This was done by sliding a window of length k over sequences in each data set. Note that if a k-gram does not appear in the data, it was not considered as a feature. The number of unique k-grams is exponential in k. However, for large values of k, many of the k-grams may not appear in the data (and, consequently, their frequency counts are zero).
Given a protein sequence x, we applied feature hashing in two settings as follows: (i) We first generated all the k-grams of a fixed length k, where k = 3. Each such k-gram was then hashed into a hash key. We refer to this setting as the fixed-length k-grams; (ii) We then generated all the k-grams of various lengths k, for values of k = 1, 2, 3, and 4. Thus, this setting uses the union of k-grams, for values of k ranging from 1 to 4. Each such k-gram (i.e., unigram, 2-gram, 3-gram, or 4-gram) was hashed into a hash key. We refer to this setting as the variable-length k-grams.
We trained Support Vector Machine (SVM) classifiers  on hash features, in both settings, fixed-length and variable-length k-grams, and investigated the influence of the hash size on the performance of the classifiers. Specifically, we trained SVM classifiers for values of the hash size (i.e., the reduced dimension) ranging from 210 to 222, in steps of 1 for the powers of 2, and compared their performance.
Furthermore, we applied feature hashing to sparse high-dimensional variable-length k-gram representations to reduce the dimensionality to a mid-size b-dimensional space, e.g., b = 216 or b = 214, and compared the performance of SVM classifiers trained using hash features with that of SVM classifiers trained using "bag of k-grams".
Specifically, the feature representations used in each case are the following:
a bag of d variable-length k-grams (where all the variable-length k-grams are used as features). This experiment is denoted by baseline.
a bag of b hash features obtained using feature hashing over all d variable-length k-grams, i.e., for each k-gram, feature hashing produces an index i such that i = h(k-gram) % b, where h represents a hash function. This experiment is denoted by FH.
In our experiments, we used the LibLinear implementation of SVM, available at http://www.csie.ntu.edu.tw/~cjlin/liblinear/. As for the hash function, we experimented with both the hashCode of the Java String class, and murmurHash. We found that the results were not significantly different from one another in terms of the number of hash collisions and classification accuracy. We also experimented with both and ξ ≡ 1 - actual counts, and found that the results were not significantly different. The results shown in the next subsection use the hashCode function and ξ ≡ 1. On all three data sets, we reported the average classification accuracy obtained in a 5-fold cross validation experiment. The results are statistically significant (p < 0.05). The classification accuracy is shown as a function of the number of features. The x axis of all figures in the next subsection shows the number of features on a log2 scale (i.e., number of bits in the hash-table).
Comparison of fixed length with variable length k-gram representations
Comparison of fixed-length with variable-length k-gram representations.
Bag of fixed or variable length k-grams
The number of variable length k-grams, for k ranging from 1 to 4, is 155,017. Feature hashing eliminates the need for storing the vocabularies in memory by implicitly encoding the mapping from strings to integers into a hash function. We conclude that feature hashing is very effective on prohibitively high-dimensional k-gram representations, which would otherwise be impractical to use. Because (1-4)-gram representation results in the highest performance, we used it for subsequent experiments.
The influence of hash sizes on classifiers' performance and the comparison of feature hashing with baseline (i.e., the "bag of k-grams" approach)
The number of variable-length k-grams and the rate of hash collisions for various hash sizes.
Value of b
As the hash size increases beyond 216, the performance of SVM classifiers does not change substantially, and, eventually, converges. For example, on the non-plant data set, with 216 hash size, SVM achieves 81.3% accuracy, whereas with 222 hash size, SVM achieves an accuracy of 82.83% (Figure 4a). On the plant data set, SVMs achieve 78.4% and 78.51% accuracy, with 216 and 222 hash sizes, respectively (Figure 4b). Furthermore, as the hash size increases beyond 216, the percentage of hash collisions decreases until no collisions occur (Table 2). For all three data sets, with 222 hash size, there are no hash collisions. The performance of SVMs trained on hash features in the 222 dimensional space is matched by that of SVMs trained on hash features in the 218 dimensional space, suggesting that the hash collisions beyond 218 does not significantly distort the data.
Because 222 (= 4,194,304) highly exceeds the number of unique features, and the rate of hash collisions becomes zero, this can be regarded as equivalent to the classifiers trained without hashing, which require storing the vocabularies in memory, referred as baseline (or the "bag of k-grams") (Figure 4). Moreover, we considered 216 as the point where the performance starts degrading. Note that the vocabulary sizes, i.e., the number of unique variable length k-grams, for non-plant, plant, and psortNeg, are 155017, 111544, and 124389, respectively.
We conclude that, if feature hashing is used to reduce dimensionality from very large dimensions, e.g., 222 to mid-size dimensions, e.g., 216, the hash collisions do not substantially hurt the classification accuracy, whereas if it is used to reduce dimensionality from mid-size dimensions to smaller dimensions, e.g., 210, the hash collision significantly distort the data, and the corresponding SVMs result in poor performance. Also, feature hashing makes it possible to train SVMs that use substantially smaller number of dimensions compared to the baseline, for a small or no drop in accuracy, for example, for a hash size of 216 = 65536 (compared to 155017 variable-length k-grams on the non-plant data set).
We presented an application of feature hashing to reduce dimensionality of very high-dimensional feature vectors to mid-size feature vectors on protein sequence data and compared it with the "bag of k-grams" approach.
The results of our experiments on three protein subcellular localization data sets show that feature hashing is an effective approach to dealing with prohibitively high-dimensional variable length k-gram representations. Feature hashing makes it possible to train SVM classifiers that use substantially smaller number of features compared to the approach which requires storing the vocabularies in memory, i.e., the "bag of k-grams" approach, while resulting in a small or no decrease in classification performance.
Because recent advances in sequencing technologies have resulted in an exponential increase in the rate at which DNA and protein sequence data are being acquired , the application of feature hashing on biological sequence data advances the current state of the art in terms of algorithms that can efficiently process high-dimensional data into low-dimensional feature vectors at runtime.
In the future, it would be interesting to investigate how the performance of hash kernels compares to that of histogram-based motif kernels for protein sequences, introduced by Ong and Zien , and the mismatch string kernels for SVM protein classification introduced by Lesli et al. . Along the lines of dimensionality reduction, it would be interesting to compare the performance of feature hashing with that of feature abstraction  on protein sequence classification tasks. Furthermore, another direction is to apply feature hashing to other types of biological sequence data, e.g., DNA data, and other tasks, e.g., protein function prediction.
This article has been published as part of Proteome Science Volume 10 Supplement 1, 2012: Selected articles from the IEEE International Conference on Bioinformatics and Biomedicine 2011: Proteome Science. The full contents of the supplement are available online at http://www.proteomesci.com/supplements/10/S1.
This research was funded in part by an NSF grant #0845487 to Prasenjit Mitra.
- Baldi P, Brunak S: Bioinformatics: the Machine Learning Approach. MIT Press; 2001.Google Scholar
- Jolliffe IT: Principal Component Analysis. Springer-Verlag, New York; 1986.View ArticleGoogle Scholar
- Blei D, Ng A, Jordan M: Latent Dirichlet Allocation. Journal of Machine Learning Research 2003, 3: 993–1022.Google Scholar
- Hofmann T: Probabilistic Latent Semantic Analysis. Proc of UAI'99, Stockholm 1999.Google Scholar
- Guyon I, Elisseeff A: An introduction to variable and feature selection. J Mach Learn Res 2003, 3: 1157–1182.Google Scholar
- Fleuret F: Fast Binary Feature Selection with Conditional Mutual Information. J Mach Learn Res 2004, 5: 1531–1555.Google Scholar
- Yang Y, Pederson JO: Feature selection in statistical learning of text categorization. In Proceedings of ICML 1997, 412–420.Google Scholar
- Shi Q, Petterson J, Dror G, Langford J, Smola A, Vishwanathan S: Hash Kernels for Structured Data. J Mach Learn Res 2009, 10: 2615–2637.Google Scholar
- Weinberger K, Dasgupta A, Attenberg J, Langford J, Smola A: Feature Hashing for Large Scale Multitask Learning. Proc of ICML 2009.Google Scholar
- Forman G, Kirshenbaum E: Extremely fast text feature extraction for classification and indexing. Proc of the 17th ACM CIKM, New York, NY, USA 2008, 1221–1230.Google Scholar
- Langford J, Li L, Strehl A: Vowpal Wabbit online learning project. 2007.Google Scholar
- Kira K, Rendell LA: The feature selection problem: Traditional methods and a new algorithm. Proc of AAAI, San Jose, CA 1992, 122–134.Google Scholar
- Saeys Y, Inza IN, Larrañaga P: A review of feature selection techniques in bioinformatics. Bioinformatics 2007, 23: 2507–2517. 10.1093/bioinformatics/btm344PubMedView ArticleGoogle Scholar
- Deerwester S, Dumais ST, Furnas GW, Landauer TK, Harshman R: Indexing by latent semantic analysis. Journal of the American Society for Information Science 1990,41(6):391–407. 10.1002/(SICI)1097-4571(199009)41:6<391::AID-ASI1>3.0.CO;2-9View ArticleGoogle Scholar
- Airoldi E, Blei D, Fienberg S, Xing E: Mixed membership stochastic block models for relational data with application to protein-protein interactions. In Proc of the Intl Biometrics Society-ENAR Annual Meeting 2006.Google Scholar
- Pan XY, Zhang YN, Shen HB: Large-Scale Prediction of Human Protein-Protein Interactions from Amino Acid Sequence Based on Latent Topic Features. Journal of Proteome Research 2010,9(10):4992–5001. 10.1021/pr100618tPubMedView ArticleGoogle Scholar
- Silvescu A, Caragea C, Honavar V: Combining Super-Structuring and Abstraction on Sequence Classification. ICDM 2009, 986–991.Google Scholar
- Ganchev K, Dredze M: Small Statistical Models by Random Feature Mixing. Proceedings of the ACL-2008 Workshop on Mobile Language Processing, Association for Computational Linguistics 2008.Google Scholar
- Wesselink JJ, delaIglesia B, James SA, Dicks JL, Roberts IN, Rayward-Smith V: Determining a unique defining DNA sequence for yeast species using hashing techniques. Bioinformatics 2002,18(2):1004–10.PubMedView ArticleGoogle Scholar
- Buhler J, Tompa M: Finding motifs using random projections. Proceedings of the Fifth Annual International Conference on Computational Biology, RECOMB '01 2001, 69–76.View ArticleGoogle Scholar
- Indyk P, Motwani R: Approximate nearest neighbors: towards removing the curse of dimensionality. Proceedings of the Thirtieth Annual ACM Symposium on Theory of Computing, STOC '98 1998, 604–613.View ArticleGoogle Scholar
- Buhler J: Efficient Large-Scale Sequence Comparison by Locality-Sensitive Hashing. Bioinformatics 2001,17(5):419–428. 10.1093/bioinformatics/17.5.419PubMedView ArticleGoogle Scholar
- Durbin R, Eddy SR, Krogh A, Mitchison G: Biological sequence analysis: Probabilistic Models of Proteins and Nucleic Acids. Cambridge University Press; 2004.Google Scholar
- Begleiter R, El-Yaniv R, Yona G: On Prediction Using Variable Order Markov Models. Journal of Artificial Intelligence Res 2004, 22: 385–421.Google Scholar
- Gardy JL, Spencer C, Wang K, Ester M, Tusnady GE, Simon I, Hua S, deFays K, Lambert C, Nakai K, Brinkman FS: PSORT-B: improving protein subcellular localization prediction for Gram-negative bacteria. Nucleic Acids Research 2003,31(13):3613–17. 10.1093/nar/gkg602PubMed CentralPubMedView ArticleGoogle Scholar
- Emanuelsson O, Nielsen H, Brunak S, von Heijne G: Predicting subcellular localization of proteins based on their N-terminal amino acid sequence. J Mol Biol 2000, 300: 1005–1016. 10.1006/jmbi.2000.3903PubMedView ArticleGoogle Scholar
- Fan RE, Chang KW, Hsieh CJ, Wang XR, Lin CJ: LIBLINEAR: A library for large linear classification. J of Machine Learning Res 2008, 9: 1871–1874.Google Scholar
- Ansorge W: Next-generation DNA sequencing techniques. New Biotechnology 2009,25(4):195–203. 10.1016/j.nbt.2008.12.009PubMedView ArticleGoogle Scholar
- Ong CS, Zien A: An Automated Combination of Kernels for Predicting Protein Subcellular Localization. Proc of WABI 2008, 186–179.Google Scholar
- Leslie C, Eskin E, Weston J, Noble WS: Mismatch String Kernels for SVM Protein Classification. Advances in Neural Information Processing Systems (NIPS 2002) 2002.Google Scholar
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.