实战:用BEAST构建带分子钟的皿蛛分子系统发育树

本文介绍了如何使用BEAST软件构建带分子钟的皿蛛分子系统发育树。首先,从Wang et.al,2015的研究中选取122个taxa的DNA序列数据。接着,通过BEAUti生成XML配置文件,其中包括设置tMRCA、分子替换模型、分子钟模型、先验分布、起始树、操作符、MCMC参数。在MCMC过程中,为了保持树拓扑一致,删除了部分操作符。最后,利用Tracer分析结果并用TreeAnnotator得到合议树。注意,对于大量树木的合议,可能需要调整jvm内存限制。" 46801067,1358276,Perl操作Cookie的示例解析,"['Perl编程', 'Web开发', 'HTTP交互', '会话管理']
摘要由CSDN通过智能技术生成

一、取样:以Wang et.al,2015中L148_r树的taxa为基础,去除重复的taxa,最终得到122个taxa。DNA序列数据L122.nex, 起始树为ML法推断出的besttree。

 二、生成XML文件

        A. 打开BEAUti,将DNA序列数据和起始树导入程序,按以下步骤生成跑BEAST所需的xml文件。。

          (1)定义tMRCA,共4个。

                    Outgroup(2taxa,2种派模蛛,非强制单系)

                     Linyphiidae(120taxa,所有传统分类中定义的皿蛛科物种,非强制单系)。

                     Linyphiidae(117taxa,except stemonyphantes亚科的3个taxa, 强制单系,用于设定化石矫正点)

                     Orsonwelles(2taxa, include Orsonwelles属两种,强制单系,用于设定地质矫正点?)

                    

           (2)核苷酸分子替换模型。Site Model:GTR。事先使用jModelTest等软件对该最适模型及其参数进行判定。此部分请参考高芳銮的博文 http://blog.sciencenet.cn/blog-460481-791370.html

            (3) 分子钟模型Clock Models。 将Molecular Clock Model选为“Relaxed Clock:Uncor

  • 10
    点赞
  • 25
    收藏
    觉得还不错? 一键收藏
  • 13
    评论
一、课程导语不少人因各式各样原因,在读研时跨行生信或从事生信相关研究,其中不少开端即自学。然漫漫自学路,迷茫困苦乃至萌生放弃做科研的念头。念此,决定在比较擅长的领域-病毒遗传进化方面做点尝试,开一门理论和实战演示并重的基础课。注:课程里数据分析演示,以病毒序列数据作为背景和例子。然大道共通,做其他生物遗传进化的同学也可以参考。 二、课程目标掌握常用的遗传进化分析方法和流程,包括序列下载和预处理、序列比对、基因突变分析、重组和重配分析、系统发育分析、基因选择压力分析、分化时间和进化速率估算、系统发生生物地理学分析(病毒时空动态分析)、蛋白质三级结构预测等。 三、适合人群(1)零基础,生物信息学入门;(2)有基础,需要提升的同学。 四、课程特色从算法原理到数据分析演示,知其所以然,避开数据分析黑箱子。思维导向,深入浅出式讲解,系统学习遗传进化相关知识。 五、课程安排【基础分析】专题一:病毒分类,序列下载,序列格式1. 病毒分类2. 病毒序列下载3. 序列格式解读4. 序列格式转换专题二:序列比对原理及应用1. 打分矩阵2. 序列比对策略3. 双序列比对4. 多序列比对专题三:基因突变分析、重组和重配分析1. 病毒基因突变分析2. 病毒重组和重配分析专题四:系统发育分析之理论基础篇1. 认识进化2. 分子假说3. p距离与泊松校准4. 核苷酸进化模型深入剖析5. 几种建方法原理详解专题五:系统发育分析之实战操作篇1. 建流程2. 估算进化模型3. 系统发育构建(NJ,ML,BI)4. 的查看和美化专题六:系统发育分析进阶1. 系统发育信号检测2. 拓扑结构比较3. 多基因联合建【高级专题】专题七:基因选择压力分析1. 自然选择VS中性进化2. 枝长模型3. 位点模型4. 枝位点模型5. 进化枝模型6. PAML进行选择压力分析专题八:分化时间和进化速率估算1. 使用Beast估算分化时间和进化速率2. 使用Treetime估算分化时间和进化速率3. 使用LSD估算分化时间和平均进化速率专题九:系统发生生物地理学1. 谱系生物地理学介绍2. 贝叶斯随机搜索量(BSSVS)分析3. 结构化溯祖分析专题十:蛋白质三级结构预测1. 蛋白质三级结构介绍2. 蛋白质结构比对3. 三级结构预测方法4. 三级结构预测实践5. 三级结构质量评估6. 三级结构查看与作图
评论 13
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值