基于基因组学的挪威两种广泛分布草地物种空间遗传结构对比：支持物种特异性种子来源策略

《Ecology and Evolution》：Contrasting Spatial Genomic Structure in Two Widespread Grassland Species Across Norway Supports Species-Specific Seed Sourcing

　　

生态修复通常依赖本地采集的种子材料，然而许多国家的种子采购策略是基于宽泛的行政区划而非经验性遗传数据。这些通用分区是否能反映具有生物学意义的种群基因组结构，尤其是在共同用于修复的同域物种之间，仍知之甚少，从而限制了基于遗传信息的修复实践。本研究利用全基因组双酶切简化基因组测序技术（double-digest RAD sequencing，ddRAD-seq）分析了挪威两种常用于生态修复的多年生草地植物——红雪轮（Silene dioica，一种专性异交且种子扩散能力有限的物种）和欧蓍草（Achillea millefolium，一个具有高基因流和扩散能力的多倍体复合体）的空间种群基因组结构。对分布于挪威广阔纬度和气候梯度上的17个红雪轮种群（355个体；19200个单核苷酸多态性位点，SNPs）和23个欧蓍草种群（428个体；21294个SNPs）的分析揭示了显著对比的遗传结构模式。红雪轮表现出中等程度的空间结构、显著的隔离距离效应（isolation-by-distance，IBD）以及多个地理上连贯的遗传群组，表明其基因流受到空间限制，并支持划分多个种子转移区。相比之下，欧蓍草表现出微弱的种群分化、广泛的混合以及较少的空间结构化特征，但检测到显著的隔离环境效应（isolation-by-environment，IBE），其中年降水量成为独立于地理距离之外预测基因组分化的最强环境因子，这为划定三个宽泛的基因组区提供了依据。由于物种特异性的生活史特征可能促成了这些差异，因此即使对于一同进行修复的同域物种而言，单一的通用种子来源框架也可能是不充分的。虽然红雪轮的结构化基因组模式支持局部种子采购，但欧蓍草的高度连通性则支持更宽泛的转移区。总体而言，本研究强调了将基因组数据纳入修复规划的重要性，即通过使种子转移距离与物种特异性连通性相匹配，以减少遗传风险并在环境变化下维持进化潜力。

研究背景与科学问题

在全球范围内，草地生态系统正因土地利用变化和气候变化而面临严重退化和生物多样性丧失。生态修复已成为缓解这些影响、促进退化生态系统恢复的关键策略。种子来源的选择是决定修复项目长期成败的核心因素之一，因为种子的遗传组成直接影响建植成功率、种群持久性和进化潜力。理论上，成功的修复需要在维持遗传多样性与适应潜力，以及避免引入不适应当地环境的基因型所导致的远交衰退风险之间取得平衡。然而，在许多国家，包括挪威在内，种子采购指南往往基于行政边界或宽泛的气候区域，而非经验性的遗传数据。这种做法隐含着一个未经检验的假设：地理或环境边界能够对应种群遗传结构，甚至同域分布的物种对环境变量会有相似的反应。但事实上，种群遗传结构受物种特有的生活史性状、繁殖系统、扩散机制和历史动态的综合塑造，不同物种间可能存在巨大差异。因此，采用通用的种子转移区可能低估了遗传结构强的物种内部的遗传分化，同时不必要地限制了基因流高度连通的物种的种子移动。挪威现行的本土野生花卉种子混合物采购策略正是基于10个行政区域，缺乏对野生植物种群内遗传分化和适应性变异的实证评估。为了解决这一知识空白，研究人员提出了三个核心问题：（1）红雪轮（Silene dioica）和欧蓍草（Achillea millefolium）在挪威范围内的基因组结构幅度和空间模式是否存在差异？（2）推断出的遗传结构模式是否与现行行政定义的种子采购区一致，还是更支持物种特异性的种子转移区？（3）如何将全基因组种群基因组数据转化为实际的种子采购区，以识别遗传相容的种子转移单元？

研究设计与技术方法

为回答上述问题，研究人员采用了全基因组双酶切简化基因组测序技术（ddRAD-seq）对两种物种进行了大规模种群基因组学调查。样本覆盖了挪威广阔的纬度和气候梯度，共采集了17个红雪轮种群（366个体）和23个欧蓍草种群（455个体）。通过流式细胞术（FCM）确认了所有红雪轮个体均为二倍体，欧蓍草个体均为六倍体，排除了倍性变异对遗传结构的干扰。经过严格的质量控制、哈代-温伯格平衡检验、缺失率过滤和连锁不平衡修剪后，最终分别获得了包含19200个SNPs（355个体）和21294个SNPs（428个体）的高质量数据集。数据分析方法包括主成分分析（PCA）、基于模型的贝叶斯聚类算法（ADMIXTURE）、分子方差分析（AMOVA）、成对遗传分化指数（F_ST）计算、以及隔离距离效应（IBD）和隔离环境效应（IBE）的Mantel检验和距离冗余分析（dbRDA）。

主要研究结果

3.1 基因组数据与SNP过滤

初始测序数据经严格过滤后，红雪轮和欧蓍草分别获得高质量、低缺失率的SNP数据集，为后续分析提供了可靠的基因组标记基础。

3.2 全基因组种群结构

3.2.1 主成分分析（PCA） PCA揭示了两种物种截然不同的空间遗传结构模式。红雪轮的前两个主成分解释了35.6%的总变异，显示出中等但空间连贯的地理结构：东南部种群沿PC1正向轴聚集，西部和中部种群位于中心，而来自特隆赫姆峡湾的TJON种群在负向轴上形成孤立簇，北部的芬马克（FINN）和特罗姆瑟（TROM）种群也表现出明显分化。相比之下，欧蓍草的遗传变异几乎是连续的，前两个主成分解释了31.4%的变异，各地区种群形成宽泛的重叠簇，无明显区域分离，仅VALU和ALST种群表现出轻微的孤立趋势。

3.2.2 混合分析（ADMIXTURE） 红雪轮的交叉验证误差持续下降未达明确最小值，基于生物学合理性确定K=7为最优分组。在该水平下，西南部、东南部和中部种群形成不同的祖先群，北部进一步分为西北极（TROM）和东北极（FINN）两个区，TJON种群则形成一个近乎纯合的孤立祖先簇。欧蓍草的交叉验证误差在K=5处达到最小值，显示弱结构和高度混合。大多数种群共享广泛的祖先成分，但VALU、ALST和北部TROM种群表现出较低的混合度，由单一祖先成分主导，反映了局部的等位基因频率偏移而非离散的空间结构。

3.3 遗传分化与隔离距离

分子方差分析（AMOVA）显示，红雪轮有7.2%的遗传变异存在于种群间（Φ_ST=0.072，p<0.001），成对F_ST范围为0.023至0.253（均值0.130），表明中等到高的遗传分化。Mantel检验和线性回归均证实了显著的隔离距离效应（IBD）（r=0.278，p=0.027）。相反，欧蓍草仅有2.4%的变异存在于种群间（Φ_ST=0.024，p<0.001），成对F_ST范围仅为0.005至0.105（均值0.041），IBD关系微弱且不显著（Mantel检验：r=0.163，p=0.145）。然而，欧蓍草检测到了显著的隔离环境效应（IBE）（Mantel检验：r=0.276，p=0.024），进一步的偏距离冗余分析（dbRDA）在控制地理距离后，发现年降水量（BIO12）是唯一能显著解释遗传分化的环境变量（F=2.46，p=0.0477）。

讨论与结论

4.1 两种广布草地物种空间基因组结构的鲜明对比

尽管两种物种在昆虫传粉、多年生生活史、广泛地理分布以及常被同域用于修复等方面具有共性，但它们在全基因组尺度上表现出截然不同的空间遗传结构模式。红雪轮呈现空间连贯的中等种群分化和显著的IBD，而欧蓍草则表现为连续的遗传变异、弱种群分化和广泛的混合。流式细胞术排除了倍性变异的影响，表明生活史性状（如繁殖系统和扩散能力）的差异是驱动这种对比的主要因素。这一发现有力地支持了种子采购应基于物种特异性种子转移区，而非通用分区的观点。

4.2 红雪轮的中等遗传结构与IBD

红雪轮中等程度的遗传结构在PCA、ADMIXTURE、AMOVA和IBD测试中均得到一致支持。其平均成对F_ST为0.13，AMOVA揭示7%的变异存在于种群间，这与瑞典北部对该物种的精细尺度研究结果相符。这种分化水平很可能源于其专性异交、有限的种子扩散能力以及挪威复杂崎岖的地形（山脉、深峡湾）对基因流的限制。此外，历史冰期后的多次回迁也对当前的谱系地理结构有所贡献。均匀为负的近交系数（F_IS）进一步印证了其异交系统维持了较高的种群内杂合度。

4.3 欧蓍草的弱遗传结构与连续遗传变异

与红雪轮相反，欧蓍草在整个挪威范围内呈现出几乎连续的遗传变异。PCA显示各区域种群形成一个单一的连续云状分布，ADMIXTURE揭示广泛的祖先共享，平均F_ST仅为0.04，IBD关系极弱。然而，显著的IBE表明环境变异，特别是年降水量，正在积极塑造其基因组景观，这与挪威的区域降水梯度相吻合。欧蓍草高效的扩散机制（风媒种子和泛化昆虫传粉）、历史上的无意识人为传播（如干草运输）以及多倍体带来的缓冲效应，共同促进了其高度连通的基因库。VALU和ALST种群的局部偏离被归因于奠基者效应或局部基因流受限等随机人口过程，而非长期的进化隔离。

4.4 从基因组数据推断的物种特异性种子转移区

基于中性基因组结构模式，研究人员分别为两种物种提出了具体的种子转移区建议。对于红雪轮，建议划分为六个基因组种子转移区：东南内陆区、西南沿海区、中挪威区、TJON本地区、西北极区和东北极区。对于欧蓍草，尽管其整体高度连通，但基于ADMIXTURE揭示的祖先梯度以及显著的IBE，建议划分为三个宽泛的基因组区：中南东部区、西部沿海区和北部区，外加VALU本地区和ALST本地区。研究人员强调，这些建议旨在描述种子来源间的遗传相容性，而非定义适应性起源边界。为了全面评估局部适应，正在进行跨越六个气候对比地点的交互移植实验，这将为完善种子采购策略提供关键的适应性证据。总体而言，本研究凸显了将基因组数据纳入修复规划的重要性，通过使种子转移距离与物种特异性连通性相匹配，可以有效降低遗传风险，并在环境变化下维持物种的进化潜力。