野生与逸生芜菁（Brassica rapa）群体局部（适应不良）遗传架构解析

《Molecular Ecology》：Insight Into the Genetic Architecture of Local (Mal)Adaptation in Wild and Feral Brassica rapa Populations

　　

摘要翻译：随着气候变化驱动快速环境变化，作物因其有限的适应性遗传多样性而面临重大挑战，威胁粮食安全。作物野生近缘种以及潜在的逸生群体可作为育种中适应性多样性的宝贵资源库。因此，探索和保护这一多样性至关重要，因为它可能涵盖广泛的适应性。在本研究中，研究人员调查了自然发生于野生状态的芜菁（Brassica rapa，field mustard）群体中气候局部适应的遗传基础，以及未来条件下的适应不良风险。共采样62个群体，显示出逸生群体（法国、斯洛文尼亚）与野生群体（阿尔及利亚、西西里）之间的强烈分化。基因组-环境关联分析（Genome–Environment Association, GEA）识别出与气候变量相关的基因组区域，更广泛的模式表明其具有多基因遗传架构。逸生群体拓宽了所探索的气候范围，但其降低的遗传多样性限制了研究人员捕获野生群体中完整适应性变异谱系的能力，表明气候适应的遗传架构在这两组之间仅部分共享。研究人员还使用基因组偏移（Genomic Offset, GO）分析评估当前与2050年气候预测下的预测最优遗传组成之间的不匹配。一项共同花园实验表明，GO捕获了适合度变异的某些组分，且逸生群体提供了超出野生群体所占据环境空间的额外预测信息。这些发现突显了芜菁气候适应遗传架构的复杂性，以及基因组工具在气候变化下指导保护和作物改良方面的潜力与局限。

论文解读：

气候变化正驱动环境快速变迁，许多物种面临适应或迁移以存续的压力。农作物因驯化瓶颈和强烈人工选择而遗传多样性较低，其适应潜力受限，进而威胁长期粮食安全。作物野生近缘种（Crop Wild Relatives, CWRs）及逸生群体（descendants of cultivated plants that have escaped cultivation）被视为适应性遗传多样性的潜在资源库。局部适应（local adaptation）概念认为，在特定环境中长期稳定存在的群体可演化出最大化其在该环境适合度的遗传适应。逸生群体虽遗传变异通常较低，但若能在野外持续存在，可能为作物改良提供适应性基因资源。芜菁（Brassica rapa）是重要蔬菜作物，也是异源四倍体甘蓝型油菜（Brassica napus napus）的祖先亲本之一，其野生群体分布相对有限。研究人员从阿尔及利亚（28个）、意大利（17个）、法国（16个）和斯洛文尼亚（1个）采集了62个自生群体，其中法国和斯洛文尼亚群体被鉴定为逸生型，其余为野生型。该研究旨在探究气候适应性遗传变异是否在野生或逸生群体中出现，评估纳入逸生群体对分析结果的影响及其作为遗传资源的价值，并利用基因组偏移（Genomic Offset, GO）方法评估未来气候变化下的适应不良风险。

研究人员对62个群体进行了基于混合样本的全基因组重测序（Pool-seq），平均深度50×，经比对和变异检测获得14,457,097个双等位SNP。气候数据取自ERA5-land再分析数据集（2010–2020年），衍生19个生物气候变量。基因组-环境关联分析（GEA）采用BayPass软件的标准协变量模型，以贝叶斯因子（Bayes Factor, BF）衡量关联强度，并运用局部得分（local score）方法界定候选基因组区域。隔离-环境（Isolation by Environment, IBE）检验采用偏Mantel检验，控制地理距离后分析候选位点遗传分化与环境距离的相关性。GO估计采用几何方法，基于10个CMIP5气候模型在RCP4.5情景下2050年（2040–2060平均）的预测，以中位数降低模型间变异。共同花园实验位于法国勒鲁（Le Rheu），种植62个群体，测量种子数与千粒重，以其乘积（生殖生物量）作为适合度代理。

结果部分包括三项主要内容。第一，群体多样性与结构：野生群体拥有约940万个私有SNP，逸生群体约130万个，两组共享约430万个；野生群体的核苷酸多样性（π=0.00916）显著高于逸生群体（π=0.00465）。随机等位基因主成分分析（PCA）显示，第一主成分分离逸生与野生群体，第二主成分进一步区分阿尔及利亚与意大利野生群体；总体遗传分化F_ST=0.328，其中组间分化F_ST=0.24，组内分化F_ST=0.122。第二，气候局部适应的遗传架构：基于全部62个群体（自发群体集）的GEA检测到3805个显著关联SNP（BF≥20），经局部得分界定29个强候选区域，候选基因包括SWEET家族糖转运蛋白（SWEET6、SWEET7、SWEET14）等；仅基于45个野生群体的分析检测到4491个显著SNP和57个候选区域，其中仅8个区域与自发分析重叠，野生群体候选基因包括热激蛋白（HSP17.4A、HSP17.6C）和叶绿体发育基因（PTAC7）。高度显著SNP在野生群体中显著富集错义变异（比值比=3.18，p=1.96×10^-13）。IBE检验显示，自发群体集识别的适应性SNP与环境距离强相关（Mantel R=0.88），但该信号在野生群体内部明显减弱（R=0.40），说明自发群体的标记未能完整捕获野生群体特有的适应性变异；进一步分析表明，中等数量标记（约4000个SNP）对应的相关性最强，且曲率检验支持适应呈多基因架构。基于17个逸生群体的单独GEA分析出现与伪变量相关的假阳性，模拟显示其检测功效不足，因而未采用。第三，基因组偏移与共同花园：GO预测显示斯洛文尼亚和法国布列塔尼群体在2050年下面临最高的潜在适应不良风险，波河及罗讷河谷和西班牙中部风险中等，沿海地区风险较低。在共同花园中，自发群体的GO值与生殖生物量呈显著负相关（R2=0.36，p=3.3×10^-7），且优于单纯环境距离的预测（R2=0.22）；仅分析野生群体时，GO与适合度代理的相关性仍显著但减弱（R2=0.11），环境距离则无预测力（R2=0）。基于野生群体集计算的GO与共同花园适合度无相关（R2=0），这归因于共同花园位置远超出野生群体的环境训练范围。

讨论部分指出，自发群体GEA结果支持气候适应具有高度多基因性，局部得分为法鉴定出少量大效应位点，但候选基因注释依赖于与远缘物种的直系同源关系，需谨慎解读。野生群体表现出更强的错义变异富集和更多私有适应性变异，与野生群体较高的遗传多样性一致；逸生群体与野生群体在气候适应的遗传架构上仅部分共享。逸生群体尽管多样性较低，仍显示出局部适应信号，可能源于长期自然选择或继承自当地适应的地方品种，但其丢失了一部分在野生群体中有利的等位基因，提示驯化和逸生瓶颈可能削弱其适应潜力。关于GO，研究强调其度量的是在新环境中群体与最优本地适应群体之间的相对适合度差异，而非群体适合度随时间的演化；GO预测应谨慎外推，但将逸生群体纳入后显著扩展了可预测的环境范围，为评估气候变化脆弱性提供了增量信息。

研究结论部分翻译如下：本研究提供了对欧洲和北非野生及逸生芜菁群体局部适应基因组基础的洞察。研究人员表明，该物种的适应高度多基因化，涉及分布在整个基因组中的大量小效应变异，以及一些大效应位点。结果证明，逸生群体尽管相对于野生群体分子多样性降低，且具有潜在复杂的种群历史，但仍可表现出清晰的局部适应信号。这与中性遗传多样性并不直接预测适应性遗传多样性的表述一致，并提示逸生谱系可能在野生近缘种缺失的环境中作为适应性等位基因的相关资源库。然而，逸生群体的适应潜力可能因驯化和逸生事件而被侵蚀，增加了未来气候变化下适应不良的风险。野生与逸生群体在气候适应遗传架构上的差异既带来挑战也带来机遇：虽然它降低了检测野生群体适应性变异的能力，但自发分析揭示了逸生群体清晰的适应信号，并显著扩展了可用于适应不良预测的环境范围。尽管预测仍不完善，纳入逸生群体使研究人员能够实现仅关注野生群体无法达成的预测。研究结果强调了在研究局部适应时，尤其是在野生类群分布有限时，将逸生群体与其野生近缘种一并考虑的重要性。