地理隔离与生态适应协同驱动东喜马拉雅—横断山区Juniperus coxii的谱系分化

《Plant Diversity》：Geographic isolation and ecological adaptation synergistically drive lineage divergence within Juniperus coxii across the Eastern Himalaya–Hengduan Mountains

　　

摘要：阐明地理隔离与生态适应如何相互作用驱动谱系分化，仍是进化生物学中的核心挑战。山地系统具有特殊的研究价值，因为陡峭的地形促进异域分布，而急剧的环境梯度则有利于局部适应。研究人员对来自东喜马拉雅—横断山脉（EHHM）19个居群的157个Juniperus coxii个体进行了靶向富集测序，解析出三个具有地理结构的基因组谱系：西部（W）、东北部（EN）和东南部（ES）。抗性表面建模和基因流分析表明，地形屏障在很大程度上限制了W与EN之间的连通性，而ES则保留了部分连通性，并通过生态位分化发生分化。群体历史推断证明了第四纪气候波动期间的持续存在，与其在局部微避难所中的长期存活相一致。基因型—环境关联分析鉴定了多个与气候因子相关的位点，表明各谱系间存在平行生态适应。基因组脆弱性评估揭示了气候变化情景下谱系特异性的适应能力差异。综上所述，这些发现表明地理隔离和生态选择共同驱动了J. coxii的种内分化和适应潜力。该研究建立了阐明山地木本植物谱系多样化的综合框架，并为环境异质性景观中的生物多样性保护提供了指导。

**东喜马拉雅—横断山区Juniperus coxii谱系分化的综合基因组学研究解读**

**一、研究背景与科学问题**

谱系分化是物种形成连续体中的早期阶段，理解其驱动机制需整合地理与生态两个经典维度。传统异域分化观点强调地理屏障对基因交流的限制作用，然而适应性生态分化在此过程中的具体贡献程度仍不明确。经典生态物种形成理论认为，分歧性环境选择通过多效效应或遗传搭车促进生殖隔离；而突变序贯与共享最适模型则提出，姊妹谱系可能在宏观选择压力大致相当的条件下起源，生殖屏障通过内在遗传不兼容等替代机制累积，显著的生态分化则在范围扩张和次生接触之后才出现。判别这些假说需要将地理与生态视为在时空尺度上相互依赖的过程，而非相互独立的驱动因子。

山地系统为检验上述相互作用提供了理想天然场所。复杂地形与深切河网促进居群隔离，而陡峭的气候与土壤梯度则创造强烈的局部适应机会。东喜马拉雅—横断山脉（EHHM）是这种环境复杂性的典型代表，构造抬升与第四纪气候波动共同塑造了隔离屏障与微生境的镶嵌格局，既细分居群又使其暴露于异质性选择之下。尽管关于山地物种多样化的研究日益增多，但许多研究仍将地理隔离与生态分化分开处理，缺乏整合群体历史与气候相关等位变异的基因组框架。Juniperus coxii A.B. Jacks.是一种主要分布于EHHM及邻近区域（约1800–3800米）的山地针叶树种，其居群常分布于孤立山脊或被河谷分隔的坡面，空间构型利于异域分化，而显著的环境梯度可能促进谱系特异适应。然而，刺柏属庞大的基因组和高度重复序列增加了全基因组研究难度，靶向富集测序为该物种的群体基因组学分析提供了高效手段。

本研究利用19个居群157个个体的靶向捕获数据，旨在回答以下三个核心问题：（1）J. coxii在EHHM范围内存在多少个基因组谱系，其在异质景观中的连通性与隔离程度如何；（2）生态分化在多大程度上与谱系边界一致，并在地形屏障影响之外解释遗传分化；（3）各谱系的群体历史与气候相关基因组特征如何，这些模式在气候变化情景下如何转化为差异化的适应能力与脆弱性。

**二、主要技术方法**

研究人员从东喜马拉雅（43个个体）和横断山脉（114个个体）采集了19个居群的J. coxii新鲜叶片样本，并以两个J. formosana个体作为外类群。利用基于柏木属（Cupressus chengiana、C. gigantea）和刺柏属（Juniperus squamata）编码区及基因间区设计的靶向探针进行序列捕获，在Illumina HiSeq X Ten平台完成测序。数据处理方面，以C. gigantea基因间区与J. saltuaria转录组编码序列组成的嵌合参考序列进行比对和变异检测，经严格质控获得高质量SNP数据集。分析方法包括：ADMIXTURE、主成分分析（PCA）和判别主成分分析（DAPC）解析遗传结构；最大似然法构建系统发育树；Stairway Plot2和fastsimcoal2推断群体历史和基因流；估计有效迁移表面（EEMS）分析空间遗传格局；MaxEnt模型和生态位重叠分析评估生态分化；Mantel检验和偏Mantel检验分别检测距离隔离（IBD）和环境隔离（IBE）；基于F_ST和核苷酸多样性的选择性清除扫描；潜在因子混合模型（LFMM）和冗余分析（RDA）进行基因型—环境关联（GEA）分析；偏冗余分析（pRDA）定量分解气候、遗传结构和地理因素的相对贡献；基于RDA的基因组偏移分析评估未来气候情景下的基因组脆弱性。

**三、研究结果**

**1. 群体基因组与系统发育分析**

ADMIXTURE分析在K=3时获得最低交叉验证误差，表明J. coxii存在三个明显遗传谱系：西部（W）、东北部（EN）和东南部（ES）。EN与ES分布于横断山南段并呈南北排列，W则局限于东喜马拉雅的雅鲁藏布大拐弯地区。最大似然系统发育分析证实了三分组结构，其中EN与W的亲缘关系较其与ES更近。PCA和DAPC进一步支持三谱系划分。群体间呈现中度遗传分化，W与EN间最大（F_ST=0.06465），而W与ES（F_ST=0.0467）和EN与ES（F_ST=0.0457）的值较低且相当。各谱系核苷酸多样性接近（W=0.00325，EN=0.00328，ES=0.00320），Tajima's D均为负值，观测杂合度高于期望杂合度，近交系数为负，表明不存在近交现象。

**2. 群体历史动态与基因流**

Stairway Plot2重建表明，三个谱系在白玉冰期（70–30千年）均经历了明显瓶颈效应，其中EN和ES的瓶颈持续时间长于W，之后所有谱系均表现出快速的有效群体大小扩张。fastsimcoal2的两阶段推断首先确定谱系分化顺序为EN最先分化，随后W与ES分离。最优模型显示EN的分化时间约为352万年，W–ES的分化约为256万年。基因流推断揭示显著的谱系间不对称性：ES到EN的迁移率（M_ES-EN=7.90×10^-5）远高于反向迁移（M_EN-ES=4.10×10^-7）；ES到W的迁移率（M_ES-W=3.26×10^-5）超过反向的50倍；EN与W之间的迁移率较低但同样不对称。ES谱系维持最大群体大小（Ne=927,511），EN最小（Ne=347,439）。EEMS分析显示一条强迁移屏障将EN与W和ES分隔，而ES与W尽管地理距离较大却表现出相对较高的迁移连通性。

**3. 物种分布模型与生态分化**

MaxEnt模型对各谱系当前生态位模拟的AUC值均超过0.95。EN和W谱系均适应于东喜马拉雅气候条件，而EN和ES的适宜栖息地位于横断山南段。生态位重叠分析显示EN与W相似性最高（D=0.3219，I=0.4017），EN与ES次之（D=0.1340，I=0.1818），ES与W最低（D=0.0215，I=0.0276）。主坐标分析（PCoA）印证了这些模式，表明ES谱系表现出明显的生态分化。

**4. 正选择基因分析**

结合F_ST和π比率扫描，共鉴定出1954个候选正选择基因。基因本体（GO）富集分析显示，EN谱系正选择基因主要涉及渗透胁迫响应、质体组织和金属肽酶活性；ES谱系基因与染色体组织、蛋白质泛素化、生长、植物型细胞壁生物合成和盐胁迫细胞响应相关；W谱系基因则与脱落酸响应、配子体和花器官发育、质体组织及丝裂原活化蛋白（MAP）激酶活性有关。京都基因与基因组百科全书（KEGG）富集进一步揭示了谱系特异通路：EN富集于ABC转运蛋白、类胡萝卜素生物合成和昼夜节律通路；ES富集于半胱氨酸和甲硫氨酸代谢、剪接体等通路；W富集于N-聚糖生物合成、卟啉代谢等通路，反映出对不同局部环境的适应。

**5. 环境适应相关基因组变异**

LFMM检测到1883个位于编码区的适应性SNP，分布于1312个基因中；RDA检测到305个与185个基因相关的离群SNP，两种方法共有29个候选基因。对1468个候选基因的功能注释显示显著富集于渗透胁迫响应、分子转导活性、细胞骨架组织、质体发育和植物器官形态发生等过程。KEGG富集涉及角质、木栓质和蜡生物合成，植物—病原互作，植物昼夜节律及苯丙烷生物合成等通路。IBD和IBE模式均显著，且控制地理距离后遗传分化仍与环境距离显著相关。偏冗余分析显示，气候因素解释了最大比例的遗传变异（43.3%，占可解释变异的45.2%），其次为遗传结构（14.7%）和地理因素（9.3%）。

**6. 未来气候情景下的基因组偏移预测**

基于RDA的基因组偏移评估在两种排放情景（SSP126和SSP585）和两个未来时段（2041–2070和2071–2100）下进行。SSP126情景下所有居群基因组偏移均低于1.0；SSP585情景下，2041–2070年间偏移低于0.8，但2071–2100年间多个区域超过1.0。所有情景下，东喜马拉雅地区的偏移值均高于横断山区，尤其以2071–2100年为著。相比之下，横断山区的EN谱系在所有情景和时段均表现出较低的基因组偏移，暗示其适应能力相对较强。

**四、讨论与结论**

本研究表明地理隔离与生态适应共同塑造了J. coxii的谱系分化。EN与W的分化主要归因于地理隔离，而ES与（EN+W）之间的分化更多与生态适应相关。从时间上看，早期山地抬升和冰川波动启动了谱系间的地理隔离，随后对异质山地环境的局部适应进一步强化了分化。EN–W边界与超过4000米的高大山系（伯舒拉岭和他念他翁山）吻合，有效阻断了花粉和鸟类传播；而ES–W边界沿海拔较低的拿加丘陵（约2000米），可能是历史花粉传播的通道。IBE的显著性和pRDA中气候因素的主导地位证实了生态隔离在塑造基因组格局中的重要性。谱系间的正选择信号和气候关联位点进一步从功能层面支持了适应性分化。基因组脆弱性评估显示W谱系在SSP585情景下表现出较高的基因组偏移和较低的适应能力，加之针叶树种世代周期长的特性，其面临较大的气候变化失配风险。

研究建议采取辅助基因流策略，将EN谱系种子引入W栖息地以增加遗传多样性和适应性等位基因，同时在雅鲁藏布大峡谷优先实施就地保护和迁地保护措施。本研究建立了将群体历史、局部适应和基因组脆弱性相整合的分析框架，不仅为保护这一分布狭窄的针叶树种提供了科学基础，也为快速环境变化下山地生态系统中特有木本植物的保护提供了可推广的研究范式。