当前位置: 首页 » NGS(二代测序了么) » Gene_prediction » 16SrRNA的预测与鉴定

16SrRNA的预测与鉴定

[文章目录] x

16S  rRNA作为微生物的特征片段,在鉴定微生物物种的过程中起着重要的作用。我们在做宏基因组分析的时候可以找一下scaffolds中的16S  rRNA,然后跟数据库比对,看里面都有什么物种;或者我们直接对环境样品进行16S  rRNA 高通量测序,然后来获知该环境样品的群落组成。总之16srRNA的分析很重要。本博文将围绕16S  rRNA的预测与鉴定两部分展开。

16s

一、预测

Rnammer、Infernal、Usearch 可以用来预测rRNA。

1.RNAmmer

网页版偶尔不好使,可以本地化该软件

RNAmmer本地化

1.1安装要求:
Unix系统,perl,hmmsearch(需要安装hmmer-2.2g版本,使用最新版本会出错)

1.2安装

安装hmmer-2.2g

安装RNAmmer
该软件需要使用edu邮箱去申请。下载地址:http://www.cbs.dtu.dk/cgi-bin/nph-sw_request?rnammer (需要学校的邮箱来认证)
下载,解压缩,然后进入rnammer修改部分地址, 一个是你安装Rnammer路径,一个是hmmsearch的路径,安装在哪里就改成哪里就可以了

检测是否正常运行

1.3. RNAmmer的运行与参数

对真核生物基因组进行rRNA预测的一个示例命令:

细菌的一条染色体,运行一下命令即可

1.4. RNAmmer的结果
rRNA.fasta,rRNA.gff2,rRNA.hmmreport,rRNA.xml。

2. infernal

官网http://infernal.janelia.org/

3.usearch寻找16S rRNA

当然,也可参照我之前的博客http://blog.sina.com.cn/s/blog_670445240101njwi.html

二.鉴定

预测出来的16S rRNA用RDP即可知道该rRNA是属于哪个物种的,由此我们推测该cluster是什么物种。

Rdp (ribosomal database project)用于对16S rRNA或者真菌LSU序列的分类地位上的鉴定。RDP naïve Bayesian rRNA Classifier的一种分层分类方法。
官网服务器:http://rdp.cme.msu.edu/classifier/classifier.jsp
在官网上将你的序列上传上去,如果是16S rRNA,就选择模式:16S rRNA training set 9,然后submit,一会的功夫就会出来结果。

两种模式:
16S rRNA training set 9
Fungal LSU

计算方式:
这个方法不需要比对。8个碱基一组(称为word),这样就有6万多种排列组合。对已知的880个属(现在数目应该大于这个)按照这个这种方式统计6万多种排列的比例情况。 上传的序列的组合数也求出来,然后随即选择word,跟参考库的来计算,(具体怎么计算,这个就是贝叶斯算法的核心罗,先不管),进行100次,100次计算都在同一个分类地位的数目,则为bootstrap cutoff。分类结果的后面都会跟上一个bootstrap cutoff。设定一个阈值,root的置信阈值一般选择的是80%。这个bootstrap cutoff又叫confidence estimate

上传序列的格式:
Fasta, GenBank and EMBL格式
长度至少为50个碱基,序列数要小于100000条。序列数超过100000条,联系
rdpstaff@msu.edu

置信阈值
通过100次抽样,分配到该分类地位的次数,则为boot cutoff ,极为最后结果分类地位后面跟着的那个百分数.
对于序列长度小于250bp,长于50bp.bootstrap cutoff 为50%则认为可以反映序列属的水平

结果:
序列后面跟着的“-”意味是根据序列反向来得出的结果。
根据置信阈值,如果小于50%,则在前面的分类地位前加上unclassfied,写了一个脚本。
allrank_rep_otu_nosingle.fasta_classified.txt为RDP反馈的结果,有需求脚本rdpformat.pl留下邮箱。

PS:
这里面还有一点我不是太确定,所有的分类地位后面都有一个bootstrap cutoff,大于50%则是一类,小于50%则是unclassfied,是这样吗? 我怎么感觉仅仅是针对属来说的呢?在这里留下一个悬念。


参考资料:
http://rdp.cme.msu.edu/classifier/class_help.jsp#model

RNAmmer的安装和使用


http://www.cbs.dtu.dk/services/doc/rnammer-1.2.readme


下一篇 :

上一篇 :

暂无评论

发表评论

电子邮件地址不会被公开。 必填项已用*标注

$(document).ready(function(){ $("#article-index").css('display','none');});