中亚平宁山脉天空岛避难所动态对比：两种特有跳甲(Alticini)的比较系统地理学(Comparative Phylogeography)

《Journal of Biogeography》：Contrasting Refugial Dynamics in Central Apennine Sky-Islands: A Comparative Phylogeography of Two Endemic Beetles

　　

摘要：第四纪(Quaternary)气候振荡深刻塑造了高山生物区系的遗传结构与地理分布。理解嗜冷(cold-adapted)类群对过去气候变迁的响应，对于揭示山地生物多样性形成过程至关重要，然而嗜冷类群可能遵循不同的避难所(refugial)历史。研究人员比较了两种同域分布的高海拔特有跳甲(flea-beetles)，以检验两种相互替代的生物地理模型：(1)冰期幸存于冰斗外围无冰山峰(nunatak，即冰川间裸露山顶)；(2)幸存于冰川外围边缘避难所(peripheral glacial refugia)。研究地点为意大利中亚平宁山脉(Central Apennines, Italy)。研究对象为高海拔(海拔1700–2700 m a.s.l.)特有跳甲 Longitarsus springeri 与 Psylliodes biondii。研究人员在21个山地样点采集112个个体，整合线粒体与核标记数据，进行系统发育最大似然(Maximum Likelihood, ML)树、单倍型网络(haplotype-network)、分歧时间估算(BEAST)、群体历史动态推演(Extended Bayesian Skyline Plot, EBSP)及空间遗传分析(Spatial Analysis of Molecular Variance, SAMOVA)。结果显示两物种表现出截然不同的分化模式、历史动态与基因流。L. springeri 呈现深层系统地理结构，具11个山体特异性谱系(lineages)及群体数量稳定性，符合长期孤立于类似冰斗(nunatak-like)避难所中持续存活；相反，P. biondii 遗传结构较浅，最高遗传多样性出现于大萨索山块(Gran Sasso massif)，群体在末次冰盛期(Last Glacial Maximum, LGM)前(~300 kya，海洋同位素阶段7，MIS 7)发生扩张，符合先存活于边缘避难所再向四周扩散的模式，此次前LGM扩张延伸了经典边缘避难所模型的时间范畴。主要结论：即使是生态相似、同域分布的特有种也可对气候与环境变化产生分化响应。L. springeri 与 P. biondii 的对比历史表明，中亚平宁山脉地形复杂性与冰川动力学如何塑造了高山生物区系的谱系存续、分布动态与遗传结构，以及同一景观内避难所动态的空间与时间维度可解耦(decouple)。中亚平宁天空岛体系既是宏观长存避难所(macro-refugium)，也是驱动谱系分化的微避难所(microrefugia)精细网络。比较系统地理学框架为理解南欧山脉天空岛体系如何产生并维持生物多样性提供有力途径。

本文发表于《Journal of Biogeography》，针对第四纪气候振荡下高山嗜冷(cold-adapted)类群避难模式存在"冰斗原地存活(nunatak refugia)"与"边缘低海拔避难后回迁(peripheral refugia)"之争，且以往欧洲系统地理学研究多聚焦温带类群，对高山协同分布近缘特有种对比研究不足。中亚平宁山脉(Central Apennines)具典型天空岛(sky-island)地形与高特有率，两种高海拔(1700–2700 m a.s.l.)特有跳甲(Alticini)——Longitarsus springeri（寄主菊科 Asteraceae 之 Senecio squalidus ssp. rupestris）与 Psylliodes biondii（寄主十字花科 Brassicaceae 之 Isatis apennina 与 Erysimum pseudorhaeticum）——同域分布却具不同物候与寄主专一性，是检验避难假说的理想对象。研究人员通过多基因标记与多类系统地理与群体遗传学分析，证实两协同分布近缘种分别符合冰斗避难与原边缘避难扩散两种对立模型，揭示同一景观内避难动态的时空解耦现象。此比较系统地理学(comparative phylogeography)研究表明中亚平宁山脉兼具宏观避难所与微避难所网络功能，对认识山脉生物多样性起源与保护具重要意义。

主要关键技术方法概述：研究人员于中亚平宁山脉21个山地位点采集 L. springeri（73个体，13个位点）与 P. biondii（39个体，8个位点），标本经形态学鉴定后提取基因组DNA，PCR扩增两个细胞色素氧化酶I亚基(cox1)线粒体片段及四个单拷贝核基因(CAD、Cv2、Rad50、Wg)并双向测序，利用PHASE算法推导配子相位；以最大似然法(IQ-TREE)重建种内系统发育关系，用统计简约法(TCS)与中位连接法(Median-Joining)构建单倍型网络并计算单倍型多样性(h)与核苷酸多样性(π)；采用BEAST软件基于cox1严格分子钟估算主要线粒体谱系分歧时间；运用离散生物地理学分析结合贝叶斯随机搜索变量选择(BSSVS)推断祖先分布区与扩散路线；以SAMOVA分析空间分子方差识别最优群体分组(K)；通过扩展贝叶斯天际线图(EBSP)基于线粒体和核基因联合重建有效群体大小(N_e)随时间变化。

研究结果：

3.1 Phylogenetic Analyses（系统发育分析）

基于cox1联合序列构建ML树与单倍型网络，L. springeri 识别出11个互为单系(monophyletic)的线粒体单倍型组，各对应一独立山块，TCS分析聚为4个子网，核标记呈弱南北结构，部分位点具较高单倍型多样性；SAMOVA显示Φ_CT随分组数K增大而升高，K=10时每个山峰为一独立群(Φ_CT=0.916)。P. biondii 识别出10个互为单系线粒体单倍型组但网络连通性更强、分化浅，核标记几乎无结构；SAMOVA最优K=6(Φ_CT=0.611)，大萨索山块内种群归为一组。BEAST时间标定显示 L. springeri 最近共同祖先(MRCA)约888 kya(95% HPD: 1.38–0.50 Mya)，最早分歧至马泰塞地块(Matese Massif)，随后向北西扩散；P. biondii MRCA约353 kya(95% HPD: 527–218 kya)，定位于大萨索山块，约326–239 kya自大萨索向外扩散。以上表明 L. springeri 具深部分化与山体隔离特征，P. biondii 具较浅分化并以大萨索为扩散中心。

3.2 Population Genetic Structure and Demographic Analyses（群体遗传结构与历史动态分析）

SAMOVA定量佐证 L. springeri 强空间结构(Φ_CT最高达0.932)与 P. biondii 较弱但显著结构(Φ_CT=0.611)。EBSP显示 L. springeri 有效群体大小在整个晚第四纪呈平稳(stationary)模式；P. biondii 约250–300 kya始显著群体扩张，约75 kya始下降。表明前者符合冰斗(nunatak)避难所长时期稳定孤立存活，后者符合边缘避难所(peripheral refugium)先保存后中更新世(MIS 7)扩张再间冰期破碎的历史。

讨论与结论翻译：

协同分布分类群未必共享相同生物地理历史。L. springeri 具深层地理结构化遗传分化、11个山块特有互为单系线粒体单倍群、核标记呈南北分异及长期群体稳定，符合冰斗(nunatak)避难所假设——各山峰在无冰期作为离散微避难所(microrefugia)保存古老谱系并通过更新世循环隔离促进分化。P. biondii 遗传结构较浅、大萨索山块具最高核苷酸与单倍型多样性、贝叶斯系统地理学与EBSP示其曾存于大萨索所在边缘避难所并于约300 kya(MIS 7)发生群体扩张随后再分布，符合边缘避难所(peripheral refugia)模型，且该扩张早于经典末次冰期(LGM)后回迁认知。两物种对比说明即使生态相似共域高山特有种亦可对气候波动产生不同进化轨迹；中亚平宁天空岛体系同时发挥长存宏观避难所与微避难所网络功能，空间与时间避难动态可在同一景观解耦(decouple)。比较系统地理学为解析南欧山脉天空岛生物多样性成因提供强有力框架。