get.coremicrobiome

The get.coremicrobiome command determines the fraction of OTUs that are found in varying numbers of samples for different minimum relative abundances.

Default Settings

The get.coremicrobiom command reads a shared or relabund file and creates a *.core.microbiome file.

mothur > get.coremicrobiome(shared=temp.job2.shared)

or

mothur > get.coremicrobiome(relabund=temp.job2.relabund)

If you open the temp.job2.unique.core.microbiome file you will see:

NumSamples Relabund-0% Relabund-1% Relabund-2% Relabund-3% Relabund-4% ... Relabund-98%    Relabund-99%    Relabund-100%   
1  1.000000    0.615385    0.461538    0.307692    0.307692 ...    0.000000    0.000000    0.000000    
2  1.000000    0.538462    0.307692    0.230769    0.230769 ...    0.000000    0.000000    0.000000    
3  1.000000    0.384615    0.307692    0.230769    0.230769 ...    0.000000    0.000000    0.000000    
4  1.000000    0.384615    0.230769    0.230769    0.230769 ...    0.000000    0.000000    0.000000    
5  1.000000    0.384615    0.230769    0.230769    0.230769 ...    0.000000    0.000000    0.000000    
6  1.000000    0.230769    0.230769    0.230769    0.230769 ...    0.000000    0.000000    0.000000    
7  1.000000    0.230769    0.230769    0.230769    0.230769 ...    0.000000    0.000000    0.000000    
8  1.000000    0.230769    0.230769    0.230769    0.230769 ...    0.000000    0.000000    0.000000    
9  1.000000    0.230769    0.230769    0.230769    0.230769 ...    0.000000    0.000000    0.000000    
10 1.000000    0.230769    0.230769    0.230769    0.230769 ...    0.000000    0.000000    0.000000    
11 1.000000    0.230769    0.230769    0.230769    0.230769 ...    0.000000    0.000000    0.000000    
12 1.000000    0.230769    0.230769    0.230769    0.230769 ...    0.000000    0.000000    0.000000    
13 1.000000    0.230769    0.230769    0.230769    0.230769 ...    0.000000    0.000000    0.000000    
14 1.000000    0.230769    0.230769    0.230769    0.153846 ...    0.000000    0.000000    0.000000    
15 1.000000    0.230769    0.230769    0.153846    0.153846 ...    0.000000    0.000000    0.000000    
16 1.000000    0.230769    0.230769    0.153846    0.153846 ...    0.000000    0.000000    0.000000    
17 1.000000    0.230769    0.230769    0.153846    0.153846 ...    0.000000    0.000000    0.000000    
18 1.000000    0.230769    0.230769    0.153846    0.153846 ...    0.000000    0.000000    0.000000    
19 1.000000    0.230769    0.153846    0.153846    0.153846 ...    0.000000    0.000000    0.000000    

The entries in the table represents the fraction of OTUs in a study that have a relative abundance at or above %X in at least %Y samples. X goes from 0 to 100, Y from 1 to the number of samples in your dataset. From the above file we can see that 0.615385 (roughly 62%) of the OTUs have a relative abundance at or above 0.01 in at least one sample, but only

0.230769 (roughly 23%) have a relative abundance at or above 0.01 in all the samples.

Options

output

The output parameter is used to specify whether you would like the fraction of OTU’s or OTU count outputted. Options are fraction or count. Default=fraction.

mothur > get.coremicrobiome(shared=temp.job2.shared, output=count)

If you open the temp.job2.unique.core.microbiome file you will see:

NumSamples Relabund-0% Relabund-1% Relabund-2% Relabund-3% Relabund-4% ... Relabund-98%    Relabund-99%    Relabund-100%   
1  13  8   6   4   4 ... 0 0   0   
2  13  7   4   3   3 ... 0 0   0   
3  13  5   4   3   3 ... 0 0   0   
4  13  5   3   3   3 ... 0 0   0   
5  13  5   3   3   3 ... 0 0   0   
6  13  3   3   3   3 ... 0 0   0   
7  13  3   3   3   3 ... 0 0   0   
8  13  3   3   3   3 ... 0 0   0   
9  13  3   3   3   3 ... 0 0   0   
10 13  3   3   3   3 ... 0 0   0   
11 13  3   3   3   3 ... 0 0   0   
12 13  3   3   3   3 ... 0 0   0   
13 13  3   3   3   3 ... 0 0   0   
14 13  3   3   3   2 ... 0 0   0   
15 13  3   3   2   2 ... 0 0   0   
16 13  3   3   2   2 ... 0 0   0   
17 13  3   3   2   2 ... 0 0   0   
18 13  3   3   2   2 ... 0 0   0   
19 13  3   2   2   2 ... 0 0   0   

From the above file we can see that 8 of the OTUs have a relative abundance at or above 0.01 in at least one sample, but only 3 have a relative abundance at or above 0.01 in all the samples.

abundance && samples

Using the abundance and samples parameters will cause mothur to create a *.core.microbiomelist file. The abundance parameter allows you to specify an abundance you would like the OTU names outputted for. Must be an integer between 0 and 100, indicating the relative abundance. The samples parameter allows you to specify the minimum number of samples you would like the OTU names outputted for. Must be an interger between 1 and number of samples in your file.

To list the OTU’s that have a relative abundance at or above 0.01 at for each number of samples, you would run:

mothur > get.coremicrobiome(shared=temp.job2.shared, abundance=1)

If you open the temp.job2.unique.core.microbiomelist file you will see:

Samples    OTUList_for_abund=1
1  Otu02,Otu04,Otu06,Otu07,Otu08,Otu09,Otu12,Otu13
2  Otu02,Otu04,Otu07,Otu08,Otu09,Otu12,Otu13
3  Otu02,Otu07,Otu08,Otu09,Otu12
4  Otu02,Otu07,Otu08,Otu09,Otu12
5  Otu02,Otu07,Otu08,Otu09,Otu12
6  Otu07,Otu08,Otu09
7  Otu07,Otu08,Otu09
8  Otu07,Otu08,Otu09
9  Otu07,Otu08,Otu09
10 Otu07,Otu08,Otu09
11 Otu07,Otu08,Otu09
12 Otu07,Otu08,Otu09
13 Otu07,Otu08,Otu09
14 Otu07,Otu08,Otu09
15 Otu07,Otu08,Otu09
16 Otu07,Otu08,Otu09
17 Otu07,Otu08,Otu09
18 Otu07,Otu08,Otu09
19 Otu07,Otu08,Otu09

To list the OTU’s that have a given relative abundance for at least 8 samples, you would run:

mothur > get.coremicrobiome(shared=temp.job2.shared, samples=8)

If you open the temp.job2.unique.core.microbiomelist file you will see:

Relabund%  OTUList_for_samples=8
0  Otu01,Otu02,Otu03,Otu04,Otu05,Otu06,Otu07,Otu08,Otu09,Otu10,Otu11,Otu12,Otu13
1  Otu07,Otu08,Otu09
2  Otu07,Otu08,Otu09
3  Otu07,Otu08,Otu09
4  Otu07,Otu08,Otu09
5  Otu07,Otu08,Otu09
6  Otu07,Otu08,Otu09
7  Otu07,Otu08,Otu09
8  Otu07,Otu08,Otu09
9  Otu07,Otu08,Otu09
10 Otu07,Otu09
11 Otu07,Otu09
12 Otu07,Otu09
...    
98 
99 
100    

To list the OTU’s that have a relative abundance at or above 0.01 at for at least 8 samples, you would run:

mothur > get.coremicrobiome(shared=temp.job2.shared, abundance=1, samples=8)

If you open the temp.job2.unique.core.microbiomelist file you will see:

Relabund%  OTUList_for_samples=8
1  Otu07,Otu08,Otu09

groups

The groups parameter allows you to specify which of the groups in your shared or relabund file you would like included. The group names are separated by dashes.

label

The label parameter allows you to select what distance levels you would like, and are also separated by dashes.

Revisions

  • 1.26.0 - First Introduced
  • 1.40.0 - Speed and memory improvements for shared files. #357 , #347
  • 1.40.3 - Bug Fix: segfault with shared option. #445