长期围封保护对草原蝗虫群落产生复合效应

《Integrative Zoology》：Decades of Enclosure Protection Exert Composite Effects on Grassland Grasshopper Communities

　　

摘要：蝗虫是评估草原生态系统健康的重要指示类群。尽管已有大量研究探讨围封保护（exclosure protection）或过度放牧（overgrazing）对蝗虫多样性的影响，但不同生境压力下——尤其是扩散能力有限的物种——蝗虫是否表现出表型可塑性（phenotypic plasticity）与遗传适应（genetic adaptation）尚不清楚。本研究系统比较了20~40年围封保护与过度放牧对蝗虫群落多样性的影响；并以低扩散能力蝗虫Chorthippus fallax为模型，比较两种管理方式下卵发育速率，并阐明其潜在遗传结构与分化模式。结果显示，围封保护显著重塑蝗虫群落结构并导致物种种群动态呈物种特异性变化，但未改变Shannon–Wiener指数或物种丰富度。来自围封样地的C. fallax卵发育历期显著长于过度放牧样地（26.34天 vs. 25.74天；p < 0.05）。种群遗传学分析显示两组间无明显遗传结构或分化（FST< 0.05）。通过BayeScan及基于FST的方法共鉴定出389个可能与环境适应相关的候选离群位点（candidate outlier loci）。后续进化模拟表明，在假定条件下（FST= 0.25；s = 0.05），这些位点在两群体间约经162 ± 19年可产生强遗传分化。本研究不仅揭示了蝗虫对围封保护的表型与遗传响应，也为制定可持续草地管理策略提供了科学依据。

《长期围封保护对草原蝗虫群落产生复合效应》——论文解读

本文发表于《Integrative Zoology》。草原覆盖全球陆地面积20%~40%，是最重要的陆地生态系统之一，放牧是塑造草原生态系统的关键人为干扰因子。过度放牧会降低植物多样性、增加裸地、升高地表温度并引发级联生态效应，而围封保护（exclosure protection，即围栏禁牧）作为常见的草地恢复措施可提升植被盖度与质量。蝗虫（grasshopper）因对环境变化敏感，被广泛用作草原生态健康的指示类群（indicator taxon）。已有研究关注放牧或围封对蝗虫多样性及肠道微生物、生态位分化的影响，但对于长期围封保护下——特别是扩散能力有限（low mobility）的蝗虫种类——是否在表型及基因组水平发生适应性响应（adaptive response），尚缺乏系统认识。本研究以内蒙古草原长期围封样地（20年和40年）与相邻过度放牧样地为对象，调查蝗虫群落多样性，以低扩散能力物种Chorthippus fallax（直翅目：蝗科）为模式种开展共同花园实验（common garden experiment）比较卵发育差异，并利用特异性位点扩增片段测序（Specific-Locus Amplified Fragment Sequencing，SLAF-seq）进行种群遗传学分析及离群位点检测，结合进化模拟预测分化时间，旨在揭示围封保护的复合生态–进化效应。

研究人员在内蒙古草原生态系统研究站选取两个围封保护样地（exA，围封始于1979年，至2019年为40年；exC，围封始于1999年，至2019年为20年）及邻近的两个过度放牧对照样地（ovA、ovC），最大样点间距8 km。于2019年6—9月采用扫网法（sweep-netting method）调查蝗虫群落，仅以成虫进行群落多样性分析，计算Shannon–Wiener指数、物种丰富度、群落丰度，并用非度量多维尺度分析（Non-metric Multidimensional Scaling，NMDS）及嵌套置换多元方差分析（Permutational Multivariate Analysis of Variance，PERMANOVA）评估群落组成差异。采集C. fallax孕雌成虫于实验室控温饲养获得卵荚，取卵开展共同花园实验，比较围封与过度放牧来源种群的卵发育历期与孵化率。同时采集四样地C. fallax个体提取基因组DNA，进行SLAF-seq文库构建与Illumina测序，将 reads比对至飞蝗（Locusta migratoria）参考基因组，双软件（GATK和SAMtools）交集call SNP并过滤得到高质量SNP。利用ADMIXTURE、主成分分析（Principal Component Analysis，PCA）和最大似然系统发育树评估种群结构，用VCFtools计算遗传多样性指标（多态位点比例Poly、观察杂合度H_O、期望杂合度H_E、近交系数F_IS、核苷酸多样性π）及种群分化指数（F_ST）。采用BayeScan与F_STtop 5%两种方法交叉鉴定受选择候选离群位点。基于候选位点，在考虑自然选择、迁移、遗传漂变和突变的框架下开展正向进化模拟，定义F_ST连续5代超过0.25时为显著分化，估算不同选择系数（s）和有效种群大小（N_e）下的预期分化时间。

3.1 Effects of Enclosure Protection on Grasshopper Community Structure（围封保护对蝗虫群落结构的影响）

通过NMDS和PERMANOVA分析发现，围封与过度放牧样地间Shannon–Wiener指数和物种丰富度无显著差异，但过度放牧样地物种丰富度方差显著更大，总蝗虫丰度显著高于围封样地（平均72.33 vs. 29.33，p < 0.01），NMDS显示群落组成明显分离（R2 = 0.3186，p < 0.01）。12个采集物种中，部分物种仅出现于一类生境；各物种对管理的响应具物种特异性——4种在过度放牧样地显著更多，2种在围封样地显著更多，其余6种无显著差异。结论：围封保护显著降低蝗虫总丰度并重塑群落结构，但不影响Shannon–Wiener指数与物种丰富度，且不同蝗虫物种对栖息地管理的响应存在明显种间差异。

3.2 Effects of Enclosure Protection on Egg Development in Chorthippus fallax（围封保护对C. fallax卵发育的影响）

共同花园实验显示，围封来源C. fallax卵发育历期（26.34天）显著长于过度放牧来源（25.74天）（p < 0.05），而孵化率无显著差异。结论：C. fallax在两种生境间种群数量无显著差异，但通过表型可塑性表现为卵发育速率的差异——较短的发育历期可能是对过度放牧草地较高地表温度的生活史适应。

3.3 SLAF Sequencing of the C. fallax（C. fallax的SLAF测序）

SLAF-seq共获得约38亿条reads（479.06 Gb），平均每个样本约2.00 Gb，Q30为95.73%，经过质控和SNP过滤后保留75 682个高质量SNP用于后续群体遗传学分析。

3.4 Population Structure, Genetic Diversity, and Differentiation（种群结构、遗传多样性和分化）

ADMIXTURE分析交叉验证误差随K值持续上升且无明显拐点，K=1~3时均无法按生境类型区分个体；PCA前三个主成分仅解释约2.8%的遗传变异；最大似然系统发育树显示不同生境个体高度混合；两群体间整体F_ST= 0.016（两两群体间F_ST范围为9.97×10?3~0.012），属极低分化。遗传多样性方面，ovC和exC的期望杂合度H_E（分别为0.253和0.252）显著低于ovA和exA（均为0.255），其余指标（Poly、H_O、F_IS、π）无显著差异。结论：经过20~40年不同管理，C. fallax围封与过度放牧种群在全基因组水平未产生明显遗传结构和分化，两群体间基因流可能抑制了全基因组分化。

3.5 Identification of Candidate Outlier Loci（候选离群位点鉴定）

以BayeScan（q ≤ 0.05）和F_STtop 5%两种方法交集鉴定出389个候选受选择离群位点。基于这389个位点构建的系统发育树在群体水平出现一定程度聚类——一支以围封个体为主，另一支以过度放牧个体为主，此模式在全基因组SNP树中未出现。结论：虽然全基因组未分化，但在与环境适应相关的局部位点可能存在初步分化信号。

3.6 Evolutionary Simulation（进化模拟）

以389个候选离群位点为输入进行正向模拟，当选择系数s = 0.05、两群体N_e均为1×10?时，达到F_ST= 0.25所需时间为162 ± 19年；s = 0.005时需约1436 ± 175年；若过度放牧群体N_e(ov) = 5×10?、N_e(ex) = 1×10?，s = 0.05时约为177 ± 29年。模拟数据构建的系统发育树显示两群体聚类更清晰。结论：受选择位点在当前（数十年尺度）尚未完成显著分化，但在持续选择压下预计需百余年至数百年可形成较强分化，分化速率取决于选择系数与有效种群大小。

讨论部分指出，本研究未发现围封与过度放牧间蝗虫Shannon–Wiener指数和物种丰富度差异，可能与研究区局地河流缓冲微气候及资源限制有关，未来需多年份分站点采样验证。围封显著改变群落结构且降低总丰度，过度放牧样地某些蝗虫种类丰度更高，提示持续重度放牧可能创造利于特定种类暴发的生境，建议采用划区轮牧（rotational grazing）或季节休牧以降低蝗灾风险。C. fallax卵发育历期差异反映了低扩散能力昆虫可通过生活史表型可塑性响应生境差异。全基因组未检测到明显分化可能与样点距离近、基因流强有关；389个候选离群位点暗示局部适应性选择信号，进化模拟提示真正显著的遗传分化需更长的时间尺度——说明当前围封保护在生态层次已实现群落重构，但遗传层次的适应性分化需更持久且稳定的选择压方能显现。

Conclusion（结论）：

长期（20~40年）围封保护显著改变了蝗虫群落结构并使物种种群动态呈物种特异性变化，但未显著影响Shannon–Wiener指数或物种丰富度。围封保护样地来源的C. fallax卵发育历期显著长于过度放牧样地（平均26.34天 vs. 25.74天）。种群遗传学分析显示围封与过度放牧种群间无明显遗传结构或分化（F_ST< 0.05）。经BayeScan和基于F_ST方法鉴定出389个候选受选择离群位点，进化模拟表明其在假设条件下（s = 0.05，F_ST= 0.25）约需162 ± 19年可产生强遗传分化。上述发现为草原生态系统保护和草地放牧管理政策制定提供了重要理论基础。