《International Journal of Molecular Sciences》：Avocado (Persea americana Mill.) Genomic Divergence Across the Victoria, Tanganyika and Rukwa Lake Basins of Tanzania
编辑推荐：
牛油果（Persea americana Mill.）是坦桑尼亚关乎粮食安全（food security）与农村收入的具有重要经济价值的木本作物。尽管牛油果产业意义重大，但此前尚未对位于坦桑尼亚湖区（Lake Zone）、西区（Western Zone）以及鲁
牛油果（Persea americana Mill.）是坦桑尼亚关乎粮食安全（food security）与农村收入的具有重要经济价值的木本作物。尽管牛油果产业意义重大，但此前尚未对位于坦桑尼亚湖区（Lake Zone）、西区（Western Zone）以及鲁夸（Rukwa）南高地的牛油果种群开展基因组分化研究。虽然南高地已被认定为重要的生产中心，但关于西裂谷（Western Rift）和维多利亚湖盆地地方品种（landraces）的基因组多样性、群体分层（population stratification）和种族祖先（racial ancestry）仍存在显著的知识空白，这些区域以独特的微气候（microclimates）和与周边大湖国家的历史跨境贸易路线为特征。因此，本研究旨在通过低覆盖全基因组重测序（low-coverage whole-genome resequencing, lcWGS）对来自五个关键区域的105棵树的牛油果种质资源的基因组多样性进行表征，这些区域包括湖区和维多利亚湖盆地的Mara、Mwanza和Kagera；西区和西裂谷坦噶尼喀湖盆地的Kigoma；以及南高地鲁夸湖盆地（与坦噶尼喀湖盆地和西裂谷相邻）的Rukwa。本研究旨在识别局部分化和基因流的模式。研究人员假设这些区域的牛油果种群并非随机交配（panmictic），而是受地理和生态屏障（ecological barriers）影响而形成的群体结构，其中（i）与阿尔贝蒂尼裂谷（Albertine Rift）接壤的区域（即Kigoma和Kagera）表现出独特的西印度（West Indian）类型，可能由跨境交换和气候异质性（climate heterogeneity）塑造；（ii）Mwanza和Mara所在的维多利亚湖东部盆地可能含有独特的中间基因型，这些基因型是低地基因型种族间杂交（inter-racial hybridization）的产物。lcWGS筛选在总计405棵研究树木中产生了超过63,577,355个单核苷酸多态性（single nucleotide polymorphisms, SNP），这些数据是在与来自坦桑尼亚Morogoro及南高地Iringa、Mbeya和Ruvuma的300棵牛油果已发表lcWGS数据（95棵）以及美洲种族对照群体（205棵）整合后获得的。系统发育（phylogenetic）、无监督主成分分析（unsupervised principal component）和ADMIXTURE分析发现，Rukwa的种质与来自Iringa、Mbeya和Ruvuma的南高地种质聚类。相反，来自Kigoma和Kagera的样本形成了独特的西裂谷谱系（Western Rift lineage），富含西印度（West Indian）祖先成分，而Mara和Mwanza的种群占据中间位置，表明在维多利亚湖东部盆地存在广泛的种族间混合（inter-racial admixture）。这些遗传集群（genetic clusters）可指导东非牛油果地方品种（landraces）的利用。
牛油果（Persea americana Mill.）作为东非重要的园艺经济树种，对坦桑尼亚的粮食安全与农村生计具有不可替代的作用。尽管南高地（Southern Highlands）已被视为该国牛油果产业的核心产区，但维多利亚湖（Lake Victoria）、坦噶尼喀湖（Lake Tanganyika）及鲁夸湖（Lake Rukwa）盆地等西部与湖区地方品种的基因组多样性、群体分层（population stratification）及种族祖先（racial ancestry）至今仍未在基因组层面得到系统解析。这些区域不仅拥有独特的湖群微气候、裂谷地形与土壤异质性，更承载着与周边大湖邻国长期跨境贸易所遗留的种质交流历史，构成了研究牛油果种族分化与生态适应的天然实验室。特别是维多利亚湖和坦噶尼喀湖作为巨大的水体，能够调节周边陆地生态系统的温湿度，全年维持高湿度环境，这与牛油果中美洲起源中心的湿润亚热带条件高度相似，为外来种质的归化与生存提供了生态基础。与此同时，西裂谷地带从湖岸约700米到高地超过2000米的剧烈海拔变化，在短距离内形成了截然不同的垂直气候带与土壤类型，这种地质与地形异质性可能通过定向选择压力促使牛油果群体产生适应性分化。然而，既有研究多集中于南部与北部高地果园，对西部裂谷（Western Rift）及湖区种质资源的遗传身份尚缺乏深入认知，这严重制约了这些珍贵地方品种的有效保护、科学利用与可持续育种。鉴于此，研究人员开展本研究，旨在利用低覆盖全基因组重测序（low-coverage whole-genome resequencing, lcWGS）技术揭示坦桑尼亚三大湖流域牛油果种质的群体遗传结构，解析地理与生态屏障如何驱动该物种的局部分化，并为东非牛油果遗传资源的科学管理与可持续利用提供关键的基因组学依据。
在具体研究设计上，研究人员利用lcWGS技术对来自坦桑尼亚Mara、Mwanza、Kagera、Kigoma和Rukwa五个区域的105份成年优树（plus trees）样本进行DNA测序，并整合已发表的95份来自Morogoro及南高地（Iringa、Mbeya、Ruvuma）的lcWGS数据，以及205份源自美洲的种族对照（涵盖中美洲墨西哥、危地马拉、西印度种族及最近报道的西北南美安第斯与加勒比群体）。测序reads经Fastp质控后，使用BWA比对至牛油果var. Hass参考基因组，通过GATK HaplotypeCaller完成SNP calling与联合基因分型。群体遗传结构分析采用FastTree构建最大似然系统发育树、GCTA进行主成分分析（principal component analysis, PCA），并运用ADMIXTURE推断祖先成分（K=2至10，依据交叉验证误差确定最优K=9）。此外，通过严格过滤（最小测序深度10×、缺失率≤20%、最小等位基因频率maf≥5%）筛选出5011个高质量SNP位点，以降低分子标记伪重复并控制连锁不平衡（linkage disequilibrium, LD）对群体推断的影响。
研究获得了一系列关键结果。在lcWGS-Based SNP Markers方面，研究人员在总计405份样本中鉴定出超过63,577,355个SNP，转换/颠换比（Ts/Tv）为2.38，符合真核生物基因组特征；SNP注释显示大部分变异集中于基因间区和内含子等非编码区域，而外显子区非同义突变虽占比较低，但为后续功能基因挖掘提供了位点资源。在Phylogenetic, Multivariate and Admixture Reconstructions方面，系统发育分析表明，Rukwa种质与南高地（Iringa、Mbeya、Ruvuma）聚为单系分支，提示共享高地遗传背景；Kigoma和Kagera则形成独特的西裂谷（Western Rift）谱系，并与西印度（West Indian）及西北南美支系互为姐妹群，表现出显著的西印度祖先富集；Mara和Mwanza位于系统发育树的中间位置，呈现出高度混合的遗传组成，符合维多利亚湖盆地东岸广泛种族间杂交与本地分化的推断。PCA结果与此一致：第一主成分解释47.86%的遗传变异，将西北南美群体与其他种族分离；第二主成分解释13.31%的变异，将西印度种族与其他中美洲种族区分，而坦桑尼亚材料则沿第二轴呈连续混合分布。ADMIXTURE分析在K=9时达到最优，进一步证实Rukwa完全嵌入南高地分支，Kigoma和Kagera保持西印度祖先主导，Mwanza和Mara维持高混合度直至K=9，而Morogoro及美洲对照则呈现纯净的种族特异祖先成分，验证了该SNP数据集解析已确立牛油果种族基因库的能力。
综合讨论与结论，本研究首次在全基因组水平证实了坦桑尼亚西部及湖区牛油果种群并非随机交配，而是呈现显著的生态地理遗传分层。Rukwa与南高地的遗传一致性可能反映了共享的高地农业生态过滤及沿历史贸易路线的种植材料交流；Kigoma和Kagera独特的西印度祖先成分可能与阿尔贝蒂尼裂谷（Albertine Rift）及大湖地区跨境历史种质流动有关，而西裂谷陡峭的海拔梯度亦可能施加了定向生态选择；Mara和Mwanza的中间遗传位置则暗示了维多利亚湖盆地内的种族间杂交与原位生态驱动分化，Lake Victoria与Serengeti平原的地理隔离亦可能促进了本地谱系固化，尽管不完全谱系分选与遗传漂变的作用仍需进一步考量。研究结论明确指出，坦桑尼亚西南部、西部及湖区的牛油果基因组分化呈现三大趋势：（i）Rukwa与南高地种质的亲缘关系由共享的高地农业生态过滤及历史与当代种植材料交换驱动；（ii）西裂谷地区的Kigoma和Kagera以West Indian遗传祖先为主，可能源于与周边大湖邻国的历史及当代跨境种质交换；（iii）Mwanza和Mara形成遗传中间型和高混合群体，提示维多利亚湖盆地内存在种族间杂交和可能的本地选择。这些发现凸显了坦桑尼亚西部与湖区作为东非牛油果多样性次级中心的重要地位，强调有必要针对这些独特的地方品种开展原位保护，并将其整合至未来的砧木预育种与选择程序中，从而增强东非牛油果产业应对未来气候变化与市场波动的韧性。
