《The Plant Journal》：Haplotype-resolved autotetraploid genome of Dioscorea nipponica provides insights into recent polyploidization, metabolic rewiring, and environmental adaptation
编辑推荐：
全基因组加倍（WGD，或多倍化）是植物进化的重要驱动力，然而与异源多倍体相比，同源多倍体的基因组、转录组和代谢后果仍了解较少。穿龙薯蓣（Dioscorea nipponica），一种产薯蓣皂苷元型甾体皂苷的药用物种，同时拥有二倍体和同源多倍体种群，为研究这些过
全基因组加倍（WGD，或多倍化）是植物进化的重要驱动力，然而与异源多倍体相比，同源多倍体的基因组、转录组和代谢后果仍了解较少。穿龙薯蓣（Dioscorea nipponica），一种产薯蓣皂苷元型甾体皂苷的药用物种，同时拥有二倍体和同源多倍体种群，为研究这些过程提供了独特体系。本研究呈现了同源四倍体穿龙薯蓣（2n = 4x = 40）的单倍型分型染色体级基因组组装。比较基因组学和群体基因组学分析表明，其起源为单一、近期的同源四倍体事件（<1.63 Ma），发生在其二倍体祖先分化（~2.5–2.9 Ma）之后，与上新世–更新世气候变迁时间吻合，后者可能促进了种群隔离与基因组加倍。四个单倍型表现出高的平均核苷酸一致性（average nucleotide identity）、广泛的共线性和相似的全局基因表达，伴随普遍的等位基因丢失（约32%的基因）和广泛的等位基因转录抑制（约60%的表达四联体），共同缓冲了剂量效应。代谢组学分析揭示了根茎中的代谢权衡：与二倍体相比，同源四倍体积累了更高的黄酮类化合物，但更低的薯蓣皂苷元型甾体皂苷，这可能增强了生殖适合度及其对边缘环境的适应能力。这一转变得到了黄酮类基因家族（如PAL、CHS）扩张的支持，其中一些表现出正剂量效应，而薯蓣皂苷元途径基因（如CYP90、CYP94）拷贝数稳定且在四倍体中大多表现出负剂量效应。选择性清除分析进一步鉴定出与基因组稳定性（如AtTOP3α、AtDMC1）和非生物胁迫响应相关的基因，这对多倍体的建立至关重要。总之，这些发现为同源多倍体如何塑造基因表达、代谢进化和生态适应提供了新见解，并对该药用物种的保护与利用具有重要意义。
研究背景与科学问题
多倍化，即全基因组加倍（Whole-Genome Duplication, WGD），是被子植物进化、物种形成与适应的主要驱动力。现存约30%-35%的开花植物为多倍体。根据染色体来源，多倍体分为由种内基因组加倍产生的同源多倍体（autopolyploid）和由种间杂交后染色体加倍形成的异源多倍体（allopolyploid）。历史上，同源多倍体被认为远不如异源多倍体普遍，且常被假定存在更大的内在劣势，导致其在进化和基因组研究中长期被忽视。然而，大规模染色体调查显示同源多倍体在被子植物中广泛分布，其对植物进化的贡献可能被系统性低估。尽管已知多倍化常伴随表型和生理变化，包括次级代谢改变及对冷、旱、盐等非生物胁迫耐受性的增强，且基因剂量效应被认为是重要机制之一，但增加基因剂量本身并不能完全解释多倍化的生态和代谢后果。多倍化常引发复杂的基因组和调控重组，包括偏向性基因保留或丢失、重复基因亚功能化或新功能化以及等位基因特异性表达，这些过程能以途径或分支特异的方式重塑代谢网络。然而，多倍体如何在基因组、转录和代谢水平上重编程代谢，以及这些变化是否赋予广泛的适合度优势或促进情境依赖的适应，仍有待探索。
穿龙薯蓣（Dioscorea nipponica）是一种多年生草本药用植物，其干燥根茎“穿山龙”是重要的传统中药，药效成分主要为薯蓣皂苷元型甾体皂苷。该物种自然种群中存在二倍体、四倍体、六倍体和八倍体等多种细胞型，且多倍体被推断为同源多倍体。值得注意的是，同源多倍体主要分布于物种分布区的边缘地带（高海拔、高纬度或寒冷气候），而二倍体则常见于中心区域，暗示同源多倍体可能增强了适应能力和生态成功。此前，薯蓣属的研究多集中于二倍体物种的基因组和薯蓣皂苷生物合成途径，而对同源四倍体与二倍体在相近遗传背景下的直接比较有限，限制了对基因组加倍直接后果的理解。特别是，尚不清楚同源多倍化是否导致次级代谢物普遍增加，还是驱动向特定化合物类别（如黄酮类或苯丙烷类）的代谢再分配，以及这种潜在代谢转变背后的基因组和转录组重组机制、自然同源多倍体群体的起源和适应机制也均不明确。
为解决上述问题，研究人员生成了同源四倍体穿龙薯蓣的染色体级、单倍型分型基因组，并整合了转录组、代谢组和群体基因组分析，旨在阐明WGD如何重塑基因表达和次级代谢，尤其是薯蓣皂苷生物合成与其他代谢物之间的潜在权衡；揭示自然同源四倍体群体的群体历史起源；以及其环境适应的遗传基础。
主要技术方法
研究人员采用流式细胞术和k-mer分布分析估计基因组大小；利用PacBio HiFi长读长测序结合高通量染色质构象捕获（Hi-C）技术，通过hifiasm和CPhasing流程完成单倍型分型的染色体级基因组组装；整合从头预测、同源比对和RNA-seq证据进行基因注释；运用OrthoFinder进行直系同源基因鉴定，RAxML重建最大似然系统发育树，MCMCTree估算分歧时间；采用AlleleFinder鉴定等位基因并进行等位基因特异性表达分析；对二倍体和四倍体根茎进行非靶向代谢组学分析；基于来自西南地区4个同源四倍体群体（31个个体）和13个二倍体群体（91个个体的全基因组SNP数据，利用fastsimcoal2推断群体历史和基因流；通过交叉群体复合似然比（XP-CLR）、遗传分化指数（FST）和核苷酸多样性（π）比值进行选择性清除分析。
主要研究结果
1. 基因组测序、倍性推断与染色体级组装
研究人员的基因组调查和多种倍性分析方法（GenomeScope、Smudgeplot、nQuire）一致支持穿龙薯蓣四倍体为同源四倍体。最终组装获得40条假染色体（10个同源染色体组，每组4个单倍型），总锚定序列1.96 Gb，BUSCO完整性达96.8%-98.1%，显示出高质量和高连续性。
2. 基因组注释
重复序列约占各单倍型基因组的66.89%-67.38%，以LTR逆转录转座子为主。各单倍型预测的蛋白编码基因数量高度相似（23,372-23,684个），基因结构特征高度保守，功能注释率达97.42%-97.51%。
3. 高度基因组共线性与近期同源四倍体化
四个单倍型之间及与二倍体穿龙薯蓣之间均展现出广泛的共线性，平均核苷酸一致性高达97.00%-98.92%。同义替代率（Ks）分析显示，四倍体单倍型间存在一个非常近期的峰值（Ks≈0.013），对应单倍型分歧时间约1.63百万年前（Ma），而四倍体单倍型A与二倍体之间的Ks峰值（~0.023）对应分歧时间约2.88 Ma，表明二倍体谱系的分歧先于随后的同源四倍体化事件。
4. 系统基因组学与基因家族进化
系统发育分析证实穿龙薯蓣与盾叶薯蓣（D. zingiberensis）为姐妹关系，分歧时间约8.34 Ma。四倍体单倍型形成单系分支，聚结年龄约1.74 Ma。在同源四倍体支系中鉴定出151个扩张基因家族，显著富集于环境适应、氧化还原酶活性、酰基转移酶活性以及黄酮类、黄酮和黄酮醇、异黄酮生物合成通路。
5. 广泛的等位基因丢失与组织特异性等位基因表达偏差
在全基因组范围内鉴定出26,190个基因，平均每个基因3.33个等位基因。约32%的基因发生了等位基因丢失。全局等位基因表达水平在四个单倍型间高度相似。在每个表达四联体中，约22%-26%表现为平衡表达，而大多数（55%-60%）表现为等位基因抑制。薯蓣皂苷生物合成途径基因存在广泛的等位基因丢失和转录沉默，表明选择性等位基因调控而非均匀的基因剂量主导着该途径。
6. 多倍化重塑了穿龙薯蓣根茎的代谢物谱与基因表达
非靶向代谢组学分析鉴定出60种差异丰度代谢物。四倍体根茎中黄酮类化合物普遍富集，而薯蓣皂苷元型甾体皂苷则一致减少，两者呈强负相关，表明存在协调的代谢再分配。转录组分析分别在叶片、茎和根茎中鉴定出大量上调（正剂量效应）和下调（负剂量效应）基因，呈现出组织特异性的转录重编程。
7. 二倍体和四倍体穿龙薯蓣中薯蓣皂苷与黄酮类生物合成基因的差异拷贝数与剂量依赖性表达
薯蓣皂苷途径关键基因（如CYP90、CYP94）的拷贝数在二倍体和四倍体间保持稳定，而黄酮类途径关键基因（如PAL、CHS、F3'H）在四倍体基因组中显著扩张。转录水平上，多个薯蓣皂苷途径基因表现出负剂量效应，构成生物合成的瓶颈；而黄酮类途径的多个基因则表现出正剂量效应，共同促进了黄酮类化合物的净积累。
8. 群体基因组学揭示同源四倍体穿龙薯蓣的起源
系统发育、ADMIXTURE和主成分分析一致将同源四倍体个体划分为独立于二倍体的单系群，且四倍体内部存在东西方向的微弱遗传分化。群体历史推断显示，四倍体有效种群规模（Ne）远小于二倍体，其与二倍体祖先的分歧时间约为2.51 Ma，东西两组四倍体随后在约0.12 Ma发生分化，并在晚更新世存在近期基因流。
9. 同源四倍体建立与适应的遗传基础
选择性清除分析的交集鉴定出155个高置信度候选基因，包括维持基因组稳定性的基因（如AtTOP3α、AtDMC1）、响应干旱、盐、冷和UV-B等非生物胁迫的基因，以及参与开花、分生组织生长等发育过程的基因。一个约1 Mb的染色体区间密集地包含了此类基因，成为选择作用于基因组稳定性、胁迫适应和发育调控的热点区域。
总结与讨论
研究人员首次呈现了薯蓣属基部分支——根茎组（Sect. Stenophora）的首个多倍体染色体级基因组。该组装质量优异，四个单倍型在序列、结构和表达上高度相似。同源四倍体通过约32%的等位基因丢失和约60%的表达四联体中等位基因抑制等机制，实现了全局表达平衡和剂量缓冲。综合Ks分析、系统发育和群体遗传学推断，该同源四倍体起源于一次相对较晚的单一事件（上限约1.63 Ma），发生在二倍体谱系分化之后，且与上新世-更新世气候转型期重叠，这一时期的地貌破碎化和严酷环境可能促进了种群隔离和基因组加倍。代谢组学揭示了根茎中显著的代谢权衡：黄酮类化合物积累增加，而薯蓣皂苷元型甾体皂苷减少。这一转变的分子机制在于黄酮类基因家族的拷贝数扩张和正剂量效应，以及薯蓣皂苷途径基因的拷贝数稳定和负剂量效应。这种代谢转向可能增强了四倍体对干旱、寒冷等边缘环境的适应能力，与其主要分布在物种分布区边缘的格局相符。选择性清除分析进一步鉴定了与基因组稳定性（确保减数分裂正常进行）和非生物胁迫响应相关的关键基因，这些基因受到的选择压力很可能促进了同源四倍体的建立和生态成功。综上，这些发现为理解同源多倍体如何塑造基因表达、代谢进化和生态适应提供了全新见解，并对这一重要药用物种的保护、利用和育种具有重要指导意义。
