岩石缝隙中的生命史诗：基因组揭示岩生银莲花属植物的物种分化、濒危机制与喀斯特适应奥秘
第一，研究的基本信息 本研究由四川大学生命科学学院，生物资源与生态环境教育部重点实验室的谢登峰（Deng-feng Xie）、张义阳（Yi-yang Zhang）、蔡竞（Jing Cai）、成瑞瑜（Rui-yu Cheng）、王渊（Yuan Wang）、谭进波（Jin-bo Tan）、余岩（Yan Yu）、何兴金（Xing-jin He）和周颂东（Song-dong Zhou）*（通讯作者）等人合作完成。研究论文题为《岩生银莲花属植物（毛茛科）在石灰岩环境中的物种形成、濒危与适应》，于2025年9月9日在线发表于期刊《植物多样性》（Plant Diversity）。论文共同作者单位还包括北京大学前沿交叉学科研究院和兰州大学生态学院。
第二，研究的学术背景 本研究的核心科学领域为进化生物学、保护基因组学和植物适应性进化。中国南方的喀斯特石灰岩地区是生物多样性的热点和特有现象中心，其环境以高钙、高镁、高pH值、低保水能力和营养贫瘠为特征，对植物构成了强大的选择压力，孕育了许多特有物种，是研究物种分化与适应性进化的“天然实验室”。然而，之前的研究多集中于比较喀斯特与非喀斯特物种，或探讨单个物种的适应机制，对于在喀斯特环境中近缘物种之间的分化过程和适应性进化仍存在认知空白。
岩生银莲花属（*Urophysa*）是毛茛科的一个特有属，仅包含两个物种：岩生银莲花（*U. rockii*）和尾囊草（*U. henryi*），二者均仅分布于中国南方的石灰岩地区。U. rockii 仅局限分布于四川盆地西部，野生个体稀少（约2000株），被列为极危（Critically Endangered， CR）物种；而 U. henryi 分布范围相对较广但也呈片段化。这两个物种在花部形态上存在关键差异（U. rockii 花瓣具距，而 U. henryi 花瓣呈囊状），且具有高度的生境特异性（专性生长于石灰岩峭壁）和异域分布格局，使其成为研究喀斯特环境中物种分化、适应性进化以及濒危机制的理想模型。
因此，本研究的主要目标是：1) 解析岩生银莲花属两个物种的进化历史与物种分化过程；2) 评估极度濒危的 U. rockii 的遗传现状并探究其濒危的遗传学原因；3) 在全基因组水平上识别并验证驱动岩生银莲花属适应高钙等恶劣喀斯特环境的关键基因和分子机制。
第三，详细的研究流程 本研究采用了一套整合了从头基因组测序、群体重测序、生理实验和转录组分析的综合性研究方案，工作流程严谨而系统。
程序一：岩生银莲花（U. rockii）染色体水平基因组的组装与注释 研究对象为一个采自四川省彭州市的野生 U. rockii 个体。 - 测序与组装：研究团队利用PacBio Sequel II平台产生了覆盖度超过70×的高保真（HiFi）长读长序列，并利用Hi-C技术获得了超过230×覆盖度的染色体构象捕获数据。基因组大小通过k-mer分析预估为约289.03 Mb。使用Hifiasm软件进行初步组装，随后使用Illumina短读长数据通过Pilon软件进行三轮纠错，得到草图基因组。最后，利用Hi-C数据，通过比对和聚类，将重叠群（contig）挂载到7条假染色体上，这与该物种已知的染色体数目（2n=14）一致。 - 基因组评估与注释：最终组装的基因组大小为303.21 Mb，Contig N50达到20.75 Mb。BUSCO评估显示完整度高达97.60%，表明组装质量很高。研究预测了32,080个蛋白质编码基因，并发现基因组中重复序列占比为48.80%。此外，还鉴定了824个 U. rockii 特有的基因家族。
程序二：系统发育、分化时间与基因家族进化分析 为了明确 U. rockii 在植物系统发育中的位置及其进化历史，研究人员从公共数据库下载了包括其他毛茛科植物（如耧斗菜属）在内的11个物种的基因组数据作为比较。 - 分析过程：使用OrthoFinder鉴定了1009个单拷贝直系同源基因（SCGs）。基于这些基因，分别采用串联法和溯祖法构建系统发育树，两者结果一致。通过分子钟模型并借助化石标定点，估算关键节点的分化时间。使用WGDI软件分析全基因组复制事件，使用CAFE软件分析各进化支系上基因家族的扩张与收缩。 - 结果目标：旨在确定 U. rockii 与近缘物种的分化时间，探究其是否经历了特殊的基因组事件，并了解其基因家族动态，为后续的适应性研究提供进化背景。
程序三：岩生银莲花属物种的群体重测序与遗传变异分析 为了研究群体层面的遗传结构、多样性和历史动态，研究对来自两个物种共14个自然种群的97个个体进行了全基因组重测序。 - 样本与测序：包括 U. rockii 的5个种群45个个体，以及 U. henryi 的9个种群52个个体。同时，以 Semiaquilegia adoxoides 和4个耧斗菜属个体作为外类群。每个个体的平均测序深度目标为25×。 - 数据处理与变异检测：将高质量的测序reads比对到新组装的 U. rockii 参考基因组上。使用GATK流程进行单核苷酸多态性（SNP） calling，并经过严格过滤（如保留基因型质量≥20、双等位SNP、过滤低深度和低频位点等），最终获得了一个高质量SNP数据集用于后续分析。
程序四：群体遗传结构、多样性及历史分析 基于高质量的SNP数据，研究团队进行了多层次的分析： - 群体结构：使用ADMIXTURE软件进行祖先成分分析，并通过交叉验证误差确定最佳K值。同时进行主成分分析（PCA）以可视化遗传结构。 - 系统发育关系：基于SNP数据，分别使用最大似然法（IQ-TREE）和溯祖法（MP-EST, ASTRAL-III）构建种群和物种水平的系统发育树。 - 遗传多样性与近交：计算了核苷酸多样性（π）、Tajima’s D值来衡量群体遗传多样性及历史动态。通过分析纯合片段（ROHs）的长度和数量，并计算基因组近交系数（FROH）和个体近交系数（FIS），来评估近交程度。 - 种群历史动态：应用PSMC模型，基于单个个体的基因组序列推断各主要支系有效群体大小（Ne）随时间的变化历史。更进一步，使用fastsimcoal2软件，基于多位点联合等位基因频谱（joint SFS），对预设的多个包含物种分化顺序、时间和基因流情况的模型进行模拟和拟合，通过AIC准则选择最优的群体历史模型，并估计分化时间、历史群体大小和基因流速率等参数。 - 生态位与隔离分析：利用MAXENT软件进行生态位建模，预测当前、历史（末次间冰期LIG、末次盛冰期LGM）和未来的潜在分布区。使用Mantel检验分析遗传分化（FST）与地理距离、环境距离之间的相关性，检验“距离隔离”和“环境隔离”效应。
程序五：喀斯特适应性基因的基因组扫描与鉴定 为了在全基因组范围内寻找与喀斯特环境适应相关的基因，研究采用了基于群体分化的扫描策略。 - 分析方法：将适应喀斯特的岩生银莲花属个体（97个）作为“喀斯特组”，将非喀斯特的外类群（耧斗菜等）作为“非喀斯特组”。在50 kb的滑动窗口（步长10 kb）中，计算两组间的遗传分化指数（FST）以及喀斯特组内相对于非喀斯特组内的核苷酸多样性比值（π比值）。将FST值最高的5%窗口与π比值最高和最低的5%窗口（分别代表喀斯特组多样性极低和极高的区域）取交集，以识别经历了强烈正选择或分化的基因组区域。对这些区域内包含的基因进行GO和KEGG功能富集分析。此外，使用PAML软件的支位点模型对这些候选基因进行正选择验证。
程序六：候选基因的钙胁迫响应与功能验证 鉴于钙离子是喀斯特环境的核心胁迫因子，研究通过水培实验和转录组测序来验证候选基因的功能。 - 实验设计：将 U. rockii 幼苗置于霍格兰营养液中，设置三个钙浓度处理：无钙（0 mM）、中钙（15 mM）和高钙（30 mM）。处理5天后，分别收集根和叶组织。 - 转录组分析：对每个处理三个生物学重复的样本进行RNA-seq。将测序reads比对到 U. rockii 基因组，进行基因表达定量和差异表达分析。鉴定在钙处理下差异表达的基因（DEGs），并与程序五中鉴定到的正选择基因（PSGs）、程序二中发现的 U. rockii 特有基因以及扩张基因家族进行整合分析，筛选出强候选的适应性基因。 - 综合选择信号分析：为了进一步确认这些强候选基因在基因组上受到的选择信号，研究计算了包含这些基因的基因组区域的多种群体遗传学统计量，包括遗传多样性（π）、Tajima’s D、复合似然比（CLR）和整合单倍型评分（iHS），最后将这些统计量整合为一个综合的DCMS统计量，以识别受选择的“热点”区域。
第四，研究的主要结果 基因组特征与进化背景：组装获得了高质量的 U. rockii 染色体水平基因组（303.21 Mb，7条染色体）。系统发育分析确认其与耧斗菜属亲缘关系最近，分化时间约为600万年前。基因组分析未检测到该属特有的全基因组复制事件，但发现 U. rockii 有754个基因家族发生扩张，这些家族显著富集于离子通道功能、钙离子结合、刺激响应等通路，暗示了其对特殊环境的潜在适应。
物种分化与群体历史：群体遗传结构分析将97个个体清晰地分为两大支，对应两个物种（K=2）。进一步细分（K=4）揭示了每个物种内部存在两个地理分化的支系：U. rockii 的西支（RW）和东支（RE），以及 U. henryi 的西支（HW）和东支（HE）。系统发育树和PCA结果支持这一结构。fastsimcoal2模拟的最优模型表明，两个物种约在595万年前分化，这与晚中新世横断山脉的剧烈隆升时期吻合。随后，在更新世气候波动的影响下，各支系的有效群体大小发生剧烈波动并经历瓶颈。模型还显示，物种分化后存在不对称的基因流，但大约在7.5万年前已完全中断。生态位建模显示，两个物种的分布区在历史上（如末次间冰期）更为广阔，后经历收缩和片段化，未来气候变化预计将导致其分布区进一步萎缩，尤其是对 U. rockii 威胁巨大。Mantel检验表明，在 U. henryi 内部及两个物种间，遗传分化与地理距离、环境距离均呈显著正相关，支持地理隔离和生态隔离共同促进了物种分化。
*U. rockii* 的濒危遗传现状：U. rockii 的整体遗传多样性（π = 0.0027）显著低于 *U. henryi*（π = 0.0055）。其西支的两个种群（CXS和PZ）遗传多样性更是低至惊人的水平（~8.7×10⁻⁵ 和 9.6×10⁻⁵），是目前报道的濒危植物中最低的之一。PSMC和fastsimcoal2分析均显示，U. rockii 的种群在历史上经历了长期、严重的瓶颈收缩。连锁不平衡衰减更慢，表明其经历了强烈的遗传漂变和/或选择。更重要的是，ROH分析显示 U. rockii 个体，尤其是西支种群，携带大量长片段纯合片段，其基因组近交系数（FROH）和个体近交系数（FIS）均极高，表明严重的近交是当前种群的普遍现象。长期的瓶颈效应、受限的基因流、高度的生境专化以及严重的近交，共同导致了 U. rockii 极低的遗传多样性和极高的灭绝风险。
喀斯特适应性机制的关键发现： 1. 细胞壁完整性是首要防线：差异表达分析发现，钙胁迫下大量与细胞壁组织、生物合成和修饰相关的基因被激活。一个编码细胞壁关联受体激酶（WAK3）的候选基因显示出强烈的正选择信号。这表明维持细胞壁的完整性是应对高钙等物理化学胁迫的第一道关键屏障。 2. 钙离子稳态的核心调控网络：基因组扫描（top 5% FST 和 π比值）鉴定出257个候选正选择基因，功能富集于钙离子结合、跨膜转运、刺激响应等通路。转录组整合分析进一步聚焦于三个关键基因：TPC1（液泡钙离子释放通道）、CAX3（液泡钙离子/质子逆向转运蛋白）和PTR/POT（肽/氮转运蛋白）。 - 表达模式：在中等钙浓度下，TPC1 和 PTR/POT 表达上调，可能促进钙离子的吸收和利用；而在高钙胁迫下，它们表达下调，以避免细胞毒性。相反，CAX3 在所有钙处理下均持续上调，其功能是将细胞质中过量的钙离子泵入液泡储存，是维持细胞质钙稳态的核心机制。 - 选择信号：对这三位基因所在基因组区域的多重选择信号扫描（高FST、高CLR、高iHS、低Tajima‘s D）均显示强烈的正选择信号，证实了它们在适应高钙环境中的关键进化作用。 3. 离子吸收、排泄与区室化的综合策略：研究还发现了其他一系列与离子调控相关的候选基因，形成了完整的适应网络： - 离子外排：如ACA11、SLAH1等基因参与将多余离子排出细胞。 - 细胞器区室化：除了CAX3介导的液隔离子储存，MRS2F基因负责将镁离子转运至线粒体，实现离子的亚细胞分隔，减轻细胞质毒性。 - 离子感应与信号传导：如KIC、CML等钙感应蛋白基因，负责监测胞内钙浓度并启动相应的信号通路。
第五，研究的结论与意义 本研究通过构建首个岩生银莲花属染色体水平基因组，并结合群体基因组学、生态学和生理转录组学等多学科手段，系统揭示了该属植物在喀斯特石灰岩环境中的进化历史、濒危机理和适应性分子机制。
科学价值： 1. 阐明了喀斯特近缘物种的分化模式：证实了横断山脉晚中新世的隆升和更新世的气候波动是驱动岩生银莲花属物种分化和异域成种的主要动力，生态位分化在其中发挥了重要作用，为理解喀斯特生物多样性热点地区的物种形成机制提供了典型案例。 2. 剖析了极度濒危物种的遗传学成因：首次在基因组层面详细揭示了 U. rockii 濒危的多重遗传因素——历史瓶颈、地理隔离、生境专化与严重近交的协同效应，为评估和保护其他具有类似生活史的濒危喀斯特特有植物提供了重要的方法论和参考依据。 3. 揭示了植物适应高钙胁迫的整合性分子网络：超越了对单个基因的认识，描绘了一个包含“细胞壁加固（第一道防线）”、“钙离子稳态核心调控（TPC1/PTR/POT/CAX3）”、“离子吸收/外排调节”和“亚细胞区室化储存”的多层次、网络化的适应性策略，深化了对植物适应钙质土壤等逆境环境的分子机制的理解。
应用价值： 1. 为 U. rockii 的保护提供直接指导：研究明确指出应对遗传多样性极度匮乏的龙门山种群（如CXS和PZ）实施优先保护，建议采取建立保护区、禁止人为干扰、进行人工辅助迁地保护和种群间人工授粉以恢复基因流、开展人工扩繁和回归引种等具体措施。 2. 为喀斯特脆弱生态系统的生物多样性保护提供借鉴：研究强调了喀斯特生境的脆弱性和物种对气候变化的敏感性，警示需加强对此类生态系统中特有种群的长期监测和前瞻性保护。
第六，研究的亮点 1. 研究对象的独特性与模型价值：首次以专性生长于喀斯特峭壁、形态特异、分布区狭窄且呈异域分布的近缘姊妹种为研究对象，完美集成了物种分化、适应性进化和濒危机制三大科学问题于一个体系内进行研究。 2. 多组学与多学科方法的深度整合：研究流程设计精妙，从基因组从头测序、群体重测序（历史动态、选择扫描），到生理胁迫实验与转录组分析（功能验证），再到生态位建模，形成了一个从“进化历史推断”到“现状评估”再到“机制解析”的完整证据链，逻辑严密，说服力强。 3. 发现了关键的适应性调控模块：不仅鉴定了大量候选基因，更通过表达模式和选择信号分析，精准定位了TPC1、CAX3和PTR/POT这三个基因构成的、在表达上受钙浓度精细调控的核心稳态模块，并阐明了其协同工作机制，是机制研究上的一个重要突破。 4. 对濒危机制的深入定量揭示：将群体历史模拟（PSMC, fastsimcoal2）、遗传多样性统计与基于ROH的近交分析相结合，定量化地揭示了历史瓶颈、基因流中断和当代近交如何逐步侵蚀 U. rockii 的遗传基础，使其濒危状态有了清晰的遗传学解释。
第七，其他有价值的内容 本研究的数据共享非常规范，组装的基因组、重测序数据和转录组数据均已存放于国家基因组科学数据中心，保证了研究的可重复性和数据资源的可利用性。此外，文中提到的通过种子进行快速组织培养和繁殖 *U. rock