I am working on a meta analysis of both prokaryotic and eukaryotic diversity in coral reefs analyzing hundreds of samples each of 16S, 18S, and CO1 amplicon data. I’m using q2 dublur denoise-16S for the 16S data and I’m confident that the default parameters for that method are good to go out of the box. I’m using q2 deblur denoise-other for the 18S and CO1 amplicon data. My collaborators and I are a bit concerned about whether the default parameters for q2 deblur denoise-other are optimal for these amplicons.
Specifically, I’m wondering if I need to worry about optimizing and changing the --p-indel-prob and --p-indel-max parameters for each specific amplicon that I’m working with. I also wonder if the default error profile (–error-dist option in the standalone version of deblur) is sub-optimal for amplicons other than 16S. I assume that all of these deblur parameters were optimized for use with 16S and that they might not be perfect for the other amplicons. It seems like the probablility of indels (and maybe the error profile) would be dependent on how highly conserved the amplicon is and the rate of evolution of the amplicon region. I want to confirm that I’m not violating any assumptions of the algorithm or the default error profile in using deblur on 18S and CO1.
If these are things that I really need to worry about in applying deblur to non-16S amplicons, could you provide some insight on the best way to approach optimizing the deblur parameters for my needs. I looked through the deblur paper and supplementary materials, but don’t see explicitly how the default settings were chosen.
I am assuming that maybe doing a multiple sequence alignment and examining the fequency of indels could help in optimizing the --p-indel-prob and --p-indel-max parameters, but I’m not sure how to go about optimizing the --error-dist parameter (or if that is even necessary). It looks like the --error-dist parameter isn’t accessible through the qiime2 interface anyway.
This question is probably best suited for @wasade or one of the other deblur devs, but I welcome any feedback from community members who have experience applying deblur to 18S and/or CO1 amplicons. Apologies for the super long post!