《Life》：Population Genomics of Chinese Melia azedarach Germplasm: Latitudinal Genetic Structure, Spatial–Climatic Associations, and Climate-Associated Candidate Genes
编辑推荐：
为了解地理如何塑造苦楝（Melia azedarach）的遗传变异，研究人员对来自中国国家种质资源库（National Germplasm Bank of China）保存的15个种源91个母本半同胞家系的91份材料进行了全基因组测序。遗传多样性从南向北增加（
为了解地理如何塑造苦楝（Melia azedarach）的遗传变异，研究人员对来自中国国家种质资源库（National Germplasm Bank of China）保存的15个种源91个母本半同胞家系的91份材料进行了全基因组测序。遗传多样性从南向北增加（期望杂合度He：0.098至0.333），而海南种群表现出较低的种群内多样性，并通过成对固定指数Fst证实其与大陆种群高度分化，形成一个独立类群。曼特尔检验（Mantel test）确认了显著的隔离距离模式（IBD，r = 0.307，p = 0.001）。冗余分析（RDA）解释了17.23%的个体水平遗传变异，但偏冗余分析（partial RDA）表明，这种关联大部分可归因于种群结构，在控制种群结构后，环境变量解释的独立方差微乎其微。全基因组关联分析（GWAS）进一步鉴定出可能与气候相关性状有关的候选位点，包括一个温度相关基因编码富含亮氨酸的受体样蛋白激酶和一个降水相关基因编码胁迫相关蛋白RAR1。这些结果共同描绘了苦楝的全基因组遗传多样性和空间-气候关联图谱。本研究为评估库藏资源、分子育种和保护这一重要树种提供了遗传基础。
研究背景、现存问题与研究目的
苦楝（Melia azedarach）是一种广泛分布于热带和亚热带地区的落叶乔木，具有重要的经济与生态价值，可用于木材生产、城市绿化和生态修复，其提取物还具备药用和农用价值。然而，此前对该物种遗传多样性的研究多基于形态学标记或低密度的分子标记如简单重复序列（SSR），这些标记基因组覆盖度有限，难以可靠地鉴定受选择位点或区分遗传漂变与适应性选择的作用。为保护与利用森林遗传资源，中国建立了国家种质资源库网络，但库藏材料为异地保存（ex situ），其遗传结构反映的是保存种质而非当代野生种群的特征。同时，由于所有材料均在统一栽培条件下生长，这些异地保存库可作为共同园平台，最大限度减少来自野生栖息地的环境噪音，从而更清晰地研究由遗传继承驱动的进化信号。因此，研究人员利用全基因组重测序（WGS）技术，系统分析了保存于国家种质资源库中的91份苦楝材料的基因组变异，旨在：（1）计算遗传多样性并确定系统发育关系；（2）描述种群结构并检验隔离距离（IBD）和隔离环境（IBE）模型；（3）利用冗余分析（RDA）鉴定与基因组变异最相关的气候变量；（4）通过全基因组关联分析（GWAS）鉴定与气候相关性状相关的候选位点，为库藏资源评估、分子育种及该树种的保护提供遗传基础。
主要关键技术方法
研究人员从位于江西省永修县的国家苦楝种质资源库中采集了91份材料，这些材料来自15个地理种源，每份材料代表一个独立的母本半同胞家系。采用改良CTAB法提取叶片DNA，构建文库后在Illumina NovaSeq 6000平台上进行双端150 bp测序，平均深度约10×。将高质量读段比对至苦楝染色体级参考基因组Melaze_AO_1.1，利用GATK进行变异检测，经严格过滤后获得6695869个高质量双等位单核苷酸多态性（SNP）。主要分析方法包括：基于邻接法（NJ）构建系统发育树；主成分分析（PCA）；ADMIXTURE软件进行祖先聚类；曼特尔检验评估IBD；冗余分析（RDA）和偏冗余分析（partial RDA）评估环境关联；利用rMVP包中的FarmCPU模型进行GWAS，鉴定与年均温（bio1）和年降水量（bio13）显著关联的SNP。
研究结果
3.1 全基因组重测序与群体遗传变异检测
研究人员对91份苦楝材料进行全基因组重测序，平均测序深度为11.67×，平均93.98%的读段成功比对至参考基因组，88.48%的基因组区域至少被覆盖一次。经筛选后获得6695869个高置信度SNP，这些SNP相对均匀地分布在14条假染色体上。功能注释显示，39.26%的SNP位于基因间区，35.24%位于内含子区，外显子区占3.87%。
3.2 群体遗传多样性分析
研究人员发现，遗传多样性在不同种源间差异明显，陕西种群普遍表现出最高值。有效等位基因数（Ne）范围为1.1666（三亚）至1.5825（勉县），期望杂合度（He）范围为0.0974（三亚）至0.3326（安康），最低值均出现在南部和沿海种群。经过稀疏化校正后，南北向递增的多样性格局依然成立，表明该趋势并非由采样不均造成。总体而言，陕西种质（安康和勉县）维持着高遗传多样性，而南部和沿海材料（金华、东方和南通）在库藏中表现出较低的基因组多样性。
3.3 群体遗传结构分析
邻接法系统发育树将91份材料分为两个主要谱系：谱系I包含海南种群和陕西种群（安康和勉县）；谱系II包含江苏、湖北及陕西其余种群（周至和眉县），呈现交织聚类模式。PCA显示前三个主成分解释了30.83%的总遗传变异，PC1分离海南岛与大陆种群，PC2分离沿海江苏与内陆陕西种群。ADMIXTURE交叉验证误差在K=3-4处达到最小平台期，依据简约性原则选取K=3为最优聚类数，其中海南种群保持高度独立且均匀的遗传背景，无近期基因渗入迹象。成对Fst计算进一步证实，海南岛与大陆之间的Fst值（0.1996-0.3670）显著高于大陆内部比较（0.0498-0.1998），表明琼州海峡作为地理屏障限制了基因流。
3.4 隔离距离分析
曼特尔检验结果显示，地理距离与遗传距离之间存在显著正相关（r=0.307，p=0.001），支持经典的IBD模型，表明即使在统一的异地保存条件下，这些材料仍保留了源于原始野生栖息地的空间遗传分化。
3.5 环境关联与冗余分析
个体水平RDA显示，空间-气候矩阵解释了17.23%的遗传变异，第一轴主要由温度相关变量驱动，分离海南与陕西种群；第二轴主要由降水相关变量驱动，分离湿润的湖北、江苏与半干旱的陕西种群。然而，在控制种群结构（前五个SNP主成分）后的偏RDA中，环境变量的独立贡献降至约4%，表明观察到的遗传-环境关联主要源于地理隔离和共同祖先，而非个体水平的局域适应响应。
3.6 气候相关性状的全基因组关联分析
研究人员采用FarmCPU模型进行GWAS，共检测到7个与年均温显著关联的SNP和131个与年降水量显著关联的SNP。通过与另一种模型取交集获得高置信位点，最终得到2个温度相关SNP和1个降水相关SNP。在温度相关SNP CM051395.1:18852443所在基因OWV82_005272内，鉴定出一个编码富含亮氨酸的受体样蛋白激酶家族蛋白的候选基因；在降水相关SNP CM051404.1:2457117所在基因OWV82_019844内，鉴定出一个编码含半胱氨酸和组氨酸丰富结构域的胁迫相关蛋白RAR1的候选基因。
讨论与结论总结
讨论部分指出，大多数SNP位于非编码区，这与其它多年生木本植物一致，暗示这些区域可能参与调控基因表达和表型变异。苦楝种群遗传多样性相对较低，可能与果实依赖鸟类传播但因含有毒三萜类化合物而降低取食效率有关。遗传多样性从南向北递增的趋势与传统纬度梯度规律相反，可能反映了北部边缘地区在冰期后扩张过程中谱系的积累与混合，以及在秦岭地区的接触带效应，但这一推断仍需古分布模拟和野生种群采样验证。三种聚类方法均一致地将海南种群识别为高度分化的独立进化单元，琼州海峡的长期隔离促进了异域分化。RDA与偏RDA的对比表明，地理隔离和共同祖先主导了遗传-环境协变，气候变量的独立贡献有限。GWAS鉴定的候选基因为后续气候适应性研究提供了有价值的目标，但需功能验证确认其实际作用。异地保存的共同园设计排除了环境塑性的混淆影响，使检测到的关联更可能反映历史适应而非当代响应。
结论翻译：本研究解码了91份苦楝材料的群体基因组景观，揭示了清晰的地理结构、北向递增的多样性梯度以及一个被琼州海峡隔离的高度独立的海南种群。这些发现，连同通过全基因组关联分析鉴定的气候相关候选位点，为保护国家种质资源库中的种质资源——特别是高多样性的陕西种源和独特的海南种群——提供了遗传基础。由于所有材料均为异地保存，研究结果主要反映保存种质而非当代野生种群。此外，将北部高多样性区域解释为避难所或接触带的推论仍缺乏古分布建模或野生种群采样的支持。未来研究应结合群体基因组扫描与独立种群的功能验证，以检验可能构成气候适应基础的候选位点。
