天空岛屿西方大力神甲虫(Dynastes grantii)种群的人口统计学轨迹差异与季风季节性多基因适应

《Molecular Ecology》：Differential Demographic Trajectories and Polygenic Adaptation to Monsoon Seasonality in Sky Island Hercules Beetle Populations of the American Southwest

　　

美国西南部的天空岛屿山脉是被干旱低地屏障包围的孤立生境片段，构成了研究种群分化、种群历史和局部适应的有力自然实验。研究人员利用西方大力神甲虫(Dynastes grantii Horn, 1870)四个主要谱系的全基因组重测序(whole-genome resequencing, WGS)数据，刻画了栖息地破碎化与气候驱动隔离的基因组后果。种群结构分析揭示了与距离隔离(isolation-by-distance, IBD)相一致的强烈层级遗传分化，而ABBA-BABA检验几乎未检测到谱系间分化后基因流的证据。溯祖建模揭示所有谱系在末次冰期(Last Glacial Period)经历了同步种群扩张，随后在全新世(Holocene)发生严重收缩；其中犹他州(Utah)种群经历了最显著的瓶颈，并持续维持较低的有效种群大小(effective population size, Ne)。犹他州比较中升高的F_ST但相当的D_XY，以及升高的纯合片段(runs of homozygosity, ROH)表明，其遗传独特性反映了加速漂变而非古老隔离。基因组-环境关联分析鉴定了1367个与降水季节性相关的高置信候选位点，其中以最湿季度降水量驱动的多基因适应信号最强。季风湿度格局是塑造基因组分化的主要选择轴。遗传偏移(genetic offset)预测将莱蒙山(Mt. Lemmon)确定为未来气候情景下最脆弱的种群；遗传拯救供体分析确定了辅助基因流(assisted gene flow)的种群特异性来源。研究结果表明，将种群历史推断、分化基因组学与景观基因组学整合全基因组数据（包括历史博物馆标本），可同时增进对进化历史的理解，并为未来的可操作保护策略提供信息。

论文解读：美国西南部天空岛屿西方大力神甲虫(Dynastes grantii Horn, 1870)种群的进化历史与气候适应

### 研究背景与科学问题

美国西南部天空岛屿(sky island)是由高海拔森林生境片段镶嵌于干旱低地中的独特生物地理系统。更新世冷期时，中生森林曾连续分布；全球变暖与干旱化使这些森林向更高海拔退缩，形成如今被低地阻隔的“岛屿”状分布。这种地理破碎化为研究种群分化、基因流与局部适应提供了天然实验。然而，以往的研究多集中于鸟类、哺乳动物和爬行动物，昆虫作为生态功能的重要支撑者，其基因组研究严重不足。西方大力神甲虫(Dynastes grantii)是分布于该区域高海拔森林的优势甲虫，幼虫依赖腐烂木材，成虫飞行能力较弱，难以跨越低地屏障，因而是研究破碎化如何影响种群遗传的合适模型。基于此，研究人员利用全基因组重测序数据，重点回答以下问题：种群遗传结构如何分布？谱系间是否存在基因流？种群历史轨迹如何？哪些环境因素驱动适应分化？未来气候情景下哪些种群最脆弱，哪些可充当基因拯救供体？

### 主要技术方法

样本队列包含2023–2024年间在美国亚利桑那州多个天空岛屿采集的当代个体，以及来自菲尔德博物馆(Field Museum)、西南生物学博物馆(Museum of Southwestern Biology)、台湾生物多样性研究博物馆及私人收藏的历史标本。最终聚焦四个主要谱系：Mogollon Rim（n=12，含Payson与Reserve两采样点）、Mt. Lemmon（n=5）、Chiricahua Mt.（n=5）和Utah（n=5）。DNA提取后利用Illumina平台进行全基因组重测序，比对至D. grantii参考基因组。分析方法包括：PCA与ADMIXTURE推断种群结构，D统计量与Mantel检验评估基因流与IBD；SMC++和GONE2分别重建深层与近期Ne轨迹，δaδi比较严格隔离与隔离加种群大小变化模型；pixy计算F_ST与D_XY，BCFtools检测ROH；RDA与LFMM2进行基因组-环境关联分析，并借助topGO进行基因本体富集；遗传偏移与供体识别基于约束RDA与未来气候预测。

### 研究结果

**3.1 种群结构与基因流**
PCA显示PC1分离犹他州种群，PC2将Mogollon Rim与南部天空岛屿区分，Mt. Lemmon与Chiricahua Mt.也明确分开。ADMIXTURE在K=4时呈现与地理分布高度一致的祖先成分。ABBA-BABA检验显示四个主要谱系间基因流信号极弱，Mantel检验表明遗传分化与地理距离显著正相关（r=0.79, p<0.0001），支持IBD格局。Mogollon Rim个体表现为Mt. Lemmon与Chiricahua Mt.的混合祖先，但主要谱系间的分化后基因流证据不足。

**3.2 种群历史与溯祖建模**
SMC++揭示所有谱系在末次冰期（约1万–2万年前）经历扩张，随后在全新世发生快速且剧烈的收缩；犹他州种群收缩最为严重并维持最低Ne。GONE2进一步显示近300年内Mogollon Rim具有最大Ne，犹他州始终保持最低。δaδi模型支持严格隔离模型，且各谱系分歧时间较浅，说明分化主要发生于异域隔离状态。

**3.3 基因组分化模式**
犹他州种群在成对比较中F_ST显著升高，但D_XY与其他谱系间无显著差异，表明该种群的遗传独特性主要源于快速漂变而非长期隔离。ROH分析显示犹他州具有最高的总ROH长度(SROH～96.0 Mb)和近交系数(F_ROH～0.140)，进一步支持其经历过严重瓶颈和近期近交。

**3.4 基因组-环境关联与未来脆弱性**
RDA模型显著且解释约6.3%的基因组变异；LFMM2鉴定出20,867个显著候选SNP，其中最湿季度降水量(precipitation of the wettest quarter, PW)关联的位点达17,448个。两种方法交集获得1367个高置信候选位点，GO富集分析发现这些位点富集于皮层肌动蛋白细胞骨架组织、N-聚糖加工、细胞连接维持等与水分胁迫相关的生物过程。遗传偏移预测在所有情景下Mt. Lemmon均具有最高偏移值，判定其为最脆弱种群。供体分析显示Mt. Lemmon的最佳遗传供体来自Chiricahua Mt.，Chiricahua Mt.的最佳供体来自Mogollon Rim，Mogollon Rim的最佳供体来自犹他州，而犹他州的最佳供体落于Mt. Lemmon附近。

### 讨论与结论

讨论部分指出，Mogollon Rim因其较连续的栖息地和较大Ne，可作为整个区域现有遗传变异的重要储备库；犹他州等小型孤立种群则面临高近交和更大脆弱性。遗传偏移预测应视为可检验的优先级假设而非绝对风险判断，辅助基因流需谨慎实施，避免因不当移殖导致远交衰退，并强调需结合野外实验与更密集采样验证。研究人员最终记载了三个关键模式：破碎化驱动的遗传漂变、对季风降水季节性的多基因适应信号，以及历史连通性对现有遗传变异的塑造作用。这些结果表明，将全基因组数据与种群历史、分化基因组学及景观基因组学相结合，能够同时揭示天空岛屿昆虫的进化历史，并为应对未来气候变化下生物多样性保护提供可操作的科学依据。