24小时热门版块排行榜    

查看: 2841  |  回复: 18
本帖产生 1 个 MicEPI ,点击这里进行查看

秋电工

新虫 (小有名气)

[求助] 请各位来看棵树

近日学习用mega4.0建系统进化树,大概会操作了,但仍然有很多疑问。坛子里关于mega的帖子看了遍,迷茫发问的朋友多,解答的大牛人也有,我的问题和很多求助的虫友相似,现系统提出,渴望知道的大侠尽早看到此贴,尽早给若干小虫指点下迷津,不甚感激!
一,附件中的小树是用mega4.0,Bootsrap Value设置1000次建成。具体步骤如下:
a、打开mega软件—Alignment 下拉—Alignment exposure/clustal—create a new alignment --然后根据实际情况(DNA or 蛋白序列)选择Yes
b、Edit—Insert sequence from file(.fastal的格式);
c、Algin by clustal W---Data—Export Alignment—Mega fomat;
d、打开上一步保存的文件—Distance—compute pairwise—compute;
e、Phylogeny—Bootstrap Test of phylogeny—Neighbor-Joining—将Boostrap(500 replicates;seed=64238修改为1000 replicates)--compute---获得进化树—设好格式---image –save as TIF file
问题来了,在坛子里看见各位大侠给出的步骤,要用Clustal X软件进行序列比对之后才能建树,请问,mega4.0中是不是嵌合了clustal么,为什么有大侠说还要下载一个clustal软件呢?上面的建树步骤C是不是算经过了clustal呢?
二,附件中小树节点处数值的含义,请各位解释下为什么HX4和AB546196.1(bootstrap值100)与HX6和AJ308316.1(bootstrap值77)两类群聚一起后又97呢?是否意味着HX4和AB546196.1亲缘关系很近,甚至是同一个东西,HX6和AJ308316.1关系较近,这两类群又聚一起,而且bootstrap值97,大于77,很费解呀很费解。
三,树下的标尺有什么用?去算枝长吗?枝长又代表什么,进化距离吗?如果枝长是进化距离,那bootstrap值又代表什么?
四,在坛子里逛的时候,看到有大侠说“遗传标度法和步长值法是用来表示两种不同的结果时用的 步长值常用在鉴定菌种的进化树中 而遗传距离法常用在表示不同菌株间进化关系的进化树中 ”。什么是遗传标度法,什么是步长值法?本人建树目的是想得到菌株间的进化关系
五,看文献过程中,发现有部分文献的进化树要么只有标尺,要么只有bootstrap值,而且也没有对数值,标尺的含义做出解释,究竟什么样的树才算规范?
昨晚码了一个求助贴,满怀希望发的时候,系统维护,伤啊伤!今天抖擞一下,再次情绪饱满的发一个,敬候各位的佳音哈!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

xp198766

铁杆木虫 (著名写手)

小木虫职业打酱油滴~~!

【答案】应助回帖

★ ★ ★
秋电工(金币+6): 很感谢,解答了我大部分问题,但是方法中有些具体的操作我还是不会,希望能更详尽 2011-04-18 11:08:26
scelab(金币+3, EPI+1): good 2011-04-18 11:54:10
1、clustalx 也是为了alignment,如果在MEGA里面alignment了,就不用clustalx了;

2、bootstrap值一般是将你的序列保留一部分,把剩下部分随机打乱,拼成不同的序列,组成1000个你的alignment文件,做树,显示的77表明,1000次做树的过程有有77%次,也就是770次得到HX6和AJ308316.1聚在一起这个结果;

3、下面的标尺是枝长,也是进化距离,bootstrap值是可信度;

4、这个我水平有限,也比较迷惑,个人觉得你现在做成的树适合的进化研究;

5、一般是要有bootstrap值的,这个表示可信度,要是值低于50,一般别人是不认同的,有的文章没有标尺,可能是只想得到树的分枝情况,没有想要计算各个物种的距离,只要得到的一个树的拓扑结果吧,但是bootstrap值一般是要的。

对了,你的alignment,排列好了之后,一般要将首尾序列截成一样齐……

只是个人意见,可能有偏差,有错误欢迎指出。
3楼2011-04-18 10:52:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xp198766

铁杆木虫 (著名写手)

小木虫职业打酱油滴~~!

引用回帖:
1266631楼: Originally posted by 红色日历 at 2012-06-05 18:06:20
请问进化树建出来以后分类地位低的物种的基因比分类地位高的物种的基因在树上的位置还要靠上,这正常吗?...

得到的系统树是根据你的基因序列来的,如果你的序列处理没有什么问题,那看下你用的是不是核心基因,有没有可能发生水平转移或者重组?如果没有的话,再来分析下原因,一般是不会出现这种情况的。。。
10楼2012-06-06 08:59:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖
2楼2011-04-18 10:04:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

秋电工

新虫 (小有名气)

引用回帖:
Originally posted by xp198766 at 2011-04-18 10:52:28:
1、clustalx 也是为了alignment,如果在MEGA里面alignment了,就不用clustalx了;

2、bootstrap值一般是将你的序列保留一部分,把剩下部分随机打乱,拼成不同的序列,组成1000个你的alignment文件,做树,显示 ...

“对了,你的alignment,排列好了之后,一般要将首尾序列截成一样齐……
”弱弱问问大侠,这个具体怎么操作?mega4.0里嵌合的clustal没有截序列的功能么?
4楼2011-04-18 11:13:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xp198766

铁杆木虫 (著名写手)

小木虫职业打酱油滴~~!

引用回帖:
Originally posted by 秋电工 at 2011-04-18 11:13:05:
“对了,你的alignment,排列好了之后,一般要将首尾序列截成一样齐……
”弱弱问问大侠,这个具体怎么操作?mega4.0里嵌合的clustal没有截序列的功能么?

MEGA里面的不知道有没有,没有怎么用MEGA,好像clustalx里面也没有,不过bioedit里面有……可以先在BIOEDIT里面截,再导入到MEGA里面,再转换成MEGA格式,做树
5楼2011-04-18 11:26:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

秋电工

新虫 (小有名气)

引用回帖:
Originally posted by xp198766 at 2011-04-18 10:52:28:
1、clustalx 也是为了alignment,如果在MEGA里面alignment了,就不用clustalx了;

2、bootstrap值一般是将你的序列保留一部分,把剩下部分随机打乱,拼成不同的序列,组成1000个你的alignment文件,做树,显示 ...

菜虫继续发问:关于标尺
1、为什么标尺有时候是1,2,有时候是0.05,0.02?如何设置标尺大小?
2、标尺用来算遗传距离,具体怎么算呢(比如,在附件的图中HX4与AB546196这一群和HX6与AJ308316这一群的距离看哪根树枝呢?或者不相邻的,再远一点儿的类群要看它们之间的进化距离怎么看呢?)
6楼2011-04-23 12:14:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xp198766

铁杆木虫 (著名写手)

小木虫职业打酱油滴~~!

★ ★
rainwander(金币+2): 鼓励交流 2011-04-23 12:54:21
标尺大小,不知道是软件看你的序列默认的还是可以设置,你可以去摸索下;
我觉得标尺是为了直观上看两个物种的距离的吧,一般是聚在一个枝里面的,说明它们相对比较接近,如果你想得到两两之间的距离的话,可以试试软件里面的dist matrix……
或者phylip软件里面的dnadist,我觉得如果不是做进化之类的分析的话,一个diversity distance 就可以吧,mage里面好像就有
7楼2011-04-23 12:20:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

红色日历

木虫 (小有名气)

引用回帖:
907076楼: Originally posted by xp198766 at 2011-04-18 10:52:28
1、clustalx 也是为了alignment,如果在MEGA里面alignment了,就不用clustalx了;

2、bootstrap值一般是将你的序列保留一部分,把剩下部分随机打乱,拼成不同的序列,组成1000个你的alignment文件,做树,显示的 ...

请问进化树建出来以后分类地位低的物种的基因比分类地位高的物种的基因在树上的位置还要靠上,这正常吗?
8楼2012-06-05 18:06:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xp198766

铁杆木虫 (著名写手)

小木虫职业打酱油滴~~!

引用回帖:
1266631楼: Originally posted by 红色日历 at 2012-06-05 18:06:20
请问进化树建出来以后分类地位低的物种的基因比分类地位高的物种的基因在树上的位置还要靠上,这正常吗?...

得到的系统树是根据你的基因序列来的,如果你的序列处理没有什么问题,那看下你用的是不是核心基因,有没有可能发生水平转移或者重组?如果没有的话,再来分析下原因,一般是不会出现这种情况的。。。
9楼2012-06-05 22:16:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 秋电工 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
[考研] 一志愿985,本科211,0817化学工程与技术319求调剂 +3 Liwangman 2026-03-15 3/150 2026-03-15 18:16 by JourneyLucky
[考研] 0703化学调剂 ,六级已过,有科研经历 +4 曦熙兮 2026-03-15 4/200 2026-03-15 18:01 by JourneyLucky
[考研] 274求调剂 +4 时间点 2026-03-13 4/200 2026-03-15 15:29 by Rambo13
[考研] 268求调剂 +5 一定有学上- 2026-03-14 6/300 2026-03-14 22:20 by 运气yunqi
[考研] 080500,材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
[考研] 一志愿安徽大学材料工程专硕313分,求调剂的学校 +8 Yu先生 2026-03-10 10/500 2026-03-14 01:04 by JourneyLucky
[考研] 调剂 +3 13853210211 2026-03-10 3/150 2026-03-14 00:47 by JourneyLucky
[考研] 304求调剂 +6 Mochaaaa 2026-03-12 7/350 2026-03-13 22:18 by 星空星月
[考研] 泣血叩求调剂恩,愿以丹心报师恩 +6 Iuruoh 2026-03-11 6/300 2026-03-13 22:06 by JourneyLucky
[考研] 求材料调剂 +5 隔壁陈先生 2026-03-12 5/250 2026-03-13 22:03 by 星空星月
[考研] 290求调剂 +9 ADT 2026-03-11 9/450 2026-03-13 21:55 by JourneyLucky
[考研] 304求调剂 +7 7712b 2026-03-13 7/350 2026-03-13 21:42 by peike
[考研] 285学硕求调剂 +6 柴郡猫_ 2026-03-12 6/300 2026-03-13 20:46 by hmn_wj
[考研] 0703一志愿211 285分求调剂 +4 ly3471z 2026-03-13 4/200 2026-03-13 13:00 by JourneyLucky
[考研] 0856化学工程280分求调剂 +4 shenzxsn 2026-03-11 4/200 2026-03-13 11:55 by ymwdoctor
[考研] 333求调剂 +3 152697 2026-03-12 4/200 2026-03-13 07:08 by Iveryant
[考研] 420求调剂 +4 莫向外求11 2026-03-10 6/300 2026-03-12 14:41 by ruiyingmiao
[考研] 293求调剂,一志愿陕师大生物学 +3 ??????.?.??? 2026-03-09 3/150 2026-03-11 10:02 by 学员8dgXkO
[考研] 0856材料与化工353分求调剂 +11 NIFFFfff 2026-03-09 11/550 2026-03-10 18:36 by suyuanhai
信息提示
请填处理意见