论文部分内容阅读
【目的】确定广西巴马小型猪内源性反转录病毒(PERV-BM)的进化年代。【方法】对先前扩增、测序获得的9个PERV-BM的LTRs,一个PERV-BM全长基因组序列,登录号为HM159246及Gen Bank上发表的所有背景清楚的猪内源性反转录病毒(PERV)的LTRs及env序列进行分析。首先利用DAMBE软件对上述序列比对分析,合并同源性较高的,筛选有差异代表性的核酸序列的数据集。然后利用MEGA5.2软件Neighbor-Joining方法计算PERV-BM遗传进化关系。利用DAMBE软件对上述筛选的序列集进行替换饱和度计算,分析核酸进化速率的稳定性;利用MEGA5.2软件,通过最大相似法对筛选数据制作的遗传进化树进行分子钟行为分析。最后对广西封闭群内的广西巴马小型猪3′-LTR序列进行替代饱和度及分子钟分析,对符合分子钟行为的LTRs序列,以假设的恒定进化速率计算PERV-BM的进化年代。【结果】经过DAMBE软件分析后,选出80多个Gen Bank上发表的有代表性的env、LTRs序列数据集。对PERV-BM(HM159246)序列全长遗传进化分析显示,PERV-BM与中国的PERV核酸序列EF133960和GU980187,荷兰的AF356697亲缘关系较近,与韩国的HQ540595和美国的AF038600和NC003059亲缘关系则较远。利用DAMBE软件对筛选数据集替换饱和度检测结果显示,S和V值随PERV序列之间进化分歧的增加呈线性的增加。LTRs序列数据集替换饱和曲线呈线性关系,没有替换饱和现象发生,LTRs的进化随时间持续且稳定。数据集env序列替换饱和曲线不完全呈线性关系,表明env序列集进化速率不稳定。利用MEGA5.2的最大相似法对LTRs和env数据集的进化树进行分子钟行为检测结果显示,LTRs序列集接受分子钟假说,env进化不符合分子钟行为,这和不饱和度检测的结果一致。因此可用LTRs进化的近分子钟行为进行PERV进化年代的计算。利用DAMBE软件对广西封闭群内的广西巴马小型猪3′-LTR序列进行替代饱和度及分子钟分析结果显示,3′-LTR的进化符合分子钟行为可用于进化年代计算。假设以灵长类ERV的积累突变率,每个核苷酸每年的替换率2.3×10-9—5.0×10-9,计算PERV-BM约进化于3.3×106—7.1×106年前。【结论】PERV-BM进化年代与欧洲小型猪PERV进化时间7.6×106年前相近。