人口历史与分化选择共同驱动濒危树种Firmiana danxiaensis在丹霞—喀斯特生境岛中的谱系特异性适应

《Plant Diversity》：Demography and divergent selection shape lineage-specific adaptation in the endangered tree Firmiana danxiaensis

　　

整合多数据揭示复杂的局部适应（local adaptation）。结合染色体尺度基因组学（chromosome-scale genomics）、重测序（resequencing）及土壤/气候数据，研究人员揭示人口历史（demographic history）与选择如何塑造濒危树种Firmiana danxiaensis的适应；该树种同时分布于丹霞（Danxia）和喀斯特（Karst）地貌。谱系特异性演化轨迹：两个地貌上的四个高度分化谱系表现出不同的人口历史和遗传负荷（genetic load），其局部适应由隔离（isolation）、漂变（drift）、基因流（gene flow）和分化选择（divergent selection）共同驱动，而非单一物种全范围轨迹。地貌特异性选择机制与共享功能主题：丹霞适应为土壤驱动（edaphic-driven），喀斯特适应与气候相关（climate-associated），二者几乎不共享适应变异；但抗旱响应（drought response）和金属离子结合（metal ion binding）是所有谱系中一致的功能主题。

研究背景：局部适应由环境选择塑造，但在破碎化生境中，人口历史可通过降低遗传多样性、增强漂变与近交、改变有害变异暴露和清除来影响选择效率；濒危土壤特化种尤其需要同时评估环境关联位点与谱系特异性人口历史、近交、遗传负荷和清除。南岭丹霞与喀斯特土壤岛形成对比选择压力，而Firmiana danxiaensis同时分布于两类地貌，既往研究缺乏整合高通量基因组、高分辨率环境和人口历史的框架。Yali Li等构建染色体尺度参考基因组，对20个野生种群225个体（15个丹霞地貌DL种群、5个喀斯特地貌KL种群，覆盖6个分布点）及11个Firmiana simplex外群重测序，测定20个种群20项土壤变量并整合19项气候变量；开展群体结构、PSMC/Fastsimcoal2人口历史、TreeMix/Dsuite基因流、ROH/遗传负荷/R_X/Y、RDA/LFMM基因组—环境关联（GEA）及XP-EHH/DCMS选择扫描。样本包括华南植物园栽培单株及野外种群。
结果：3.1 染色体尺度参考基因组：获得1.51 Gb组装，contig N50为69.08 Mb，97.37%序列锚定于20条假染色体，BUSCO完整度98.8%，预测31,938个蛋白编码基因。3.2 与环境适应可能相关的基因组特征：比较基因组显示与Heritiera littoralis共享约133和20 Ma的WGD（全基因组复制，whole-genome duplication），分化约20.5 Ma；鉴定682个扩张基因家族、2,238个转录因子（TF，transcription factor），LTR-RT（长末端重复反转录转座子）占67.12%，近期插入峰约1.3–1.4 Ma，bHLH、MYB、C2H2、NAC、bZIP、WRKY等胁迫相关家族扩张。3.3 群体结构揭示四个高度分化谱系：ADMIXTURE、系统树和PCA支持BZY、DX、SXNX、BWYD四个谱系，与地理和地貌相关；保留5,471,157个SNP（单核苷酸多态性，single nucleotide polymorphism）。3.4 分化人口历史与非对称基因流：PSMC显示BZY长期低N_e（有效种群大小，effective population size），SXNX/DX约0.9 Ma达峰后收缩；TreeMix/Dsuite支持非对称基因流，SXNX基因流最高、BZY最低；谱系间F_ST为0.41–0.75，π在BZY最低（4.94×10^-5）。3.5 谱系特异性近交、遗传负荷与R_X/Y清除：BZY的F_IS最高（0.96）、F_ROH最高（0.50）；BZY潜在LoF负荷最高（0.04440），实现负荷谱系间差异小；R_X/Y提示LoF分化大于missense，BZY相对BWYD/DX的LoF R_X/Y<1，但相对SXNX>1。3.6 地貌环境异质性：20个土壤变量中14个显著差异（13个KL更高），19个气候变量中14个显著差异（10个KL更高）。3.7 地貌特异性适应变异：RDA与LFMM交集获DL 1,031个核心变异/685基因、KL 750个/531基因，二者无共享（Fisher检验p=0.87）；DL土壤变量解释74.91%，关键为ACd、TK、AZn，KL气候变量解释67.10%，关键为Bio9、Bio11；DL有130个变异（12.61%）重叠XP-EHH，KL仅6个（0.8%），DL候选包括NAC090、MAPKKK20、BHLH42、ITPK2、GES、GLR1.4。3.8 谱系特异性适应架构：BZY仅与总铜关联，DX/SXNX多因素关联分别占96.21%/98.37%；各谱系均以抗旱响应最富集，金属离子结合次之；约10%核心变异位于编码区，多数位于非编码区，上游最丰富；BZY、DX、SXNX分别有1,977、132、91个适应变异具正选择信号（31.09%、7.98%、2.77%）。
讨论总结：LTR-RT、WGD和胁迫相关TF扩张可能提供适应底物，但原文强调需功能验证；BZY人口侵蚀最强、潜在负荷最高，R_X/Y提示清除不足以完全抵消有害等位基因积累；南部谱系多样性较高、近交较低且历史连通性较强；DL适应主要由土壤驱动，KL更受气候驱动，谱系与地貌分布混杂及分化时间未知限制候选变异解释；干旱响应和金属离子结合为跨地貌/谱系反复主题，未检测TPC1，原文认为更符合多基因架构，非编码变异可能具调控作用但未证明；保护上BZY应优先，辅助基因流需风险评估，SXNX可作同DL候选供体，DL重土壤稳定与污染控制，KL重气候韧性监测。研究结论（译）：通过整合染色体尺度基因组、全基因组重测序、野外土壤、气候和GEA，研究表明人口历史与异质选择塑造F. danxiaensis适应演化；四个高度分化谱系具有不同人口历史、多样性、近交、遗传负荷和适应特征，局部适应非单一物种全范围轨迹，而由历史隔离、漂变、基因流和分化环境选择共同作用；DL与KL适应变异几乎不重叠，DL以土壤驱动、KL以气候关联为主，但抗旱响应和金属离子结合跨谱系一致。