- Open Access
RAFTS3G: an efficient and versatile clustering software to analyses in large protein datasets
BMC Bioinformaticsvolume 20, Article number: 392 (2019)
Clustering methods are essential to partitioning biological samples being useful to minimize the information complexity in large datasets. Tools in this context usually generates data with greed algorithms that solves some Data Mining difficulties which can degrade biological relevant information during the clustering process. The lack of standardization of metrics and consistent bases also raises questions about the clustering efficiency of some methods. Benchmarks are needed to explore the full potential of clustering methods - in which alignment-free methods stand out - and the good choice of dataset makes it essentials.
Here we present a new approach to Data Mining in large protein sequences datasets, the Rapid Alignment Free Tool for Sequences Similarity Search to Groups (RAFTS3G), a method to clustering aiming of losing less biological information in the processes of generation groups. The strategy developed in our algorithm is optimized to be more astringent which reflects increase in accuracy and sensitivity in the generation of clusters in a wide range of similarity. RAFTS3G is the better choice compared to three main methods when the user wants more reliable result even ignoring the ideal threshold to clustering.
In general, RAFTS3G is able to group up to millions of biological sequences into large datasets, which is a remarkable option of efficiency in clustering. RAFTS3G compared to other “standard-gold” methods in the clustering of large biological data maintains the balance between the reduction of biological information redundancy and the creation of consistent groups. We bring the binary search concept applied to grouped sequences which shows maintaining sensitivity/accuracy relation and up to minimize the time of data generated with RAFTS3G process.
Since the emergence of large-scale genomic sequencing
, in 2002, the analyses of genomes and proteomes begun to be used and have strength, mainly in recent years. However, it was noticed that there was an exponential increase of more sequences to be deposited resulting in the need to create large databases to store such information which we call Big Data . Currently works highlight the importance of the study of large clusters: as in the prediction of structural families, identifying biologically relevant molecular features in large-scale omics experiments with variable measurements at multiple conditions and to detect in the expansion of the network of interaction between groups and subgroups of biological sequences [2,3,4]. Clustering methods are essentials for partitioning biological samples and are useful in minimizing the complexity of needed information in extensive datasets  and in bioinformatics is the first strategy to search information in biological datasets. In addition, as the size of large biological databases is extensively larger - billions of sequences are currently available for analysis - clustering algorithms generate large number of clusters and superclusters which makes manual curation of these impracticable  – i.e. UniRef consortium contains clusters with more than 302,000,000 clusters . Most methods apply the same approach: First, the similarity is calculated and then used to group objects - e.g., experimental samples or biological sequences - into clusters, however the clustering output is useful only if the clusters correspond to the biologically relevant data features that were not used to define the grouping . Currently, two tools are considered as “golds standards” in the clustering sequences to minimize redundancy in large proteins dataset: CD-HIT  and UCLUST . CD-HIT is one of the most popular tools and is the state-of-art method . UCLUST is a tool used by thousands of users around the world as high-performance clustering considered faster than the CD-HIT algorithm . However, those tools use greedy strategies for clustering. Furthermore CD-HIT does not support values lower than 40% of similarity and in lower identities whereas UCLUST degrades the quality of alignment . It is also worth pointing out that both the CD-HIT and UCLUST tools require a manual preprocessing step in which the data to be rotated by the algorithms must be organized in order of sequence size, because both algorithms select the largest to minor sequences to choose the representative sequence to the group and align the others from them, not being a random process. Therefore, both CD-HIT and UCLUST are not reliable choices for clustering in large datasets with values less than 30% of similarity so trivial to search sequences with homologies in remotely structures . The most efficient techniques for this prediction use as gold standard the Basic Local Alignment Search Tool (BLAST) ‘all-against-all’ or, in another cases, Markov Clustering (MCL) method adaptations . However, these tools are dependents on alignment metrics requiring a lot of processing and time to generate results mainly in large datasets [16,17,18].
Alignment-free methods are strong alternatives to alignment-dependent techniques and are also efficient in minimizing the redundancy of biological data its computationally fast and use less memory compared to alignment-based methods . A method that has been highlighting among the clustering techniques of large databases to solve the main time and memory bottlenecks of existing clustering the algorithms is MMSeqs2-Linclust, a deep clustering approach . This method explores the alignment-free analyses and apply two main steps to clustering: the global Hamming distance and the gapless local alignment extending the k-mer match. Sequence pairs are generated under the conditions that satisfying the clustering criteria - e.g., on the E-value, sequence similarity, and sequence coverage- and are linked by an edge. In the end, the greedy incremental algorithm locates a cluster so that each input sequence has an edge to the representative sequence of its cluster . Ultimately, alignment-free methods have been applied to problems ranging from whole-genome and are particularly useful for processing and analyzing Next-Generation Sequencing (NGS) data. However, the benchmark data sets are required to explore the full potential of alignment-free methods .
The validity of the clusters is challenging: information from external clusters are needed because they are not known in advance. At this point, the lack of a priori knowledge about the number of clusters underlying in the dataset makes it indispensable and an efficient metric is necessary to compare clustering solutions with different number of clusters . Validity is constantly being questioned because there is a need for standardization of metrics, besides the application of internal and external metrics and the use of consistent bases of biological value . Another point is the application of a high level of programming skills on the part of researchers to analyze large volumes of data : generally, each tool uses a different output and makes difficult the manipulation of data which hinders the fluidity of the researches .
To explore the potential of the alignment-free method associated with a strategy that combines hashes and BCOM matrices to reduce the need for the slow sequence alignments, we have developed the RAFTS3G. We incorporated the binary search as an option cluster input criterion to align the best n candidates, a new alternative proposal for clustering analyses in proteins sequences data. We compared RAFTS3G with three main clustering methods exploring standard metrics applied to database “gold standard” of enzymes family adopting as criterion the default parameters of all methods.
To minimize time and maintaining consistency in data analysis with proteins, we developed Rapid Alignment Free Tool for Sequences Similarity Search to Groups (RAFTS3G) tool. RAFTS3G was written in MATLAB v2017a explores the RAFTS3 engineer (Additional file 1: Figure S1) and uses integrates functions, the Bioinformatics Toolbox and an in-house library.
The RAFTS3G overview
RAFTS3G applies as search engine RAFTS3  tool, which purpose is to perform faster by minimizing disk access storing sequences information in RAM and in addition to reducing the need for slow sequence alignments. RAFTS3 has a hashing strategy based on k-mers to directly access sequence data – the sequence itself and the Co-Occurrence Matrix of amino acid residues (BCOM). BCOM are sets of 50 bytes containing a binary matrix within amino-acid sequential co-occurrence data for a given sequence. The comparison between BCOM of two sequences is faster than to alignment them to get similarity metric. When RAFTS3 searches for sequence similarities, however, it allows the user to choose to align a set of the top n selected candidates within some k-mer match against to a query sequence. The metric provided by BCOM  is effective to sort a set of sequences according to their similarity, the similarity measure based on identities, enabled when alignment is performed, is desirable when the intention is to hold clusters and it is often selected as cut-off criterion . Once aligning every subject candidate would be impeditive to a rapid approach sequence grouping algorithm, we studied ways to minimize the need of alignment in RAFTS3G; it will be discussed forward, while we present the algorithm.
From a set of input sequences in a FASTA format - variable or file -, for each sequence not grouped yet, RAFTS3G exploits a formatted RAFTS3 data base searching for similar sequences. Candidates are ordered by higher BCOM similarity to the query. To select which from candidates should be in the same cluster of the query sequence, given a cut-off value (RAFTS3 self-score), the user can choose:
Align the query with up to a limited n number of the BCOM ordered candidates, living behind the rest.
Make a binary search aligning candidates/query to find the cutting point where all sequences of lower order should be as similar or more than the sequence in this point. Sequences of higher order are likely less similar then the stipulated by the cutoff criterion and are left.
The step in ii) is the only change we made in original RAFTS3 approach in order to program RAFTS3G. The main gain of the binary search approach is to allow the constrution of a cluster within less steps, since it finds most sequences related to a query in a single search, aligning only a relatively small number of candidates (O(log2(n)).
In both cases we have a list of sequences to group that are supposed to be at least as similar to the query as the measure defined in cut-off.
It remains now review the assembled groups based on the sequences to group:
if the query found already grouped sequences, all the groups found are joined in a single one and all other sequences to group are added in this group;
if none of the sequence to group is member of a previously created group then a new group is built and these sequences are added to it.
While there are sequences to be analyzed these steps will be repeated for each of them. See (Fig. 1). The RAFT3G output is easier to be manipulated by the end user because it is in FASTA format with an extra log is generated with clusters information (Additional file 1: Figure S2).
RAFTS3G clustering in large dataset
We performed RAFTS3G using the Ref-Seq Non-Redundant protein from NCBI database (NCBI/NR)  - with 78,002,046 sequences deposited at this release. We generated 12,594,179 Total clusters of which 4,127,885 are non-unique clusters and 8,466,294 are unique clusters. Twenty-one clusters have more than 100,000 grouped protein sequences and in nine of them exceed 200,000 sequences clustered. In Fig. 2 the 30 largest clusters are represented, according to the number of sequences in each cluster. Therefore, with these results, RAFTS3G it is possible to generate clusters in a higher set of data. Due to this large set of data we are evaluating the results obtained allow us to bring more information about the developed clustering techniques in future works.
Benchmark standardization with F1-score
The choice of a good basis is essential for the reliability of the metrics, so we chose the GOLD/Brown base from ASTRAL/SCOPe . For the validation of clusters, we used F1-Score, an external metric that provides the balance between the accuracy and sensitivity measures [31, 32]. The GOLD database - a collection “gold standard” of enzymes families experimentally validated  totalizing 866 sequences - to evaluation of clusters generated for RAFTS3G compared to three highlighted methods. The Brown database is a collection of experimentally classified enzymes with extreme remote similarities and this database is a challenge to be correclty grouped because extreme remote similarities sequences have low identity which generates many false positives in the clustering process . In comparison with CD-HIT we exemplifying this difficult evaluated the F1-Score, accuracy and sensibility metrics (Additional file 1: Table S5) and we are improving the RAFTS3G to obtain more hits with these data sets. We analysed RAFTS3G in 0.5 of similarity threshold in 3 representative clusters from Swissprot/UniProtKB with remote similarity: Apolipoprotein C-IV, Period circadian protein and Ribulose bisphosphate carboxylase/oxygenase activase. We generated the distance matrix calculing the sequences alignments to each cluster and we found that RAFTS3G had grouped sequences with great distances and no false positives (Additional file 1: Figure S4). These suggests that RAFTS3G was able to group distance sequences with low similarities.
According to the results obtained with GOLD database, in low similarities, between 0.2–0.4 intervals of threshold, RAFTS3G presents sensitivity above the other compared tools but without significance. We noticed that all tools seem to have similar performance in similarity of 0.3 - excepts CD-HIT because does not generate groups with this threshold. From the cut-off lines between 0.4 and 0.9 of similarity, we observed the ability of RATS3G to group consistent sequences compared to MMSeqs2 (Linclust algorithm) - method which stands out in relation the others two tools Usearch (Uclust algorithm) and CD-HIT. As all the methods compared are developed to reduce redundancy, in the higher similarities between the values of 0.8–0.9 of similarity we observed an equity between the results obtained between MMSeqs, USEARCH and CD-HIT. In this range RAFTS3G has a 10% gain of F1-Score in relation to the others. (Comparison with CD-HIT and UCLUST performed against Astral/SCOPe of proteins database in 20 to 90% of similarity is available at Additional file 1: Table S2 and S4).
Analyzing these points, RAFT3G is the best choice optimized to be more permissible to members inclusion when the clusters increase (Fig. 3). This is interesting when the user wants to “guess” or to “risk” a data set when the similarity does is not known by user. Other methods generate more restricted clusters and choose to lose these informations. In metagenome data, for example, where the collected material is very heterogeneous and abundant, using a strategy which increases sensitivity or probability of clustering sequences mainly at an early stage of data mining is crucial to the success of the experimentation and analysis.
Binary search input criteria
In the RAFTS3G overview, we bring the proposal of a binary search to the assembly of the clusters after the selection of the candidates obtained by the RAFTS3 engineering, instead of the cut-off for the groups to be based on the alignment of the sequences by the selection of n candidates. Results of clusters generated with the GOLD base (Astral / SCOPe) suggest that this type of strategy maintains the sensitivity / accuracy ratio (Fig. 4). In addition to being significantly high - around 91% of F1-Score for RAFTS3G in relation to 0.87 in MMSeqs, 0.73 of USEARCH and 0.72 of CD-HIT - another observable advantage is in reducing time - binary search reduced by up to 73% of the overall execution time of RAFTS3G - maintaining the quality of the data generated.
The goal of this study is to provide an alternative to clustering analyses with reduced losses of biological data information improving the alignment-free concept. RAFTS3G is able to group up to millions of sequences. Furthermore, we brought a benchmark analysis using the F1-score as an external metric to evaluate the performance of the main clustering methods by exploring a wide range of similarity and found that the RAFTS3G strategy is the best optimized - to be more permissive - which reflects in greater accuracy and sensitivity in generating clusters with consistent biological content. The binary search input criteria for creating groups demonstrates to be efficient to create or to integrate candidate groups as the overall alignment of n candidates.
We hope the RAFTS3G algorithm will prove helpful to assist the researcher to explore the widest range of available data and to make them more consistent.
Data and RAFTS3G availability
Project name: RAFTS3G.
Project Home Page: https://sourceforge.net/projects/rafts-g/
Operating System: Windows and Linux (× 86 and × 64 versions).
Programming Language: Designed in Matlab® v2012.
Other requirements: MCR runtime (v7.17) is required to runs.
License: the software is under licensed by Matlab® v 2012.
Any restrictions to use by non-academics: none.
Availability of data and materials
The RAFT3G is freely accessible and can be downloaded without user registration at: https://sourceforge.net/projects/rafts-g/ and additional informations in supplementary material.
Emms DM, Kelly S. OrthoFinder: solving fundamental biases in whole genome comparisons dramatically improves orthogroup inference accuracy. Genome Biol. 2015;16(1):157. https://doi.org/10.1186/s13059-015-0721-2.
Schwämmle V, Jensen ON. VSClust: feature-based variance-sensitive clustering of omics data. Bioinformatics. 2018;34(17):2965–72. https://doi.org/10.1093/bioinformatics/bty224.
Adams J, Mansfield MJ, Richard DJ, Doxey AC. Lineage-specific mutational clustering in protein structures predicts evolutionary shifts in function. Bioinformatics. 2017;33(9):1338–45. https://doi.org/10.1093/bioinformatics/btw815.
Städler N, Dondelinger F, Hill SM, Akbani R, Lu Y, Mills GB, Mukherjee S. Molecular heterogeneity at the network level: high-dimensional testing, clustering and a TCGA case study. Oxf J Bioinforma. 2017;33(18):2890–6. https://doi.org/10.1093/bioinformatics/btx322.
Bursteinas B, Britto R, Bely B, Auchincloss A, Rivoire C, Redaschi N, et al. Minimizing proteome redundancy in the UniProt knowledgebase. Database. 2016;2016:1–9. https://doi.org/10.1093/database/baw139.
Chen C, Huang H, Wu CH. Protein Bioinformatics Databases and Resources. In: Wu C, Arighi C, Ross K, editors. Protein Bioinformatics. Methods in Molecular Biology, vol. 1558. New York: Humana Press; 2017.
Consortium TU. UniProt: a worldwide hub of protein knowledge. Nucleic Acids Res. 2019;47(D1):D506–15. https://doi.org/10.1093/nar/gky1049.
Altman N, Krzywinski M. Points of significance: clustering. Nat Methods. 2017;14(6):545–6. https://doi.org/10.1038/nmeth.4299.
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–9.
Apweiler R, Bairoch A, Wu CH, Barker WC, Boeckmann B, Ferro S, Yeh L-SL. UniProt: the universal protein knowledgebase. Nucleic Acids Res. 2004;32(Database issue):D115–9.
Chen Q, Wan Y, Lei Y, Zobel J, Verspoor K. Evaluation of CD-HIT for constructing non-redundant databases, Proceedings - 2016 IEEE International Conference on Bioinformatics and Biomedicine, BIBM 2016, 703–706; 2017. https://doi.org/10.1109/BIBM.2016.7822604.
Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010;26(19):2460–1.
Kopylova E, Navas-Molina JA, Mercier C, Xu ZZ, Mahé F, He Y, et al. Open-source sequence clustering methods improve the state of the art. MSystems. 2016;1(1):e00003–15. https://doi.org/10.1128/mSystems.00003-15.
Chen J, Long R, Wang XL, Liu B, Chou KC. DRHP-PseRA: detecting remote homology proteins using profile-based pseudo protein sequence and rank aggregation. Sci Rep. 2016;6(August):1–7. https://doi.org/10.1038/srep32333.
Azad A, Pavlopoulos GA, Ouzounis CA, Kyrpides NC, Buluç A. HipMCL: a high-performance parallel implementation of the Markov clustering algorithm for large-scale networks. Nucleic Acids Res. 2018;46(6):1–11. https://doi.org/10.1093/nar/gkx1313.
Vinga S, Almeida J. Alignment-free sequence comparison--a review. Bioinformatics. 2003;19:513–23. https://doi.org/10.1093/bioinformatics/btg005.
Mahmood K, Webb GI, Song J, Whisstock JC, Konagurthu AS. Efficient large-scale protein sequence comparison and gene matching to identify orthologs and co- orthologs. Nucleic Acids Res. 2012;40:e44. https://doi.org/10.1093/nar/gkr1261.
Tabari E, Su Z. PorthoMCL: parallel orthology prediction using MCL for the realm of massive genome availability. Big Data Anal. 2017;2(1):4. https://doi.org/10.1186/s41044-016-0019-8.
Zielezinski A, Vinga S, Almeida J, Karlowski WM. Alignment-free sequence comparison: benefits, applications, and tools. Genome Biol. 2017. https://doi.org/10.1186/s13059-017-1319-7.
Steinegger M, Söding J. MMseqs2 enables sensitive protein sequence searching for the analysis of massive data sets. Nat Biotechnol. 2017;35(11):2–4. https://doi.org/10.1038/nbt.3988.
Steinegger M, Söding J. Clustering huge protein sequence sets in linear time. Nat Commun. 2018;9(1). https://doi.org/10.1038/s41467-018-04964-5.
Ren J, Bai X, Lu YY, Tang K, Wang Y, Reinert G, Sun F. Alignment-free sequence analysis and applications. Annu Rev Biomed Data Sci. 2018;1:93-114. https://doi.org/10.1146/annurev-biodatasci-080917-013431.
Srivastava A, Baranwal M, Salapaka S. On the persistence of clustering solutions and true number of clusters in a dataset. Retrieved from arXiv 2018. http://arxiv.org/abs/1811.00102.
Wiwie C, Baumbach J, Röttger R. Comparing the performance of biomedical clustering methods. Nat Methods. 2015;12(11):1033–8. https://doi.org/10.1038/nmeth.3583.
Nichio BTL, Marchaukoski JN, Raittz RT. New tools in orthology analysis: a brief review of promising perspectives. Front Genet. 2017;8(OCT):1–12. https://doi.org/10.3389/fgene.2017.00165.
Pavlopoulos GA. How to cluster protein sequences: tools, tips and commands. MOJ Proteomics Bioinform. 2017;5(5). https://doi.org/10.15406/mojpb.2017.05.00174.
Vialle RA, Pedrosa FO, Weiss VA, Guizelini D, Tibaes JH, Marchaukoski JN, Raittz RT. RAFTS3: rapid alignment-free tool for sequence similarity search. bioRxiv. 2016;55269. https://doi.org/10.1101/055269.
Krawetz S. Bioinformatics for systems biology. Cap. 27 Clustering algorithms, vol. 9781597454407: Humana Press; 2009. ISBN 978-1-59745-440-7. https://doi.org/10.1007/978-1-59745-440-7.
Marchler-Bauer A, Schoch CL, Canese K, Schneider VA, Hefferon T, Bolton EE, Kimchi A. Database resources of the National Center for biotechnology information. Nucleic Acids Res. 2018;47(D1):D23–8. https://doi.org/10.1093/nar/gky1069.
Brown SD, Gerlt JA, Seffernick JL, Babbitt P. A gold standard set of mechanistically diverse enzyme superfamilies. Genome Biol. 2006;7(1):1–15. https://doi.org/10.1186/gb-2006-7-1-r8.
Xiong M, Liu X, Hao M, Li Y, Shugart YY, Qiao C, et al. Nuclear norm clustering: a promising alternative method for clustering tasks. Sci Rep. 2018;8(1):1–7. https://doi.org/10.1038/s41598-018-29246-4.
Bernardes JS, Vieira FRJ, Costa LMM, Zaverucha G. Evaluation and improvements of clustering algorithms for detecting remote homologous protein families. BMC Bioinformatics. 2015;16(1). https://doi.org/10.1186/s12859-014-0445-4.
Fox NK, Brenner SE, Chandonia JM. SCOPe: structural classification of proteins - extended, integrating SCOP and ASTRAL data and classification of new structures. Nucleic Acids Res. 2014;42(D1). https://doi.org/10.1093/nar/gkt1240.
Federal University of Paraná (UFPR), CAPES (Coordination for the Improvement of Higher Education Personnel) & Araucária Fundations to support this work.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Support material - system requirements, extra information about RAFTS3 engineering, methodology overflow, tests, additional links and literatures. (DOCX 808 kb)