PICRUSt2

PICRUSt2(Phylogenetic Investigation of Communities by Reconstruction of Unobserved States)是一款基于标记基因序列来预测功能丰度的软件。“功能”通常指的是基因家族,如KEGG同源基因和酶分类号,但可以预测任何一个任意的特性。同样,预测通常基于16S rRNA基因测序数据,但也可以使用其他标记基因.

PICRUSt2相关使用方法来啦!!!

一、安装

#使用 bioconda 安装 PICRUSt2 环境

conda create -n picrust2 -c bioconda -c conda-forge picrust2

#激活环境

source activate picrust2

#退出环境

source deactivate picrust2

二、准备数据

PICRUSt2无需再以GreenGene注释的OTU表为输入,可以直接读取OTU的代表序列自动完成物种注释,并进一步根据物种丰度组成预测群落功能。文件otu_table.txt为OTU丰度表格,仅包含丰度组成信息即可,无需添加注释列。文件otu.fasta中包含了OTU丰度表中各OTU的代表序列。

三、使用

picrust2_pipeline.py -s otu.fasta -i otu_table.txt -o picrust2_result -p 4

四、结果文件

1)out.tre,所有OTU代表序列构建的系统发育树文件。
2)KO_metagenome_out/,该结果路径中记录了细菌群落KO(KEGG Orthology)功能的丰度预测结果。
3)KO_metagenome_out/seqtab_norm.tsv.gz,对于很多细菌而言,一个个体可能包含多条16S(多拷贝16S),因此在原始OTU 16S rRNA丰度表的基础上,根据物种所含16S rRNA拷贝数对物种丰度进行标准化,得到校正16S rRNA拷贝数后的OTU丰度表。
4)KO_metagenome_out/pred_metagenome_unstrat.tsv.gz,该文件中即为预测得到的细菌群落功能丰度表,记录了各样本中所包含KO功能的丰度,丰度计算由上述校正16S rRNA拷贝数标准化后的OTU丰度表推断得到。功能以KO ID为名称,代表了特定的功能基因。
5)KO_metagenome_out/weighted_nsti.tsv.gz,各样本预测功能的加权NSTI值,由OTU的NSTI值通过标准化后的丰度加权所得。
6)EC_metagenome_out/,该结果路径中记录了细菌群落酶(EC)功能的丰度预测结果。文件结构同上述KO_metagenome_out/,不再展示。
7)pathways_out/path_abun_unstrat.tsv.gz,上述为预测得到的以KO ID为名称的KO功能,实则代表了特定的功能基因,将这些功能基因映射到具体的KEGG代谢途径(KEGG pathway)中,并统计各途径在各样本中的丰度,获得该表。
8)KO_predicted.tsv.gz和EC_predicted.tsv.gz,两个矩阵文件中记录了OTU对预测功能丰度的贡献,即可以理解为每个OTU所代表的物种个体基因组中,分别有多少数量的基因与对应的KO功能或酶功能有关。如果期望关注哪些OTU是否对群落功能是重要的,这些表格(该表仅代表了单个物种个体基因组的特征,可能还需结合OTU的丰度信息)可以提供参考。
9)marker_predicted_and_nsti.tsv.gz,记录了OTU代表物种基因组中,16SrRNA拷贝数以及功能预测的NSTI值信息。
10)Intermediate/,一些中间文件。

END

PICRUSt2软件相关推荐

  1. QIIME2得到PICRUSt2结果后如何分析

    作者:聋言瞎面 审核:Listenlii 本文主要解决的问题: 1.QIIME2做完PICRUSt2后,只输出pathway id,如何比对得到pathway description? 2.PICRU ...

  2. PICRUSt2分析实战:16S扩增子OTU或ASV预测宏基因组、新增KEGG层级

    PICRUSt2分析实战:16S扩增子OTU或ASV预测宏基因组.新增KEGG层级 更新时间:2021年7月8日 PICRUSt推出了近8年,引用5000余次. 现推出PICRUSt2,202年再次霸 ...

  3. PICRUSt2分析实战:16S扩增子OTU或ASV预测宏基因组EC、通路、KO(200806更新)

    PICRUSt2分析实战:16S扩增子OTU或ASV预测宏基因组 更新时间:2020年8月6日 PICRUSt推出了近7年,引用4000余次. 现推出PICRUSt2,再次霸气发表于顶级期刊Natur ...

  4. NBT:PICRUSt2预测宏基因组功能

    PICRUSt2预测宏基因组功能 PICRUSt2 for prediction of metagenome functions Nature Biotechnology [IF:31.864] 20 ...

  5. PICRUSt2:OTU或ASV等16S随便预测宏基因组且数据库增加10倍

    文章目录 简介 工作流程 引用 安装 一条命令完成分析 核心输出结果: 额外输出文件: 参数详解 猜你喜欢 写在后面 PICRUSt推出了近6年,引用2500余次. 现推出PICRUSt2 https ...

  6. 16S的细菌群落功能预测工具PICRUSt2学习

    2018年推出了全新版本的PICRUSt,即PICRUSt2(https://github.com/picrust/picrust2) PICRUSt2 (Phylogenetic Investiga ...

  7. picrust2功能预测-从qiime2安装到数据分析

    记录一下跌跌撞撞地摸索过程 picrust2功能预测 参照公众号<宏基因组>刘永鑫的原创:<PICRUSt2:OTU或ASV等16S随便预测宏基因组且数据库增加10倍>http ...

  8. PICRUSt2:OTU/ASV等16S序列随意预测宏基因组,参考数据库增大10倍

    PICRUSt推出了近6年,引用2500余次. 现推出PICRUSt2 https://github.com/picrust/picrust2 具有以下三大优势: 任何OTU/ASV直接预测功能: 数 ...

  9. picrust2基于ubuntu的安装

    安装源方式安装 参考于Installation · picrust/picrust2 Wiki · GitHub wget https://github.com/picrust/picrust2/ar ...

最新文章

  1. 写文章 TEE技术分析【转】
  2. 可以作为python字符串常量的是_Python 复习
  3. java内存溢出怎么排查_【转】Java学习---内存溢出的排查经历
  4. SQlite-数据库的访问实例(转)
  5. 一看就懂,图解 Spring 循环依赖,写得老好了!
  6. Ubuntu 下连接Github
  7. 深究angularJS——(上传)FileUploader中文翻译
  8. shell脚本执行php文件_分享两个shell脚本实例--批量生成随机字符文件名和批量改名...
  9. clob字段怎么导出_Oracle 11g及12c+版本下为啥有些表不能exp导出?
  10. java 提取url参数_Java提取URL某个参数的值
  11. java接口自动化(四) - 企业级代码管理工具Git的应用
  12. 九度oj题目1153:括号匹配问题
  13. 刚刚看到 PNaCl, 这才是我一直期待的跨平台的好东西!
  14. 服装计算机辅助设计论文,服装设计计算机辅助设计应用探索-服装设计论文-设计论文.docx...
  15. Nim理论初探——编程之美1.12
  16. 微信小程序 开发工具 项目创建步骤
  17. 数据库中查出的树形,获取其所有子孙(递归方式)
  18. Excel使用---excel2016___一般操作(搬,侵删)
  19. 诚立业科技面试经历.
  20. 微信服务号开发之jssdk-微信充值、支付见解

热门文章

  1. 基于Hadoop的58同城离线计算平台设计与实践
  2. 别再问我 new 字符串创建了几个对象了!我来证明给你看!
  3. 自己动手写一个能操作redis的客户端
  4. 创业公司崛起靠加班?别闹了,早死才靠加班
  5. 为什么大公司还在采用过时的技术?
  6. SpringBootH ttpInvoker接口调用
  7. 分布式存储系统考虑因素-一致性
  8. 嵌入式系统学习笔记之五-uboot常用命令之补充
  9. Ubuntu16.04下安装mysql5.7数据库
  10. 数字加密c语言程序_大厂程序员整理的 C++ 资源大全,不私藏了,都送给你们吧...