《Journal of Integrative Plant Biology》：Structural variation drives rhizome innovation and adaptive divergence in sister Medicago species
编辑推荐：
野生多年生姐妹物种Medicago archiducis-nicolai（具根茎/高山）和M. ruthenica（无根茎/旱生）构成牧草改良的重要遗传资源。为解码其对比性状和生境适应的基因组基础，研究人员构建了染色体规模的基因组组装，对128个个体进行了重测
野生多年生姐妹物种Medicago archiducis-nicolai（具根茎/高山）和M. ruthenica（无根茎/旱生）构成牧草改良的重要遗传资源。为解码其对比性状和生境适应的基因组基础，研究人员构建了染色体规模的基因组组装，对128个个体进行了重测序，在冷/热胁迫下分析了转录组，并功能验证了因果等位基因。研究人员证明，结构变异（SVs）——特别是基因重复——是根茎形成和高山/旱生适应的主要驱动力。此外，非编码调控区域中普遍存在的存在缺失变异（PAVs）支撑了控制根茎发育和胁迫反应的等位基因特异性表达分化。关键的是，这些调控性PAVs在性状发育和胁迫适应过程中诱导了对比性的表达模式。研究人员的发现揭示了一种双重机制，即编码和调控SVs协同调控姐妹物种的表型创新和生态特化，为豆科植物进化研究和苜蓿育种提供了宝贵的基因组资源。
**研究背景与问题**
物种适应性分化驱动生物多样性和生态系统韧性，其遗传基础涉及基因重复、同源基因功能分化及谱系特异性创新等多个机制。然而，现有研究常孤立分析这些机制，缺乏对它们在生态特化和物种形成中协同作用的整体理解。姐妹物种由于共享近期共同祖先和相似遗传背景，是解析这些机制相对贡献的理想模型。本文聚焦于两个姐妹物种——高山具根茎的Medicago archiducis-nicolai和旱生无根茎的M. ruthenica，它们分布于青藏高原高海拔和半干旱西伯利亚-蒙古-华北地区，具有显著的形态和适应性差异（如根茎有无、耐寒/耐热）。尽管先前研究涉及单系群或特定性状（如形态变异、花期、胁迫响应），但不同遗传机制对它们适应性分化的集体贡献仍不清楚。此外，两者作为栽培苜蓿（M. sativa）的野生近缘种，具有重要农艺价值，可用于牧草改良和抗逆育种。为此，本研究通过高质量基因组组装、群体重测序、转录组分析和功能验证，旨在揭示驱动这两个物种性状和适应性分化的基因组因素。
**研究方法**
研究人员主要采用以下关键技术方法：
（1）染色体水平基因组组装：利用PacBio HiFi测序和Hi-C染色质构象捕获技术，分别对M. archiducis-nicolai（来自甘肃兴隆山）和M. ruthenica（来自同一地区低海拔）进行基因组组装，并辅以Illumina短读长抛光。
（2）比较基因组与进化分析：基于16个物种（包括14个豆科物种和拟南芥）的单拷贝直系同源基因构建系统发育树，估算分化时间；通过同义替换率（Ks）和四倍简并位点（4DTv）分析全基因组重复事件；利用CAFé5分析基因家族扩张/收缩。
（3）转录组分析：对两个物种不同组织（根、根茎分节、茎芽）以及冷（4°C）和热（40°C）胁迫下1、3、6、12小时叶片样本进行RNA-seq，鉴定差异表达基因（DEGs），并采用加权基因共表达网络分析（WGCNA）识别根茎发育和热胁迫响应模块。
（4）结构变异（SV）鉴定：结合基于组装（MUMmer/SyRI）和基于HiFi读长映射（NGLMR/Sniffles2）的方法，在全基因组范围内鉴定SV；再对128个个体（89个M. archiducis-nicolai和39个M. ruthenica，包括3个同域种群）进行重测序，利用DELLY和LUMPY检测群体水平的SV，最终筛选物种固定SV。
（5）功能验证：对候选基因OSBP和TLP所在的物种固定PAV进行RT-qPCR表达验证和双荧光素酶报告基因（LUC）分析，以评估PAV对转录活性的调控效应。
**研究结果**
**基因组测序、组装与注释**
M. archiducis-nicolai基因组大小为940.03 Mb，contig N50达105.98 Mb，99.78%序列锚定到8条假染色体；M. ruthenica基因组经优化后大小为915.09 Mb，contig N50提升至26.26 Mb（较先前版本提高41.7倍，缺口数减少96.1%）。两个基因组质量评估均达标（BUSCO >93%，LAI>29，QV>32）。TE注释显示LTR反转录转座子占主导（约59%），且均在最近1百万年内扩张。
**系统发育关系与物种分化**
基于16个物种的系统发育树表明，M. archiducis-nicolai与M. ruthenica为姐妹物种，分化时间约2.23百万年前（更新世初期）。两个物种的基因家族扩张/收缩富集于不同通路：高原适应物种富集于“种子成熟调控”、“根发育调控”、“DNA损伤修复”等；旱生物种富集于“生长素运输”、“根毛起始”、“次生代谢”等。
**影响基因组大小的因素**
无近期全基因组重复（WGD）事件发生，但LTR-RT（特别是Gypsy）长度与基因组大小高度相关（R2=0.91），是基因组大小差异的主要驱动力。此外，检测到M. archiducis-nicolai第4号染色体发生TE介导的不平衡易位事件，该区域在M. ruthenica中富含物种特异性基因并富集于胁迫响应通路。
**基因重复对环境适应的贡献**
两个物种的重复基因类型相似（WGD约26%，TE介导的重复约45%-49%），但正选择基因中约43%为物种特异性。共享正选择重复基因富集于通用胁迫应答（如“UV/冷/热响应”），而物种特异性重复基因则显示分化适应：M. archiducis-nicolai富集于缺氧响应、DNA修复、磷饥饿；M. ruthenica富集于盐/渗透胁迫响应、热响应、黄酮醇生物合成、先天免疫。
**根茎生长的遗传基础**
M. archiducis-nicolai的根茎特异性基因（534个）中，27%为物种特有，富集于细胞周期调控、发育可塑性、DNA损伤修复和胁迫响应。WGCNA分析鉴定到一个蓝色模块（532个DEGs）与根茎尖端表达高度相关，功能涉及向性、根形态建成、环境感知和细胞增殖。该模块的20个枢纽基因（如bZIP29、PELPK1、XTH8等）在根茎尖端特异性高表达，经RT-qPCR验证。
**热胁迫的遗传分化**
冷胁迫诱导两个物种相似的DEG时间动态（12小时内递增），但热胁迫下M. archiducis-nicolai早期（1小时）快速响应（1,445个特有DEGs），而M. ruthenica晚期（12小时）响应强烈（1,053个特有DEGs）。尽管多数DEGs（86.9%-93.1%）在另一物种中存在直系同源基因，但仅8.8%-27.7%具有保守表达模式，表明直系同源基因功能广泛分化。物种特异性基因（占6.9%-13.1%）富集于温度相关功能。共表达网络分析进一步揭示两个物种应激响应枢纽基因和网络拓扑结构完全分化。
**全基因组结构变异**
共检测到98,439个SVs，其中90.7%为存在缺失变异（PAVs），且约99%位于非编码区。5,485个高效应PAV与3,849个基因相关，富集于根形态建成、DNA损伤修复、水稳态、免疫响应等。约38%的热胁迫差异表达直系同源基因和44%的根茎特异性直系同源DEGs含有PAVs。
**PAV驱动直系同源基因功能分化**
通过群体分析（128个个体）鉴定出5,071个物种固定PAV。以OSBP基因（编码氧甾醇结合蛋白）为例，M. archiducis-nicolai在第一内含子中含1,728 bp序列（冷胁迫下增强转录活性），而M. ruthenica中该序列缺失，导致冷诱导表达消失。以TLP基因（类甜蛋白）为例，M. ruthenica下游71 bp缺失（热胁迫下增强转录活性），而M. archiducis-nicolai中该缺失固定，导致热诱导表达缺失。双荧光素酶实验证实这两种PAV分别增强冷/热胁迫下的转录活性，且RT-qPCR验证了表达模式。
**讨论与结论**
本研究表明，结构变异（包括基因重复、物种特异性基因、PAV介导的直系同源基因调控分化）共同驱动了姐妹物种的适应性分化。基因重复（尤其是TE介导的重复）提供了适应性进化的原材料，物种特异性基因在根茎发育和热胁迫响应中发挥关键作用，而PAV通过改变顺式调控元件重塑直系同源基因的表达模式，从而产生环境特异性适应。例如，OSBP内含子缺失导致冷诱导表达丧失，可能影响高山物种的耐寒性；TLP下游缺失赋予旱生物种热诱导表达。同域种群中PAV的杂合状态暗示可能通过杂交进行适应性渐渗。这些发现重新定义了近缘物种的适应性进化机制，强调了编码和调控SVs的协同作用。研究为苜蓿属作为进化基因组学模型提供了宝贵资源，并为利用野生遗传资源进行苜蓿抗逆分子育种提供了实践途径。论文结论：结构变异，包括基因重复、物种特异性基因以及PAV与直系同源基因调控分化的相互作用，共同驱动了两个姐妹物种的适应性分化。这些发现重新定义了近缘物种的适应性进化，强调了这些遗传机制在快速生态位分割中的协同作用。作为苜蓿的近缘种，本研究产生的基因组资源和机制见解不仅增强了苜蓿属作为进化基因组学模型的价值，而且为通过靶向分子育种提高苜蓿胁迫耐受性提供了实用途径。
