全基因组测序揭示的马其顿本地甜樱桃（Prunus avium L.）的遗传多样性与群体结构

《Horticulturae》：Genetic Diversity and Population Structure of Autochthonous Macedonian Sweet Cherry (Prunus avium L.) Revealed by Whole-Genome Sequencing

　　

甜樱桃（Prunus avium L.）是一种经济上重要的水果物种，具有相当大的遗传变异，特别是在本地品种中。研究人员旨在利用全基因组重测序评估来自北马其顿奥赫里德地区的六个本地甜樱桃品种的遗传多样性和群体结构。共鉴定出约227万个单核苷酸多态性（SNP）和26.3万个插入缺失（InDel），表明基因组具有高度变异。群体结构和系统发育分析揭示了所研究品种中的两个不同遗传簇。Ohridska bela品种表现出最高的遗传分化程度，凸显了其作为宝贵遗传资源的重要性。遗传变异的功能注释显示，与开花时间、休眠和胁迫响应相关的基因具有显著变异，表明其对当地环境条件的适应性，而与果实成熟相关的基因则高度保守。此外，快速连锁不平衡衰减证实了群体内的高遗传多样性。这些发现强调了马其顿本地甜樱桃种质资源在育种计划、保护工作以及对未来环境变化的适应性方面的重要性。

### 论文解读：马其顿本地甜樱桃遗传多样性与群体结构的全基因组分析

#### 研究背景、问题与意义
甜樱桃（*Prunus avium* L.）是全球重要的鲜食水果，具有显著的遗传变异，尤其在本地品种中。马其顿奥赫里德地区保留了多个本地甜樱桃品种，这些品种长期适应局部土壤气候条件，表现出丰富的表型差异（如果实大小、颜色、成熟期）和较高的营养品质。然而，这些本地种质资源在现代生产系统中利用不足，其基因组层面的遗传特征尚不清楚。已有研究多依赖于简单序列重复（SSR）标记，而全基因组测序能提供更高分辨率的变异信息。本研究旨在通过全基因组重测序，系统评估六个马其顿本地甜樱桃品种的遗传多样性、群体结构，并鉴定与重要农艺性状相关基因的变异模式，为本地种质的保护、育种和适应气候变化提供基因组资源。论文发表在《Horticulturae》。

#### 主要关键技术方法
研究人员从北马其顿奥赫里德地区特雷贝尼什塔村的私人种质园采集了六个本地甜樱桃品种（Ohridska rana、Ohridska crna、Ohridska bela、Dolga siska、Krcka和Dalbazlija），每品种选取单株健康幼叶，送至Macrogen公司进行DNA提取、350 bp插入片段的文库构建及Illumina NovaSeqX平台150 bp双端测序。原始数据经质量过滤后，使用minimap2比对至甜樱桃参考基因组（*P. avium* cv. Satonishiki v1.0.a1），通过bcftools检测单核苷酸多态性（SNP）和插入缺失（InDel），并用SNPeff进行功能注释。群体结构分析采用最大似然系统发育树（Kimura 2-参数模型）、主成分分析（PCA）和ADMIXTURE模型聚类。连锁不平衡（LD）衰减通过PLINK计算r²与物理距离的关系。

#### 研究结果
**3.1 基因组测序与变异鉴定**
通过全基因组测序共获得33.4 Gb高质量数据，平均测序深度37.91×。鉴定出2,270,912个SNP和263,218个InDel，表明基因组高度变异。

**3.2 变异的全基因组分布与功能注释**
约65.62%的SNP位于基因间区，11.64%位于编码区。非同义突变与同义突变比值（dN/dS）为1.74，高于蔷薇科其他物种（如杏1.1、梨1.044），提示强烈正选择。高影响（high-impact）SNP占0.5%，包括移码和提前终止密码子突变；中度影响占7.098%。InDel中以移码突变为主。转换/颠换比（Ts/Tv）为1.403–1.417，杂合率平均38.68%。

**3.3 连锁不平衡（LD）分析**
LD在50 kb内快速衰减，符合高度多样、异交物种特征，表明关联作图需要高密度标记。

**3.4 品种间遗传关系与群体结构**
系统发育和PCA将六个品种分为两个簇：簇1包括Dalbazlija、Krcka和Ohridska crna；簇2包括Ohridska rana、Dolga siska和Ohridska bela。Ohridska bela分支最长，遗传分化最大。ADMIXTURE分析（K=2）支持上述分组；K=3时Ohridska bela形成独特谱系。表型变异（如果实大小、颜色）与遗传聚类一致。

**3.5 病原防御相关基因的变异**
在抗性相关基因中发现44个高影响SNP（如提前终止密码子），包括*RPM1*-样和*RGA2*-样基因（属NB-LRR家族）。Krcka和Dalbazlija携带高影响等位基因最多，Ohridska bela具有独特高影响变异（如*Pav_sc0000387.1_g020.1.br*纯合移码）。抗性相关基因存在大量高影响InDel（85%为移码），可能使部分基因成为假基因。但发病机制相关基因高度保守，缺乏高影响变异。

**3.6 开花时间与休眠相关基因的变异**
开花和休眠相关基因高度变异：在5511个SNP中，939个为高中影响变异，包括59个高影响突变（如终止获得、剪接破坏）。Ohridska bela拥有96个品种特异性变异，包括与休眠相关的MADS-box基因*DAM6*的提前终止突变（p.Gln152*）。Dalbazlija中*FAR1*基因纯合高影响突变（终止丢失和剪接区），解释其最早开花表型（4月10日）。其他基因如组蛋白甲基转移酶、*FRIGIDA*-样、*SOC1*-样、*HOS1*、*ELF3*等也发现高影响变异。InDel分析显示40个高影响InDel，Ohridska bela独特突变最多。

**3.7 果实裂果与成熟期相关基因的变异**
裂果相关基因中，所有品种共享一个高影响终止获得突变（*AGAMOUS*-样转录因子基因p.Gln141*），位于已知GWAS标记上游36 bp。Ohridska rana中独特的错义突变p.Leu122Gln可能改变裂果响应。剪接供体位点变异在Dalbazlija中纯合参考等位基因，使其保持功能剪接位点。成熟期基因高度保守：唯一的变异是*NAC5*转录因子基因中固定纯合终止获得突变（p.Gly349*），无品种特异性高中影响变异，表明成熟期受到稳定选择。

#### 总结讨论与结论
讨论部分指出研究局限性：每个品种仅选单株，缺乏生物学重复；六份材料的群体结构和LD分析初步，外推需谨慎。结论翻译如下：
本研究首次对马其顿本地甜樱桃种质池进行了全面基因组特征描述，揭示了关键物候和品质性状中遗传多样性与选择的复杂景观。分析识别出分层突变图谱：开花与休眠途径高度多样化（尤其Ohridska bela），而果实发育相关性状显著保守。所有品种中成熟基因固定终止获得突变表明区域成熟模式为经早期选择固定的祖先性状。裂果响应则受靶向破坏驱动，Dalbazlija因保留功能剪接位点而成为遗传异常值。这些发现凸显了马其顿本地甜樱桃种质作为独特等位基因库的价值，可用于提升气候适应性和果实品质的育种计划。最终结果为零星品种的保护及其融入现代全球樱桃产业奠定了关键基础。