《Ecology and Evolution》：Do Abies religiosa Seedlings That Survived Hot/Dry Range-Edge Conditions Exhibit Genetic Evidence of Selection?
编辑推荐：
圣冷杉（Abies religiosa）是墨西哥越冬君主斑蝶（monarch butterfly）种群的唯一夜宿树，但其地理分布范围预计会因气候变化而缩小。2019–2022年开展的同质园实验（common garden experiment）显示，在3400
圣冷杉（Abies religiosa）是墨西哥越冬君主斑蝶（monarch butterfly）种群的唯一夜宿树，但其地理分布范围预计会因气候变化而缩小。2019–2022年开展的同质园实验（common garden experiment）显示，在3400和3000 m处幼苗存活率高，但在2600 m处存活率仅为6%（代表未来可能更暖、更干的条件）。研究人员使用PCAdapt检验这些幸存者是否与起始种群（由高存活同质园中的幼苗代表）存在遗传差异。研究人员还利用源种群的气候信息开展基因型-环境关联（Genotype-Environment Association, GEA）检验，以识别与气候梯度相关的单核苷酸多态性（single nucleotide polymorphism, SNP）变异。PCAdapt识别出2202个异常SNP（outlier SNPs），其等位基因频率可能受到选择影响。LFMM2识别出33个与5个源种群气候梯度相关的SNP，尤其是冬季降水。两类SNP的次等位基因频率（Minor allele frequencies, MAF）在幸存者中均表现出更大波动。这表明局地气候条件可能在短期和长期时间尺度上对Abies religiosa的遗传差异产生选择作用，但要充分理解这些响应，需要更完整的基因组或转录组数据，以及进一步野外试验验证。
论文解读
圣冷杉（Abies religiosa）分布于Trans-Mexican Volcanic Belt，主要在海拔2800–3500 m，少数可低至2400 m。它是墨西哥越冬君主斑蝶（Danaus plexippus L.）的重要夜宿树，其密集树冠可为斑蝶提供适度温度。气候变化导致热干旱和树皮甲虫爆发，增加圣冷杉死亡；预测到2090年全国适宜面积可能缩小96%，Monarch Butterfly Biosphere Reserve（MBBR）是分布区严重收缩的区域之一。树木自然迁移受世代长、种子扩散有限和生境破碎限制，可能难以跟上气候变化，因此辅助迁移被作为适应性管理策略考虑。Cruzado-Vargas等和Sáenz-Romero曾将11个海拔3000–3550 m种源幼苗植入三个同质园，发现上移400 m可改善幼苗表现。2019–2022年同质园中，3400和3000 m园存活率高，2600 m园存活率仅0–13.3%（总6%）。源海拔仅弱解释存活，研究人员因此假设幸存幼苗可能具有与耐热或耐旱相关的遗传变异，这些变异可跨海拔存在，但可能在低海拔更常见。为此，研究人员比较低海拔园幸存者与中高海拔园初始群体，利用PCAdapt和GEA检测潜在选择信号。结果表明，低海拔园幸存者存在细微选择特征；PCAdapt识别2202个异常SNP，LFMM2识别33个气候相关SNP，尤其与冬季降水相关；幸存者中两类SNP的MAF波动更大。该研究提示局地气候可在短期和长期尺度上选择圣冷杉遗传差异，现存遗传变异可能通过辅助迁移和自然选择帮助该物种应对气候变化，但仍需更完整基因组或转录组资源及野外验证。
关键技术方法概括（不超过250字）：研究人员在MBBR利用11个种源（海拔3003–3552 m）的Abies religiosa幼苗，植入2600、3000和3400 m三个同质园（2019–2022）；低海拔园幸存19株作为幸存组，中高海拔园216株作为初始组。针叶样品经基因分型测序（Genotyping-By-Sequencing, GBS）建库并在NovaSeqX测序；用STACKS-2.59进行de novo组装和SNP calling，VCFtools质控；开展主成分分析（principal component analysis, PCA）和F
ST分析；用PCAdapt（K=3）检测异常SNP；以源种群1961–1990年气候样条数据筛选5个气候变量，用潜因子混合模型2（latent factor mixed model 2, LFMM2；K=6和11）开展GEA，并以基因组膨胀因子（Genomic Inflation Factor, GIF）校正；比较SNP类型和等位基因频率。
研究结果
3.1 迁移距离分析：研究人员用二项分布广义线性模型（generalized linear model, GLM）检验迁移距离与存活关系。低海拔园中迁移距离与存活呈正相关但不显著（p=0.180）；3000 m园（p=0.537）和3400 m园（p=0.311）均无显著关系。向低海拔迁移显著降低存活，但在同一园内迁移距离不能显著预测存活，说明种源表现差异涉及更复杂的气候、环境或遗传因素，故进一步开展基因组分析。
3.2 SNP发现：GBS数据从236株幼苗中识别1,511,227个变异位点；初始组346,114个和幸存组38,509个SNP显著偏离Hardy–Weinberg平衡。经VCFtools质量和深度过滤后去除98.76%位点，保留18,816个高质量SNP；其中转换11,780个（62.61%），颠换7,036个（37.39%）；C/T为5937个，A/G为5843个，A/C为1901个，G/T为1867个，C/G为1736个，A/T为1532个。
3.3 种群结构与遗传分化：基于18,816个SNP的PCA显示，前三个主成分分别解释3.02%、1.92%和1.85%的遗传变异。初始组与幸存组间F
ST为0.00256；部分种源分化较强，但初始组与幸存组在全SNP数据中整体分化不强，符合幸存者为初始群体子集的预期。
3.4 PCAdapt的SNP异常检测：PCAdapt（K=3）识别2202个异常SNP，其中转换1374个（62.40%），颠换828个（37.60%）；A/G为663个，C/T为711个，A/T为144个，A/C为219个，C/G为225个，G/T为240个，比例与全部过滤SNP相似。
3.5 与环境变量相关的SNP：LEA admixture分析显示K=6和K=11分别识别29个和31个位点，合并后33个SNP与气候变量相关；其中转换23个（69.70%），颠换10个（30.30%）；A/G占42.42%（其中11/14为G到A），C/T占27.27%（其中6/9为C到T）。33个SNP中，1个与全部5个气候变量相关，6个与多个变量相关，26个仅与1个变量相关；冬季降水相关SNP最多（15个，14个唯一），其次为dd0（degree-days < 0°C；11个，7个唯一）、Sday（Julian date of the last spring freeze；11个，4个唯一）、Smrp（summer precipitation；9个，1个唯一）和Smrpb（summer precipitation balance；2个，0个唯一）。
3.6 多种方法识别的幸存者与初始幼苗群体遗传差异：PCAdapt识别的2202个SNP与GEA识别的33个SNP有4个重叠，其中3个与冬季降水相关，1个与dd0相关。PCAdapt SNP的MAF均值±标准差在初始组、幸存组和全部幼苗中分别为0.124±0.089、0.127±0.104和0.118±0.087；LFMM SNP分别为0.137±0.100、0.133±0.139和0.127±0.099，幸存组MAF变异更大。F检验未发现两组方差有统计学显著差异；Cohen检验显示少数位点（locus20、locus33、locus21）效应量为0.4–0.46，接近中等效应，提示多位点累计分化可能具有自适应选种潜力，但单个位点差异细微。
讨论部分总结：种群遗传结构方面，花粉广泛传播可模糊邻域遗传差异，但种源聚类仍可见，说明基因流未完全均质化；PCAdapt和GEA可在低遗传分化群体中识别潜在选择位点。气候适应方面，栖息地适宜性受最暖月均温和夏降水与冬季温度交互影响；冬季降水相关SNP比例较高，可能反映真实关系，因为该地区冬季总体降水较少。MAF变异方面，LFMM显著SNP少于PCAdapt，源于方法目标不同；幸存者LFMM SNP的MAF标准差较高可能受采样影响，但等位基因频率变化明显。SNP类型方面，转换/颠换比在总SNP、PCAdapt和LFMM中分别为1.67、1.66和2.3，符合转换更常见、颠换更易受选择的一般格局；LFMM中A/G尤其G到A偏高。局限包括GBS和de novo SNP calling错误率及缺失率高、针叶树连锁不平衡低可能漏检、采样地点未覆盖全分布、幸存者数量少、家族效应或亲缘关系可能影响结果、LFMM主要检测线性关联。未来研究建议开发更多Abies基因组或转录组资源，采用megagametophyte测序、与Abies alba比对、转录组序列捕获，并在有注释基因组时开展GO和KEGG通路分析；还需种植实验验证SNP对田间表现的预测价值。
研究结论翻译：研究人员分析揭示了暴露于强热和干旱胁迫的Abies religiosa幼苗中细微的选择特征。研究人员还识别出33个与源种群气候变化相关的SNP，尤其是冬季降水，并确认这些SNP及PCAdapt识别SNP的次等位基因频率在低海拔幸存者中更易波动，提示相对于初始种群，不同等位基因受到偏爱。与既往研究一致，在分布区内种植并向上迁移的幼苗表现略好，而低海拔种源种植在分布区下缘以下时至少略优于高海拔种源。总体而言，这表明A. religiosa种群内存在现存遗传变异，通过辅助迁移和自然选择的结合，可能有助于该物种在气候变化下存活。然而，需要进一步工作——尤其是基因组或转录组资源的开发——以更好理解为何某些遗传变异与长期气候变化或热干条件下近期选择相关。
