基于全基因组SNP数据解析墨西哥东部云雾林特有树种Ocotea psychotrioides的精细尺度遗传结构与谱系地理历史

《American Journal of Botany》：Fine-scale genetic structure of the cloud forest tree Ocotea psychotrioides revealed by genome-wide SNP data

　　

前提：基因组技术的最新进展使研究人员能够在精细空间尺度上重建进化历史，揭示塑造自然不连续、群岛状生态系统（如云雾林）中生物多样性格局的复杂动态。
方法：研究人员利用基因分型测序（GBS）数据，对墨西哥东部狭域特有树种Ocotea psychotrioides（Lauraceae）的进化历史与谱系地理进行研究；数据包括来自10个群体101个个体的8579个单核苷酸多态性（SNPs）。
结果：研究人员的分析在Ocotea psychotrioides中鉴定出两个主要遗传谱系。除这种地理结构外，研究人员还揭示了群体间遗传多样性与连通性的复杂格局，突显了墨西哥东部云雾林的生态与进化异质性。更新世（Pleistocene）期间的历史气候波动，叠加地形屏障，可能促成了所观察到的遗传分化。值得注意的是，形态变异与遗传格局一致，为未来关于物种形成、局域适应以及气候驱动栖息地变化的进化后果研究提供了稳健框架。
结论：研究人员的发现突显了基因组方法在揭示云雾林分类群精细尺度进化格局方面的效用，揭示了复杂山地系统中的谱系分化、群体连通性和历史生物地理结构。

中美洲云雾林是生物多样性极高的生物群系，分布在较高大气湿度与温和温度区域，但自然破碎、分布有限，呈“群岛状”森林斑块。其北部中美洲区系由热带与温带亲缘物种组合，并受地理位置、海拔、坡度和邻近植物群落影响，斑块间物种周转高、植被结构差异明显。该地区经历长期隔离、过去气候波动和局部造山事件，特万特佩克地峡、东马德雷山脉和跨墨西哥火山带等地理屏障被认为影响云雾林生物区系的进化轨迹；已有干燥模型、湿润模型、更新世避难所模型和屏障模型等，但实证证据有限，分子数据常显示弱或可忽略的人口扩张信号。尤其在同一种孤立斑块内部群体间分化（within-patch model）仍少被检验。Ocotea psychotrioides是墨西哥东部云雾林狭域特有小乔木，叶大小和毛被由南向北增加，且叶形态与当前气候不相关，而可能受历史气候波动影响；其群体遗传分化程度及火山活动形成的物理屏障对遗传连通性的影响尚不清楚，因此研究人员开展全基因组SNP研究，检验地理距离、当前与历史气候及物理屏障等分化驱动因素，并比较遗传与叶形态分化，以理解云雾林精细尺度物种形成与保护意义。
研究人员于2022年3月在墨西哥东部10个群体采集130个个体，最终保留101个基因型、8579个SNP，结合群体遗传结构、空间遗传变异、人口历史、形态距离和古分布模型，识别出Sierra de Chiconquiaco南北两个主要遗传谱系及中北部亚结构，发现过去气候而非当前气候更能解释遗传分化，群体迁移低且不对称，形态分化与遗传分化相关，并重建了过去120,000年高度不稳定的适宜分布。该研究强调基因组数据可揭示狭域特有云雾林树种的隐蔽精细进化历史，并提示评估这些谱系在气候变化下的适应与保护意义。
关键技术方法概括（不超过250字）：研究人员采用基因分型测序（GBS）获取全基因组单核苷酸多态性（SNP），样本来自墨西哥东部10个Ocotea psychotrioides群体（2022年3月采集，最终101个体、8579个SNP）；用iPyRAD流程进行序列组装与变异检测，并以Persea americana参考基因组比对；开展STRUCTURE、主成分分析（PCA）、F_ST与每代迁移者数（N_m）估计及分子方差分析（AMOVA）；用广义相加模型（GAM）和分段回归比较地理距离、当前与历史气候距离对F_ST的贡献；用位点频率谱（SFS）和dadi比较群体分化与迁移模型；用Maxent生态位模型重建古分布；并以叶形态距离与F_ST相关及残差分析评估遗传—形态一致性。
研究结果：
（1）遗传多样性与遗传分化（Genetic diversity and genetic differentiation）：通过GBS最终获得8579个SNP和101个基因型，10个群体各含6—13个基因型；物种水平遗传多样性低，观测杂合度（H_o）=0.1571、群体内遗传多样性（H_s）=0.1730、近交系数（F_IS）=0.0919。Tlapacoyan是唯一F_IS为负值（?0.034）的群体，Xalapa和Xicotepec的F_IS最高。STRUCTURE分析支持K=2，南部Huatusco、Xalapa、Naolinco与中北部Teziutlán、Tlapacoyan、Xicotepec、El Texme、Huayacocotla、Tlanchinol、Xilitla分化；Teziutlán有19.58%南部祖源成分，Huatusco有3.75%中北部成分。中北部子集分析识别7个遗传簇，对应7个采样群体，Tlanchinol混合最高达27.63%。F_ST、N_m和PCA与结构分析一致，PCA第一轴解释7.1%、第二轴解释3.7%，并显示Tlanchinol与Xilitla进一步分化。
（2）空间遗传变异格局（Patterns of spatial genetic variation in O. psychotrioides）：通过提取23,000年以来每1000年的19个生物气候变量并与成对F_ST相关，发现13,000年前（R²=0.2096，AIC=?169.23）和17,000年前（R²=0.2314，AIC=?170.49）气候条件拟合最佳，近期气候与遗传结构无显著关联；分段回归检测到约9743年前的结构断点。GAM显示地理距离（F=11.474，P=0.00159）和历史气候距离（F=5.913，P=0.01960）显著，当前气候距离和气候不稳定性不显著；简化模型解释43.6%偏差、调整R²=0.391、AIC=37.99。AMOVA显示Sierra de Chiconquiaco屏障分组解释的组间变异最高（12.36%），其次为水文盆地与地质特征模型（11.93%）和云雾林区域模型（7.30%）；无预设分组时65.52%变异位于群体内、34.48%位于群体间。
（3）有效群体大小变化（Changes in effective population size）：通过SFS和dadi比较三种人口模型，支持分化伴随不对称迁移模型（split_asym_mig；AIC=536.99，Akaike权重=0.9994）。参数估计显示祖先有效群体大小（N_ref）约417，南部谱系（N₁）约1925、北部谱系（N₂）约4873，分化时间约17,966年前；迁移率低且高度不对称，北到南m₂₁=0.00035（N_{m NorthSouth}≈0.37），南到北m₁₂=0.00008（N_{m SouthNorth}≈0.67），符合长期隔离和精细群体结构。
（4）遗传变异与形态分化的关系（Relationship between genetic variation and morphological differentiation in O. psychotrioides）：通过叶形态的非度量多维标度前两轴均值与Gower距离构建形态距离（Distmorph），其与F_ST的线性相关最强（调整R²=0.425，F=33.6，P=7.23E?07，AIC=?183.59），强于地理距离（调整R²=0.219）和历史距离（调整R²=0.193）；残差分析显示控制地理距离后形态距离仍显著（调整R²=0.274，P=0.0001311）。
（5）古分布与气候范围变化（Changes in the paleodistribution and climatic range of O. psychotrioides）：通过Maxent生态位模型重建过去120,000年分布，模型拟合良好（平均测试AUC=0.748±0.177；partial ROC平均AUC比=1.30，P<0.001）。末次间冰期（约120 ka）适宜条件集中于当前分布南缘；末次冰盛期（约21 ka）适宜栖息地大幅缩小并限于南部；Heinrich Stadial 1（17.0—14.7 ka）适宜条件向北扩展至当前范围；B?lling–Aller?d（14.7—12.9 ka）和Younger Dryas（12.9—11.7 ka）期间栖息地连续性破坏并形成两个主要气候隔离区；全新世至今适宜区扩大并重新连接先前破碎区域。
讨论部分总结：讨论指出，该研究揭示了以往中美洲云雾林谱系地理研究较少关注的精细地理尺度分化，Sierra de Chiconquiaco作为先前未被植物谱系地理识别的屏障，南北谱系分化明显；基于Wright公式和dadi的基因流估计均显示每代迁移者远少于1且不对称，可解释Teziutlán等群体的混合格局。低遗传多样性和轻微杂合子缺失与长期隔离、精细群体结构和有限重组一致；遗传与形态距离相关，但讨论明确指出这不意味着直接因果，可能由适应机制或漂变、历史隔离等中性过程共同造成。云雾林地质、土壤、海拔和气候的短距离剧烈差异会强化隔离，昆虫传粉限制花粉流、鸟类传播种子未能克服历史和 contemporary 屏障，进一步维持精细结构和低多样性。讨论还强调历史气候振荡、人口扩张和云雾林破碎化塑造了该物种进化历史，并使其对持续栖息地丧失具有脆弱性；在更广尺度上气候波动和地理屏障驱动谱系形成与区域分化，在更细尺度上相同过程维持遗传与生态多样性动态镶嵌。
研究结论翻译：尽管分布范围有限，Ocotea psychotrioides具有清晰的谱系地理结构，在Sierra de Chiconquiaco南北存在不同进化谱系。群体间遗传距离最好由末次冰期期间的气候条件解释，突显历史气候在塑造谱系分化中的作用。这些基于SNP的分析揭示了以往研究在很大程度上忽视的精细尺度遗传分化格局，表明即使狭域分布物种也可能拥有意外复杂的进化历史。通过揭示这种隐藏结构，研究人员的发现强调亟需评估这些谱系在当前和未来气候变化下的适应与保护意义。