Reference database for clustering

[SoilRotifer (Mike Robeson)]

Hi @Rakaya,

RE 1:

@Rakaya

By this I would have the dereplicated database, is this correct?

Yep.

RE 2

@Rakaya

Can I use dereplicated database instead of clustered one for clustering?

There is no need to cluster your reference sequences to match the clustering of your data, for reasons previously mentioned within these threads:

@SoilRotifer

Anyway, the old way of pre-clustering reference reads to 97% or 94% sequence similarity (as done with Greengenes and SILVA) was simply a logistical way to reduce the size of the reference database even further, so that it'll run on computers with limited memory and cpu power.

@SoilRotifer

If possible, it is best to try and use reference reads that have been minimally clustered or simply dereplicated, i.e. at 99% or 100% similarity. As you'll be able to more accurately classify your reads, that is you have more reference data to use.

@SoilRotifer

So, if you are using reference reads that are 99% - 100% then you can use them to classify any reads that are clustered at different similarities, e.g. 97%, 94%, etc...

In fact, some consider classifying 97% clustered OTUs against a 97% clustered reference database not a great idea. That is, as your nearest reference sequence can be up to 6% away, further reducing classification.

Does this help?

Additionally:

@SoilRotifer

In general, we do not recommend clustering reads for the purposes of making a reference database as your ability to correctly assign taxonomy to your reads declines. This is covered in our RESCRIPt manuscript. Usually, clustering reference sequences is performed to reduce the file and memory size of the reference database when computational resources are limited.

RE 3:

@Rakaya

If I want to cluster my database at 85% identity threshold, what would be the command to use?

Simply add the paramater --p-perc-identity 0.85 to that command. : I would not recommend using 85%. This is just provided as an example to drastically reduce the database size.

RE 4:

@Rakaya

Can I use 99% silva to cluster at 85%?

Yes you can classify your 85% clustered reads against the 99% SILVA database. Again, I'd avoid 85%. Also, the content I linked for RE 2, applies here too.