沿后冰期迁移路线解析中欧—东南欧夏栎（Quercus robur）核质种群遗传结构与遗传保护单元

《Diversity and Distributions》：Along Ancient Migration Routes: Population Genetic Structure of Pedunculate Oak (Quercus robur) Across Central and Southeastern Europe

　　

引言与目的：理解树种种内遗传变异的空间分布，对森林遗传资源保护和经营至关重要。夏栎（Quercus robur）是欧洲具有重要经济和生态价值的树种，其南部和东南部种源被认为是暖干气候局部适应的假定储备库，但在全分布区种群遗传研究中仍代表性不足。研究人员分析了中欧和东南欧Q. robur的核与叶绿体种群遗传结构，包括亚种robur与pedunculiflora接触带，并评估后冰期再殖民、地理屏障、精细地形和人类影响对当代多样性的贡献。
地点：从阿尔卑斯山以北的西欧中部至黑海的东西向样带。
方法：研究人员对来自56个种群的1500个个体，使用16个核微卫星（nuclear microsatellites, nSSRs）和8个叶绿体微卫星（chloroplast microsatellites, cpSSRs）进行基因分型。采用贝叶斯和空间多变量分析评估核遗传结构与多样性，量化分化（分子方差分析（analysis of molecular variance, AMOVA）、分化指数（F_ST）），并推断叶绿体谱系/单倍型分化及系统地理结构。
结果：nSSRs揭示种群间分化低但显著，并呈现三个集群模式：西欧中部、迪纳拉阿尔卑斯—潘诺尼亚盆地和东南巴尔干。遗传多样性自西向东南欧增加。叶绿体单倍型变异显示与冰期避难所和迁移路线一致的空间结构。两个亚种分化弱，在接触带共享叶绿体谱系，表明存在实质基因流。
主要结论：后冰期再殖民和景观特征塑造了Q. robur遗传变异。东南部种群保留深层谱系结构和较高区域等位基因丰富度，并含许多私有单倍型/等位基因，符合后冰期扩张的后缘位置。北阿尔卑斯接触带多样性升高表明再殖民路线间混合，与阿尔卑斯—喀尔巴阡走廊附近隐蔽晚冰期微避难所相容。整合核与叶绿体模式，建议保护东南部基因池、维持接触带连通性，并使遗传保护单元（Genetic Conservation Units, GCUs）与所识别集群对齐。

研究背景与问题：理解种内遗传变异空间分布是森林遗传资源保护与经营的基础，尤其在气候变化背景下。夏栎（Quercus robur）是欧洲阔叶林基础树种，具重要经济与生态价值，分布广且生态可塑性强；其南部和东南部种源可能适应暖干气候，可为气候适应性造林提供种源，但全分布区核标记研究仍少，多数为区域尺度或聚焦白栎杂交，东南后缘种群代表性不足。其遗传结构受末次冰盛期（Last Glacial Maximum, LGM）避难所、后冰期再殖民、阿尔卑斯—迪纳拉—喀尔巴阡地形、人类土地利用和种子调运影响。该研究拟解析西欧中部—东南欧梯度上的核与叶绿体结构、多样性和基因流，评估后缘/前缘框架下避难所与接触带贡献，并检验亚种robur与pedunculiflora分化，为遗传保护单元（Genetic Conservation Units, GCUs）提供操作框架。

研究人员对56个种群1500个个体跨样带取样，开展nSSR与cpSSR分型及空间遗传分析。结论为：核分化低但显著并呈三集群，东南部多样性高且私有等位多，叶绿体空间结构强，亚种分化弱；研究强调保护东南基因池、维持接触带连通性并将GCUs与集群对齐，对气候适应型森林遗传资源管理具有重要意义。

关键技术方法（不超过250字）：样本队列来自中欧至东南欧56个Q. robur种群1500个个体（平均27株/种群），并纳入ACORN项目Q. petraea和Q. pubescens各6个参考种群核验物种；cpSSR对每第二株共729个个体分型。关键方法包括16个nSSR与8个cpSSR分型；STRUCTURE/CLUMPAK贝叶斯聚类；COLONY推断同胞关系；主坐标分析（PCoA）、空间主成分分析（sPCA）和Alleles in Space（AIS）分析空间结构；AMOVA、F_ST和Mantel检验隔离距离（IBD）；HP-RARE估计等位基因丰富度与私有等位基因丰富度；PermutCpSSR计算h_S、h_T、G_ST、R_ST及系统地理检验；partial RDA与dbMEM评估亚种效应和空间协变量。

研究结果：
3.1 物种鉴定：多物种nSSR贝叶斯聚类在K=2区分Q. robur与Q. petraea/Q. pubescens；以Q≥0.85保留1484个“纯”Q. robur个体、56个种群，剔除16个个体。
3.2 种群内家族结构：COLONY显示多数种群家族结构弱，但380、381（迪纳拉）及300、303、305（奥地利）全同胞对多；去除近缘个体后51个种群1321个个体重跑，K≤3与全数据高度一致，故后续用全数据。
3.3 种内种群遗传结构与分化：56个种群核分化低但显著（AMOVA种群间约3.32%，F_ST=0.035）；STRUCTURE K=2—3与sPCA/PCoA一致，分西欧中部、迪纳拉—潘诺尼亚、东南巴尔干三集群，界限渐变；sPCA全局显著、局部不显著；AIS表面西部较均质，东南和部分迪纳拉异质性高。
3.4 区域种群结构：西欧中部内部呈西—东分裂（法国/瑞士/德国与奥地利），AMOVA种群间约1.94%，IBD显著；迪纳拉—潘诺尼亚无额外区域亚结构，380/381/383分化最强，无显著IBD；东南巴尔干分化最强，K=2分希腊/西黑海/保加利亚—罗马尼亚边界与西北保加利亚/东喀尔巴阡，395最分化，IBD显著；主要遗传集群大体对应pedunculiflora，但罗马尼亚397例外；亚种间AMOVA小但显著（1.87%，p=0.0008，总体3.7%），partial RDA控制空间后亚种效应不显著，空间解释大部分。
3.5 遗传多样性：全数据Na均9.16、Ne 4.78、Ho 0.692、He 0.697、F_IS低；AR₁₇约7.98、PAR₁₇约0.06；He、AR₁₇东南显著高于西欧和迪纳拉—潘诺尼亚，PAR₁₇东南最高（边缘显著），奥地利高多样且有私有等位；最高F_IS见于367及瑞士342/343。
3.6 叶绿体单倍型地理格局：729个个体检出17个cpSSR单倍型、4个谱系A/B/C/E；A最常见（65%），单倍型6约32%（对应RFLP 7）沿北阿尔卑斯、迪纳拉和黑海，单倍型4约13%；B稀有，C在北阿尔卑斯，E在东南和迪纳拉/奥地利；分布支持巴尔干避难所、伊比利亚残留及多路径再殖民。
3.7 cpSSR遗传多样性：总体h_T=0.854、h_S=0.247、G_ST=0.711、R_ST=0.697，全局系统地理结构不显著；西欧R_ST=0.756显著高于G_ST（p=0.023），提示区域系统地理结构；迪纳拉—潘诺尼亚h_S最高（0.311）；东南G_ST最高（0.725）。

讨论总结：讨论将三集群归因于阿尔卑斯、迪纳拉、喀尔巴阡等屏障与后冰期动态；东南后缘长期近避难所，保留深层谱系、高区域多样性和私有等位，局部低境内多样与隔离一致；迪纳拉—潘诺尼亚呈意大利与巴尔干来源混合，cpDNA保留双路径而核标记被花粉流均质化，崎岖地形和人类利用导致局部分化及380/381低多样强家族结构；西欧核分化低，北阿尔卑斯为核质接触带，奥地利高多样可能反映迁移路线交汇及隐蔽晚冰期微避难所贡献，但作者称尚需加密采样验证；cpSSR显示母系、种子限制导致强空间结构，西欧有显著系统地理信号，迪纳拉混合，东南马赛克；亚种robur与pedunculiflora核分化弱，地理和局部基因流解释大部分结构，cpDNA不能区分亚种，但pedunculiflora含私有核等位。

研究结论（译）：尽管中欧和东南欧低地夏栎林长期人类利用（如皆伐、临时农地化、再造林），研究观察到与已知cpSSR谱系广泛一致及nSSR空间结构相符，符合后缘持续和前缘扩张动态，未发现远距离种子移动或区域基因池广泛改变证据，符合中世纪和近代早期运输储存受限的本地播种实践；私有等位在若干后缘和迪纳拉—潘诺尼亚种群持续，提示长期局部连续性。东南巴尔干边缘避难和残遗种群具独特基因池和高私有等位，保护价值高，但面临水文改变（河流整治、排水、洪泛区断连）和暖干化累积压力，增加夏栎衰退风险。覆盖稀疏处应优先纳入国家/欧盟框架，特别是EUFORGEN及其EUFGIS登记原地GCUs；样带内最南种群（希腊360、保加利亚395）和迪纳拉—潘诺尼亚点（塞尔维亚311、312）因私有等位和区域独特性为GCU强候选，而EUFGIS目前仅列两个巴尔干GCU（西黑海沿岸和塞尔维亚），希腊、保加利亚、黑山等未代表。西欧接触带形成高多样、私有等位的cytonuclear mosaic，建议巩固西南GCU（含瑞士）并承认奥地利走廊，以维持谱系多样性和连通性；迪纳拉部分种群（380、381）低多样强家族结构，需原地保护和迁地备份，虽已用于生长试验和繁殖材料，但低多样限制广泛部署；操作上种子区和GCU可与三集群对齐，接触带作为单元内多样性库，并保护地形限制基因流的走廊（如北阿尔卑斯）。