当前位置： 首页 > news >正文

外国人做的购物网站,文章网站模板,在线设计平台用户分析,河南郑州房价Tajimas D Test 已经开发了几种中性检验#xff0c;用于识别模型假设的潜在偏差。在这里#xff0c;我们将说明一种有影响力的中性检验#xff0c;即Tajimas D#xff08;Tajima 1989#xff09;。Tajimas D通过比较数据集中的两个#x1d703; 4N#x1d707;估计值来…Tajimas D Test 已经开发了几种中性检验用于识别模型假设的潜在偏差。在这里我们将说明一种有影响力的中性检验即Tajimas DTajima 1989。Tajimas D通过比较数据集中的两个 4N估计值来工作。我们已经推导出了is等于平均成对杂合性average pairwise heterozygosity当我们讨论共祖时也称为Tajima的估计器。当考虑DNA序列集合中的等位基因或SNP总数以及将它们联合到共同祖先的共祖树内包含的预期世代数时可以得出另一种推导。这被称为Watterson的估计器或WYong 2019。正如我们在共祖中所展示的n个谱系的预期共祖时间是 一组n个初始谱系到一个单一祖先的所有共祖时间的总和是 在每次共祖事件之间的步骤中有i 1个谱系可能发生突变。因此当考虑一段时间内所有谱系可能产生的等位基因数量时我们乘以i 1 我们将树上所有谱系的总时间乘以每代突变率得到我们期望在n个DNA序列样本中的等位基因总数而4N是一个常数所以我们可以将其放在求和之外 其中S是序列集合中SNP的数量。这可以重新排列以从SNP的数量估计W4N 请注意Watterson的估计器需要了解一组谱系的共祖但这是在1975年发表的当时还没有发表超过两个谱系的共祖Kingman 1982。Tajimas D是平均成对杂合性估计值与从样本中SNP数量估计的之间的差异除以该差异的预期方差的平方根 而 这看起来有点乱但在各种位置反复出现两个不同的n求和只需要计算一次然后填入。这可以用以下代码计算平均成对差异和S也可以从数据集中计算但为了简洁起见我们在这里省略了。  Tajimas D是用来评估一个种群中中性突变即没有自然选择影响的突变的假设是否成立。它通过比较两个不同的估计器来衡量种群的遗传多样性和种群规模的变化我们通过R语言实现  

Calculates Tajima’s D# 平均成对差异用于估计theta_IS

theta_IS - 2.8# 数据集中的SNP数量 S - 16# 采样的等位基因拷贝数 n - 20# 初始化求和变量 i1_sum - 0.0# 循环计算i1的和这是Wattersons theta的一部分 for(i in 2:n-1){i1_sum - i1_sum 1/i }# 计算Wattersons theta它是基于序列多态性的一个种群规模的估计器 theta_W - S / i1_sum# 初始化第二个求和变量 i2_sum - 0.0# 循环计算i2的和用于后续计算 for(i in 2:n-1){i2_sum - i2_sum 1/i^2 }# 计算期望值e1它是Tajimas D公式中的项 e1 - ((n1)/(3(n-1)) - 1/i1_sum) / i1_sum# 计算期望值e2它也是Tajimas D公式中的项 e2 - (2(n^2n3)/(9n(n-1)) - (n2)/(n*i1_sum) i2_sum/i1_sum^2) / (i1_sum^2 i2_sum)# 计算Tajimas D值它衡量的是theta_IS和theta_W之间的标准化差异 (D - (theta_IS - theta_W) / sqrt(e1*S e2S(S-1))) 前三个变量将根据您的数据进行调整。在这个例子中返回的D -1.409。大于或小于2的D被认为是显著的然而实际的p值是通过模拟确定的。D的正值表示中间频率等位基因过多这可能是由于人口减少或平衡选择因为这两种情况都会延长人口历史较老部分的共祖事件时间。在更大的种群中有更多的祖先可供选择共祖是一种罕见的事件并会膨胀IS因为相对于W较老的谱系在后代中以更高的频率共享。这种负D表示稀有频率等位基因过多共祖的最近尖端被放大它们对S和W的贡献比对平均成对差异的贡献更多因为它们很稀有这表明人口扩张、选择性清除或对有害等位基因的低效净化选择。关键是要计算多个位点的D值并寻找异常值以标记假定选择候选者。人口统计学效应如人口规模的变化应该影响基因组中的所有位点而选择通常被认为在其影响上是位点特异性的。还有许多其他的中性检验如HKA检验Hudson等人1987、McDonald-Kreitman检验McDonald和Kreitman1991、Fay和Wu的HFay和Wu2000以及dN/dS比率Yang和Bielawski2000。其中许多也利用了物种之间发生的遗传变化它们都有各自的优点和缺点。 linkage disequilibrium (LD) 群体遗传学的一个独特性质在进化博弈论等类似领域中并未发现即不同位点甚至不同染色体上的等位基因可以“链接”尽管并不总是与经典遗传学中的重组图谱同义并且比随机预期更频繁地一起遗传。 连锁不平衡LD的程度由量化不要与Tajima的D混淆。考虑两个位点一个具有A/a多态性另一个具有B/b多态性。我们对于跨位点一起遗传的等位基因之间的关联感到好奇。使用概率的乘法规则我们期望AB单倍型pAB的频率是两个等位基因频率pApB的乘积如果它们是独立遗传的话。这两者之间的差异由量化作为连锁不平衡的一种度量。 根据AB单倍型是过量还是不足可以是正数或负数或者如果你从ab单倍型任意计算符号会改变。也可以从所有单倍型频率计算得出。 为了说明假设我们有一个包含两个SNP的单倍型频率的小型数据集。一个是A/G多态性另一个具有C/T等位基因 让我们关注A-C单倍型。A等位基因的频率是0.57C等位基因的频率是0.65 遗传漂变、种群结构和强选择是推动偏离零的力量。在大种群中预测随时间呈指数衰减轨迹返回到零就像在小种群中遗传漂变下的杂合性一样 其中r是感兴趣的位点对之间预期的重组分数。这可以用来估计单倍型的年龄。最后即使对于在不同染色体上独立分配的位点也需要时间衰减。哈代-温伯格基因型可以在一代中恢复但过去的事件对LD有持续影响这可以用来推断更远的过去的过程如种群中不再存在的种群结构。 当从实际数据集计算时双重杂合子是不明确的。假设我们有一个个体的C/TA/G SNP集合。C等位基因与第二个位置的A还是G等位基因相关联通常我们不知道。但是不明确的单倍型频率为我们提供了关于解决双重杂合子可能方式的信息。如果C-G单倍型非常常见而C-A单倍型很少见那么这表明C/TA/G个体可能具有C-G/T-A单倍型。使用这种方法计算太繁琐无法手工完成。 幸运的是这正是EM算法发挥作用的地方。Kalinowski和Hedrick2001使用大角羊Ovis canadensis数据集Boyce等人1997来估计LD。这个物种很罕见样本量很小所以我们需要从可用的数据中获得尽可能多的信息。 以下R代码实现了Kalinowski和Hedrick2001给出的方程式。它从猜测相等的单倍型频率和 0开始。然后它更新这个猜测并迅速达到最大似然解 ≈0.0779和na-B单倍型频率基本为零。 # 定义各个复合基因型的频率 AABB - 2 # 两个位点都是纯合子AABB的个体数量 AaBB - 0 # 一个位点是杂合子另一个是纯合子AaBB的个体数量 aaBB - 0 # 一个位点是纯合子aa另一个是纯合子BB的个体数量 AABb - 0 # 第一个位点是纯合子AA第二个位点是杂合子Bb的个体数量 AaBb - 1 # 两个位点都是杂合子AaBb的个体数量双杂合子 aaBb - 0 # 第一个位点是纯合子aa第二个位点是杂合子Bb的个体数量 AAbb - 1 # 第一个位点是纯合子AA第二个位点是纯合子bb的个体数量 Aabb - 0 # 第一个位点是杂合子Aa第二个位点是纯合子bb的个体数量 aabb - 0 # 两个位点都是纯合子aabb的个体数量# 使用上述输入运行函数‘Dcalc’ Dcalc(AABB, AaBB, aaBB, AABb, AaBb, aaBb, AAbb, Aabb, aabb) 与用于衡量线性相关性的统计相关系数皮尔逊“r”有关。让我们使用ℛ表示相关系数以避免将其与重组分数r混淆。2除以所有等位基因频率的乘积等于ℛ^2。 此外奇怪的是如果我们将ℛ2乘以采样的染色体总数如果我们观察n个二倍体个体通常为2n那么我们就会得到一个具有一个自由度的2统计量 然而这并不令人惊讶。在如此小的样本量下即使LD非常强检测偏差的能力也非常有限。最后我们希望指出EM算法是一种“爬山”算法它找到一个局部最大似然峰值。可能存在其他峰值可以使用MCMC方法来处理这个问题并更全面地探索复杂的似然表面。