蛇纹岩栖息地独立殖民化塑造谱系分化的异质基因组景观

《Journal of Systematics and Evolution》：Independent colonizations of serpentine habitats shape heterogeneous genomic landscapes of lineage diversification

　　

蛇纹岩土壤以高含量重金属、低养分可利用性和水资源匮乏为特征，对植物物种构成显著的生态挑战。与蛇纹岩露头相关的土壤变异在短地理距离内产生剧烈的环境对比，为研究异质性地土壤景观中的谱系分化和种群结构提供了机会。在此，研究人员检验了三种常同域分布的植物物种——Lavandula stoechas L.、Phlomis purpurea L.（均为唇形科Lamiaceae）以及Halimium atriplicifolium (Lam) Spach.（半日花科Cistaceae）——的进化历史，这些物种栖息于伊比利亚半岛南部的相邻蛇纹岩与非蛇纹岩生境中。利用比较种群基因组学框架，研究人员检验了蛇纹岩种群是代表独立的进化谱系还是对这一土壤生境的平行殖民，以及基因组分化在何种程度上同时反映土壤条件和地理隔离。高度支持的最大似然种群系统发育表明，三个物种均对蛇纹岩栖息地发生了重复的、平行的殖民事件。隔离距离（isolation-by-distance, IBD）的程度以及种群中祖先群体成员身份的地理分布在这三个同域物种间存在显著差异。人口统计模拟揭示了低至中等程度的近期不对称基因流，其特征在于从蛇纹岩斑块向邻近非蛇纹岩种群的定向遗传贡献，而基因型-环境关联分析确定了土壤选择下的潜在候选位点。不同土壤上种群间株高的显著变异代表了与环境相关的表型分化，并提示生态型差异，而该关联的强度在不同物种间有所变化。综合而言，研究结果表明，平行辐射至蛇纹岩栖息地、种群祖先、隔离距离以及持续的基因流对这三个物种的种群基因组结构产生了不同程度的贡献。

**论文解读：蛇纹岩栖息地独立殖民化与谱系分化的异质基因组景观**

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

长期困扰进化生物学的一个核心问题是：相似的选择环境是否会导致独立谱系产生可预测且平行的进化响应。当同域分布的物种在地理空间上遭遇相似的环境挑战，尤其是在新栖息地殖民的过程中，进化响应的可预测性仍是一个悬而未决的问题。此类动态过程通常涉及对新环境的殖民以及在严苛生态条件下的持续生存，而相应的基因组、地理和人口统计模式能够为物种如何响应过去的环境挑战提供重要线索。

蛇纹岩土壤以低养分可利用性、高重金属含量和水资源匮乏为特征，对植物构成极端生理胁迫。尽管蛇纹岩生境与非蛇纹岩生境往往近在咫尺，蛇纹岩耐性种群通常维持独特的基因组身份，并由开花时间差异和生殖不兼容等基因流障碍所强化。蛇纹岩露头的碎片化、隔离性和重复性分布，为判断蛇纹岩土壤上的种群究竟源自单一进化事件还是由景观中多个谱系重复独立殖民提供了天然实验平台。

已有研究大多来自北美和中欧两个地区，且各自仅研究了焦点群落中的一个物种，地理范围有限且缺乏重复，制约了对同域物种间平行进化过程普遍性的认识。本研究在伊比利亚半岛南部选取三种常同域分布的植物物种——Halimium atriplicifolium (Lam) Spach.、Lavandula stoechas L.和Phlomis purpurea L.——开展重复检验，旨在区分蛇纹岩殖民的两种竞争模型：单一栖息地转变与平行进化。

**研究内容与核心发现**

研究人员在西班牙马拉加省采集了20个种群（构成10对，每对包含一个蛇纹岩种群和一个非蛇纹岩种群），每种种群各取15个个体，共885个个体。通过双酶切限制性位点相关DNA测序（double digest restriction-site associated DNA sequencing, ddRAD）获得全基因组单核苷酸多态性（single nucleotide polymorphism, SNP）数据，结合种群结构推断、系统发育重建、人口统计历史模拟和基因型-环境关联分析，系统解析了三个物种在蛇纹岩与非蛇纹岩生境间的进化历史。

研究结果表明，三个物种均发生了多次独立的蛇纹岩栖息地殖民事件，而非单一蛇纹岩关联谱系在景观中的扩张。隔离距离信号在三个物种中呈现不同程度：Halimium atriplicifolium表现出强烈的空间遗传结构，而Lavandula stoechas和Phlomis purpurea则呈现更均匀的基因组聚类模式。人口统计模拟显示，所有分析的生态对均支持含有不对称迁移的模型优于严格隔离模型，基因流方向为从蛇纹岩斑块向非蛇纹岩种群扩散，而反向迁移极为有限。基因型-环境关联分析分别识别出21、7和45个高置信候选SNP位点，并通过BLAST比对发现了涉及细胞壁成分、蛋白激酶/磷酸酶、转录因子、激素响应基因和初级代谢酶等多种功能类别的候选基因。

**主要关键技术方法**

研究使用了以下核心分析方法：ddRAD测序技术用于全基因组SNP发掘；ADMIXTURE v1.3.0进行祖先群体聚类分析；TreeMix v1.13推断种群历史分裂和基因流事件；fastsimcoal2 v2.8基于二维联合位点频谱（2D-JSFS）进行人口统计模拟；BayPass v3.1与偏冗余分析（partial Redundancy Analysis, pRDA）交叉验证候选适应位点；线性混合效应模型检验功能性状差异。样本涵盖20个种群，来自西班牙马拉加省，蛇纹岩土壤性质通过土壤化学分析（Eurofins Spain）确认。

**研究结果**

**3.1 功能性状**：线性混合效应模型检测到Lavandula stoechas和Phlomis purpurea蛇纹岩种群的株高分别平均矮13.56 cm ± 5.65 s.e.和16.35 cm ± 4.84 s.e.，而Halimium atriplicifolium蛇纹岩种群株高虽有降低趋势（6.86 cm ± 8.99 s.e.）但未达显著水平。三个物种的比叶面积（specific leaf area, SLA）在蛇纹岩与非蛇纹岩种群间均无显著差异。

**3.2 基因组结构与聚类**：主成分分析（PCA）和ADMIXTURE分析显示，H. atriplicifolium具有强地理结构，种群间几乎无混合；L. stoechas和P. purpurea基因组结构较松散，中央聚类涵盖多数种群，少数外围种群分化明显。隔离距离检验显示H. atriplicifolium的Fst与地理距离高度相关（r = 0.763, P = 0.001），L. stoechas（r = 0.404, P = 0.001）和P. purpurea（r = 0.347, P = 0.005）呈中度相关。L. stoechas和P. purpurea的蛇纹岩种群分别表现出更高比例的祖先群1和祖先群2成分。

**3.3 系统发育关系与基因流**：TreeMix分析为每个物种构建了分辨良好的系统发育树。H. atriplicifolium推断出2个历史基因流事件，均为非蛇纹岩种群间的长距离连接；L. stoechas识别出3个基因流事件，蛇纹岩种群散布于树中但呈现一定聚类；P. purpurea推断出5个基因流事件，其中仅有一例是蛇纹岩与非蛇纹岩邻近种群间的基因流（PAV75-PAV76）。fastsimcoal2模型选择在所有生态对中均支持含有迁移的模型。所有配对分析中，迁移率均极低，2Nm值均低于1，且严格定向为从蛇纹岩到非蛇纹岩种群（2Nm SN），而反向迁移（2Nm NS）可忽略不计。最高2Nm均值出现在地理距离最近的种群对中。

**3.4 基因型-环境关联**：BayPass和pRDA交叉验证分别鉴定出H. atriplicifolium 21个、L. stoechas 7个、P. purpurea 45个高置信候选SNP。一致性分析验证了这些等位基因在10个生态重复中的平行频率变化方向。BLAST注释识别出一系列功能基因，包括果胶乙酰酯酶、五肽重复（pentatricopeptide repeat, PPR）蛋白、蛋白磷酸酶2C 52、糖原磷酸化酶1类似蛋白、CBL互作丝氨酸/苏氨酸蛋白激酶5类似蛋白、生长素诱导蛋白X10A类似蛋白等。

**讨论总结与研究结论**

研究结果表明，三个物种均经历了多次独立蛇纹岩栖息地殖民，符合研究人员的预期。不同物种的基因组聚类模式差异显著，反映了物种特异性生命史特征、扩散机制和种群历史的影响。H. atriplicifolium的地理隔离格局与其种群稀少、分布分散和花粉限制一致；L. stoechas和P. purpurea的均质化格局则与其连续分布和移动昆虫授粉的生物学特性相符。不对称基因流模式表明，蛇纹岩生境对不适应的迁移者构成强烈的选择屏障，但蛇纹岩种群的遗传变异性仍可向外扩散，为邻近非蛇纹岩种群提供潜在适应性遗传资源。候选基因涉及ABA信号通路调控、渗透胁迫耐受、根系结构重塑等多个生理层面，反映了蛇纹岩适应需要同时应对干旱、金属毒性和养分缺乏等多重胁迫的复杂性。

在进化与保护意义上，蛇纹岩种群作为定向进化的遗传来源，可能对非蛇纹岩种群在气候变化下的适应性潜力具有重要价值。维持蛇纹岩与非蛇纹岩种群的连通性对保护完整的环境梯度基因组多样性至关重要。研究同时指出，ddRAD测序仅覆盖基因组的一小部分，可能遗漏其他适应位点；未来需结合交互移植和生态位建模等整合方法，进一步阐明环境、人口统计和选择在塑造异质景观进化响应中的作用。