# denoising vs dereplication

**URL:** https://forum.qiime2.org/t/denoising-vs-dereplication/18990
**Category:** User Support
**Created:** [March 26, 2021, 9:39am UTC](https://forum.qiime2.org/t/denoising-vs-dereplication/18990 "2021-03-26T09:39:06Z")
**Posts on this page:** 6
**Page:** 1

<div class="post-metadata">

### Author: ![taf1g17](https://forum.qiime2.org/letter_avatar_proxy/v4/letter/t/e99b99/32.png) [@taf1g17](https://forum.qiime2.org/u/taf1g17)
#### Post date: [March 26, 2021, 9:39am UTC](https://forum.qiime2.org/t/denoising-vs-dereplication/18990/1 "2021-03-26T09:39:06Z")

</div>

Hello,

I would just like to ask what the difference is between denoising and dereplication when used inside work-flows. I have been using vsearch and I have joined, quality filtered and dereplicated (as described in the overview), however I have NOT performed any denoising specifically, either using DADA2 or DeBlur. I know that DADA2 is an independent plugin therefore performed denoising, dereplication, filtering etc internally, however if I have dereplicated using vsearch dereplicate-sequences, do I still need to run additional denoising using deblur? Basically are the vsearch dereplicate-sequences and deblur denoise-16S commands homologous, and if not, what is the difference and should I therefore use both to achieve accurate dereplication and denoising?

Thanks and have a great day,

Tom

---

<div class="post-metadata">

### Author: ![llenzi](https://forum.qiime2.org/letter_avatar_proxy/v4/letter/l/9fc29f/32.png) [@llenzi](https://forum.qiime2.org/u/llenzi)
#### Post date: [March 26, 2021, 10:30am UTC](https://forum.qiime2.org/t/denoising-vs-dereplication/18990/2 "2021-03-26T10:30:47Z")

</div>

Hi @taf1g17

In my view, denoising is a step where you try to infer the set of original amplicon sequences by analysing the error-containing sequences obtained by the combination of PCR/sequencing production steps. Dereplicating a sequence set is merely collapsing sequences at 100% of similarity, which does retain any sequences containing errors.

By using vsearch to join - quality filter - dereplicate your sequences, you obtained a non-redundant set of sequences common for all samples, which you now need to use as reference to align all the sequences on, in order to obtain an abundance table ("to see how many times each sequences is present in any sample"). You may probably use vsearch again for this step! At this point you may notice that there are plenty of sequences with very low total abundances, which may be due to sequences containing error, and you need to discard these before proceeding with the analysis.

Deblur and dada2 are two alternative tools designed to denoising the dataset as well as report how many times each original amplicon (feature in the final table) is present in each sample. Hence they produce two main outputs: the final feature table and the set of identified amplicons (you may see them as error-correcting the seuences -\> dereplicating the error free sequences -\> counting the abundance of the sequences in each sample).

Hope it helps  
Luca

---

<div class="post-metadata">

### Author: ![taf1g17](https://forum.qiime2.org/letter_avatar_proxy/v4/letter/t/e99b99/32.png) [@taf1g17](https://forum.qiime2.org/u/taf1g17)
#### Post date: [March 26, 2021, 11:21am UTC](https://forum.qiime2.org/t/denoising-vs-dereplication/18990/3 "2021-03-26T11:21:15Z")

</div>

Hello Luca,

Thanks a lot for getting back to me so soon!

So what does denoising with deblur achieve that I cannot achieve independently using Vsearch? Basically do I need to include additional denoising steps IF I have already quality filtered, dereplicated etc?

Thanks,

Tom

---

<div class="post-metadata">

### Author: ![llenzi](https://forum.qiime2.org/letter_avatar_proxy/v4/letter/l/9fc29f/32.png) [@llenzi](https://forum.qiime2.org/u/llenzi)
#### Post date: [March 26, 2021, 11:46am UTC](https://forum.qiime2.org/t/denoising-vs-dereplication/18990/4 "2021-03-26T11:46:59Z")

</div>

Hi @taf1g17

Both deblur and dada2 should return error-free sequences, which will be beneficial for the taxonomy assignment step. Also a smaller data set to handle after the denoising, because sequences with errors are not discarded but reverted to the original state. You could compensate this via filtering the low abundance OTUs, but you always taking the risk to discard real low abundant clusters as well by applying a filtering by abundance.

There are few interesting discussion on comparing OTUs (clusters obtained by vsearch, as in old qiime1) and ASVs/ESVs (amplicon sequence variants obtained by denoisers as dada2 or deblur), one is:

> [@ESV vs OTU deflation (QIIME1 vs 2)](https://forum.qiime2.org/t/esv-vs-otu-deflation-qiime1-vs-2/14867):
>
> This is a more general question. Now that we are switching completely to QIIME2, i did one comparison between 97% OTU clustering in QIIME 1 and 99% ESV clustering in QIIME 2/deblur. The dataset is a fairly standard MiSeq 2x250 bp run with \>100 bp overlap. I compared two experiments by reads surviving the pipeline raw OTUs/ESVs identified by sample. I found that QIIME2 accepts more reads, while identifying much less ESVs than OTUs. This was surprising because i would expect that 99% clusteri…

Another good one is:

> [@Major differences in alph diversity between q1 and q2](https://forum.qiime2.org/t/major-differences-in-alph-diversity-between-q1-and-q2/2419/3):
>
> Hi @uria, The short answer: yes, these differences are expected, sensible, and acceptable. @colinbrislawn's answer on the qiime1 forum does a great job of summarizing the cause of this discrepancy — the methods used by that user (and yourself) in qiime1 and qiime2 are very different — but does not go into the nitty gritty details. Let's discuss those here. The long answer: The differences that you are observing are not differences between qiime1 and qiime2, strictly speaking. Each of these i…

Which remind me about chimera filtering steps embedded into deblur and dada2, that seems missing in your pipeline from the step you mention.

One question, what are you going to do next with the dereplicated sequences?

Keep in mind, working with ASVs or OTUs are both accepted way to get your results, so as long you are doing correctly they are both valid! Mostly is a matter of preference and knowing what you are working with!

Luca

---

<div class="post-metadata">

### Author: ![Mehrbod\_Estaki](https://forum.qiime2.org/user_avatar/forum.qiime2.org/mehrbod_estaki/32/4001_2.png) [@Mehrbod\_Estaki](https://forum.qiime2.org/u/Mehrbod_Estaki)
#### Post date: [March 26, 2021, 9:05pm UTC](https://forum.qiime2.org/t/denoising-vs-dereplication/18990/5 "2021-03-26T21:05:20Z")

</div>

Hi @taf1g17 ,  
To add to @llenzi 's already excellent and comprehensive answer.

> [@taf1g17](#):
>
> Basically do I need to include additional denoising steps IF I have already quality filtered, dereplicated etc?

Denoising with Deblur/DADA2 requires you to operate on the raw fastq files. If you have already used vsearch for clustering/dereplication you can no longer denoise even if you wanted to. If you want (and I personally recommend this) is to use either Deblur/DADA2 to first denoise your raw reads, utilizing their superior quality controls and error correction, and then if you are interested in OTU clustering then you can also cluster the output from Deblur/DADA2 afterwards using vsearch as shown [here](https://docs.qiime2.org/2021.2/tutorials/otu-clustering/). There's a video I made on this topic [here](https://www.youtube.com/watch?v=PmtqSa4Z1TQ) which also goes a bit more in detail.

---

<div class="post-metadata">

### Author: ![system](https://forum-qiime2-org.s3.dualstack.us-west-2.amazonaws.com/original/3X/2/1/21af5fe23cb6f4579467c66a9ed94e55274ca7bd.svg) [@system](https://forum.qiime2.org/u/system)
#### Post date: [April 27, 2021, 3:06am UTC](https://forum.qiime2.org/t/denoising-vs-dereplication/18990/6 "2021-04-27T03:06:07Z")

</div>

This topic was automatically closed 31 days after the last reply. New replies are no longer allowed.
