人类影响下山茱萸（Cornus officinalis）遗传分化与气候适应：药用经济林可持续管理与保护

《Trees, Forests and People》：Understanding Genetic Divergence and Climate Adaptation in Cornus officinalis for the Sustainable Management and Conservation of Medicinal Economic Forests under Human Impact

　　

气候变化正严重威胁药用植物的多样性和药效，对药用植物经济林的可持续发展产生负面影响，因而需要利用景观基因组学（landscape genomics）等工具分析气候适应性、管理种质资源并促进气候适应型育种（climate-resilient breeding）。山茱萸（Cornus officinalis）是一种重要的木本传统中药植物，主要分布于中国秦岭-伏牛山和天目山，是重要的经济林资源。本研究利用中国核心分布区原生种群的双酶切限制性位点相关DNA测序（ddRAD-seq）数据，研究C. officinalis的种群结构、局部适应及对气候变化的基因组脆弱性（genomic vulnerability）。遗传结构分析揭示三个与主要山脉相对应的不同谱系：秦岭、伏牛和天目山。广泛的基因流（gene flow），尤其是涉及伏牛谱系的基因流，提示显著的人为介导引种或自然基因流模糊了遗传边界。检测到强的距离隔离（IBD）和环境隔离（IBE）。共有1113个受选择的候选单核苷酸多态性（SNP）被注释到参与胁迫响应、次生代谢和发育的基因。鉴定出谱系特异性环境驱动因子：伏牛谱系的温度季节性（Bio4）和等温性（Bio3），以及天目谱系的年降水量（Bio12）。未来气候情景下的预测显示显著的遗传偏移（genetic offset），尤其是当使用基因型-环境关联SNP时，秦岭谱系种群表现出更高的基因组脆弱性。这些发现强调地理屏障、局部适应和人为活动在塑造C. officinalis当前遗传结构中的作用，并突出需要针对C. officinalis的谱系特异性保护策略和气候适应型育种指南。

研究背景：气候变化由人类活动驱动，改变温度、降水、极端事件和生态系统，威胁药用植物多样性、药效、传统知识和保护。快速环境变化可能超过自然进化速度，导致适合度下降和局部灭绝，因此药用植物种质资源管理和保护需考虑遗传多样性和进化潜力，尤其长寿木本药用经济林。识别高适应种群、预测未来气候下栖息地适宜性、推进气候适应型育种（climate-resilient breeding，筛选有益等位基因并整合适应性状的育种策略）很重要。景观基因组学（landscape genomics，整合全基因组与空间环境变量以揭示适应性遗传基础并预测基因组偏移的交叉学科）可量化遗传偏移（genetic offset，当前遗传组成与未来气候需求之间的错配程度），识别脆弱种群并指导辅助迁移和标记辅助选择。山茱萸（Cornus officinalis）是重要木本传统中药，果实干制品“Yurou”/“Zaopi”富含环烯醚萜苷、三萜和多糖，有抗炎、抗氧化、抗凋亡、抗糖尿病、神经保护和心血管保护等活性；核心分布于秦岭-伏牛和天目山，河南、陕西、浙江为重要经济林产区，长期栽培和贸易导致无完全野生种群，种质交换可能影响遗传边界。以往RAPD、ISSR、SSR等中性标记揭示高多样性和部分地理结构，但无法连接遗传变异与环境驱动，限制适应理解和保护策略设计。因此研究人员用景观基因组学识别与关键环境梯度相关的适应性位点，揭示局部适应机制并服务气候适应型品种选育。研究人员以核心分布区原生种群为对象，基于ddRAD-seq获取全基因组SNP，检验异质环境促进秦岭、伏牛、天目山种群遗传分化与局部适应，以及长期贸易栽培导致种质交换、自然扩散和人工引种贡献谱系间基因流两个假设。研究得出三谱系结构、广泛基因流、强IBD和IBE、1113个候选受选择SNP、谱系特异性环境驱动和未来遗传偏移，尤其GEA SNP下秦岭谱系脆弱性高。其意义在于为谱系特异性保护、种质资源管理和气候适应型育种提供基因组基础。

关键技术方法：样本队列来自陕西（5）、河南（12）、安徽（1）和浙江（4）的22个山茱萸种群，每种群8株、采样间距>15 m，并取陕西5株C. chinensis为外类群。关键技术包括ddRAD-seq（双酶切限制性位点相关DNA测序，一种简化基因组测序方法）和STACKS流程获取SNP（单核苷酸多态性，DNA序列中单个碱基的变异）；ADMIXTURE、PCA（主成分分析）和ML（最大似然）树分析种群结构；TreeMix检测基因流；BayeScan、PCAdapt及LFMM/RDA识别受选择与GEA（基因型-环境关联，检测与环境变量显著相关的遗传位点）位点；Mantel检验和偏RDA解析地理与环境贡献；GF（梯度森林，非参数机器学习回归树方法）和MaxEnt（最大熵模型，物种分布模型）预测SSP126/SSP585下遗传偏移。

结果：3.1 种群结构与遗传多样性：对176株C. officinalis的9348个高质量SNP进行ADMIXTURE、PCA和ML树分析，K=3三个谱系对应秦岭、伏牛、天目；安徽JC群体为伏牛与天目混合，后续景观分析剔除。遗传多样性伏牛最高（π=0.00041，观测杂合度H_O=0.00026），天目次之（π=0.00039，H_O=0.00019），秦岭最低（π=0.00032，H_O=0.00018）；所有种群近交系数F_IS>0，天目近交最高；谱系间遗传分化指数F_ST为0.075-0.150，秦岭-天目最高，秦岭-伏牛最低；中性SNP结果一致。3.2 基因流分析：TreeMix选m=13，拓扑与种群关系一致；秦岭谱系遗传漂变高；检测到频繁迁移事件，秦岭与伏牛间基因流较高，并见祖先节点到现存种群迁移；说明广泛基因流/混合。3.3 受选择候选SNP：BayeScan和PCAdapt分别检测415和697个离群位点，交集81，任一1031；LFMM和RDA分别检测22和177个环境关联位点，交集17，任一182；合计1113个候选SNP。注释涉及胁迫响应（热激蛋白、水通道蛋白、热胁迫转录因子、受体激酶、丝氨酸/苏氨酸激酶）、次生代谢（花青素、纤维素、苯丙烷）、发育（细胞壁、生长素、开花时间、表观修饰）及转座元件和防御，提示多基因适应。3.4 遗传变异的地理与环境关联：Mantel检验显示强IBD（距离隔离，遗传距离随地理距离增加而增加；r=0.7742，p=0.001）和IBE（环境隔离，遗传距离随环境差异增加而增加；r=0.7735，p=0.001），地理与环境距离相关r=0.9059（p=0.001）。RDA（冗余分析）显示环境解释22.7%-41.6%、地理22.6%-32.1%、种群结构25.1%-40.2%；GEA SNP中环境更高；纯种群结构显著（3.3%-5.0%），纯环境仅GEA显著（4.9%），纯地理不显著。谱系特异性：伏牛受等温性（Bio3）、温度季节性（Bio4）和纬度影响，天目受年降水量（Bio12）影响；局部适应主要作用于环境关联位点。3.5 梯度森林分析与未来气候脆弱性：GF显示等温性（Bio3）最重要，其次降水季节性（Bio15）、年降水量（Bio12）、年均温（Bio1）、温度季节性（Bio4）；阈值约Bio3 26-27、Bio4 820、Bio12 1200 mm、Bio15 57-60 mm。MaxEnt预测精度高（平均测试AUC=0.996）。未来高排放情景（SSP585）下遗传偏移大于低排放情景（SSP126）；GEA SNP偏移高于全SNP；全SNP下伏牛最高、天目次之、秦岭最低，GEA SNP下秦岭最高、伏牛次之、天目最低，提示谱系脆弱性不同且GEA位点更敏感。

讨论部分总结：三谱系与地理对应，山地地形和长江中下游平原等屏障塑造结构；广泛混合和基因流可能来自人为引种/种苗交换或自然二次接触。秦岭谱系位于ML树基部、多样性低且F_IS>0，可能受人工选择、近交和片段化影响；伏牛多样性最高可能与混合效应有关；天目中等多样性反映基因流、局部适应和漂变。IBD和IBE显著但环境、地理和结构共线性强；纯结构显著而纯环境仅GEA显著，说明背景变异主要受种群结构影响，局部适应集中于环境关联位点。伏牛温度因子与其气候过渡带春季温度波动、晚霜和夏季高温相关；天目降水与其湿润亚热带气候一致；同质园实验支持性状和药用化合物与温度降水相关。候选基因注释支持多基因适应。GF阈值可能是生态临界点，Bio3阈值可能关联温度波动耐受和冷耐受/物候可塑性，Bio12阈值可能关联水分可用性和耐旱。全SNP与GEA SNP的遗传偏移反映整体背景与适应性位点两类证据，应互补评估。管理上建议谱系特异性就地保护：秦岭保护核心区并辅助授粉/混合采种缓解近交；伏牛设气候适应优先区保护温度敏感种群；天目保护栖息地和水文环境。迁地保护建立种质资源圃并按谱系、遗传背景和性状编目构建核心种质。谨慎实施辅助基因流（assisted gene flow，人为将适应未来气候的等位基因引入目标种群），以伏牛为桥梁亲本，通过控制杂交筛选多气候适应等位基因，显著分化谱系采用逐步基因流（秦岭→伏牛→天目）并监测物候同步和适合度。遗传偏移预测不足独立指导管理，需交互移植实验（reciprocal transplant experiments）和同质园实验（common garden experiments）及功能基因组验证。

研究结论翻译：本研究基于全基因组数据阐明了C. officinalis的种群结构和适应进化机制。该物种当前遗传架构可能源于历史进化过程（尤其是地理屏障驱动的地理隔离（vicariance）和适应进化）塑造的三大谱系（秦岭、伏牛和天目山）分化，以及随后可能由人为引种和自然过程驱动的基因流。这种基因流丰富了遗传多样性（在伏牛谱系尤为突出），同时造成遗传混合（如JC群体）。地理距离（IBD）和环境异质性（IBE）共同驱动遗传变异，自然选择通过与环境胁迫耐受和代谢功能相关的候选SNP调控局部适应。存在显著的谱系特异性环境敏感性：伏牛谱系主要响应温度变异（Bio3、Bio4），天目谱系受降水（Bio12）调控。梯度森林建模预测在高排放情景（SSP585）和使用GEA SNP时存在显著遗传偏移，且GEA SNP偏移高于背景基因组，秦岭谱系种群风险最高，提示该谱系存在严重适应不良风险。为确保C. officinalis在气候变化和遗传侵蚀下的存续和药用价值，研究人员提出综合保护策略：实施谱系特异性就地保护、在气候稳定地区建立种质库，并通过辅助基因流增强进化潜力。