有限群体中基因频率概率分布的模拟实验 A Simulation Study on Probability Distribution of Gene.pdfVIP

遗传 HEREDITAS(Beijing)15(3):12-151993 有限群体中基因频率概率分布的模拟实验 杨长锁 吴常信 青（海畜牧兽医学院畜牧系，西宁810003) 北（京农业大学畜牧系，北京100094) 本文采用计算机模拟的方法，对有限群体中基因频率的概率分布进行了研究． 关锐词 计算机模拟，基因频率，概率分布 ASimulationStudyonProbabilityDistribution ofGeneFrequencyinFinitePopulation YangChangsuo (Dept.ofAnirnalScience,QinghaiAnimalHusbandryandVeterinaryMedicineCollege,Xining810003) WuCbangxin (Dept.ofAnimalScience,BeijingAgriculturalUniversity,Beijing100094) Theprobabilitydistributionofgenefrequencywasstudiedbysimulativemethodinfinitepopulation. Keywords Genefrequency,Probabilitydistribution,Simulation 基因频率是研究群体遗传组成的一项重要指标，在种群的繁衍过程中，因有限繁殖群体的随 机留种使群体的基因频率随世代不断变化，这种随机变化的过程已有较多的研究 ’一〔，,5,7,9)，而对有 限群体基因频率概率分布的研究却相对较少。Kimura(6）用数学方法推算出有限群体中基因频 率概率的理论分布，在数学推导时作了许多假定条件，且很难得知某世代时基因可能固定或丢失 的概率；Buric4)用普通果蝇的105个世系每（个世系包含8只雄虫和8只雌虫）的棕·眼位·点上的 两个等位基因bw75和bw，研究了基因频率的概率分布，但世系内个体数较少，且实验进行的世代 数有限。由于采用生物实验来研究这个问题比较困难，因此本文用计算机模拟方法对基因频率概 率分布进行了一些讨论，也为研究有限保种群体基因频率的概率分布提供了一种参考方法。 一、模拟过程描述 一（）参数组合 假定有限群体中某性状受一个位点两个等位基因的作用，其基因型模式为： 基因型 AIA, A,A2 A,A, 基因型频率 产 2Pq 扩 其中P为 A1基因的频率，q为A：基因的频率，并假定零世代基础群体很大，且处于 Herdy-Weiberg遗传平衡状态，A,和A2基因的初始频率均为0.5，有限群体的有效大小N（e)分别 12 · 为20和40，亚群体数为100,模拟连续50个世代． 二（）模拟过程 模拟模型的基本思想与陈贵才 ’〔〕所采用的模型基本一致，但在模拟模型的程序设计上采用 了模块结构．主控程序控制各模块的执行，其主要模块为： 模块I：由乘同余法产生伪随机数模块。 模块II：计算各亚群体基因频率，并统计具有相同基因频率范围内亚群体数以及A,基因固定 或丢失的亚群体数目． 模块ul：绘制直方图及打印输出模块。 模拟程序用Basic语言编写，各世代亚群体大小始终保持恒定，整个模拟过程在长城0520CH 微机上完成。 二、结果与分析 为了说明有限群体基因频率概率分布的过程，仅选择了相隔5代的模拟结果和典型的概率分 布图． 一（）基因频率概率分布 由表1可知，有限群体基因频率的概率分布随世代逐渐分散，即基因频率分布的标准差逐渐 加大。表2的数值也具有此特点，但因有效群体大小增加了一倍，使其概率分布的分散性减慢，这 一特征可从图1和图2清楚看出．当Ne二20时，分布