《Molecular Phylogenetics and Evolution》：Genome-wide SNP data reveal geographic structure and landscape-associated genomic differentiation in a widespread lizard in arid Eastern Central Asia
编辑推荐：
干旱景观为研究地理结构和环境异质性如何塑造基因组分化提供了重要系统。然而，在地形复杂的荒漠地区，确定种群结构主要反映景观阻力、地理距离还是当代环境变异仍具挑战性。本研究利用全基因组SNP数据，对分布于干旱中亚东部（Arid Eastern Central As
干旱景观为研究地理结构和环境异质性如何塑造基因组分化提供了重要系统。然而，在地形复杂的荒漠地区，确定种群结构主要反映景观阻力、地理距离还是当代环境变异仍具挑战性。本研究利用全基因组SNP数据，对分布于干旱中亚东部（Arid Eastern Central Asia, AECA）的广布蜥蜴——变色麻蜥（Eremias vermiculata）的种群结构、系统发育关系、历史基因流、种群历史动态及基因组分化的景观关联因素进行了研究。对164个个体分析恢复了六个与主要荒漠盆地和山地限界区域相关的地理结构核基因组聚类。核基因组系统发育树解析出与物种分布范围东北部和西南部相对应的两个广泛区域支系，而主成分分析（PCA）和ADMIXTURE分析恢复了六个更细尺度的遗传聚类。基于来自相同个体的NCBI来源细胞色素b（Cyt b）和细胞色素氧化酶I（COI）序列联合构建的线粒体系统发育树恢复了四个更深的母系谱系。这些结果表明核基因组与线粒体数据集在系统地理格局上总体一致，而全基因组SNP提供了更细尺度的种群结构分辨率。种群历史重建进一步揭示了各聚类间晚更新世历史的区域性异质性，包括扩张、稳定和衰退信号。景观基因组分析表明，基因组分化主要与景观阻力（特别是海拔和土地覆盖）以及地理距离相关，而在控制空间结构后，当代环境变量解释的变异相对较少。综上，研究结果表明变色麻蜥的基因组分化反映了干旱中亚东部持续存在的景观格局、历史连通性和区域特异性种群历史之间的相互作用。更广泛而言，本研究凸显了整合系统地理学与景观基因组学方法对理解地形异质性荒漠生态系统中种群分化与进化历史的价值。
论文解读：全基因组SNP揭示干旱中亚东部广布蜥蜴的景观关联基因组分化
一、研究背景与科学问题
干旱地区覆盖地球陆地表面约41%，但其地质历史、地形复杂性和气候动态在不同区域间存在显著差异。在干旱中亚东部（Arid Eastern Central Asia, AECA），广布的山系、山间盆地、高原和极端干旱荒漠共同构成了高度异质的景观格局。该区域的地质演化深受印度–欧亚板块碰撞相关构造过程的影响，这一过程推动了中新世和上新世期间青藏高原及周边山系的隆升，进而通过雨影效应加剧了大陆性干旱，促进了内流盆地和广阔荒漠系统的发展。
在此复杂背景下，一个核心科学问题尚未得到充分解答：景观结构（如山脉阻隔、盆地隔离）与当代环境变异在塑造基因组分化中的相对贡献究竟如何？尽管此前已有研究对干旱适应类群的系统地理格局进行了探讨，但鲜有研究在统一基因组框架内同时评估种群结构、景观阻力、环境变异和种群历史动态。对于地形复杂的荒漠区域而言，物理屏障限制扩散并促进空间遗传结构的形成，而气候波动则改变栖息地连续性、种群规模和基因流，二者往往协同作用，难以区分各自贡献。
变色麻蜥（Eremias vermiculata）是一种广泛分布于中亚东部荒漠和半荒漠地区的蜥蜴科（Lacertidae）物种，其分布范围横跨多个荒漠盆地和山地限界区域，是研究干旱景观如何影响基因组结构的理想模型。此前基于线粒体DNA的系统地理学研究已识别出与主要荒漠盆地大致对应的地理结构化谱系，但线粒体DNA仅代表母系遗传历史，可能无法完整反映全基因组水平的种群结构、混合和连通性模式。
基于此，研究人员提出三个核心问题：（i）核基因组SNP数据是否能恢复与主要荒漠盆地和山地限界区域对应的地理结构化种群；（ii）核基因组模式是否与相同个体的线粒体谱系大体一致，或提供比线粒体数据更细尺度的种群结构分辨率；（iii）基因组分化与景观阻力（isolation-by-resistance, IBR）、地理距离（isolation-by-distance, IBD）或当代环境变异（isolation-by-environment, IBE）的关联强度如何。
二、主要技术方法
研究人员利用基因分型测序（Genotyping-by-Sequencing, GBS）技术，对来自中国和哈萨克斯坦共77个采样点的164个变色麻蜥个体进行全基因组SNP检测，经严格质控后获得13,906个高质量SNP位点。分析方法涵盖：主成分分析（PCA）、判别主成分分析（DAPC）和ADMIXTURE模型聚类解析种群结构；最大似然法（ML）和邻接法（NJ）构建核基因组系统发育树，并与相同个体的线粒体Cyt b和COI序列系统发育进行比较；TreeMix推断历史基因流；Stairway Plot重建各遗传聚类的种群历史动态；ResistanceGA优化景观阻力面，结合Mantel检验比较IBR、IBD和IBE；冗余分析（RDA）划分环境、空间和祖先相关基因组变异；BayeScan、PCAdapt和潜在因子混合模型（LFMM）扫描候选选择位点。
三、研究结果
3.1 遗传结构分析
PCA和ADMIXTURE分析一致揭示了六个地理结构清晰的遗传聚类（Cluster I–VI），分别对应准噶尔盆地、蒙古高原—阿拉善高原—河西走廊、吐鲁番—哈密盆地、敦煌盆地、塔里木盆地东部和塔里木盆地西部。ADMIXTURE交叉验证确定K=6为最优聚类数。PCA显示Cluster I与II、Cluster V与VI在PC空间中存在部分重叠，表明这些组群虽然可作为独立ADMIXTURE聚类区分，但共享大量等位基因频率变异，与历史基因流或近期分化一致。DAPC较PCA改善了多数聚类间的分离度。
3.2 核基因组与线粒体系统发育模式
ML系统发育树恢复了两个高支持度的区域支系：东北支系（UFBoot=100%）包含准噶尔、阿拉善和吐鲁番—哈密盆地个体（Cluster I–III），西南支系（UFBoot=100%）包含敦煌和塔里木盆地个体（Cluster IV–VI）。然而，ML树中Cluster I、III和IV呈现非单系状态，而NJ树则恢复这些组群为单系，这种不一致可能反映了近期分化种群中的不完全谱系分选（incomplete lineage sorting, ILS）。线粒体系统发育恢复了四个地理结构母系谱系，将核Cluster I+II和IV+V分别合并，而核Cluster III和IV为单系。总体而言，核基因组与线粒体数据集在AECA范围内呈现大致一致的系统地理格局，但核SNP提供了更细尺度的分辨率，包括塔里木盆地东西部种群的分化。
3.3 历史基因流
TreeMix分析识别出两个显著迁移事件。最强迁移边指示从包含Cluster III、IV、V和VI支系的祖先节点向Cluster I的基因流，提示西南和东北区域组合之间的历史连通性。次级较弱迁移边连接同一祖先节点至Cluster II。值得注意的是，TreeMix未检测到Cluster V和VI之间的迁移，与ADMIXTURE结果不同，可能反映TreeMix更关注深层历史事件而非近期混合。
3.4 种群历史动态
Stairway Plot重建结果显示六个遗传聚类呈现异质性种群历史：Cluster I在约5000年前经历种群增长；Cluster II在约50,000年前经历早期扩张，随后在约20,000–15,000年前衰退；Cluster III和IV在约15,000年前扩张，其中Cluster IV在扩张后进入衰退；Cluster V和VI在整个晚更新世保持相对稳定的有效种群大小。这些结果表明种群对晚更新世气候波动的响应具有区域特异性，而非全分布区同步的种群动态变化。
3.5 景观与环境关联分析
ResistanceGA分析中，土地覆盖与海拔组合阻力面获得最强支持（AICc=-3728.95，模型权重w=0.91）。Mantel检验显示三类预测因子均与基因组分化显著相关，但相对强度不同：地理距离相关性最强（DPS: r=0.470；FST/(1?FST): r=0.518），海拔阻力次之（r=0.467），环境距离相关性较弱（r=0.129）。地理距离与海拔阻力高度相关（r=0.921），表明二者存在显著空间非独立性。EEMS分析显示天山和北山等主要山系区域有效迁移降低，而低地荒漠盆地和廊道区域连通性较高。
3.6 环境、空间与祖先相关基因组变异
未约束RDA中，六个生物气候变量解释了较大比例的SNP变异，其中bio9（最干季度平均温度）贡献最大（28.67%）。然而，在控制空间自相关（dbMEMs）后，环境变量解释的变异比例骤降至2.1%；在控制ADMIXTURE祖先系数后进一步降至0.86%。dbMEM变量单独解释了66.9%的SNP变异。偏RDA分析表明，控制种群结构后，bio14（最干月降水量）解释最大比例残余变异（25.78%），而bio9贡献降至10.87%。
3.7 候选选择位点扫描
BayeScan检测到110个候选SNP，PCAdapt检测到216个，LFMM检测到501个。三种方法间重叠有限：PCAdapt与LFMM共享14个候选位点，PCAdapt与BayeScan共享9个，无位点被三种方法同时检出。这种有限重叠反映了不同方法假设和敏感性的差异，以及强种群结构和空间自相关的影响。
四、讨论与结论
研究人员认为，变色麻蜥的基因组分化主要受AECA持续存在的景观格局驱动。核基因组与线粒体数据集在系统地理格局上总体一致，但核SNP提供了更细尺度的种群结构分辨率，反映了不同遗传标记系统在继承模式和溯祖特性上的差异。种群对晚更新世气候波动的响应具有区域异质性，表明局部地理和环境背景塑造了不同的种群历史轨迹。
景观基因组分析的核心结论是：基因组分化主要与景观阻力（特别是海拔和土地覆盖）和地理距离相关，而当代环境变量在控制空间结构后解释的变异相对较少。这支持了持续景观格局、历史连通性和区域特异性种群历史共同塑造基因组分化的观点，而非仅由近期环境梯度驱动。该研究强调了整合系统地理学和景观基因组学方法对理解地形异质荒漠生态系统中种群分化和进化历史的价值，也为干旱生态系统保护规划中识别维持长期连通性和进化多样性的关键区域提供了科学依据。论文发表于《Molecular Phylogenetics and Evolution》。
