| 查看: 1364 | 回复: 0 | |||
[求助]
构建进化树相关问题，求助，谢谢大家
|
|
本人刚接触分子生物信息学这一块，最近刚开始学习做系统发育进化树，然后在做的过程中出现了很多的问题，现在说出来，希望做过的大神能够给一些意见和建议。 我是测的高通量测序，做功能菌AOAAOB，因为功能菌基因的注释结果不是很好，所以想要做一个系统发育的进化树，来看看自己得到的OUT到底是属于哪一个属的，我的操作就是拿代表OUT的序列在BLAST上比对，然后勾选uncultured，会得到一些有具体菌种名称的序列，只是identity没有100%，可能就是92%，94%之类的，然后QUERY COVER也是有可能不会是100%，就99%，98%，97%这样子的，就拿AOB举例子来说，我取丰度最高的前十代表OTU来按照上面的操作搜索参比序列，发现前十代表OTU的序列本身就比较相近，可能是土壤里面本身的菌落结构比较单一，最后得到了一个比较粗糙的进化树。这个树太粗糙以至于有些序列放了几遍我也开始没发现。只是举个例子而已。 现在存在的问题如下： 1.我的代表OTU序列本身就比较相近，所以在进化树里面OTU自己都聚在了一起，而已经的参比序列聚到了一起，而且尽管相似度都很高，但是自展值有的却很低，有的也就是三十几，四十几，这完全没有达到我的目的，看不出哪个OTU跟哪个菌比较近，不知道这种情况应该怎么办？怎么调整？ 2.我用blast找参比序列的时候是勾选了unclutured的，就是不知道在进化树能不能放入uncultured的菌属，像我上面描述这样子的操作到底可不可行？ 3.高通量序列测出来的有很多代表OTU序列，不知道是应该放丰度前1%，还是放丰度排名前30位的OTU序列？ 4.一般都是会用亲缘关系比较远的序列作为根，亲缘关系比较远的序列如何找，如何确定？感觉自己找亲缘关系远的或者近的都比较盲目。 5.有的文章中会在后面加上是哪些类，比如是nitrosospira,nitrosomonas等，甚至有的是Group Ⅰ.1之类的，就像知道如果序列是uncultured，那如何知道它是属于哪个属的？ 自己搞了一周了， 还不解，希望各位帮帮忙  屏幕快照 2017-11-24 下午7.27.25.png  屏幕快照 2017-11-30 下午2.47.22.png  屏幕快照 2017-11-30 下午2.47.32.png  屏幕快照 2017-11-30 下午2.48.00.png  1D8A9823-7B19-44C0-A928-44E047C560F0.png |
回复此楼» 猜你喜欢
含氟糖皮质激素的外用抗炎利器醋酸氟轻松
已经有0人回复
-
阿莫奈韦（Amenamevir）：解旋酶-引物酶靶向分子的结构特征与合成路径全解析
已经有0人回复
-
化学工程及工业化学论文润色/翻译怎么收费?
已经有265人回复
-
利瑞鲁单抗：靶向KIR的NK细胞检查点工具抗体，结构基础与科研应用全解析
已经有0人回复
-
Wnt-C59：PORCN靶向化学探针，Wnt信号通路阻断机制与科研应用全解析
已经有0人回复
-
Fexaramine：肠道限制性FXR化学探针，结构特征与信号通路研究全解析
已经有0人回复
-
NVS-STG2：靶向STING跨膜域的分子胶类激动剂，结构特征与寡聚化机制研究全解析
已经有0人回复
-
DOTA-GPC3肽（CAS 3066439-94-6）：双功能分子探针前体的结构特征与标记应用解析
已经有0人回复
-
科研用 TLR7/8 双激动剂怎么选？现货与质控要点
已经有0人回复
-
SEL120-34A HCl：CDK8/CDK19高选择性激酶抑制剂
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
