用Phylomatic和PhyloCom进行.doc

用phylomatic和phylocom进行群落系统进化分析张金龙j（中国科学院植物研究所）phylomatic和phylocom软件，是哈佛大学的cam webb博士、俄勒冈大学的steve kembel博士、加州大学伯克利分校的david ackerly教授编写的一套用于群落物种组成系统发育关系的软件。phylomatic是在线软件，可以利用植物名录，按照apgiii的被子植物科的拓扑结构，生成进化树。由于当前物种测序还不够充分，某一个群落中或某一个地区所有物种的某几个基因的序列还不能全部拿到，因此可以用phylomatic建立基于apgiii骨架的进化树。phylocom软件是用来进行群落系统发育与进化分析的。phylocom可以为phylomatic软件得到的进化树拓扑结构，按照一定的规则拟合一定的枝长，其功能还包括：计算群落的系统发育多样性（pd），计算群落的系统发育结构(community structure)，计算群落的系统发育距离(community phylogenetic distance)，分析群落的性状进化(aot)。下面就介绍如何使用phylomatic和phylocom进行相应的分析。一phylomatic建立进化树1物种名录的准备复制拉丁文属名，打开/phylomatic/kewlookup.html查询每个属所在的科等信息图1查询属名所在的apg科输出结果如下：图2属名查询结果新建一个excel空白文档，将结果粘贴到excel中选中粘贴过来的列，点击excel的菜单，数据分列分隔符号空格从中选择apg fam列图3选取apg科名将没有查到的属，手动添加相应的科名。整理种名，去掉去空格和括号，报名命名人等信息将种名中所有的空格用“.”或者“_”代替。按照科/属/种的顺序，将各列粘贴到一个新的excel表格中。图4删除物种的命名人和括号图5将物种内的空格用“_”替换，并粘贴到对应的科属删除其中的非被子植物。将非被子植物删除后，选中全部内容，复制后粘贴到记事本中。将其中的制表符，用“/”来替换。所得数据，形式如下：图6将excel表格中的数据粘贴到记事本中，并将tab制表符替换为“/”将所有的科名首字母，改为小写字母。将所得数据，粘贴到phylomatic在线对话框中。2phylomatic建树图7将科的起始字母替换为小写字母，并粘贴到phylomatic的对话框中。数据输出形式，选择默认的newick即可。点击提交，即可得到无枝长的，newick格式的进化树。(castanopsis_eyrei,castanopsis_carlesii,castanopsis_fargesii,castanopsis_tibetana)castanopsis,(cyclobalanopsis_gracilis,cyclobalanopsis_glauca)cyclobalanopsis,lithocarpus_glaber)fagaceae,myrica_rubra),elaeocarpus_decipiens),syzygium_buxifolium),(daphniphyllum_oldhamii,loropetalum_chinense),(rhododendron_latoucheae,rhododendron_ovatum)rhododendron,vaccinium_carlesii)ericaceae,(schima_superba,adinandra_millettii,eurya_muricata,camellia_fraterna,ternstroemia_gymnanthera)theaceae),(machilus_thunbergii,neolitsea.aurata.var.chekiangensis)lauraceae,chimonanthus_salicifolius)laurales)magnoliales_to_asterales)euphyllophyte;图8输出的进化树将输出网页中的进化树,复制到记事本中，并另存为phylo文件，该文件不要有任何扩展名。二phylocom软件的使用本部分内容包括：用phylocom的bladj模块为进化树添加枝长,.利用phylocom中的练习数据，计算群落的系统发育多样性，系统发育结构等。1下载phylocom软件/phylocom/图9填写信息，并下载phylocom首先提交相应的信息，之后，便可下载phylocom软件将下载的phylocom4.1.zip文件解压缩，图10将下载zip文件解压缩example_data中是练习数据mac中是苹果机的运行软件r，是驱动phylocom的r脚本src是phylocom的源程序，phylocom是c语言写的。w32是windows平台下可以运行的exe程序phylocom_manual是说明书readme注意事项2创建工作路径在c盘根目录下，创建一个名为phylocom的文件夹。将w32中的phylocom.exe文件，拷贝到该文件夹c:phylocom。再将做好的包含进化树的phylo文件，拷贝到该文件夹c:phylocom将phylocom-4.1example_databladj_example文件夹下的wikstrom.ages文件，拷贝到c:phylocom，并将wikstrom.ages改名为ages，也不留任何扩展名。图12将ages文件和phylo文件phylocom,exe程序拷贝到c:phylocom文件夹3用bladj为进化树拟合枝长运行phylocom开始运行输入cmdcd c:phylcom图13用cd c:phylocom切换路径。图14用phylocom的bladj拟合枝长，并存储到out文件中输入phylocom bladj out.tre则生成的out.tre就是含有枝长的进化树该进化树可以用treeview软件，figtree软件查看。图15用figtree软件绘制的进化树4 phylomatic软件生成的进化树在r中的操作在r软件的ape程序包中，计算out.tre,即所得有枝长进化树中，物种两两之间的进化距离（分化时间）图16在r软件的ape程序包中，读取phylomatic树出现的问题phylomatic建立的进化树，在ape程序包中，不能正常读取，提示“there is apparently two root edges in your file: cannot read tree file.reading newick file aborted at tree no.1”图17用记事本打开拟合好枝长的进化树,删除最外面的euphyllophyte及相应括号，保留“;”删除进化树中最外面的一层,euphyllophyte，包括前面的一对括号，即紧挨着euphyllophyte的“)”和最前面的“（”。最后的分号保留。之后，即可在ape中读取。在r中读取进化树，并计算物种之间的进化距离脚本:library(ape)setwd(c:/phylocom/)tr - read.tree(out.tre)phylodist cmd cd c:phylocom输入：phylocom comstruct comstruct.output计算群路的系统发育结构。图18用phyloco