《The Plant Journal》：Natural selection and recombination interact to structure genome-wide variation in pines
编辑推荐：
为解析进化与基因组特征如何塑造裸子植物（尤其松属）的“基因组分化景观”，研究人员对3种松树（P. banksiana、P. contorta、P. nigra）100个个体进行全基因组测序分析，揭示了重组率与基因长度的负相关关系（尤其内含子驱动），并明确了低分化“基因组谷”主要通过平衡、背景或平行选择维持，而“基因组岛”主要源于异域选择。该研究为理解长寿、非开花植物谱系中基因组分化与进化的相互作用提供了新视角。
生命的蓝图——基因组，并非一张平整无瑕的图纸，而是一片被自然之力与历史尘埃雕琢得沟壑纵横的“景观”。在这片景观中，不同种群之间遗传分化的程度在基因组各处差异巨大，那些分化异常剧烈的区域被称为“基因组岛”（Genomic Islands），而分化程度异常平缓的区域则被称为“基因组谷”（Genomic Valleys）。长期以来，科学家们试图理解塑造这片复杂景观的力量，特别是自然选择与遗传物质交换（重组）的相互作用。在开花植物中，相关研究已取得丰硕成果，但作为地球上另一大类重要植物群，拥有巨大基因组和漫长寿命的裸子植物（如松树），其基因组景观的成因却仍是一个谜团。是什么力量塑造了松树基因组中那些独特的高地与低谷？它们是如何形成和维持的？对这些问题的探索，不仅关乎我们对植物进化基本规律的理解，也可能揭示这些重要树种如何适应环境变化、维持遗传多样性。
为了解答这些问题，一个研究团队聚焦于三种具有重要生态和经济价值的松树：北美短叶松（Pinus banksiana）、扭叶松（Pinus contorta）和欧洲黑松（Pinus nigra）。他们利用中国油松（Pinus tabuliformis）的染色体级别基因组作为参考，对来自这三个物种总计100个个体的叶片样本进行了高深度的全基因组测序。通过一系列严谨的生物信息学分析，研究人员系统描绘了这三种松树的种群结构、演化历史，并深入探究了基因组范围内的重组率、遗传多样性、分化模式以及不同选择机制的信号。这项成果最终发表在了国际知名植物学期刊《The Plant Journal》上。
研究人员主要运用了以下几项关键技术：首先，对三个物种的100个个体进行全基因组重测序，样本来自捷克共和国阿布瑞特姆·索夫隆卡（Arboretum Sofronka）收集的北美和欧洲地理种源。其次，通过一系列比对、变异检测和质量控制流程获得高质量单核苷酸多态性位点数据。接着，利用软件（如pixy）估算核苷酸多样性、遗传分化、绝对遗传分歧等参数，并使用主成分分析和ADMIXTURE软件分析种群结构。然后，基于连锁不平衡（LD）利用LDhelmet软件估算种群尺度重组率。最后，通过整合遗传分化、多样性、重组率数据和Tajima‘s D、iHS等统计量，识别并分类“基因组岛”和“基因组谷”区域，推断其背后可能的选择机制。
遗传多样性与种群结构
分析显示，三个物种的核苷酸多样性存在差异，其中欧洲黑松最高，可能与亚种分化或与参考基因组的亲缘关系更近有关。功能区域中，四重简并位点多样性最高，其次是基因间区。主成分分析和群体聚类分析清晰地分开了三个物种，并部分支持了欧洲黑松的亚种划分，但未能完全支持扭叶松基于形态的亚种分类，表明其可能分化时间较近或存在持续的基因流。研究未发现近期种间杂交个体。
演化历史
通过SMC++模型推演的种群动态历史表明，三种松树在末次冰盛期（约2.5-5万年前）都经历了种群规模的收缩，随后在约1-1.25万年前（冰期后）发生了扩张。这印证了气候波动对松树种群历史的深远影响。对北美短叶松和扭叶松这对姐妹种的模拟分析表明，它们大约在1440万年前分化，且在分化早期存在不对称的基因流，即扭叶松向北美短叶松的基因流动更强。基因流速率的变化与北美大草原扩张等地质气候事件在时间上吻合。
全基因组连锁不平衡与重组率估计
与以往基于短基因片段的研究一致，松树全基因组范围内的连锁不平衡衰减很快。然而，基因组不同区域的重组率（以种群尺度重组率ρ衡量）差异巨大。研究发现，基因间区的重组率高于基因区，且重组率随与基因距离的增加而降低。一个重要的发现是，重组率与基因长度呈显著负相关，特别是与内含子累计长度负相关，而与外显子累计长度关系不明确。这意味着，长基因（尤其是富含长内含子的基因）倾向于位于低重组区域。此外，重组率还与转座元件（TE）密度负相关。这暗示，富含转座元件的长内含子可能通过抑制局部重组，促进了超长基因在松树中的保留。
分化、多样性景观及其相关性
研究人员比较了物种间的相对遗传分化、绝对遗传分歧和物种内的核苷酸多样性。总体而言，基因组大部分区域的多样性低于物种间的差异，表明存在异质性的基因组分化。尤其对于分化时间较短的北美短叶松和扭叶松对，基因组窗口的分布模式揭示了分化、多样性、重组等多个基因组特征之间存在复杂的相互关联。
“基因组岛”与“基因组谷”：不同选择机制的证据
基于极端值的分位数阈值，研究人员从基因组中识别出符合五种演化情景的“岛”和“谷”，并比较了它们的特征。结果表明，在“基因组岛”中，异域选择是主导机制，占主导地位的“岛”面积显著大于随机预期。而“基因组谷”主要通过平衡选择、背景选择或平行选择来维持。与背景水平相比，平衡选择区域的重组率显著更高。Tajima‘s D和iHS（综合单倍型分数）统计量的分析结果也支持了这些推断：平衡选择区域呈现正Tajima‘s D值（表明中等频率等位基因富集），而异域选择、背景/平行选择等区域呈现负值（表明稀有等位基因富集）。
在讨论和结论部分，研究者强调，他们的工作首次在裸子植物中系统揭示了重组率与基因长度（尤其内含子长度）的负相关关系，为理解松树等物种中异常长基因的进化与维持提供了新视角。长内含子中富集的转座元件可能通过抑制局部重组，在功能上“保护”了长基因。此外，研究明确了在已充分分化的松树物种中，异域选择是形成基因组分化高地（岛）的主要力量，而低分化区域（谷）则由平衡、背景或平行选择共同维持。这表明即使在长期分化的谱系中，不同类型的自然选择仍在持续塑造着基因组景观的细节。特别重要的是，研究发现重组率的局部变异与多种选择形式强烈相关，突显了重组在塑造基因组分化模式中的关键作用。本研究不仅弥补了裸子植物，特别是松属物种在比较基因组学和“基因组景观”研究方面的空白，也为理解长寿、非开花植物谱系中基因组结构的长期进化动态提供了宝贵的框架。这些发现对预测物种对环境变化的适应潜力、理解生物多样性形成机制以及指导林木遗传资源保护与育种均具有重要科学意义。
