• 扩增子分析神器USEARCH 简介 v11新功能 v11命令大全
• 扩增子分析神器VSEARCH 分析流程 2.8.1中文帮助文档

otutab_core

http://www.drive5.com/usearch/manual/cmd_otutab_core.html  

鉴定核心微生物组——大多数样品中出现的OTUs，这也是Usearch11新增的功能。

本质上是统计每个OTUs在大量样品中出现的频率。比如在所有样本中都出现，即100%。特别大量样本时，如几千个，可能很少有OTU存在于所有样本，可选择95%，或90%的样本中出现的OTUs作为核心OTUs。

Identifies a possible “core microbiome” of OTUs which are present in more samples than others.

输入文件为经典的QIIME格式OTU表

Input is an OTU table in QIIME classic format.

在一些样品或许多样本中出现的OTU可能是假的，原因可能是串道(cross-talk)或扩增、测序错误的假OTU。为方便检查，otutab_core命令产生分析报告，标明哪些OTU可能是串道、哪些可能是测序错误。

The presence of an OTU in some or many samples can be spurious because of cross-talk or because the OTU itself is spurious. To enable manual review, the otutab_core command generates a report indicating cases where the presence of an OTU may be spurious due to cross-talk, and where an OTU may be spurious due to sequence errors.

如果使用-sintaxin参数指定物种注释，报导中也会包括物种信息。

If a sintax tabbed file is provided using the -sintaxin option, then the taxonomy of the core OTUs is included in the report.

如果使用distmxin选项提供距离矩阵，可用于鉴定主导的OTUs，例如在报告中高丰度OTUs与低丰度OTUs相似，如果存在主导OTUs，那这些低丰主OTUs是假的。

If a distance matrix is provided using the distmxin option, this is used to identify possible dominant OTUs, i.e. high-abundance OTUs which are similar to a low-abundance OTUs in the report. If there is a dominant OTU, this may indicate that the low-abundance OTU is spurious.

-tabbedout参数指定输出文件。OTUs按样品中出现频率排序，包括如下12个字段。

The -tabbedout option specifies the output file. OTUs are sorted in order of decreasing number of samples where they are present. Fields are:

1. OTU = name of the OTU.
2. Samples = number of samples where the OTU has a non-zero count.
3. Size = total number of reads assigned to this OTU.
4. DomOTU = high-abundance “dominant” OTU which is very similar to this OTU, if any.
5. DomSize = total number of reads assigned to the dominant OTU.
6. DomId = identity of the dominant OTU with this OTU.
7. Min = minimum count for this OTU.
8. LoQ = low quartile count for this OTU.
9. Med = median count for this OTU.
10. HiQ = high quartile count for this OTU.
11. Max = maximum count for this OTU.
12. Taxonomy = condensed taxonomy prediction.

If the minimum or LoQ count is much smaller than the maximum count, this suggests that the smaller counts may be due to cross-talk.

If the size of an OTU is much smaller than a neighboring “dominant” OTU, then the OTU itself may be spurious due to sequence error.

使用实例

Example

1. 基于OTU序列计算距离矩阵

usearch -calc_distmx otus.fa -tabbedout distmx.txt \
  -sparsemx_minid 0.9 -termid 0.8

1. 物种注释(己完成可跳过)

usearch -sintax otus.fa -strand both -db ref16s.txt \
  -tabbedout sintax.txt

1. 鉴定核心OTUs

usearch -otutab_core otutab.txt -distmxin distmx.txt \
  -sintaxin sintax.txt -tabbedout core.txt

在使用中，我碰到了报错。可以把-distmxin distmx.txt去掉。可以正常获得结果。

结果文件如下：

OTUID   Samples Size    Freq    DomOTU  DomSize DomId   Min     LoQ     Med     HiQ     Max     Taxonomy        Core
OTU_2   1000    5079019 0.131   .       .       .       162     1915    3270    5470    23217   d:Bacteria,p:"Proteobacteria",c:Betaproteobacteria,o:Burkholderiales,f:Burkholderiaceae,g:Ralstonia,s:Ralstonia_mannitolilytica 100
OTU_34  999    180434  0.00466 .       .       .       1       40      83      174     2484    d:Bacteria,p:"Proteobacteria",c:Gammaproteobacteria,o:Clostridiales,f:Chloroplast,g:Streptophyta,s:Porticoccus_litoralis        99.9154

具体每列的意义见上方帮助文档。最主要的结果是Samples列，即该OTU在多少个样本中检测到。我们还需要将此数值除以总样本量，才能获得Core OTU的比例，方便筛选核心OTUs。

猜你喜欢

• 10000+: 菌群分析
  宝宝与猫狗 提DNA发Nature 实验分析谁对结果影响大  Cell微生物专刊 肠道指挥大脑
• 系列教程：微生物组入门 Biostar 微生物组  宏基因组
• 专业技能：生信宝典 学术图表 高分文章 不可或缺的人
• 一文读懂：宏基因组 寄生虫益处 进化树
• 必备技能：提问 搜索  Endnote
• 文献阅读 热心肠 SemanticScholar Geenmedical
• 扩增子分析：图表解读 分析流程 统计绘图
• 16S功能预测   PICRUSt  FAPROTAX  Bugbase Tax4Fun
• 在线工具：16S预测培养基 生信绘图
• 科研经验：云笔记  云协作 公众号
• 编程模板: Shell  R Perl
• 生物科普: 肠道细菌 人体上的生命 生命大跃进  细胞暗战 人体奥秘  

写在后面

为鼓励读者交流、快速解决科研困难，我们建立了“宏基因组”专业讨论群，目前己有国内外1800+ 一线科研人员加入。参与讨论，获得专业解答，欢迎分享此文至朋友圈，并扫码加主编好友带你入群，务必备注“姓名-单位-研究方向-职称/年级”。技术问题寻求帮助，首先阅读《如何优雅的提问》学习解决问题思路，仍末解决群内讨论，问题不私聊，帮助同行。

学习扩增子、宏基因组科研思路和分析实战，关注“宏基因组”

点击阅读原文，跳转最新文章目录阅读
https://mp.weixin.qq.com/s/5jQspEvH5_4Xmart22gjMA