农业景观及植被结构的时空变化影响新热带林下灌木Palicourea prunifolia基因组多样性的研究

《Perspectives in Plant Ecology, Evolution and Systematics》：Spatial and temporal changes in agricultural landscapes and vegetation structure affect genomic diversity of a Neotropical understory shrub

　　

摘要：集约化农业扩张导致的景观变化影响植物群落结构并最终影响遗传多样性。由于生活周期长，植物对持续景观变化的响应可能存在时滞(time lag)，特别是在集约化耕作方式近期才被引入的地区。研究人员采用景观基因组学(landscape genomics)方法探讨了景观变化及其时滞效应对Palicourea prunifolia（新热带森林林下专性物种）的空间和时间影响。研究人员利用26,349个SNP位点数据估算遗传参数，并采用溯祖分析(coalescent analyses)估算有效迁移面(effective migration surfaces)。研究人员利用地面激光雷达(LiDAR)扫描分析植被结构，并提取1985–2022年的景观组成与配置指标。结果表明：基因组多样性在森林覆盖较低、林冠较开阔的残存片段中更高；近交系数(inbreeding coefficient, f)在农业用地占比较高、林冠闭合度较高的残存片段中更高；所有种群间基因流极低；种群间基因组分化(genetic differentiation, FST)较高，且最好由早期景观配置（主要是牧场比例）解释。遗传多样性和近交系数未表现出对景观变化的时滞响应，但遗传分化已对历史景观变化产生响应。研究结果表明高度的生境丧失与片段化及大量农业生态系统已对P. prunifolia种群产生影响。恢复森林特别是提高景观连通性的修复措施，应作为促进P. prunifolia种群间基因流与遗传多样性、降低种群分化的优先策略。

论文解读：农业景观时空变化与新热带林下灌木Palicourea prunifolia基因组多样性——景观基因组学研究**

研究背景与意义

全球近半数宜居土地已被农业生态系统(agroecosystem)覆盖，集约化农业快速扩张是生物多样性危机的主要驱动因子之一，通过生境丧失、片段化、孤立化和景观简化导致物种丰富度下降及遗传多样性降低。对于生活周期较长的植物而言，其对景观变化的遗传响应可能存在时滞效应(time-lag effect)，即当前遗传结构反映的是历史上而非最近的景观格局。巴西塞拉多(Cerrado)是世界物种最丰富的热带稀树草原及生物多样性热点地区，但因农业扩张已丧失逾50%原生植被，然而农业景观变化如何影响塞拉多植物遗传多样性尚知之甚少。已有研究显示不同传粉/种子传播系统的树种响应各异，而林下灌木作为生境特化类群，其基因流受限且对片段化敏感，是探讨景观基因组学(landscape genomics)的理想对象。本研究以新热带林下灌木Palicourea prunifolia（茜草科Rubiaceae，河岸林片段特有种，花被大型蜂类和蝶蛾类 pollinated，果实为鸟播小型核果）为模式种，结合37年景观变化遥感数据与地面LiDAR植被结构测定、ddRADseq基因组测序，旨在阐明景观组成(composition)、配置(configuration)及植被结构对新热带林下灌木种群基因组变异的空间与时间效应，检验时滞假设，并为塞拉多破碎化生境保护提供遗传学依据。论文发表于《Perspectives in Plant Ecology, Evolution and Systematics》。

主要关键技术方法

研究人员自塞拉多农业景观长期生态研究(LAND LTER)区14个P. prunifolia种群采集150个成年个体叶片，经Qiagen DNeasy Plant Mini Kit提取DNA后行双酶切(EcoRI+MseI) ddRADseq建库并Illumina HiSeq 2000单端测序；用ipyrad进行去接头、解复用、聚类和denovo组装，VCFtools过滤保留双等位SNP、缺失率≤40%、每群体选9个低缺失个体得26,349个SNP；PLINK剔除连锁不平衡(LD)；用hierfstat算期望杂合度(expected heterozygosity, He)、观察杂合度(observed heterozygosity, Ho)、等位基因丰富度(allelic richness, Ar)和近交系数(inbreeding coefficient, f)，用Weir & Cockerham法算F_ST，sNMF做群体结构聚类；基于MapBiomas 2022土地覆被图设250 m–7000 m六级缓冲区算景观组成（森林%、农田%、牧场%、稀树草原%）与配置指标（斑块密度patch density PD、聚集度clumpy、边缘密度edge density ED、香农多样性指数Shannon Diversity Index、面积加权平均斑块大小MESH）；用地面LiDAR(FARO Laser Focus 70)扫10×10 m样方，Fusion/LDV算冠层起伏比(canopy relief ratio, CRR)及各层植被高度/密度；用AIC选最优空间尺度，线性回归分析景观/植被结构对基因组参数的影响，GLMM分析37年(1985–2022)景观变化时序β系数与遗传参数关系；用MMRR（多元矩阵随机化回归）分析历史各年份景观变量与成对F_ST的关系；用Estimated Effective Migration Surfaces(EEMS)可视化基因流屏障与廊道。

研究结果

3.1. Genomic structure（基因组结构）

研究人员基于26,349个SNP分析发现，所有种群期望杂合度He＝0.008–0.095、观察杂合度Ho＝0.010–0.078，近交系数f多数显著为正，整体F_ST＝0.396(p＝0.001)表明种群间基因组分化高；sNMF显示K＝4个基因簇且各群体存在混合，中部景观内某些种群对F_ST较低。结论：P. prunifolia具低基因组多样性与高种群分化特征。

3.2. Spatial and temporal effects of landscape changes and vegetation structure on genomic diversity（景观变化时空效应及植被结构对基因组多样性的影响）

空间效应分析显示：等位基因丰富度Ar和期望杂合度He随景观中森林覆盖比例(% forest)升高而降低（1000 m尺度，p＜0.05），观察杂合度Ho同理（5000 m尺度）；近交系数f随农业用地比例(% agriculture)升高而增加（500 m尺度，p＜0.05）。时序分析(1985–2022景观变化β系数)对基因组多样性及近交系数无显著影响。植被结构分析表明：Ar、He、Ho均随冠层起伏比CRR增大（林冠闭合度高、隙少）及林下层平均高度(strata mean elevation)增大（植被密度高）而降低；交互项显示开放林冠（低CRR）缓解森林覆盖率高对遗传多样性的负作用；f随CRR增大和% agriculture升高交互增加。结论：当代景观组成（低森林覆盖、低农业占比）与残存片段内较开阔林冠利于维持较高基因组多样性与较低近交；遗传多样性对景观变化无时滞响应。

3.3. Spatial and temporal effects of landscape changes on population genomic differentiation（景观变化时空效应对种群基因组分化的影响）

MMRR分析表明：种群间F_ST与历史（1985–2005年间）牧场比例(% pasture)变化呈正相关，2000年农业镶嵌比例(% mosaic)及2005年clumpy亦正向关联F_ST；地理距离及其他景观指标不显著。EEMS显示种群间有效迁移率极低（周边种群＜中心种群P03/P06/P08），有效多样性低，偏离隔离-by-距离模式。结论：历史景观组成变化（尤牧场扩张）削弱连通性致种群分化升高；基因流总体极低且无隔离-by-距离。

讨论与结论翻译

讨论部分总结：景观组成空间变异是P. prunifolia基因组多样性的主要驱动因子，与预期相反，低森林覆盖景观具更高基因组多样性，这与物种偏好具林窗/较开阔林冠河岸残存片段有关；农业比例升高伴较高近交系数，暗示传粉与食果鸟受基质阻隔致基因流降；历史牧场比例增致连通性失并提升F_ST，反映历史景观变化影响分化；未发现遗传多样性/近交的时滞效应，但分化响应历史格局；全群遗传多样性极低，危及长期存活力。作者建议保护并恢复原生植被以提高塞拉多景观连通性，促进P. prunifolia基因流与遗传多样性并降低分化。

研究结论（Conclusion）：景观组成显著影响P. prunifolia基因组多样性与种群分化；遗传多样性与近交系数对景观变化无时滞响应，但遗传分化受历史景观变化驱动；LAND景观高强度农业 mosaic下高度生境丧失与片段化已影响种群；较低农业占比景观近交系数为负指示异交频率高；需优先恢复森林并提高景观连通性以促进基因流、增遗传多样性并减种群分化。