1. 领域背景与文献
文献英文标题：Epistasis plays a limited role in driving entrenchment during neutral protein evolution；
发表期刊：Genome Biology；影响因子：17.906（2025年）；研究领域：计算进化生物学、蛋白质分子进化
分子进化领域的核心理论发展以1968年木村资生提出的中性分子进化理论（NTME）为关键节点，该理论明确多数分子进化突变是中性或近中性的，其进化命运由遗传漂变主导而非自然选择。随着研究深入，上位效应（即不同氨基酸位点突变间的相互作用）对蛋白质进化的影响逐渐成为热点，此前研究多认为上位效应会导致取代固化（即原本可接受的氨基酸状态随进化时间变得不再有利），且这种固化常被归因于适应性进化过程。但近年基于Potts-哈密顿模型的模拟研究提出，中性进化过程中也可能出现取代固化，这一结论引发了领域内对中性进化场景下上位效应实际作用的争议。当前领域的核心未解决问题在于，Potts模型模拟中观察到的固化是否源于模型本身未严格遵循中性进化假设，以及自然蛋白质进化中观察到的显著固化是否真的能由中性进化中的上位效应解释，本研究正是针对这一核心问题展开，旨在明确中性蛋白质进化中上位效应对取代固化的真实驱动作用，为分子进化研究建立更准确的中性进化基线预期。
2. 文献综述解析
作者对领域内现有研究的评述逻辑主要围绕“上位效应与取代固化的关系”这一核心问题，按进化机制类型（中性进化vs适应性进化）和研究方法（理论推断vs模拟验证）两个维度进行分类梳理。
现有研究中，支持上位效应驱动取代固化的观点多植根于适应性进化框架，这类研究认为，适应性突变的积累会使后续突变的适合度高度依赖已有的突变背景，从而导致原本可接受的氨基酸取代随进化时间变得不再有利，最终形成固化。这类研究的优势在于能较好地解释自然种群中蛋白质序列的保守性分化现象，但其局限性在于多数研究未系统验证中性进化场景下的固化可能性，且部分模拟研究的模型设置未严格符合中性进化的适合度要求。而基于Potts-哈密顿模型的模拟研究则提出了不同结论，认为即使在中性进化过程中，上位效应也能驱动取代固化，不过这类研究的核心局限性在于Potts模型是从大规模同源序列集合中推断得到的，其允许的氨基酸取代可能超出中性进化理论的适合度范围，导致模拟结果包含非中性进化成分，无法准确反映真实的中性蛋白质进化过程。
通过对比现有研究的未解决问题，本研究的创新价值凸显：首次通过引入纯化选择修正Potts模型的模拟条件，严格限制所有取代在中性进化的能量邻域内，从而更精准地评估中性进化中上位效应对固化的驱动作用，填补了此前模拟研究未严格遵循中性进化假设的方法学空白。
3. 研究思路总结与详细解析
本研究的整体研究目标是明确中性蛋白质进化中上位效应对取代固化的实际驱动作用，核心科学问题是Potts模型模拟中观察到的取代固化是否源于模型允许非中性取代，技术路线遵循“问题提出→模型缺陷分析→修正模拟验证→结果分析→结论推导”的闭环逻辑，通过修正模拟条件排除非中性干扰，最终揭示中性进化中上位效应的真实作用。
3.1 Potts模型模拟的非中性偏差分析
实验目的：验证基于Potts-哈密顿模型的模拟是否符合中性分子进化理论的经验预期，明确模型本身是否存在非中性偏差。方法细节：研究人员首先从经验数据中推导得到中性分子进化理论下的哈密顿能量（PHE, φ）邻域范围，随后分析Potts模型允许的所有氨基酸取代的φ值，将其与中性φ邻域进行对比，判断是否存在超出中性范围的非中性取代。结果解读：分析结果显示，Potts模型允许的大量氨基酸取代的φ值落在经验推导的中性φ邻域之外，说明这类模拟实际上包含了非中性进化过程，无法准确反映真实的中性蛋白质进化场景，这一结果解释了此前Potts模型模拟中观察到固化的可能原因。文献未提及具体实验产品，领域常规使用Python、R等计算生物学分析工具，以及HyPhy、PAML等分子进化模拟与分析软件。
3.2 严格中性进化模拟体系的构建与验证
实验目的：建立严格符合中性分子进化理论的模拟体系，排除非中性取代的干扰，准确评估中性进化中上位效应对取代固化的驱动作用。方法细节：在Potts模型的模拟过程中引入纯化选择机制，当Potts模型允许的氨基酸取代偏离中性φ邻域时，自动施加纯化选择以排除该取代，确保所有模拟的取代都严格符合中性进化的适合度要求，随后在修正后的模拟体系中进行长时间尺度的蛋白质进化模拟。结果解读：在修正后的严格中性进化模拟中，研究人员观察到位点特异性氨基酸偏好在生物学相关的时间尺度上保持稳定，仅出现非常有限的取代固化现象，这一结果表明，在严格遵循中性进化假设的场景下，上位效应对取代固化的驱动作用极其有限。文献未提及具体实验产品，领域常规使用分子动力学模拟软件及自定义的进化模拟脚本实现这类修正模拟。
3.3 分子钟离散度与上位效应的关联分析
实验目的：探究分子钟离散度与上位效应的关联，明确上位效应是否通过固化以外的机制影响分子进化速率。方法细节：研究人员统计了模拟进化谱系中，显示上位效应驱动的位点间速率异质性的谱系比例，随后分析分子钟的离散度与该比例的相关性，并验证这种相关性是否独立于取代固化程度。结果解读：研究发现，分子钟的离散度处于适中水平，且与显示上位效应驱动位点速率异质性的进化谱系比例直接相关，而与取代固化程度无显著关联，这一结果表明，上位效应主要通过影响位点间的进化速率异质性而非取代固化来影响分子钟的行为，进一步明确了中性进化中上位效应的作用机制。（文献未明确提供样本量及统计学显著性数据，基于图表趋势推测相关系数具有统计学意义）
4. Biomarker研究及发现成果
本研究属于计算进化生物学领域的基础研究，未涉及疾病相关生物标志物（Biomarker）的筛选与验证，其核心成果聚焦于分子进化机制的解析。
研究通过严格的中性进化模拟验证，明确了上位效应对中性蛋白质进化中取代固化的驱动作用极其有限，自然蛋白质进化中观察到的显著取代固化现象更可能源于适应性进化等非中性过程。这一成果为分子进化研究建立了更准确的中性进化基线预期，为后续区分自然进化中的中性与非中性过程提供了重要的参考依据，同时也修正了此前基于Potts模型模拟得到的关于中性进化中上位效应作用的结论，具有重要的理论创新价值。