韩国孑遗侧柏的基因组结构揭示隔离谱系和石灰岩走廊多样性热点

《Conservation Genetics》：Genomic structure of relict Platycladus orientalis in Korea reveals isolated lineages and a limestone-corridor diversity hotspot

　　

孑遗和边缘树木种群可能结合较低的种群内变异和较强的种群间分化，使保护决策复杂化。研究人员结合了一个从头组装的侧柏（Platycladus orientalis）基因组草图（约8.73 Gbp；N50=248,484 bp）与基因分型测序（GBS）数据，包含来自韩国7个林分的95个个体的38,820个单核苷酸多态性（SNP）。丹阳-堤川石灰岩走廊的核苷酸多样性（π=0.245–0.274）高于孤立的安东、英阳和大邱种群（π=0.156–0.216）。孤立种群之间相互高度分化（成对F_ST=0.396–0.431），而走廊内的分化较低但异质（F_ST=0.117–0.223）。分子方差层次分析将29.26%的变异归因于种群间差异（Φ_ST=0.293）。未检测到重复或近乎相同的多位点基因型。常见变异中稀疏私有等位基因丰富度在走廊群体水平高于孤立群体（每座位0.1214 vs. 0.0415），尽管稀有变异被低估。F_IS在严格深度和缺失数据过滤下强烈下降，表明基线杂合子缺陷存在实质性的技术成分。当代有效种群大小（N_e）估计高度异质且不精确，范围从丹阳石门的2.0到大邱的202.8，而丹阳哈溪里无界。探索性Mantel检验提示名义上的地理和气候关联，但在多重检验校正后均不显著。三个焦点林分的成苗F_ST值≤0.035。这些结果支持将孤立林分视为候选保护单元，将石灰岩走廊视为管理相关的多样性储存库，而基因流或种子转移应通过风险评估框架进行评估。

### 论文解读：韩国孑遗侧柏的基因组结构揭示隔离谱系与石灰岩走廊多样性热点

**研究背景与问题**

孑遗和边缘树木种群通常面临遗传漂变和基因流受限的双重压力，导致种群内多样性降低而种群间分化增强，这给保护决策带来复杂性。侧柏（Platycladus orientalis）是一种耐逆性强的东亚针叶树，在韩国以少量孤立林分形式分布于石灰岩悬崖、岩石斜坡和河流陡岸。此前关于侧柏的系统地理学、景观基因组学及SSR研究主要聚焦中国材料，而韩国自然林分的遗传状态——是遗传侵蚀的孤立片段、相互连接的种群，还是独特的边缘谱系——尚未澄清。本研究旨在利用物种特异性参考基因组和全基因组GBS标记，回答五个核心问题：遗传多样性如何在丹阳-堤川石灰岩走廊与孤立悬崖边缘林分间分布；多位点基因型重复和SNP过滤敏感性是否影响克隆性和杂合子缺陷推断；种群分化程度及SNP聚类是否支持管理相关分组；基于连锁不平衡（LD）的当代有效种群大小（N_e）异质性如何；成树与幼苗群体是否存在遗传连续性。

**研究内容与结论**

研究人员结合de novo侧柏基因组草图（约8.73 Gbp，N50=248,484 bp）与来自韩国7个自然林分95个个体的38,820个SNP的GBS数据，进行了全面的群体遗传学分析。结果表明，丹阳-堤川石灰岩走廊种群的核苷酸多样性（π=0.245–0.274）显著高于孤立种群（安东、英阳、大邱，π=0.156–0.216）。孤立种群间高度分化（成对F_ST=0.396–0.431），而走廊内分化较低但异质（F_ST=0.117–0.223），其中丹阳永川里与丹阳哈溪里分化最大（F_ST=0.223）。分子方差分析（Φ_ST=0.293）表明29.26%的变异来自种群间。未检测到重复基因型。稀疏私有等位基因丰富度在走廊群体水平（0.1214/座位）是孤立群体（0.0415/座位）的约2.93倍，但稀有变异被低估。F_IS在严格深度与缺失过滤后显著下降，表明基线杂合子缺陷存在技术成分。N_e估计高度异质且不精确，丹阳石门低至2.0，大邱达202.8，部分种群无穷上限。成苗F_ST≤0.035，显示阶段连续性。这些结果支持将孤立林分作为候选保护单元，并将石灰岩走廊视为管理相关的多样性储存库，同时建议以风险评估框架评估基因流或种子转移。

**重要意义**

本研究发表于《Conservation Genetics》，首次为韩国7个自然侧柏林分提供了基因组水平基线，揭示了隔离谱系与石灰岩走廊多样性热点之间的鲜明对比，为侧柏保护提供了基于风险评估的管理框架，区分了即时保护优先事项与未来遗传拯救或辅助基因流的决策。

**关键技术方法**

研究人员采用PacBio Revio平台进行de novo基因组组装（HiFi reads，hifiasm软件），获得8.73 Gbp的草图。GBS文库采用双酶切（MspI和PstI）方案，在Illumina NovaSeq X平台测序（151 bp双端），产生约9.7亿条reads。通过BWA-MEM2比对到参考基因组，DeepVariant进行SNP calling，经质量过滤（深度≥3，缺失率≤30%，MAF>0.05）获得38,820个SNP。主要分析包括：fastSTRUCTURE（K=1–10）进行贝叶斯聚类；主成分分析（PCA）和判别主成分分析（DAPC）进行多变量排序；Weir-Cockerham F_ST和分子方差分析（AMOVA）评估分化；基于LD的N_e估计（NeEstimator v2.1，P_crit=0.02/0.05，随机交配模型）；稀疏私有等位基因丰富度（超几何稀有化法）；成苗比较（F_ST和G检验）；以及探索性Mantel检验（地理距离和气候PCA距离）。样本来源为韩国7个自然林分：丹阳永川里（N=22）、丹阳哈溪里（N=8）、丹阳石门（N=10）、堤川苏山里（N=12）、安东九里（N=12）、英阳甘川里（N=15）和大邱都洞里（N=16），其中部分样本为幼苗。

**研究结果**

- **参考基因组与SNP数据集质量**：通过PacBio Revio测序和hifiasm组装，获得57,500条序列的8.73 Gbp草图，N50为248,484 bp。经SNP calling和过滤后，最终数据集包含38,820个双等位SNP，覆盖95个个体。平均每个样本约72,533个SNP观测值，整体纯合率64.43%，杂合率12.09%。

- **遗传多样性**：七个种群的平均期望杂合度（He）为0.243，平均观测杂合度（Ho）为0.190，所有种群Ho低于He，导致正基线F_IS（均值0.191）。但严格过滤后F_IS大幅降低，因此不作为直接近交证据。核苷酸多样性（π）在走廊种群中较高（丹阳永川里0.245，丹阳哈溪里0.274，堤川0.267，丹阳石门0.248），孤立种群较低（安东0.156，英阳0.175，大邱0.216）。

- **多位点基因型筛选与F_IS稳健性**：在4,465对个体中，未发现重复或近乎相同的多位点基因型（最低错配率0.06994，远高于阈值0.005）。过滤敏感性分析显示，12种深度与缺失组合（SNP数3,973–31,970）下，F_IS从DP≥3/MISS30%的0.067–0.217降至DP≥10/MISS10%的-0.071–0.117，表明基线F_IS包含显著技术成分。丹阳石门和堤川在严格过滤后仍保留正F_IS，提示可能仍有生物学因素。

- **群体结构**：fastSTRUCTURE支持K=7为最佳聚类。孤立种群（安东、英阳、大邱）显示种群特异性祖先成分，而四个丹阳-堤川种群共享成分，但各自仍有分化。PCA前三个轴解释24.88%变异，孤立种群在排序空间中分离。DAPC（20个PC）平均分配成功率为79.5%，进一步确认了分离模式。种群级UPGMA树也显示孤立与走廊种群间的大区分。

- **遗传分化**：AMOVA将29.26%变异归因于种群间（Φ_ST=0.293，P=0.001）。孤立种群间F_ST为0.396–0.431，孤立-走廊对比为0.261–0.339。走廊内F_ST为0.117–0.223，丹阳永川里与丹阳哈溪里最高（0.223）。

- **地理与气候关联**：探索性Mantel检验显示遗传距离与地理距离（r=0.507，P=0.038）及气候PC1-PC2距离（r=0.602，P=0.047）存在名义正相关，但Bonferroni校正后均不显著（α=0.0167）。仅7个种群且预测变量共变，故不作因果解释。

- **当代有效种群大小**：LD法N_e在P_crit=0.02/0.05下高度异质：丹阳永川里64.5/115.9，丹阳哈溪里无穷，丹阳石门2.0/2.0，堤川4.0/3.2，安东9.3/8.8，英阳19.1/17.8，大邱202.8/202.8。多数置信上限无穷。低估计值表明部分种群面临强烈漂变，但需谨慎解释。

- **成苗遗传连续性**：在丹阳永川里、大邱和英阳三个焦点林分，成树与幼苗间F_ST≤0.035，显示种群内遗传连续性。He和等位基因丰富度相似，但幼苗核苷酸多样性较低（如英阳0.113 vs. 0.225），可能源于抽样方差或家族结构。F_IS在成苗间方向不一致，未证实随机交配。

**讨论与结论**

讨论部分指出，孤立悬崖边缘林分（安东、英阳、大邱）的强分化与低多样性模式与长期小规模栖息地隔离一致，这些林分可能作为生态避难所。丹阳-堤川石灰岩走廊的高多样性和共享祖先成分支持其为多样性热点，但堤川林分虽未受法律保护，却具有第二高核苷酸多样性，值得纳入管理。F_IS对过滤的敏感性表明基线杂合子缺陷主要源于技术因素，而非近交；部分种群（丹阳石门、堤川）残留正F_IS需进一步空间与亲缘分析。成苗连续性不代表遗传安全，低N_e和低多样性仍是关切点。孤立林分建议作为候选保护单元，避免未筛选的转移；走廊应保持多个本地种子来源，而非视为单一泛交区。风险评估框架强调，基因流或种子转移需在评估效益和风险后实施。

**研究结论**：本研究利用38,820个GBS衍生SNP为韩国7个自然侧柏林分建立了基因组基线。核心结果是三个强分化、低多样性的孤立悬崖边缘林分与一个更丰富的丹阳-堤川石灰岩走廊之间的对比。孤立种群间F_ST=0.396–0.431，最好作为候选保护单元；走廊种群整体保持更高核苷酸多样性和组水平私有等位基因丰富度，但内部分化异质（F_ST=0.117–0.223），需维持多个本地种群。额外分析显示，无克隆基因型，但基线F_IS受技术影响；N_e高度异质且不精确；Mantel关联不稳健；成苗连续性存在但未排除观赏树基因渗入。保护应结合隔离谱系保护、走廊栖息地维持、多本地种子源保留、直接种群监测，并谨慎评估未来遗传拯救或辅助基因流。最终决策需更广泛的种植与自然种群采样、稀有私有变异评估、染色体级连锁控制及适应性分化证据。