种群历史限制下猕猴桃重复遗传适应由祖先 standing genetic variation 主导

《The Plant Journal》：Repeated genetic adaptation under demographic limits in kiwifruit

　　

群体遗传理论(population genetic theory)预测，与广布近缘物种相比，地理分布狭窄的物种经历更强的遗传漂变(genetic drift)、更高的近交(inbreeding)以及更大的突变负荷(mutational burden)，然而种群历史(demographic history)如何塑造适应的基因组基础(genomic basis of adaptation)仍知之甚少。研究人员假设：(i) 分布受限(range-restricted)物种由于长期有效种群规模(effective population size, N_e)较小，表现出遗传多样性降低、近交升高和突变负荷(mutation load)增大；(ii) 尽管存在这些种群约束(demographic constraints)，共同环境压力仍可驱动重复适应(repeated adaptation)；(iii) 重复适应可能通过祖先 standing genetic variation、适应性 introgression(adaptive introgression)或独立的从头突变(de novo mutations)产生。研究人员利用两个猕猴桃物种 Actinidia hemsleyana（窄域）和 Actinidia eriantha（广布）检验这些假设，二者在种群历史和分布范围上差异显著，但在中国亚热带占据重叠环境。研究人员表明，种群收缩(demographic contraction)和近交已侵蚀分布受限物种的全基因组多样性，并伴随显著更高的有害突变负荷(deleterious mutation load)。尽管种群历史遗留(demographic legacies)对比鲜明，两个物种均对共享气候梯度表现出重复基因组响应，适应靶向相似基因。这种重复适应主要由祖先 standing variation 驱动，适应性 introgression 贡献较小，而独立 de novo mutations 的证据有限。这些结果揭示，种群约束塑造全基因组多样性和突变负荷，但当共享遗传变异仍可获得时，并不排除重复适应。

论文解读：种群历史限制下猕猴桃的重复遗传适应

研究背景：群体遗传理论预测，地理范围狭窄或有效种群大小较小的物种会受到更强遗传漂变、更高近交和更大突变负荷影响。然而，种群历史如何塑造适应的基因组基础仍不清楚，尤其小且隔离的种群能否通过清除有害等位基因维持或恢复适应潜力，还是累积突变负荷会长期限制进化，仍是进化与保护生物学的核心未解问题。尽管种群限制可侵蚀遗传多样性并提高突变负荷，当物种经历相似选择压力时，适应可能仍沿可重复基因组路径发生。重复适应可由独立突变、祖先 standing genetic variation（祖先种群中已存在的遗传变异）或适应性 introgression（通过基因流传播适应性等位基因）产生，但这些机制可留下相似基因组特征，区分困难，尤其对分布重叠且种群历史不对称的近缘物种。比较群体基因组可检验这些理论预测。研究人员比较两个同域猕猴桃物种，窄域 Actinidia hemsleyana 和广布 Actinidia eriantha，二者种群历史和分布范围差异显著，但在中国亚热带占据重叠环境。

研究人员开展了比较群体基因组研究，结论是窄域 A. hemsleyana 遗传多样性降低、近交升高、有效种群较小且有害突变负荷更高；但两物种仍对共享气候梯度表现重复适应，靶向相似基因，且主要由祖先 standing genetic variation 驱动，基因流贡献较小，独立 de novo mutations 证据有限。这表明种群约束塑造全基因组多样性和突变负荷，但当共享遗传变异存在时并不排除重复适应，对理解进化可重复性、种群限制下的适应架构和保护具有意义。

关键技术方法概括：研究人员采集中国浙江、福建、江西11个同域自然分布区的205个 A. hemsleyana 和 A. eriantha 个体，并采集10个 Actinidia lanceolata 个体作为外群；进行全基因组重测序。主要分析包括：读取比对与SNP calling；种群结构推断（最大似然树、ADMIXTURE、PCA）；遗传多样性与近交估算（θ_π、θ_W、N_e、ROH与F_ROH）；遗传负荷预测（est-sfs、SnpEff、Grantham、SIFT4G、R_X/Y）；环境关联分析（GF、LFMM、RDA）；选择清除检测（SweeD、OmegaPlus）；DMC模型区分重复适应机制；GO富集分析。

结果：
1. 种群历史预测对比鲜明的多样性与近交（Demographic history predicts contrasting diversity and inbreeding）：通过种群结构和多样性分析，两个物种形成分离遗传簇，确认尽管同域但长期分化。窄域 A. hemsleyana 的全基因组多样性低于广布 A. eriantha，θ_π为2.54×10^?3对3.76×10^?3，θ_W为1.43×10^?3对2.35×10^?3，N_e为27 908对45 863。多个100-kb滑窗支持该差异。A. hemsleyana 杂合度降低，F_ROH近两倍，且杂合度随近交增加而下降。

2. 种群侵蚀转化为不对称突变负荷（Demographic erosion translates into asymmetric mutation load）：通过蛋白编码基因SNP及多种有害突变预测，A. hemsleyana 携带更高比例的衍生有害等位基因。标准化R_X/Y分析显示，A. hemsleyana 强有害变异过量，而轻度有害变异未富集，表现为长期近交下部分清除。

3. 种群约束下的重复适应（Repeated adaptation despite demographic constraint）：通过GF和GEA分析，降水与太阳辐射是两物种主导选择轴。两物种共享2092个气候相关候选基因，富集于胁迫响应、激素信号和生殖发育，显示广泛基因重用。SNP水平仅2个适应性SNP共享，说明重复适应主要通过基因重用而非相同核苷酸位点。SweeD和OmegaPlus识别近期正选择，458个基因位于跨物种重叠选择清除区，富集于生物或外部刺激响应、防御和向性。与GEA交集得到41个环境适应候选基因，富集于ABA激活信号、细胞应激响应和固有免疫调节；这些候选基于统计关联，需实验验证。

4. 重复适应主要源于祖先变异（Repeated adaptation arises primarily from ancestral variation）：DMC模型比较显示，近90%候选基因支持祖先 standing genetic variation，迁移模型解释9.8%，独立 de novo mutations仅2.4%。表明重复适应主要通过预存变异重用，而非独立产生新适应变异。

讨论总结：窄域 A. hemsleyana 低多样性、高近交符合小N_e和长期种群下降，与其他窄域与广布近缘植物类似；因两物种采样环境重叠，差异不太可能反映局部生态，而凸显种群历史。突变负荷方面，A. hemsleyana 强有害变异负担更高，符合长期种群收缩下纯化选择效能降低；近交增加纯合性可能暴露弱有害等位基因并导致部分清除，但强有害变异仍通过漂变累积。保护上，该窄域物种更易受栖息地丧失、随机事件和未来气候变化影响，建议人工杂交以降低近交衰退风险并防止遗传多样性继续下降，同时需基因组脆弱性评估。重复适应方面，尽管种群背景差异显著，两物种仍表现适应靶标重叠；窄域 A. hemsleyana 环境关联基因更多（11 950对4832），广布 A. eriantha 共享候选比例更高（43.3%对17.5%），符合基因流影响适应架构的理论预期。降水与太阳辐射一致重要，41个候选包括钙依赖蛋白激酶、蛋白磷酸酶2Cs和MYB转录因子，参与耐旱和氧化应激，但功能需实验验证。祖先 standing variation 与基因流促进重复适应：重复适应主要由祖先 standing variation 驱动，基因流次要；两物种分化约9.1 Ma（95% CI: 9.0–10.7 Ma），基因流约0.98 Ma（95% CI: 0.48–0.99 Ma），较晚基因流提供祖先多态性持续并独立选择的时间窗。近缘杂交和遗传连通性增加共享适应性等位基因可能性。预存等位基因可立即供选择且可能已经过滤，适应更快更可重复。总体，猕猴桃重复适应反映共享祖先和持续基因流维持的预存变异被选择性重用。

结论翻译：通过对比两个具有分化种群历史但环境重叠的近缘猕猴桃物种，本研究揭示种群约束如何塑造进化潜力而不消除适应重复性。与群体遗传理论一致，窄域 A. hemsleyana 较广布 A. eriantha 表现出遗传多样性降低、近交升高和有害突变负担更重，反映长期有效种群大小和纯化选择效能差异。尽管存在这些约束，两物种对共享气候梯度表现重复基因组响应，适应涉及与降水和太阳辐射相关的相似基因和通路。因此，种群限制改变适应的基因组架构，而非完全阻止适应。结果进一步表明，重复适应主要由祖先 standing variation 驱动，基因流有额外贡献，突出共享祖先和进化连通性的重要性。总体，这些发现显示漂变、突变负荷和共享变异如何相互作用塑造种群约束下的适应结果。然而，这些发现仅基于单一物种对，未来研究应纳入更多物种对，尤其涵盖更广泛种群对比，以确定所观察模式是否广泛普遍。