Genera shared between groups

47 views
Skip to first unread message

Camila Fontoura

unread,
Mar 20, 2019, 4:33:40 PM3/20/19
to Qiime 1 Forum


Hi

I would like to get genera results shared beteen groups, rather than OTUs, and then create a venn diagram. Can anyone help me on what scripts to use to get this data?

Colin Brislawn

unread,
Mar 20, 2019, 7:24:35 PM3/20/19
to Qiime 1 Forum
Hello Camila,

You can groups your OTUs into genera using this script:

Then, to find out how many genera are shared, you can use this script:

Qiime 1 does not include a way to make Venn diagrams, but you can make them using other programs. Maybe this one?

Colin

Camila Fontoura

unread,
Mar 21, 2019, 10:18:18 AM3/21/19
to Qiime 1 Forum
Hello

Thank you for your help! I like this Venn diagram maker, I usually use venny, but it does have some limitations, for example it does not have to choose the color etc.

I've grouped the otus for pick_closed_reference.py and to know which OTUs are shared between the groups I used the group_significance.py script. Is correct?

But I would like to cite genera, not OTUs. Do you have a script that does this?

One thing I doubted is: when I make this comparison of shared OTUs between groups, many OTUs are unique. However, when I do alpha diversity statistics, there is no significant difference between the groups. Do you know any explanation?

Colin Brislawn

unread,
Mar 21, 2019, 12:18:06 PM3/21/19
to Qiime 1 Forum
Hello again,

The group_significance.py script compares the abundances of OTUs that found in both groups, while shared_phylotypes.py counts the number of OTUs found in both groups.

But I would like to cite genera, not OTUs. Do you have a script that does this?
All these scripts count OTUs... unless you use summarize_taxa.py to turn your OTU table into a Genera table. Then all these scripts count genera! :-)  

One thing I doubted is: when I make this comparison of shared OTUs between groups, many OTUs are unique. However, when I do alpha diversity statistics, there is no significant difference between the groups.Do you know any explanation?
Yes! What you are seeing is that the alpha diversity between groups is the same, but the beta diversity is different. 

Great questions Camila!

Colin

Camila Fontoura

unread,
Mar 21, 2019, 2:55:12 PM3/21/19
to Qiime 1 Forum
Thank you for the explanation. Now I understand the differences.

But for example: in the shannon index, the higher the index the more uniformity of species? or are there more dominant species than others?

And the other indexes. Chao1, simpsom, observed_otus? Do they all show the same thing?

I ran shared_phylotypes.py but it generates results per sample. I would like to get this result by group. How do I do? Can you help me? And in the exit file it does not generate the taxonomic classification, should it not contain this information?

Colin Brislawn

unread,
Mar 22, 2019, 12:48:59 PM3/22/19
to Qiime 1 Forum
Hello Camila,

All the alpha diversity metrics measure slightly different things, in slightly different ways. Read more about that here: 

I would like to get this result by group.
OK. Then you first have to merge your samples based on group.

Colin

Camila Fontoura

unread,
Mar 26, 2019, 9:21:48 AM3/26/19
to Qiime 1 Forum
Hi Colin,

Thanks!

I merged the samples based on the groups and then make the shared_phylotypes.py. But I'm having trouble filling the diagram. In the diagram, there is a field
Set 2 ^ Set 3 ^ Set 4 but I do not know how to put this information. My output file was:


Column1 SDE FDS FDE SDS
SDE 1459 1119 891 1163
FDS 1119 1392 892 1104
FDE 891 892 1051 867
SDS 1163 1104 867 1448

Can you help me?

Colin Brislawn

unread,
Mar 26, 2019, 3:08:58 PM3/26/19
to Qiime 1 Forum
Good afternoon,

I'm glad you got shared_phylotypes.py working! That results looks promising. 

I'm not sure how best to calculate these shared groups you described. It looks like they are not included in the shared_phylotypes.py script. :-(

Maybe you could color-code this table and present it as a heatmap? I know it's not the Venn diagram you were hoping for, but I think it's still a valid way to present this result. 

Colin

hg.maryo

unread,
Jul 10, 2019, 11:46:17 PM7/10/19
to Qiime 1 Forum
Collin,

Thanks for the thread.
Is this the valid valid way you meant to represent the Camila's shared genera numbers data instead a Venn diagram?

What do you think?


Camila's_data.jpeg


Mario Hg.

Colin Brislawn

unread,
Jul 11, 2019, 11:40:35 AM7/11/19
to Qiime 1 Forum
Hello,

Yes, that's what I had in mind. This looks good! 
(My only advice is to make sure the labels are in the same order on both sides that way the heatmap is symmetrical.) 

Colin

Camila Fontoura

unread,
Aug 5, 2019, 3:05:31 PM8/5/19
to Qiime 1 Forum
Hi Collin and Mario,

Thanks for your help!

How did you make this chart? How would I generate symmetrically?

Another thing....I would like to generate a heatmap with several variables. As the attached picture. Can you help me please?
heatmap_ok.png

Colin Brislawn

unread,
Aug 5, 2019, 3:20:06 PM8/5/19
to Qiime 1 Forum
Hello again,

I think that plot was made outside of Qiime 1, maybe using a program like R. 

Here is a tutorial on how to make heatmaps using R and phyloseq:

Colin

hg.maryo

unread,
Nov 13, 2019, 12:02:52 AM11/13/19
to Qiime 1 Forum
Thanks Collin,

I was able to make it better now.




Dear Camila,


To plot the later figure I used the pheatmap package in R.
It has a lot of documentation to follow and you can custom your heatmap almost as you wish.





Saludos.
Reply all
Reply to author
Forward
0 new messages