重组染色体上竞争性有益突变间的干扰

《Biology》：Interference of Competing Beneficial Mutations on Recombining Chromosomes

　　

在基因组中寻找选择性清除（selective sweep）的信号是当前群体基因组学（population genomics）的主要目标，因为这能够估算正在走向固定的有益突变速率，并识别参与选择的基因。重复选择性清除模型传统上假设，在正常重组率的染色体区域中，最多只有一个有益等位基因正在走向固定。研究人员在此回顾并扩展了关于紧密连锁有益突变间干扰的理论研究，表明该假设可能被违背。研究人员证明，即使在正常重组率的染色体区域，有益突变间的干扰也可能导致固定时间显著延长。此外，研究人员讨论了如何通过分析选择性清除的基因组足迹，在群体基因组学研究中检测干扰，并在已发表的数据集中寻找干扰的经验证据。

发表于《Biology》的这篇论文针对传统选择性清除模型的一个关键假设——正常重组率区域中每次只有一个有益等位基因走向固定——提出质疑。研究背景方面，有益突变在自然群体中走向固定的速率和时间是群体遗传学的核心问题；然而，两个同时出现并在固定过程中相互作用的有利突变如何共同固定，研究相对较少。Otto和Barton研究过第二个突变有利度低于第一个的情形，Cuthbertson等、Bossert和Pfaffelhuber则考虑了第二个突变更有利、且两个突变可通过重组形成最高适应度重组型的情况。Stephan用确定性常微分方程模型分析了两个竞争性有益突变，本文在此基础上进行扩展，并将固定时间与选择性清除的基因组特征联系起来，以探讨如何在真实数据中检测干扰。研究人员通过理论分析和经验数据重分析得出，两个紧密连锁的有益突变即使在正常重组率区域也可能彼此干扰，显著延长重组型的固定时间；这一发现对解释基因组中选择性清除信号、估计有益突变固定速率具有重要意义。

研究方法上，研究人员采用理论建模与实证数据分析相结合的方式。理论部分建立双座位选择-重组常微分方程（ordinary differential equations, ODEs）模型，利用非归一化变量近似求解两个竞争性有益突变A和B的固定时间；随后将理论预测与Chevin等的Monte Carlo模拟（基于Wright-Fisher模型）进行比较，分析选择性清除的基因组足迹。实证部分重新分析了欧洲黑腹果蝇（Drosophila melanogaster）群体polyhomeotic座位的多态数据，以及内布拉斯加州deer mice（鹿鼠）Agouti座位的数据；后者包括沿330 km南-北样带采集的11个群体样本，其中Sand Hills上（ON）群体271个个体、Sand Hills外（OFF）群体266个个体，并使用了复合似然比（composite likelihood ratio, CLR）检验、HapFLK分化统计、Tajima‘s D等群体遗传学方法。

研究结果按原文小标题概括如下：

2.1 模型：通过双座位模型，假设两个座位上的有益等位基因A和B具有加性适应度，忽略上位效应，并考虑座位间重组率r。利用ODE描述等位基因频率变化，区分“相斥”（repulsion）与“相引”（association）两种突变起源背景。结论是，该模型能够近似描述两个竞争性有益突变在有限群体中的频率动态。

2.2 固定时间：通过求解固定时间方程，得到相斥情形下重组型AB的固定时间T_rep和相引情形下的T_assoc。结果表明，固定时间依赖于B引入时A的拷贝数或频率；当该值较小时，两个有益突变间的干扰最强，固定时间显著延长；重组率越低，固定时间越长。相引情形下的固定时间较短，且随引入频率增加趋于渐近值。

3 利用选择性清除检测竞争性有益突变之间的干扰：将理论预测与Chevin等的模拟结果比较后发现，在两个被选择座位之间的中性座位上，当A的引入频率较低时，Tajima‘s D呈强正值，表明中间频率变异过多，这是干扰导致的特征性基因组足迹。因此，通过分析选择性清除的基因组足迹可以检测竞争性有益突变间的干扰。

4 寻找有益突变间干扰的证据：在黑腹果蝇欧洲群体的polyhomeotic区域，检测到五个紧密连锁的选择性替换，但这些位点附近的Tajima’s D并未升高，提示它们可能不是独立、依次固定的，而是以单倍型块形式从非洲祖先群体共同迁入，是一个有价值的阴性结果。在deer mice的Agouti座位中，CLR检验在Exon 1A附近检出两个紧密连锁的选择性清除，但Tajima‘s D未升高；而在Intron 1内的Region 2，HapFLK分化的最高峰与强连锁不平衡（linkage disequilibrium, LD）和高Tajima’s D共定位，且存在两个由紧密连锁位点组成的簇，在ON与South-OFF群体间表现出较大的频率差异，提示可能存在干扰；Region 1中频率差异位点较少，Tajima‘s D升高不明显。

讨论部分指出，模型假设加性适应度并忽略上位效应，该假设是否限制结论的普适性尚不清楚；确定性模型在强选择条件下有效，与Bossert和Pfaffelhuber的随机模型结果部分一致，但重组参数仅出现在确定性结果中。利用黑腹果蝇的选择参数估算，每条染色体上平均约有15.6个强选择替换可在第一个突变固定期间出现，但需要满足紧密连锁且第一个突变仍处于低频的条件。两个经验案例都涉及群体结构，基因流可能是干扰之外的替代解释，未来研究应予以考虑。

结论部分翻译如下：对于两个紧密连锁的有益突变A和B，研究人员推导了重组单倍型AB固定时间的公式，假设B比A具有更高适应度，并且B在A固定之前出现。两个突变间的干扰即使在正常重组率区域也可能显著延长AB的固定时间。研究人员解释了如何利用选择性清除研究中估计的选择系数和选择性替换频率，在DNA多态性数据中检测干扰证据。研究人员讨论了两个案例：欧洲黑腹果蝇的polyhomeotic座位虽有五个紧密连锁的选择性替换，但周围Tajima‘s D较低，表明这些突变并非独立演化，因此是有价值的阴性结果；deer mice的Agouti座位也检出两个紧密连锁的选择性清除，但周围Tajima’s D未升高，同样未显示干扰特征。另一方面，Sand Hills上与Sand Hills外群体间的遗传分化分析提示干扰证据：在Intron 1的一个区域中，HapFLK的最高峰与强连锁不平衡和高Tajima‘s D共定位。