实施方案
[0025] 下面结合附图具体实施例来详细描述本发明的技术方案。
[0026] 参照图1,一种问号赖型钩端螺旋体复制基因挑选方法,包括以下步骤:
[0027] 1)获得问号赖型钩端螺旋体
[0028] 问号赖型钩端螺旋体的全基因组可在NCBI下载,包括其各段基因信息;
[0029] 2)挑选复制基因
[0030] 通过Blast软件进行比对挑出复制基因,首先,我们筛除了可能的没有功能的序列,然后,用blast对比所有的问号赖型钩端螺旋体复制基因,选取对比分数在10-8以下的,通过BLAST软件基于二参数模型,MATLAB,VC++软件进行挑选、比对,选择复制基因,再用Kestim61软件分别计算出每对复制基因的每个同义位点和非同义替代位点的替代数目,总共获得519对复制基因;
[0031] 3)进化距离信息提取
[0032] a)通过Blast软件进行比对挑出复制基因.首先,我们筛除了可能没有功能的序列,然后,用blast软件对比所有的问号赖型钩端螺旋体基因,选取对比分数在10-8以下的.[0033] b)通过BLAST软件基于二参数模型,MATLAB,VC++软件进行挑选、比对,选择复制基因.并同时对进化距离参数信息进行提取.
[0034] 首先,采用基于差分方程建立二参数模型,建立方程分别计算每对复制基因的每个同义位点和每个非同义位点的替代数目.设分歧时间为T,α代表每单位时间单位位点的转换替代率,β代表每单位时间单位位点的颠换替代率,P指同源位点转换替代率,Q指同源位点颠换替代率:
[0035] ΔP(T)/ΔT=2α-4(α+β)P(T)-2(α-β)Q(T)
[0036] ΔQ(T)/ΔT=4β-8βQ(T)
[0037] P(0)=Q(0)=0.
[0038] 同源位点转换替代计算公式:
[0039]
[0040] 同源位点颠换替代计算公式:
[0041]
[0042] 获得的进化距离计算公式:
[0043]
[0044] 然后,采用Kestim61软件分别计算每对复制基因的每个同义位点和每个非同义位点的替代数目,
[0045] 最后,采用C++软件对Blast软件的结果进行筛选、精炼和识别,并采用Matlab和origin获得效果图.最后获得问号赖型钩端螺旋体519对复制基因.同样的复制基因筛选方法也可以用在其他的生物的复制基因的挑选。
[0046] 同时,本方法也可以推广到其他基因或基因组。方法如下:可以在步骤(1)下载相关基因,顺序按照步骤(2)(3)进行即可。
[0047] 以上所述,仅为本发明较佳的具体实施方式,本发明的保护范围不限于此,任何熟悉本技术领域的技术人员在本发明披露的技术范围内,可显而易见地得到的技术方案的简单变化或等效替换均落入本发明的保护范围内。