在物种水平上的宏基因组比对分析流程

文章链接:http://biorxiv.org/content/early/2016/10/15/081141

作者:Yee Voan Teo, Nicola Neretti

时间:2016.10.15

摘要:

许多宏基因组分类工具在宏基因组学领域增长速度飞速。然而,在这个领域相近物种的分类仍然是一个挑战。这里,我们用两个宏基因组数据集,人类宏基因组数据和环境宏基因组数据对比MetaPhlAn2, kallisto 和 Kraken的性能。研究表明kallisto比MetaPhlAn2 和Kraken有更高的敏感度,并且在物种水平上的分类精确度也比kraken高。我们的研究还表明在运行全基因组的分类工具会将不存在的物种分类到存在的物种上。为了减少这种FP,在我们的流程中MetaPhlAn2我们介绍marker基因,它使用Kallisto为分类步骤,作为用于物种检测一个附加的过滤步骤。

时间: 2024-10-05 07:40:06

在物种水平上的宏基因组比对分析流程的相关文章

宏基因组测序及分析

宏基因组测序: 浓度>=50ng/ul OD260/280:1.8-2.0 DNA两次需要量>=3ug 宏基因组测序需要提供什么样品要求? (1)提供环境微生物的基因组DNA或者扩增产物,OD值在1.8-2.0 之间:样品浓度大于30 ng/ul:每次样品制备需要10μg样品,如果需要多次制备样品,则需要样品总量=制备样品次数*10 ug. (2)DNA样品请置于-20℃保存:请提供DNA样品具体浓度.体积.制备时间.溶剂名称.请同时附上QC数据,包括电泳胶图.分光光度或Nanodrop仪器检

一个全基因组重测序分析实战

Original 2017-06-08 曾健明 生信技能树 这里选取的是 GATK best practice 是目前认可度最高的全基因组重测序分析流程,尤其适用于 人类研究. PS:其实本文应该属于直播我的基因组系列,有两个原因把它单独拿出来, 首先,直播我的基因组阅读量太低了,可能是大家觉得错过了前面的,后面的看起来没有必要,这里我可以肯定的告诉大家,这一讲是独立的,而且是全流程,你学好了这个,整个直播我的基因组就可以不用看了. 其次,最近有一些朋友写了一些GATK的教程,但是大多不合我意,

宏基因组扩增子图表解读2散点图:组间整体差异分析(Beta多样性)

散点图 数据点在直角坐标系平面上的分布图.在宏基因组领域,散点图常用于展示样品组间的Beta多样性,常用的分析方法有主成分分析(PCA),主坐标轴分析(PCoA/MDS)和限制条件的主坐标轴分析(CPCoA/CCA/RDA). Beta多样性 Beat多样性是生态学概念,专指不同组或生态位间物种组成的差异. 分析方法 在读文章中经常可以看到PCA分析.PCoA分析,NMDS分析,CCA分析,RDA分析.它们在本质上是排序(ordination)分析.排序的过程就是在一个可视化的低维空间(通常是二

zabbix上的宏(macro)介绍

宏:macro,预设的文本替换模式: 宏是一种抽象概念(Abstraction),它根据一些列预定义的规则替换一定的文本模式,而解释或编译器在遇到宏时会自动进行这一模式替换.类似地,zabbix基于宏保存预设文本模式,并且在调用时将其替换为其中的文本. 内置{MACRO}    不需要加$就可以引用: zabbix有许多内置的宏,如{HOST.NAME}.{HOST.IP}.{TRIGGER.DESCRIPTION}.{TRIGGER.NAME}.{TRIGGER.EVENTS.ACK}等. 关

线代笔记 #01# 几何水平上的理解 VS. 数值水平上的理解

源: 线性代数的本质 内容来自字幕. Let me quote without further comment from Dieudonné's "Foundations of Modern Analysis, Vol. 1" There is hardly any theory which is more elementary than linear algebra, in spite of the fact that generations of professors and te

Genome-wide Complex Trait Analysis(GCTA)-全基因组复杂性状分析

GCTA(全基因组复杂性状分析)工具开发目的是针对复杂性状的全基因组关联分析,评估SNP解释的表型方差所占的比例(该网站地址:http://cnsgenomics.com/software/gcta/).目前GCTA工具可实现以下功能: 1 评估全基因组SNP的亲缘关系(遗传关系) 2 评估全基因组SNP的近交系数 3 评估所有的常染色体SNP对于变异的解释度 4 评估遗传方差与X-染色体的关联 5 检测遗传方差对X-染色体的剂量补偿效应 6 预测单个个体和单个SNP的全基因组加性遗传效应 7

vsftpd开启日志记录上传、下载、删除,分析xferlog日志

在web服务器上搭建了vsftpd用来上传代码程序,开启日志记录客户端的上传.下载删除等操作: #修改/etc/vsftpd/vsftpd.conf [[email protected] ~]# vim /etc/vsftpd/vsftpd.conf anonymous_enable=NO anon_root=/ftpanon local_enable=YES write_enable=YES local_umask=022 dirmessage_enable=YES dual_log_enab

【宏基因组】MEGAN4,MEGAN5和MEGAN6的Linux安装和使用

MEGAN(Metagenome Analyzer)是宏基因组学进行物种和功能研究的常用软件,实际上现在的Diamond+MEGAN6已经是一套比较完整的物种和功能注释流程了. 但是由于各种原因,我们在流程中使用的并非最新版.不同版本的MEGAN使用方法差别较大,尤其在命令行模式下.网上的关于这方面的资料也寥寥无几,这里简单记录下使用方法,主要是针对Linux平台的使用. MEGAN的GUI版相对友好,如果你在Windows平台使用过该软件,那么在Linux上使用和理解起来相对容易些. MEGA

怎样在Android.mk上加宏定义【转】

本文转载自:http://blog.csdn.net/ttxgz/article/details/7591282 很简单, LOCAL_CFLAGS += -DWHATEVERDEFINE 就可以了