How to cluster peptide/protein sequences using cd-hit software?

/

Cd-hit is one of the most widely used programs to cluster biological sequences [1]. It helps in removing the redundant sequences and provides better results in the sequence analyses. Cd-hit is used for sequence-based clustering by making clusters of a particular cut off provided as an input. It uses a single linkage clustering and finds a representative sequence for each cluster. In this article, we will learn how to cluster a set of protein sequences using cd-hit software.

Cd-hit package has many programs for clustering different kinds of sequences. For example, the cd-hit program is used to cluster peptide sequences, cd-hit-est is used to cluster nucleotide sequences, and even this package can compare two different databases using cd-hit-2d and cd-hit-est-2d to compare peptide and nucleotide databases respectively [1]. In this tutorial, we are using the cd-hit program which is used to cluster a group of peptide sequences. The complete package of cd-hit can be downloaded from here.

Prepare input file

The input file consists of all the peptide or protein sequences in FASTA format. There is no need to format the FASTA header of the sequences. The software manages it on its own.

Basic commands

$ cd-hit -i input.fasta -o db100 -c 1.00 -n 5 -M 2000

where,

-i = input

-o = output

-c = cut-off

-n = word size:

n=5 for thresholds 0.7 ~ 1.0

n=4 for thresholds 0.6 ~ 0.7

n=3 for thresholds 0.5 ~ 0.6

n=2 for thresholds 0.4 ~ 0.5

-M = maximum available memory

T0 cluster the sequences at 97% similarity cut-off

$ cd-hit -i input.fasta -o db90 -c 0.97 -n 5 -M 2000

Output

The output of cd-hit provides two different files:

1. A FASTA file of the representative sequences of all the clustered sequences.

2. A text file listing all the clusters showing a representative sequence signified with a ‘*’ at the end of the header of the sequence.

There are many other options which you can define in the command line including -G to use global sequence identity, -t to set tolerance for redundancy, -l to set length of throw_away_sequences, and -d to adjust the description of sequences in the .clstr output file. You can read about the command-line options either in the user guide provided at the cd-hit website (http://www.bioinformatics.org/cd-hit/cd-hit-user-guide.pdf) or by entering the help command ($ cd-hit --help).

References

  1. Li, W., & Godzik, A. (2006). Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics22(13), 1658-1659.
Tariq is founder of Bioinformatics Review and a professional Software Developer at IQL Technologies. His areas of expertise include algorithm design, phylogenetics, MicroArray, Plant Systematics, and genome data analysis. If you have questions, reach out to him via his homepage.

Leave a Reply

HOW TO CITE THIS ARTICLE Tariq Abdullah (2019). How to cluster peptide/protein sequences using cd-hit software?. Bioinformatics Review, 5 (09)
Previous Story

How to search motif pattern in FASTA sequences using Perl hash?

Next Story

How to blast against a particular set of local sequences (local database)?

Latest from Software

Willing to stay updated?

By investing less than 30 seconds you can start recieving all our new articles in your mailbox. Stay updated with latest Bioinformatics Research, trends and tools of trade.

 

Lost your password? Please enter your email address. You will receive mail with link to set new password.

0 $0.00