「错误阈值」:基因组长度的天花板,是复制保真度写死的——以及为什么高突变率下「更差」的序列反而会赢
【锚点】这一篇只钉一个判断:Eigen 的【错误阈值】不是「突变太多会出错」这种废话,而是一条可算的【信息容量预算线】——基因组长度 L 与每位点突变率 u 之间存在硬约束 L·u ≲ ln σ。而它的反面比它本身更反直觉:当 u 逼近阈值时,赢的不是适应度最高的序列,而是【适应度景观最平】的那一个。
【显式】Eigen 1971(Naturwissenschaften 58:465)的准种模型。设主序列每位点复制保真度 q,基因组长度 L,则整基因组保真度 Q = qL;设主序列相对其突变谱的平均适应度优势为 σ。选择-突变平衡下,主序列在稳态分布中的占比在 Q·σ = 1 处发生【离域】(delocalization):分布铺满整个序列空间,主序列携带的信息被抹平。令 u = 1-q 且 u 很小,临界条件就是
u_crit ≈ ln σ / L ⇔ L_max ≈ ln σ / u
【隐式】这条式子一读就出结论:σ 只进对数,几乎不动;所以【基因组长度上限基本被 u 反向钉死】。代入真实量级:RNA 病毒每位点 u ~ 10⁻⁴(Qβ 噬菌体 4220 nt、脊髓灰质炎病毒 7500 nt 都在这档),ln σ 取 1–10,给出 L_max 10⁴——而 RNA 病毒基因组确实卡在 3×10³ 到 3×10⁴ 之间。DNA 生物 u ~ 10⁻⁹–10⁻¹⁰,同一个不等式给出 L_max ~ 10⁹–10¹⁰——细菌与真核基因组正落在这一档。这不是巧合,是同一个不等式在两端各钉了一次。
【流程节点链路】起点:RNA 依赖的 RNA 聚合酶缺少 3’→5’ 外切酶校对活性 → 中间节点①:每位点错误率 u 被钉在 10⁻⁴ 量级 → 中间节点②:整基因组保真度 Q = qL 随 L 指数衰减 → 中间节点③:当 Q·σ < 1 时选择-突变平衡解离域,主序列信息无法维持 → 终点:L 被反向约束在 ln σ / u,RNA 病毒基因组无法越过 10⁴。整条链上没有任何一步是「因为病毒小所以基因组小」——因果方向是反的。
【细微】真正反直觉的是【存活最平者】(survival of the flattest)。Wilke、Wang、Ofria、Lenski、Adami 2001(Nature 412:331)用数字演化实验证明:高突变率下,一个【峰值适应度更低、但突变邻域更宽】的序列能击败峰值更高、邻域陡峭的序列;Bull、Meyers、Lachmann 2005 的准种综述把这一点系统化。机制上这几乎是【率失真】的翻版:适应度峰的高度与宽度是一次权衡,而 u 决定了你该站在这条权衡曲线的哪一点。低 u 时买高度,高 u 时买宽度。
【极难发现】由此暴露一个盲点:错误阈值之所以「锐利」,是两个建模假设的产物——【无限种群】与【单峰适应度景观】。有限种群下,离域不是相变而是平滑漂移;空间结构、重组、区室化都会把阈值推走甚至拆掉。所以「错误阈值禁止 X」这类断言,通常是拿一个无限种群极限去断言有限种群的事实。可证伪的用法只有一个方向:把突变率推到阈值以上应当导致【种群崩溃】(致死性突变发生)——Crotty、Cameron、Andino 2001(PNAS 98:6895)用利巴韦林在脊髓灰质炎病毒上直接测到了 error catastrophe。但反方向——「因为没崩溃所以阈值不存在」——不成立,因为有限种群本来就会把锐利阈值磨平。
【结论】错误阈值给的不是一条禁令,而是一条【预算线】:你能维持多少比特的遗传信息,取决于复制机器每复制一个比特要付多少错误。DNA 的出现不是「更高级」,而是一次把 u 从 10⁻⁴ 压到 10⁻⁹ 的【信道升级】——它买来的不是适应度,是【基因组长度】。
【待验证点】如果错误阈值本质是率失真约束,它应当能推广到非遗传的信息维持系统:语言、代码库、制度。一个系统的「突变率」(每次传递的失真概率)是否也钉死了它能维持的规则条数上限?这个类比最先在哪里断掉——是在「选择」这一环,还是在「复制」这一环?







0 Comments