低覆盖度全基因组重测序揭示 Quercus robur 复合体内 Q. pedunculiflora 的弱分化及其对中欧-东欧气候适应性森林管理的启示

《Trees, Forests and People》：Genome-wide evidence for differentiation of Quercus pedunculiflora within the Quercus robur (Fagaceae) species complex—implications for forest management in Central-Eastern Europe

气候变化预计将加剧中欧大部分地区的干旱胁迫，这引出了南部橡树群体是否可能成为未来气候适应性森林管理中有价值遗传资源的问题。在这项探索性研究中，研究人员利用低覆盖度全基因组重测序（low-coverage whole-genome resequencing）分析了采自保加利亚、罗马尼亚和匈牙利的 Quercus robur（pedunculate oak）、Q. pedunculiflora（greyish oak）以及形态过渡个体之间的基因组关系。全基因组分析揭示了 Q. robur 复合体内微弱但一致的核基因组分化，这与将 Q. pedunculiflora 识别为一个分化谱系相一致。在本研究采样设计的限制内，该模式与亚种等级相容，而不支持将其处理为独立物种或 Q. robur 的简单生态变型。在采样的匈牙利材料中，未发现 Q. pedunculiflora 存在的基因组证据，包括此前基于形态学被归入该分类群的个体。系统发育基因组学分析显示染色体水平信号具有异质性，而形态过渡个体未形成单独的基因组谱系。对21个先前鉴定的干旱响应候选基因（drought-responsive candidate genes）的分析表明，分类群和地理组之间的归一化映射剂量（normalised mapped dosage）存在大量重叠，且无统计学支持的差异。总之，这些结果未提供证据表明所检测的干旱响应候选基因中的系统性剂量差异是观察到的基因组分化的基础。从应用角度看，巴尔干橡树群体是后续在种源试验和辅助基因流试验中进一步评估的候选对象，而不是已获遗传证明、可用于立即业务化部署的干旱适应性变异来源。

研究背景与问题：橡树属约435种，是北半球生态和经济上最重要的森林树种群体之一。其物种界定困难，原因包括表型可塑性、频繁杂交和持续基因流，导致近缘分类群间遗传分化弱，常被视为典型的橡树syngameon，即由有限但反复基因交换连接的谱系集合。Q. robur 是欧洲低地关键种，具有重要生态经济价值，但在气候变化下，中欧-东欧低地及干旱暴露林分长期存续不确定。东南欧，尤其巴尔干半岛群体因适应暖干条件而受到关注，但引种前需遗传评估。Q. pedunculiflora 与 Q. robur 形态相似，分类地位长期争议，曾被处理为生态型、亚种或独立物种；遗传证据显示两者分化很低、缺乏诊断等位基因并广泛共享叶绿体单倍型。其是否存在于匈牙利尚不确定。因此，研究人员使用低覆盖度全基因组重测序（low-coverage whole-genome resequencing）表征匈牙利和北巴尔干 Q. robur、Q. pedunculiflora 及形态过渡个体，评估复合体内分化、匈牙利样本的基因组身份，并分析干旱响应候选基因变异，以支持森林管理建议。

关键技术方法：样本来自保加利亚、罗马尼亚和匈牙利，共15株，包括14株先前形态计量研究个体及匈牙利 Bikács 1株，采样点包括孤立树、路旁、城市公园和局部小林分，非代表性群体调查。采集休眠芽并提取DNA后，进行 Illumina 双端150 bp低覆盖度全基因组重测序，平均约10×覆盖度。质体基因组用 GetOrganelle 组装，经 MAFFT 比对、BMGE 曲率和 IQ-TREE 最大似然法构建质体系统发育树。核基因组 reads 比对至 Q. robur 参考基因组，FreeBayes 进行全基因组单核苷酸多态性（single nucleotide polymorphism, SNP） calling，bcftools 过滤，PLINK 进行连锁不平衡修剪。随后用 IQ-TREE 构建核系统发育树并计算基因一致因子（gene concordance factor, gCF），进行主成分分析（principal component analysis, PCA）、ADMIXTURE 聚类及 50-kb 滑窗固定指数（F_ST）分析。对21个干旱响应候选基因构建 mini-reference，估算归一化映射剂量（normalised mapped dosage），并用 Wilcoxon 秩和检验和线性混合效应模型比较。

研究结果：3.1 质体系统发育。基于全质体比对的最大似然树显示，Q. robur、Q. pedunculiflora 和形态过渡个体未形成清晰分离的互单系分支，而是分类学混合分布。匈牙利 Q. robur 形成完全支持亚支，但 Bikács 个体也聚于其中；该个体先前被归为 Q. pedunculiflora，却与匈牙利 Q. robur 聚类，而非巴尔干 Q. pedunculiflora。巴尔干 Q. pedunculiflora 和过渡个体分散于多个位置。结果表明质体变异不支持 Q. robur 与 Q. pedunculiflora 的清晰系统发育分离，提示共享母系。

3.2 核系统发育与全基因组群体结构。基于全基因组 SNP 的核系统发育显示，染色体水平树总体分支支持高但基因一致度低，拓扑 discordance 明显。忽略两个异常个体后，Q. robur 和 Q. pedunculiflora 形成互单系分支；Bikács 个体与匈牙利 Q. robur 聚类，Muselievo 两个 Q. robur 为离群。Q. robur 内部存在匈牙利与巴尔干两组，Q. pedunculiflora 较内聚。形态过渡个体未形成独立簇，而持续与 Q. pedunculiflora 聚类。PCA 显示三组，即匈牙利 Q. robur、巴尔干 Q. robur 和 Q. pedunculiflora，另有两个 Muselievo 离群；PC1 和 PC2 分别解释8.8%和8.2%方差。ADMIXTURE 分析中，K=2 最佳，两个祖先成分大致对应先验分类；K=3 进一步区分匈牙利与巴尔干 Q. robur。Bikács 个体在 K=2–5 均与匈牙利 Q. robur 同一成分，无 Q. pedunculiflora 贡献。全基因组分析未发现匈牙利个体形成对应 Q. pedunculiflora 的独特簇。窗口化 F_ST 显示大部分基因组分化低，通常为0.01–0.06，但散布局部高分化峰，部分超过0.15–0.20，稀有个体接近0.23，染色体4、9及10–11峰值较明显。

3.3 候选基因剂量模式。对21个干旱响应候选基因的归一化映射剂量分析显示，Q. robur、Q. pedunculiflora 和过渡样本间广泛重叠，基因间绝对剂量差异大于组间差异。逐基因 Wilcoxon 检验经多重校正后均无显著差异，校正后 P≥0.72。全局混合效应模型中，分类群对归一化映射剂量无显著总体效应，χ²=1.24，d.f.=1，P=0.266；Q. pedunculiflora 系数为正，β=0.119±0.104 SE，对应约1.13倍增加，但置信区间重叠无效应。基因内标准化后分析仍无显著分类群效应，χ²=1.05，d.f.=1，P=0.305。候选基因 PCA 结构有限，PC1 和 PC2 分别解释9.5%和8.6%，匈牙利与巴尔干样本沿 PC1 分离，一个 Q. pedunculiflora 个体为离群。地理分组匈牙利与巴尔干比较亦无显著差异，χ²=1.66，P=0.197；标准化后 χ²=2.07，P=0.150。21个候选位点均未与最高 F_ST 前1%窗口重叠。总体无统计学支持的干旱相关候选基因剂量差异，仅见弱趋势。

讨论部分总结：核基因组分析支持 Q. pedunculiflora 为 Q. robur 复合体内弱分化谱系，质体与核基因组不一致可能反映历史渐渗和叶绿体捕获。研究人员倾向于将其视为 Q. robur subsp. pedunculiflora，而非独立物种或简单生态变型。Muselievo 样本可能为人为引入，故排除；Bikács 个体不支持匈牙利存在基因组上独特的 Q. pedunculiflora。巴尔干同域 Q. robur 个体显示部分混合祖先成分和形态中间性，提示局部渐渗；形态过渡个体反而与 Q. pedunculiflora 聚类，说明毛状体形态中间性不一定对应全基因组中间性。候选基因分析为探索性，未发现差异，可能与干旱适应多基因性、样本量有限及候选位点未捕获适应性变异有关。橡树 syngameon 表现为低分化背景上散布高分化峰，核基因组保留谱系一致性而质体共享。管理上，结果不支持立即大规模将北巴尔干繁殖材料转移至匈牙利，但无强基因组屏障，支持通过种源试验和小规模辅助基因流进一步评估。

结论翻译：本研究提供了中欧和东南欧 Q. robur 复合体分类地位与基因组模式的全基因组视角。考虑到橡树 syngameon，结果支持将 Q. pedunculiflora 解释为一个弱分化谱系，与亚种等级一致，即 Q. robur subsp. pedunculiflora，而非独立物种或 Q. robur 的简单生态变型。显然，在 Q. robur 复合体内使用高鉴别力 SNPs 进行全分布区系统地理学研究，可对研究人员的分类处理提供坚实检验。尽管候选基因分析未检测到干旱相关位点的显著差异，全基因组数据揭示了区域间结构化变异，表明适应性分化可能细微且分布于全基因组。鉴于 Q. robur 作为关键种的生态重要性，将 Q. pedunculiflora 认定为 Q. robur subsp. pedunculiflora，突显其与气候适应性森林管理的潜在相关性，尤其对识别适合匈牙利未来气候条件的树种具有重要意义。特别是来自保加利亚北部和罗马尼亚南部的树木，可能是在变化环境条件下森林更新或林分改造中值得考虑的宝贵遗传变异来源。然而，样本数量有限意味着当前结果应被解释为首次基因组评估，而非整个区域基因库的确定性表征。尽管如此，它们为未来研究提供了首个稳健基因组框架，并为在中欧-东欧气候适应性森林管理中整合进化与应用视角提供了有价值基础。