# Taxonomy assignement issues

**URL:** https://forum.qiime2.org/t/taxonomy-assignement-issues/4365
**Category:** User Support
**Tags:** feature-classifier
**Created:** [May 30, 2018, 11:40am UTC](https://forum.qiime2.org/t/taxonomy-assignement-issues/4365 "2018-05-30T11:40:55Z")
**Posts on this page:** 8
**Page:** 1

<div class="post-metadata">

### Author: ![BiOMan](https://forum.qiime2.org/letter_avatar_proxy/v4/letter/b/0ea827/32.png) [@BiOMan](https://forum.qiime2.org/u/BiOMan)
#### Post date: [May 30, 2018, 11:40am UTC](https://forum.qiime2.org/t/taxonomy-assignement-issues/4365/1 "2018-05-30T11:40:55Z")

</div>

Hello Everyone,

I am new qiime2 user, I tried to assign taxonomy to my feature so I have tried:  
Sklearn but at certains point it crashes ( problem of memory I have only 8GB on my machine and I alocate 4GB for my VM).  
I have tried also BLAST and vserach for the first one after one hour I was oligied to quite the process (too long) the same for Vsearch.

Is there any other way to use to have my toxonomy assignement ?? could I use Qiime1 (because the last year when i used qiime1 it works).

Thanks in advance.

Manil

---

<div class="post-metadata">

### Author: ![Nicholas\_Bokulich](https://forum.qiime2.org/user_avatar/forum.qiime2.org/nicholas_bokulich/32/19937_2.png) [@Nicholas\_Bokulich](https://forum.qiime2.org/u/Nicholas_Bokulich)
#### Post date: [May 30, 2018, 12:20pm UTC](https://forum.qiime2.org/t/taxonomy-assignement-issues/4365/2 "2018-05-30T12:20:10Z")

</div>



---

<div class="post-metadata">

### Author: ![Nicholas\_Bokulich](https://forum.qiime2.org/user_avatar/forum.qiime2.org/nicholas_bokulich/32/19937_2.png) [@Nicholas\_Bokulich](https://forum.qiime2.org/u/Nicholas_Bokulich)
#### Post date: [May 30, 2018, 3:06pm UTC](https://forum.qiime2.org/t/taxonomy-assignement-issues/4365/3 "2018-05-30T15:06:25Z")

</div>

Hi @BiOMan,

> [@BiOMan](#):
>
> I have tried also BLAST and vserach for the first one after one hour I was oligied to quite the process (too long) the same for Vsearch.

Runtime is really going to depend on:

1. the number of query sequences (features)
2. the length of query/reference sequences
3. the number of reference sequences

1 hour is really not much time for many datasets, depending on the upstream analysis steps. I would recommend letting those jobs run a little longer.

You could also try trimming your reference sequences to the same length/primer sites as the query sequences by using `qiime feature-classifier extract-sequences`. That would speed up blast and vsearch.

Using a smaller reference set, e.g., greengenes instead of silva, would also help.

> [@BiOMan](#):
>
> could I use Qiime1 (because the last year when i used qiime1 it works).

Again, runtime is going to depend on several factors; when you did this a year ago you may have had fewer or shorter query sequences.

The uclust-based taxonomy classifier in qiime1 is much faster than other methods; you could give that a try if 1 hr is really too long to wait.

Good luck!

---

<div class="post-metadata">

### Author: ![Nicholas\_Bokulich](https://forum.qiime2.org/user_avatar/forum.qiime2.org/nicholas_bokulich/32/19937_2.png) [@Nicholas\_Bokulich](https://forum.qiime2.org/u/Nicholas_Bokulich)
#### Post date: [May 30, 2018, 3:06pm UTC](https://forum.qiime2.org/t/taxonomy-assignement-issues/4365/4 "2018-05-30T15:06:33Z")

</div>



---

<div class="post-metadata">

### Author: ![BiOMan](https://forum.qiime2.org/letter_avatar_proxy/v4/letter/b/0ea827/32.png) [@BiOMan](https://forum.qiime2.org/u/BiOMan)
#### Post date: [May 31, 2018, 2:00pm UTC](https://forum.qiime2.org/t/taxonomy-assignement-issues/4365/5 "2018-05-31T14:00:35Z")

</div>

Hi @Nicholas_Bokulich

Thank you for your response, sorry I did not give you the details.

\*For my data set it is only 9 samples (number of query sequences is 380 759).

\*For the reference sequences I am using SILVA-199-99-nb-classifier.qza (my mentors insists that I use the SILVA database because we are working on Crohn's patients gut micofolra and for them greengenes is not appropriate as a database because it is used more for ecology (can you confirm this statement?!).

\*For the third point I did not really get it : do you mean that for example I am working on the V3-V3 region so I need to trim the ref seq only on this region instead of using the hole sequences in the SILVA-199-99-nb-classifier.qza??

- For the runtime at the beginning I thought that was a bug but now it is clear for me after your explanation.

- we figure out a solution in my lab and I will Run the SKearn algo from a cluster (so I hope no memory issues anymore). Because after that I will have a data set of 50 patients.

Thanks a lot @Nicholas_Bokulich

Manil

---

<div class="post-metadata">

### Author: ![Nicholas\_Bokulich](https://forum.qiime2.org/user_avatar/forum.qiime2.org/nicholas_bokulich/32/19937_2.png) [@Nicholas\_Bokulich](https://forum.qiime2.org/u/Nicholas_Bokulich)
#### Post date: [June 1, 2018, 11:12am UTC](https://forum.qiime2.org/t/taxonomy-assignement-issues/4365/6 "2018-06-01T11:12:21Z")

</div>

> [@BiOMan](#):
>
> my mentors insists that I use the SILVA database because we are working on Crohn’s patients gut micofolra and for them greengenes is not appropriate as a database because it is used more for ecology

Greengenes was not designed specifically for non-host-associated samples, nor was SILVA designed specifically for host-associated samples. Both are quite general, so I would disagree with that statement on the grounds of "they are used more for X". However, there are other reasons folks might choose one over the other — SILVA has been updated much more recently, so taxonomy IDs may be better updated for your system of interest.

> [@BiOMan](#):
>
> \*For the third point I did not really get it : do you mean that for example I am working on the V3-V3 region so I need to trim the ref seq only on this region instead of using the hole sequences in the SILVA-199-99-nb-classifier.qza??

You do not need to, but it will reduce runtime, and will (slightly) increase classification accuracy. See [here](https://docs.qiime2.org/2018.4/tutorials/feature-classifier/) for more details (we also have pre-trained V4 classifiers [here](https://docs.qiime2.org/2018.4/data-resources/) if those are useful to you).

> [@BiOMan](#):
>
> - we figure out a solution in my lab and I will Run the SKearn algo from a cluster (so I hope no memory issues anymore). Because after that I will have a data set of 50 patients.

Excellent! Note that by using a cluster you can also take advantage of multiprocessing with the `--p-n-jobs` parameter. This will speed up time if you can harness multiple CPU on your cluster.

---

<div class="post-metadata">

### Author: ![antgonza](https://forum.qiime2.org/user_avatar/forum.qiime2.org/antgonza/32/1236_2.png) [@antgonza](https://forum.qiime2.org/u/antgonza)
#### Post date: [June 1, 2018, 2:29pm UTC](https://forum.qiime2.org/t/taxonomy-assignement-issues/4365/7 "2018-06-01T14:29:24Z")

</div>

Hi @BiOMan,

> [@BiOMan](#):
>
> greengenes is not appropriate as a database because it is used more for ecology

Just to add to @Nicholas_Bokulich reply. If you check [Extended Data Figure 2](https://www.nature.com/articles/nature24621/figures/6) from [A communal catalogue reveals Earth’s multiscale microbial diversity](https://www.nature.com/articles/nature24621); you will see in panel b a comparison of the 2 databases on how well they recover/match sequences when doing close reference; which could be used as a guideline to which database might be better for which kind of environment. Now, if you see panel c, you will see that Silva recruits more observed OTUs, while Greengenes more phylogenetic diversity; this is due to how these references are created; if you are interested, suggest reading more on how Silva and Greengenes are actually created.

---

<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: [July 2, 2018, 8:29pm UTC](https://forum.qiime2.org/t/taxonomy-assignement-issues/4365/8 "2018-07-02T20:29:27Z")

</div>

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