关于《古代杂交与组合机制：铁线莲属灌木铁线莲组适应性进化的驱动力》研究的学术报告
本研究由北京林业大学森林资源高效生产全国重点实验室的何健、李明阳、肖佳敏、罗一轲、李文赫、吴桓宇、罗月欣和谢磊共同完成。研究成果于2025年9月18日被接受，发表于植物学领域知名期刊《New Phytologist》（2026年，第250卷，第1884-1904页）。论文题为“Ancient hybridization and combinatorial mechanism: drivers of adaptive evolution in Clematis sect. Fruticella (Ranunculaceae)”。这项研究综合运用系统基因组学、生态位建模和群体遗传学方法，深入探究了铁线莲属中适应干旱的灌木铁线莲组及其近缘类群的进化起源与快速多样化机制，为理解杂交在植物适应性辐射和物种形成中的作用提供了关键案例。
一、 学术背景与研究目的
本研究的科学领域主要涉及进化生物学、系统基因组学、植物生态学和生物地理学。传统物种形成理论通常认为物种形成是一个缓慢、渐进的过程，依赖于新突变的积累。然而，越来越多的证据表明，种间杂交（interspecific hybridization）可以作为一种“组合机制”（combinatorial mechanism），通过重组已有的亲本遗传变异，快速产生新的基因型和表型组合，从而驱动快速的适应性进化和物种形成。尽管这一理论框架日益受到重视，但在植物类群中，尤其是基于大规模基因组数据验证该机制驱动适应性进化的实证研究仍然相对匮乏。
铁线莲属的灌木铁线莲组（Clematis sect. *Fruticella*）为此提供了一个理想的研究系统。该组植物为适应干旱的灌木，分布于中国北部和青藏高原东北缘的干旱、半干旱地区，其内部物种在形态上高度相似，但在短短约70万年内迅速分化并适应了不同的干旱生境，暗示了快速的适应性进化。此前基于叶绿体基因组和核糖体DNA的研究揭示了该组与近缘类群（如黄花铁线莲组 sect. *Meclatis*、准噶尔铁线莲C. songorica 以及厚萼铁线莲类群 the Otophora group）之间存在显著的细胞核-细胞器系统发育冲突（cyto-nuclear discordance），暗示了可能存在复杂的杂交历史。
因此，本研究旨在：（1）利用大规模核基因数据厘清灌木铁线莲组及其近缘类群的系统发育关系；（2）探究基因树与物种树之间不一致性的原因，特别是检测古代杂交事件；（3）验证“组合机制”是否以及如何驱动了灌木铁线莲组和黄花铁线莲组的适应性进化与物种形成。
二、 详细研究流程
本研究是一个多步骤、整合性的系统基因组学研究，其工作流程严谨而复杂，主要包括以下环节：
1. 样品采集与基因组数据收集： * 研究对象与样本量： 研究共采集了64份铁线莲属样品，涵盖10个组的40个物种。重点包括了灌木铁线莲组所有已知物种、准噶尔铁线莲、黄花铁线莲组以及厚萼铁线莲类群等关键类群，每个物种包含多个个体。以独叶草属（*Anemoclema*）作为外类群。 * 数据生成方法： 采用两种高通量测序技术。通过RNA测序（RNA-seq）获取核基因转录组数据（平均每个样本约6GB原始数据），用于核基因系统发育和杂交分析；通过基因组浅层测序（genome skimming）获取完整的叶绿体基因组数据，用于构建质体系统发育树并与核基因树进行比较。
2. 数据组装与正交同源基因识别： * 转录组组装： 使用Trinity软件对每个样本的RNA-seq原始数据进行从头组装，预测开放阅读框，并通过BUSCO评估组装完整性，确保数据质量可用于后续分析。 * 叶绿体基因组组装与注释： 使用GetOrganelle或Geneious Prime软件从头组装叶绿体基因组，并进行注释。 * 核单拷贝正交同源基因筛选： 这是构建可靠系统发育树和进行杂交分析的基础。研究团队使用ProteinOrtho软件对所有样本的蛋白编码序列进行聚类分析，筛选出在所有样本中至少存在70%的单拷贝正交同源基因。随后，通过比对去除可能来源于细胞器的基因，并利用TreeShrink软件检测并剔除基因树中具有异常长分支的序列，以排除因错误组装或深度溯祖等因素造成的干扰。最终，获得了一个包含6536个核单拷贝正交同源基因的高质量数据集。
3. 系统发育重建与分化时间估算： * 分析方法： 为了获得可靠的物种树，研究采用了多种方法进行对比和验证。包括基于串联法的最大似然法（使用RAxML软件），以及基于溯祖模型的物种树推断方法，如两步汇总法（使用ASTRAL软件）和全溯祖模型法（使用StarBEAST3软件）。鉴于串联法和两步汇总法在存在强烈基因树不一致时可能产生偏差，研究最终选择StarBEAST3构建的全溯祖物种树作为后续分析的框架。 * 分化时间估算： 在StarBEAST3框架内，采用两步化石标定策略，共估算了铁线莲属主要分支的分化时间。
4. 基因树不一致性与杂交检测： * 不一致性量化： 使用PhyParts软件量化每个物种树节点上一致与不一致基因树的比例，并通过模拟分析（使用Phybase软件）评估不完全谱系分选（incomplete lineage sorting, ILS）能否单独解释观察到的基因树冲突以及细胞核-质体不一致性。 * 杂交检测与可视化——核心创新方法： 这是本研究的方法学亮点。研究团队主要使用HyDe软件包检测杂交信号。为了解决HyDe在推断古老杂交事件和量化亲本贡献方面的挑战，研究团队自主开发了一个名为“VisualHyDe”的Python可视化脚本。该脚本包含两种模式： * “叶节点模式”（Leaf Mode）： 为每个疑似杂交个体生成一个热图，展示其与所有潜在亲本对之间的杂交显著性及遗传贡献比例（γ值）。 * “节点模式”（Node Mode）： 这是用于推断祖先杂交事件的关键创新。 该算法通过分析一个杂交支系后代所有个体的“叶节点模式”热图，寻找跨后代谱系一致且稳定的杂交信号。只有当对于某一特定亲本对，其后代内部节点间显示出高度一致的γ值差异（即“差异共识”）时，才将该杂交信号归因于它们最近共同祖先（MRCA）发生的一次古代杂交事件，并估算祖先的遗传贡献比例。这极大地增强了对古老杂交推断的稳健性。 * 假阳性过滤流程： 为了排除因谱系速率异质性（如长枝吸引）导致的假阳性信号，研究设计了一套严谨的过滤流程，包括零假设模拟、自举重抽样、内部拓扑结构随机化以及使用不同外类群进行敏感性测试。 * 补充验证： 同时使用PhyloNet软件推断系统发育网络，作为HyDe分析结果的交叉验证。
5. 干旱相关基因的拓扑偏向性与分化时间分析： * 目的： 为了验证杂交是否特异性地引入了适应性基因，研究团队从6536个基因中筛选出与干旱胁迫响应相关的基因集。 * 拓扑偏向性分析： 比较干旱相关基因集与非干旱相关基因集在支持物种树拓扑（(（Otophora群，目标组），*C. songorica*））与支持杂交拓扑（(（目标组，*C. songorica*），Otophora群））之间的比例差异。显著更高的杂交拓扑比例表明干旱相关基因更可能来源于*C. songorica*。 * 分化时间分布分析： 分别估算干旱相关和非干旱相关基因集中，目标组（如sect. *Meclatis*）与*C. songorica*之间的分化时间。如果存在从*C. songorica*的基因渗入，预计干旱相关基因的分化时间分布会出现一个更近的次要峰值。
6. 生态位比较分析： * 方法： 收集所有研究类群的分布点数据（共647条记录）和19个生物气候变量。通过主成分分析（PCA）在环境空间中可视化各类群的生态位，并计算Schoener’s D生态位重叠指数，量化各类群与其推测亲本以及内部各支系之间的生态位差异。
7. 染色体倍性确定： * 方法： 使用流式细胞术（FCM）测定关键类群的DNA相对含量，并结合传统的染色体计数法（对准噶尔铁线莲等样本进行染色体计数），以确定所检测到的杂交事件是否为同倍体杂交（homoploid hybridization）。
三、 主要研究结果
1. 系统发育关系与分化时间： 基于6536个核基因构建的全溯祖物种树高度支持灌木铁线莲组（sect. *Fruticella*）和黄花铁线莲组（sect. *Meclatis*）各自为单系群。Clematis songorica 是灌木铁线莲组的姐妹群。黄花铁线莲组与厚萼铁线莲类群（Otophora group）构成姐妹群。分化时间估算表明，灌木铁线莲组的冠部年龄约为130万年（更新世早期），其与*C. songorica*的分化时间约为310万年（上新世-第四纪交界）。黄花铁线莲组与*Otophora*群的分化时间约为370万年（上新世中期），其冠部多样化发生在约230万年前（第四纪早期）。
2. 广泛的基因树不一致性与杂交证据： PhyParts分析显示，在包含sect. Fruticella, sect. Meclatis, C. songorica 和 Otophora 群的大支系基部，存在极高的基因树冲突。细胞核-质体系统发育树之间存在显著的不一致性（例如，C. songorica 在质体树中部分样本嵌入sect. *Meclatis*或sect. *Fruticella*）。模拟分析表明，仅凭不完全谱系分选无法解释这些冲突（P值显著），强烈暗示了历史杂交事件的存在。
3. 古代杂交起源的检测与量化： 应用VisualHyDe的“节点模式”分析，研究获得了强有力的证据： * 灌木铁线莲组（sect. *Fruticella*）具有古代杂交起源，其祖先是干旱适应的*C. songorica*（贡献了约73%的基因组）与中生适应的*Otophora*群祖先（贡献了约27%的基因组）杂交的产物。 * 黄花铁线莲组（sect. *Meclatis*）同样具有古代杂交起源，亲本与前者相同，但遗传贡献比例相反：主要来源于*Otophora*群祖先（约80%），少量来源于*C. songorica*（约20%）。 * 这一结果得到了PhyloNet网络分析、严格的假阳性过滤以及多项敏感性测试（自举分析、拓扑随机化、更换外类群）的稳健支持。此外，研究还检测到组内更近期的杂交事件，如C. canescens 可能是C. fruticosa 和 C. tomentella 的杂交后代。
4. 干旱相关基因的偏向性遗传： 对黄花铁线莲组的分析为适应性杂交提供了关键证据： * 拓扑偏向性： 干旱相关基因集显著更多地支持黄花铁线莲组与C. songorica 聚类的杂交拓扑结构（与非干旱基因集相比，P值极显著），表明干旱适应性基因有偏向地从干旱适应的亲本C. songorica 渗入到黄花铁线莲组。 * 分化时间分布： 干旱相关基因集中，黄花铁线莲组与C. songorica 的分化时间呈现明显的双峰分布，其中一个较近的峰值与杂交发生的时间相符；而非干旱相关基因集则主要呈现一个较古老的分化峰值。这进一步证实了干旱相关基因是通过较近的杂交事件而非古老的共同祖先遗传获得的。
5. 生态位分化与遗传贡献的相关性： 生态位分析显示，两个杂交起源的组（sect. Fruticella 和 sect. *Meclatis*）与其任何一个亲本类群的生态位重叠度都很低（D值分别为0.24和0.17，与另一亲本无重叠），表明它们已占据新的生态位。更重要的是，在组内，各支系的生态偏好与其从干旱亲本*C. songorica*继承的遗传比例呈显著相关。例如，在灌木铁线莲组中，分布于更干旱的中国北部地区（支系5和6）的物种，比分布于相对湿润的横断山区（支系4）的物种，从*C. songorica*获得的遗传贡献高出约5%。黄花铁线莲组内也观察到类似模式。这提示，杂交后的谱系在多样化过程中，其生态适应性的细微差异可能继续受到亲本遗传贡献比例调控。
6. 倍性水平： 流式细胞术和染色体计数结果表明，所有检测的关键类群（包括亲本和杂交后代）均为二倍体，证实了所检测到的杂交事件为同倍体杂交。
四、 研究结论与意义
本研究得出结论：灌木铁线莲组和黄花铁线莲组均起源于一次古老的同倍体杂交事件，杂交亲本为适应极端干旱的准噶尔铁线莲（*C. songorica*）和适应中生环境的厚萼铁线莲类群（Otophora group）的祖先。 这一杂交事件通过“组合机制”将两个亲本已有的、经过自然选择检验的遗传变异（如干旱适应性和特定的花部形态）重新组合，快速产生了兼具亲本特征但又占据新生态位的新谱系。
科学价值： 1. 提供了“组合机制”驱动植物适应性进化的强有力实证案例： 研究不仅检测到古代杂交信号，还通过关联遗传贡献比例、特定功能基因（干旱相关基因）的偏向性渗入以及生态位分化，完整地勾勒出杂交如何通过重组现有变异促进快速适应和物种形成的路径。 2. 阐明了铁线莲属复杂系统发育关系的成因： 研究揭示了铁线莲属内广泛的细胞核-质体不一致性主要由频繁的杂交和质体捕获导致，解决了该属系统学中的长期争议。 3. 展示了杂交在植物快速辐射演化中的重要作用： 研究表明，同倍体杂交在铁线莲属从晚中新世到第四纪的快速多样化中扮演了关键角色，为理解其他植物类群的快速辐射提供了新视角。 4. 方法学贡献： 开发的VisualHyDe工具及其配套的假阳性过滤流程，为从基因组数据中稳健推断古老杂交事件和量化亲本贡献提供了新的方法学框架。
五、 研究亮点