分子实验学习总结——郁娟娟.doc

DNAstar：序列拼接MEGA：分析碱基组成CLUSTAL：序列比对PAUP：跑树DABE：饱和性分析总思路：序列拼接-比对-四、下载序列的方法打开.Fasta格式的序列-File export-sequance alignment-保存在NCBI上找到相应的序列方法：Search Nucleotide for 文章的genibank序列号-reports-复制序列到txt文件中，删除没用的东西，只留下序列和学名。保存成txt文件后，打开Editseq-将,txt序列复制到其中-save保存为seq格式五、将拼接序列翻译成蛋白质方法打开拼接的序列-全选中-Goodies-translate DNA2、如何计算碱基含量比值、答：MEGA-Nucleotide Composition六、DAMBE序列饱和性分析序列的饱和性分析在DAMBE程序中，以转换、颠换数为纵轴，以TN93模型（Tamura and Nei, 1993）校正的距离为横轴做散点图。如果这些点随序列间分歧增大呈线性分布就说明序列间替换没有达到饱和，序列可以用于后续的系统发育分析。DAMBE还长用在单倍型的归纳、基因频率和碱基组成分析、遗传距离计算、序列排序、5.2.2碱基替换饱和效应的检验当核苷酸序列分歧较小时，序列之间被观察到的核苷酸差异数接近实际替换数。如果序列分歧较大，DNA序列的变异可能会受到饱和效应的影响，这时观察到的碱基差异可能包含了部分进化杂音。一般说来，碱基转换发生的频率高于颠换，但是随物种分歧程度的增加，转换/颠换的比率接近或小于1，说明转换可能已经趋于饱和并可能成为进化杂音，所以转换和颠换发生的频率与序列间分歧度的关系就可以反映出碱基替换是否受饱和效应的影响。在着手构建分子系统树前，为了排除这种进化杂音获取正确的进化信息，必须对目的片段的碱基替换是否受到饱和效应的影响进行检验。用 DAMBE（Xia，2000）分别对四个mtDNA片段的碱基替换模式进行检测，以此估计所研究mtDNA片段的碱基替换是否受到饱和效应的影响。打开DAMBE-FILE-Open standard sequence file-类型unknown-打开proteincoding nuc.seq-invMtdna-trains-table go-graphics-trainsition and tranversion versis advance tn93-g0-graphstyle-curve only-保存.bmp问题1、 如何计算密码子的不同位上碱基变化，转换颠换比值2、 如何进行序列的饱和性分析序列的饱和性分析在DAMBE程序中，以转换、颠换数为纵轴，以TN93模型（Tamura & Nei, 1993）校正的距离为横轴做散点图。如果这些点随序列间分歧增大呈线性分布就说明序列间替换没有达到饱和，序列可以用于后续的系统发育分析。3、 如何看用Editseq检验拼接后的序列是否能正确翻译为蛋白质4、 如何计算遗传距离（DAMBE）一、DNAstar：序列拼接在序列的拼接时，修改红色和小写字母的碱基具体步骤：1、 ABI文件包括图形和序列、SEQ文件只有序列没有图。2、 DNAstarseqmanFileNEW，把序列都放近去，有两个正向，俩个反向的序列Time endsLOWSEAN ALL-ASSEMBLE 出现contig 打开-双击contig_核对-contig-save consenus-single file -保存位置-sequce-add-选择保存到的位置-点序列名字-add-双击-拼接完关闭。3、 原始序列和拼接的序列各保存在一个文件夹内。4、 EDIT-GOODIES-TRALS 翻译成蛋白质。5、 从ncbi上下载相关的序列和测得的序列保存在一个文件夹内，把前面的内容都删除，保留种名。加符号-保存二、序列的比对6、 先切齐序列：用Seqman-sequence-add-skip-把所有的序列加进去进行删除，使序列整齐。-contig双击,将序列切齐，再复制到一个新建的文本文档。7、 切齐后保存方式：contig-export sequences-multiple files=set location-file-保存切齐8、 当在切齐的序列中发现反向序列要改正：用editseq打开序列-select all goodies-reserse complement-全选粘贴到文本文档中，目的是比对。三、构建NJ树方法(MEGA)1、将TXT格式保存为nexus和clustal格式（生成alndndnxs）CLUSTAR-进行比对-FILE-LOAD SEQUENCE-打开新建的文本文档txt-加入进去后-aligment-do complete Alignment-File-Save sequence as format-custal-ok-File- Save sequence as format-NEXUS最后保存为nexus和clustal格式，NEXUS是PAUP格式-关闭clustal.2、点击MEGA-File-open data-.aln文件-打开-ok-点击向下的绿色图标-ok-ok-删除最后没用的东西-保存-切齐序列,meg-3、点击MEGA-File-open data-切齐序列,meg-ok-yes-inverttebrate Mitochondrial-ok-4、phylogency-construct phylogency - NJ树-none改成bootstrap-k 2p-compute-保存（image-EMF）5、使用Mega 2.0分析序列间的p-distance, 即两两序列间的核酸差异比例，分析序列变异情况方法1：Distance-compute pairwise-distance only改为std err-compute方法2：pattern-compute composition distance三、构建MP树方法1、MrBays，目的是转化为.nex文件步骤为：用bioedit将.txt文件比对转换成.nex文件格式方法：打开BioEdit，将.txt文件打开-Accessory application-clustal multiple alignment-File-export-sequence Aligment-选paup/nex命令.nex-保存-把.nex begin 前面的都删除。 File-open-切齐序列.txt-accessory appliacation-clustal multiple aligment-file-expore-sequence aligment-选paup/nenu命名.nex-保存2、在.nex中加命令模板;end;begin paup;log file=coimp;set maxtrees=1000;set criterion=parsimony;hsearch addseq=random nreps=1000;showtrees;describetrees 1/plot=phylogram brlens=yes;savetrees file=coi(mp).tre root=yes brlens=yes;contree all/file=coi.tre;pscore all/CI=yes RI=yes RC=yes;bootstrap nreps=1000 keepall=yes search=heuristic/addseq=random nreps=100;savetrees file=coi(mpb).tre root=yes;END;2、打开paup-打开1.nex-open 开始跑树3、用treeview 打开.tree树。MrBays建树1、 用bioedit将.txt文件比对转换成.nex文件格式方法：打开BioEdit，将.txt文件打开-Accessory application-clustal multiple alignment-File-export-sequence Aligment-选paup/nex命令.nex-保存-把.nex begin 前面的都删除。 File-open-切齐序列.txt-accessory appliacation-clustal multiple aligment-file-expore-sequence aligment-选paup/nenu命名.nex-保存2、 用paup运行.nex文件3、 接着继续用paup打开modeltest3.7 文件夹中的model paupblock命令-运行结束后，在modeltest 的modelfolder中生成一文件名为”model.scores”文件。4、 运行modeltest批处理文件，生成一文件名为”mt.txt”文件，则为最适合的模型文件 5、 将最适模型文件中的参数(hLRTs)拷贝到.nex文件后，加入贝叶斯运行命令，并改正.ne格式等。删除前后没用的，改正外群名字等。只能用一个外群。begin mrbayes;log start filename=coibayes.log;outgroup Cinara;Prset statefreqpr=dirichlet(0.3321,0.1058,0.1271,0.435);Lset nst=6 rates=gamma;mcmc ngen=1000000 printfreq=100 nchains=4 savebrlens=yes startingtree=random;end;6、 将加好命令的.nex文件拷贝到MrBayes文件夹所在路径内，运行人头的MrBayes , 输入execute. 文件名.net,回车。7、 输入 sump burnin=1000回车8、 输入 sumt burnin=1000回车9、运算完毕后点击close 弹出一对话框，输入命令语句“ save tree file文件名.tre；”回车。菜单栏Tree中的Show internal edge labels 即可显示各结点的bootstrap值。9、 生成的.con，用treeview打开注意：贝叶斯树只能用一个外群。建设一个新文件夹，将.scores、mt.txt、txt、nex文件都放到一个文件夹内。ML建树方法同MrBays相似。预测模型利用贝叶斯得到的模型。把ML命令模板拷贝到.nex，把log file改成自己的名(=coiML)，外群也改(Daktulosphaira_vitifoliae Cinara_edulis C_arizonica);-把最适模型lset base-likhood下面savetree file=coi savetrees file=coin(mb) Lset Base=(0.3365 0.0976 0.1280) Nst=6 Rmat=(26794342.0000 8028280.5000 10365153.0000 1.4084 255543472.0000) Rates=equal Pinvar=0.7355;生成的.nex.con文件-treeview打开。注意：外群要写成Cinara_edulis形式，可以有多个外群。打开方式和mp树一致。戴传银笔记PAUP NJ treebegin paup;outgroup outgroupname;setcriterion=distance;bootstrap search=nj nreps=1000 keepall;savetrees file=nj.tre;end;PAUP MP TREEbegin paup;outgroup outgroupname;setcriterion=parsimony;hsearch addseq=simple(random) bootstrap nreps=1000;describetrees 1/plot=phylogram brlens=yes;savetrees from=1 to=1000;end;PAUP ML TREEbegin paup;outgroup outgroupname;setcriterion=likelihood;gettrees file=xxx.tre;(xxx modeltest文件保存的树的名称)hsearch addseq=asis bootstrap nreps=1000;savetrees file=ml.tre;end;MEGA3.1分析其核酸组成1、打开MEGA，点击“File”中的“Open Data”打开“COI.aln”文件，点击绿色箭头“Convert to Mega format”快捷键，出一把COI.aln转换成Mega格式的对话框，点击OK，即转换成COI的Mega file，保存。2、用Mega 打开该文件。选Data Type。Proteincoding nucleotide sequence data？Yes。Select genetic Code。程序自行计算。3、 点击TA图标，出现Sequence Data ExplorerC：Conserved sites 291/660,分母613是总位点数V：Variable sites 352/660P: Parsimony-informative sites 26/660N: Nonparsimony-informative sites 352-26326。4、 Mega-File-open date-运算-sequence Data Explorer-write data to file-title-写名字.nex注意：序列一定要正确，不能是反向的！5、 打开paup-打开1.nex-open 开始跑树6、 Bootstrap 值检验及树的保存输入命令“bootstrap nreps1000 keepallyes searchheuristic/addseqrandom nreps100；”点击回车键，程序开始运算。运算完毕后点击close 弹出一对话框，输入命令语句“ save tree file文件名.tre；”回车。菜单栏Tree中的Show internal edge labels 即可显示各结点的bootstrap值。7、 用treeview 打开.tree树。阿征笔记begin paup; set autoclose=yes notifybeep=yes; log file=log.txt; outgroup A00; set root=outgroup outroot=monophyl; set criterion=parsimony; hsearch addseq=random nreps=500 nchuck=5 chuckscore=1; savetrees format=nexus brlens=yes append=yes file=mp.tre root=yes; contree all/file=mp(contree).tre strict=no majrule=yes percent=50 le50=yes; bootstrap nreps=1000 conlevel=50 keepall=no search=heuristic/addseq=random nreps=100; savetrees file=mp(bootstrap).tre root=yes brlens=yes savebootp=both from=1 to=1000;end; 征征 20:31:57我每次都是先写到contree那一步，然后在里面敲命令，后两条分别征征 20:32:13begin paup; set autoclose=yes notifybeep=yes; log file=log.txt; outgroup A00; set root=outgroup outroot=monophyl; set criterion=parsimony; hsearch addseq=random nreps=500 nchuck=5 chuckscore=1; savetrees format=nexus brlens=yes append=yes file=mp.tre root=yes; contree all/file=mp(contree).tre strict=no majrule=yes percent=50 le50=yes;end; 征征 20:32:25然后bootstrap nreps=1000 conlevel=50 keepall=no search=heuristic/addseq=random nreps=100;征征 20:32:36跑完后 savetrees file=mp(bootstrap).tre root=yes brlens=yes savebootp=both from=1 to=1000;征征 20:32:54保存一个有bootstrap的树征征 20:33:35再给你个阿荣给我的征征 20:33:42set autoclose=yes criterion=parsimony notifybeep=yes;log file=.;Hsearch addseq=random nreps=100 start=stepwise savereps=yes randomize=addseq rstatus=yes hold=1 swap=tbr multrees=yes;savetrees format=nexus brlens=yes append=yes file=.;contree all/strict=no majrule=yes percent=50 le50=yes;bootstrap nreps=1000 conlevel=50 keepall=yes search=heuristic/addseq=random nreps=10;征征 20:34:16对了，我的outroot=monophyl;你们分子不是monophylPAUP软件使用简要说明 1.数据输入格式 将需要分析的一组DNA数据经Clustal软件比