宏基因组生信分析方法介绍

news/2024/11/8 0:51:42/

随着高通量技术的发展,宏基因组学(metagenomics)已经成为研究微生物群落物种及功能的前沿科学,在肠道微生物、环境微生物等研究领域具有广泛应用。宏基因组学通过对微生物群落全部DNA进行高通量测序,将测序序列与公共数据库进行比对或从头组装出微生物基因组,从而识别微生物群落的物种和功能基因

目前主流的宏基因组数据分析方法包括三种:基于测序结果进行组装的分析方法;基于reads直接和已知数据库进行比对的方法;基于分箱bin的方法。

常用到的分析工具见下图:

 

一、基于组装

数据分析从下机原始序列开始,首先对原始序列进行去接头、 质量剪切以及去除污染等优化处理。然后使用优质序列进行拼接组装和基因预测,将各样本预测得到的基因集合并在一起去冗余,得到非冗余基因集;对得到的非冗余基因集与NCBI的nr数据库进行比对,得到物种信息,与各个功能数据库(抗性基因ARDB/CARD、碳水化合物CAZy数据库、KEGG数据库,eggnog数据库、VFDB数据库)进行比对,得到不同功能的注释结果。并使用BWA软件将优化序列比对到非冗余基因集,计算得到各基因在各样品中的丰度信息(RPKM); 最后对物种和功能注释结果进行统计及后继分析。

具体分析的流程见下图:

 

二、基于Reads

基于Reads的宏基因组分析,使用质控好的优化序列与已知物种和功能数据库进行比对,从而得到每个样本的物种、功能注释信息和相对丰度。此方法最大的优势是不用进行基因组的组装,从而快速的得到群落的物种与功能组成信息。

物种注释:从质控好的优化序列出发,使用Kranken2软件和微生物数据库(包括RefSeq最近的99版本基因组数据,涵盖细菌、真菌、病毒、古菌、原生动物)鉴别样本中所含的物种,再用Bracken对Kraken2得到的分类结果进性分类后贝叶斯重新估算丰度来估算宏基因组样本的种级别相对丰度;

功能注释:从质控好的优化序列出发,使用HUMAnN3软件和蛋⽩质数据库(UniRef90)进⾏⽐对,根据UniRef90 的ID 和各个功能数据库ID的对应关系,统计各个功能数据库对应功能相对丰度;

其他功能数据库:使⽤FMAP软件将各样本质控好的优化序列与参考数据库进⾏⽐对(基于DIAMOND),根据⽐对结果,统计出每个样本⽐对到各参考序列的reads数,从⽽计算相对丰度;

具体分析流程见下图:

三、基于Bin分析

宏基因组分箱(Binning)是将序列组装得到的Contigs按物种分开归类的过程。基于Binning的宏基因组分析流程,数据分析从下机原始序列开始,首先对原始序列进行去接头、 质量剪切以及去除污染等优化处理。然后使用优质序列进行拼接组装得到Contigs;使用metabat2和maxbin2软件分别对每个样本的Contigs进行分箱;不同软件分箱得到的bin进行合并(binning_refiner)、提纯(MAGpurify);然后将所有样优化后的bin进行去冗余(dRep);而后分别从物种和功能方面进行信息统计。

具体分析流程见下图:

四、三种方法各自的优缺点:

1、基于组装的方法可以得到各个ORF的丰度信息,不管该ORF是否有对应的功能或是物种信息。如果是环境样本,微生物非常复杂,里面会得到的ORF信息量大,导致分析需要的时间延长。

2、基于Reads与已知序列Mapping的方法,没办法表征大量尚未分离和测序的微生物群体。现有的比对算法只能发现近缘关系,且仅能用于分析数据库中与基因组最密切相关的环境微生物。

但是基于基于Reads与已知序列Mapping的方法可以快速的得到物种和功能分析的结果。

3、Binning:Metagenome 组装完成后,我们得到的是成千上万的 contigs,我们需要知道哪些 contigs 来自哪一个基因组,或者都有哪些微生物的基因组。所以需要将 contigs 按照物种水平进行分组归类,称为 "bining"。主要依据是:来自同一菌株的序列,其核酸组成是相似的。优势:通过宏基因组测序尽可能完整的组装出样本中菌株的基因组。


http://www.ppmy.cn/news/301324.html

相关文章

宏基因组分析步骤Linux,宏基因组分析专题研讨班

计算中心生物部一直以来都以最前沿的课程、理论与实践相结合和高品质的服务为广大科研工作者提供全面的生物信息课程。 为了解决科研工作者在宏基因组样品制备、数据分析、功能基因挖掘等的困扰,生物部特推出《宏基因组分析专题研讨班》。专业的讲课老师团队与您一起…

宏基笔记本linux系统安装图解,宏基笔记本重装系统详细图文教程

大家都知道宏基笔记本重装系统是比较受欢迎的,因为它的用户很广泛,那么宏基笔记本重装系统到底要如何去操作呢?这里面就很有技巧了,接下来小编就要给大家带来宏基笔记本重装系统详细图文教程,手把手教大家如何给自己的宏碁笔记本…

多少年前现在依旧爆火的游戏——星际战甲

《星际战甲》是Digital Extremes研发的一款科幻题材的第三人称射击网游,于2015年9月25日在中国公测发行 。 游戏采用以战甲代替职业的设定,玩家选择不同的战甲,每个战甲的攻击方式和技能都不同。 《星际战甲》Switch版2018年11月21日正式发…

微软xp壁纸rgb

微软xp壁纸rgb,打开截图就可以使用了哦&#xff01; <!DOCTYPE html> <html><head><meta charset"UTF-8"><title></title></head><body><div id"test1" style"background-color:rgb(0,78,152);w…

Ubuntu系统下玩星际争霸

我在windows首先安装了星际争霸&#xff0c;再安装一个免CD的补丁&#xff0c;然后上battle,net升级到1.14 在Ubuntu系统下的新利得软件管理器搜索到wine这个软件&#xff0c;选择安装&#xff0c;执行以后就会自动安装好了 在media目录下找到星际争霸的目录&#xff0c;并且…

关于星际争霸的

神族真正的进攻主力: ①金甲(攻击效果惊人,但是速度太慢,须配以运输机(为战时操作简便,一般用小叮当跟住运输机,这样的话运输机就承担了辅助进攻和反隐形的双重任务))、 ②白球(进攻的时候将电兵混编在大部队之中,撒电将魔法用完以后就合成白球.白球相当于虫族的大象,攻高血厚…

星际争霸:人族(terran)魔法与科技

兴奋剂消耗&#xff1a; 10 生命 机枪兵和火焰兵可以使用。兴奋剂是一种合成肾上腺素和有止痛作用的药丸&#xff0c;它可以让火焰兵和机枪兵移动更快&#xff0c;开火速度也更快。其缺点就是&#xff0c;使用了兴奋剂后就会对部队立即造成10点的伤害&#xff0c;这个伤害无法…

030 星际争霸网页版!!

星际争霸网页版&#xff01;&#xff01; 2007-08-18 09:10:45| 分类&#xff1a; 月夜沉思处——感 题记&#xff1a;比较好玩的一款WEB游戏&#xff0c;希望自己在玩游戏之余&#xff0c;还能收获些别的东西。 黑客反病毒论坛 http://bbs.hackav.com 这是冰点的大哥们写的一…