大西洋森林山地特有凤梨科Pitcairnia flammea沿海拔梯度的种群基因组分化与种群历史

《Heredity》：Population genomic differentiation and demographic history across an elevational gradient in an endemic bromeliad of the Atlantic Forest mountains

　　

山地海拔梯度为研究环境异质性、地理隔离和历史气候动态如何塑造种群历史、种群分化和多样化提供了天然实验室。地形复杂性和历史气候振荡是全球山地生物多样性格局的主要驱动因素。巴西大西洋森林（Brazilian Atlantic Forest, BAF）是热带生物多样性热点，以天然片段化的山地生境为特征，为研究这些过程提供了合适系统。此处，研究人员利用分布于海拔梯度上的Pitcairnia flammea Lindl.（Bromeliaceae）种群的1000多个单核苷酸多态性（single nucleotide polymorphisms, SNPs），研究遗传多样性、种群结构和种群历史格局。贝叶斯（Bayesian）和基于溯祖（coalescent）的分析揭示了强种群结构，两个主要遗传类群大致表现为高海拔种群与低海拔和中海拔种群的对比，尽管海拔本身并未被直接支持为遗传分化的驱动因素。种群历史建模支持两个主要遗传类群之间早期不对称基因流，从类群2向类群1的迁移较高，随后发生隔离。分化时间估计和有效种群大小（effective population size, Ne）的异质性变化进一步提示中更新世和晚更新世（Middle and Late Pleistocene）期间存在不同的种群历史轨迹。总之，结果揭示了单一特有物种内的分化进化轨迹，并提示大西洋森林山地的种群分化反映了生境片段化、历史连通性和隔离、遗传漂变（genetic drift）以及与更新世气候振荡潜在相关的种群变化的联合效应。

研究背景：山地海拔梯度是研究环境异质性、地理隔离和历史气候动态如何塑造种群分化与多样化的天然系统。地形复杂性、环境梯度和生境片段化可通过空间隔离和降低基因流促进种群分化，历史气候振荡亦被认为是热带山地多样性格局的重要驱动因素。巴西大西洋森林（BAF）是热带生物多样性热点，其岩石露头和孤山等自然片段化生境具有高度地形与环境异质性，海拔跨度大且温降水变化显著。Pitcairnia flammea Lindl.（Bromeliaceae）广泛分布于BAF自然隔离的岩石露头和孤山，具有较高形态生理变异、强遗传结构和有限基因流，是研究海拔梯度及更新世气候振荡如何影响种群历史与遗传分化的合适模型。已有研究报道该种沿海拔梯度存在分化生态策略，但对其种群遗传结构、主要遗传类群分化时间、历史基因流及有效种群大小（Ne）变化仍缺乏整合分析。研究人员因此采集BAF沿海拔梯度的8个P. flammea种群，利用1000多个单核苷酸多态性（SNPs）研究遗传多样性、种群结构、基因流和种群历史，旨在回答遗传变异如何沿海拔梯度结构化、主要遗传类群何时分化及何种种群过程表征其分化、推断的种群变化是否与更新世气候波动一致、以及遗传变异和种群历史如何与已报道生态策略相关。研究结论显示强种群结构，两个主要遗传类群大致对比高海拔与低中海拔种群，但海拔本身未被直接支持为遗传分化驱动因素；早期不对称基因流后发生隔离，Ne呈异质性变化；生境片段化、历史连通性与隔离、遗传漂变及与更新世气候振荡潜在相关的种群变化共同作用。其意义在于揭示单一特有物种内不同种群可沿不同进化轨迹分化，并深化对BAF山地特有种多样性形成过程的理解。

关键技术方法（不超过250字）：样本队列来源为BAF沿海拔梯度（海平面至2140 m）的8个P. flammea种群，共159个个体（每种群约20个）。研究人员采用限制性位点相关DNA测序（restriction site-associated DNA sequencing, RAD-seq）获取SNPs，构建四套过滤数据集，并进行异常位点检测。种群结构分析使用STRUCTURE、ADMIXTURE、主成分判别分析（discriminant analysis of principal components, DAPC）、SplitsTree邻接网络和有效迁移表面估计（EEMS），并以Mantel检验评估隔离距离（isolation by distance, IBD）和隔离海拔。种群树用SNAPP推断。基因流动态与Ne变化用FASTSIMCOAL基于位点频谱（site frequency spectrum, SFS）建模，Ne历史另用STAIRWAY PLOT 2推断；分化时间校准采用3×10^?9替换/位点/世代及5年世代时间。

研究结果：
Genetic diversity and structure（遗传多样性与结构）：通过STACKS计算的多样性指标显示，高海拔种群MAR和PAP遗传多样性最高，私有等位基因数、多态位点比例、核苷酸多样性及期望和观察杂合度均较高；中海拔MAC多样性模式与高海拔种群相似；各种群近交系数（F_IS）均低，且高海拔趋于更低。通过STRUCTURE（ΔK=2）、DAPC（K=4）、ADMIXTURE、SplitsTree网络、EEMS和Mantel检验，研究发现强遗传结构：STRUCTURE将种群分为Group 1（PAP、MAR、MAC）与Group 2（PES、ITA、COR、RAN、UBA）；DAPC进一步区分MAC及MAR、PAP；ADMIXTURE显示除地理邻近的低地RAN和COR未分化外，各种群多形成独立聚类；SplitsTree显示两个主要分裂且各种群多为独立遗传单元；EEMS支持当前基因流有限；Mantel检验未发现地理距离或海拔距离与遗传距离显著相关，控制地理距离后海拔亦不显著。

Population relationships（种群关系）：通过SNAPP对8个种群构建贝叶斯种群树，研究恢复4个主要分支，与结构分析大体一致。结果显示高海拔MAR和PAP与较低海拔种群约在275,000年前发生初始分裂（95% CI：226,000–325,000）；随后高海拔种群间约在150,000年前分裂，低中海拔种群约在60,000年前分化。MAC（959 m）显示连接两个主要类群的替代拓扑，约在250,000年前从Group 2祖先谱系分化（95% CI：180,000–297,000），但其在STRUCTURE中归入Group 1，反映其遗传多样性和结构与MAR、PAP相似。除MAC外，SNAPP关系与结构分析一致；分化估计被原文视为近似，因依赖Puya raimondii的突变率和Pitcairnia的假定世代时间。

Gene flow dynamics and demographic modeling（基因流动态与种群历史建模）：通过FASTSIMCOAL比较无基因流、早期基因流、持续基因流和近期基因流模型，研究发现早期基因流模型对Group 1与Group 2之间基因交换支持最强（AIC=55,035.09，最大似然=422.83，ΔAIC=71.65）。参数估计显示两组间基因流约在67,000年前停止（CI=2705–76,970），且明显不对称，Group 2向Group 1的迁移率（7.3×10^?5）高于反向迁移率（1.3×10^?6）。通过FASTSIMCOAL对ADMIXTURE识别的7个遗传簇进行Ne变化建模，研究发现多数簇中收缩或扩张模型优于稳定模型；MAR、ITA、UBA、PAP和RAN/COR以收缩模型最佳，PES和MAC以扩张模型最佳，但UBA和PAP支持较弱。当前Ne估计为75,685–495,326个二倍体个体，祖先Ne为5,615–46,737,128个二倍体个体，种群事件时间为18,275–937,133年前。通过STAIRWAY PLOT 2，研究发现多数种群在最近100,000年内发生Ne变化并多呈近期下降：MAR约750,000年前开始扩张后趋稳并在最近20,000年内略降；UBA约500,000年前开始扩张；RAN和COR约100,000年前开始下降并持续至现今；ITA约500,000年前Ne降低后逐渐下降并趋稳；PES在20,000年前后显著下降，约与末次盛冰期吻合；PAP类似但近期下降更强；MAC约300,000–150,000年前Ne增加，约150,000年前下降，并从约50,000年前至现今趋稳。

讨论部分总结：讨论指出，P. flammea存在强遗传分化和异质种群历史轨迹，两个主要遗传类群大致对比高海拔与中低海拔种群，但海拔本身未被直接支持为遗传分化驱动因素；Group 1遗传多样性较高，分化时间较早，种群历史建模支持Group 2向Group 1的早期不对称基因流，随后隔离，Ne变化具种群特异性。遗传多样性与结构方面，强结构与此前基于质体和核标记报道的有限基因流一致，岩石露头和孤山作为“天空岛”可促进隔离与分化；未检测到IBD或隔离海拔，可能反映采样分辨率限制及历史种群过程、生境片段化、扩散动态等作用；Pitcairnia种子扩散受限，偶发长距离扩散和水媒传播可能参与分布格局；Group 1较高多样性可能与较大Ne、长期种群存续或高地局部避难所有关，也可能受传粉者、花期和交配系统影响，但原文指出这些机制尚未检验。种群历史方面，SNAPP提示Group 2内种群关系更近，MAC位置存在替代拓扑；Group 1与Group 2约275,000年前分化，MAC约250,000年前从Group 2祖先分化，MAR和PAP约150,000年前分化，Group 2约60,000年前多样化，发生在中晚更新世；与Mota et al. (2020)报道的约2百万年前分化不同，可能源于采样空间尺度、基因组区域和方法差异；早期基因流模型最佳且约67,000年前停止，不对称基因流可能贡献Group 1较高多样性，并与“闪烁连通假说”一致；间断基因流可能类似岛屿“物种泵”过程；FASTSIMCOAL和STAIRWAY PLOT 2均支持Ne随时间变化，多数种群近期下降，但冷期下坡扩张不一定增加Ne，地形复杂山地中气候转变可导致局部扩张、收缩或片段化；MAC扩张获强支持且历史Ne较高、遗传多样性最高，但其系统位置不确定，值得进一步研究。遗传与生态分化方面，强遗传分化与同一海拔梯度上已报道的耐性—回避权衡并存：低海拔种群以节水性状避热，高海拔种群以耐热冷生理机制应对；热耐受更与群落多样性和系统发育背景相关而非海拔本身；综合生态生理、种群遗传结构和种群历史证据，沿海拔梯度的进化格局不能仅归因于海拔，而可能反映气候、生物和历史过程交互作用。

结论翻译：在测试框架下，结果揭示了大西洋森林山地沿海拔梯度分布的P. flammea种群存在强遗传分化和不同种群历史。种群历史建模支持两个主要遗传类群之间早期不对称基因流，随后发生隔离，而有效种群大小的异质性变化表明种群随时间发生种群特异性历史响应。结合先前记录的沿海拔梯度生态生理变异，研究提示P. flammea中观察到的进化格局不能仅归因于海拔，而可能反映气候、生物和历史过程的相互作用。特别是，生境片段化、受限连通性、遗传漂变以及与更新世气候振荡潜在相关的连通性时间变化，可能共同促成了观察到的种群分化和种群历史格局。总之，这些发现强调单一特有物种的种群可在环境异质和片段化的山地景观中遵循不同的进化轨迹。