巴拿马两大可可主产区种质组成与繁殖实践的对比研究

《Genetic Resources and Crop Evolution》：Contrasting germplasm composition and propagation practices in the two major cacao-growing areas in Panama

　　

巴拿马的可可农场主要由土著社区管理，记录可追溯至公元前1500年。然而，所使用的具体可可种质及相关管理实践在很大程度上仍未有文献记载。研究人员的主要目标是：（i）表征巴拿马两个主要栽培区（博卡斯德尔托罗省Bocas del Toro与埃雷拉省Herrera）农场中的可可种质；（ii）探讨可可种质类型与其管理及环境之间的关系。每农场选取5棵树（共16个农场）的根和叶样品，使用96-单核苷酸多态性(Single Nucleotide Polymorphism, SNP)阵列进行分析，同时收集环境与土壤数据。在博卡斯德尔托罗省，可可基因型表现为所有主要遗传组的连续混合，而埃雷拉省的 genotype 主要为阿梅隆纳多(Amelonado)和克里奥洛(Criollo)，包括具有纯克里奥洛祖先的个体。区域间有限的种质交换导致了截然不同的遗传组成。虽然农民主要使用种子建立种植园，但博卡斯德尔托罗省的嫁接（55%）显著多于埃雷拉省（16%）。回归分析确定全氮(Total Nitrogen)和pH是与遗传差异最相关的土壤变量。巴拿马可可的高遗传多样性对国家产业发展具有潜在价值。可可农场主要由土著社区管理，记录可追溯至公元前1500年。然而，所使用的具体可可种质及相关管理实践在很大程度上仍未有文献记载。研究人员的主要目标是（i）表征巴拿马两个主要栽培区（博卡斯德尔托罗省与埃雷拉省）农场中的可可种质，以及（ii）探讨可可种质类型与其管理及环境之间的关系。每农场选取5棵树（共16个农场）的根和叶样品，使用96-单核苷酸多态性(SNP)阵列进行分析，同时收集环境与土壤数据。在博卡斯德尔托罗省，可可基因型表现为所有主要遗传组的连续混合，而埃雷拉省的基因型主要为阿梅隆纳多和克里奥洛，包括具有纯克里奥洛祖先的个体。区域间有限的种质交换导致了截然不同的遗传组成。虽然农民主要使用种子建立种植园，但博卡斯德尔托罗省的嫁接（55%）显著多于埃雷拉省（16%）。回归分析确定全氮和pH是与遗传差异最相关的土壤变量。巴拿马可可的高遗传多样性对国家产业发展具有潜在价值。

论文解读：《Genetic Resources and Crop Evolution》发表研究解析巴拿马两大可可产区遗传与农艺差异

研究背景与问题提出

可可（Theobroma cacao L.，2n = 2x = 20）是具有重大经济意义的多年生作物，其遗传多样性是成功栽培的基础，可为遗传改良、抗逆及提升产量品质提供必要的变异来源。可可起源中心位于亚马逊上游盆地，向北传播至中美洲由玛雅和奥尔梅克人半驯化，向东进一步分化为其他遗传组。传统上分为克里奥洛（Criollo）、外来种（Forastero）和特立尼达（Trinitario），分子技术现已将其重新定义为十个不同的聚类群，包括Criollo、阿梅隆纳多（Amelonado）、纳纳西（Nanay）、孔塔马纳（Contamana）等。尽管全球开展了大量育种工作，但可可育种面临世代时间长、种子顽拗性（Recalcitrance）及嫁接不亲和等挑战。巴拿马作为假设的可可半驯化地峡，其可可栽培可追溯至公元前1500年，主要由Ng?be-Buglé等土著社区管理。博卡斯德尔托罗（Bocas del Toro，大西洋沿岸）自1920年代引入联合果品公司（United Fruit Company）克隆及CATIE基因型，较早采用嫁接技术；而埃雷拉省（Herrera，太平洋山区拉斯米纳斯Las Minas）则主要依赖历史悠久的种子繁殖低产Criollo类型。目前巴拿马全国平均产量仅103 kg ha^?1，远低于中美洲均值473 kg ha^?1，且西部巴拿马栽培的种质尚未有研究分析。鉴于种质交流有限及环境差异，明确两地遗传组成、繁殖实践与环境关联对于育种计划、就地保护及气候适应性策略至关重要。因此，研究人员旨在表征两区种质、探讨其与管理和环境的关系、验证嫁接普及度差异及土壤气候与遗传背景的关联。

主要关键技术方法

研究人员在巴拿马两大可可主产区博卡斯德尔托罗（11个农场，含IDIAP实验场）和埃雷拉拉斯米纳斯（5个农场）共16个农场，随机选取80棵可可树（原计划81棵因样本不足为80），分别采集每棵树的叶（接穗Scion）和根（砧木Rootstock）组织于硅胶干燥。同时收集农场管理与环境数据及土壤理化性质。DNA提取后在美国农业部（USDA）可持续多年生作物实验室采用96-SNP（Single Nucleotide Polymorphism）阵列在Fluidigm EP1^TM系统上进行基因分型，并纳入150份已知参考种质。数据分析使用GenAlEx计算个体杂合度（H_i），以≤6%差异定义遗传重复组（GDG）；采用STRUCTURE v2.3.4进行模型贝叶斯聚类（K=10），R语言包进行PCA、岭回归（Ridge Regression）及PERMANOVA环境差异检验。

研究结果

遗传结构与多样性

对80棵树根叶组织及参考群体SNP分析识别出10个稳健聚类，其中8个对应参考群体（Guiana、IMC Iquitos、Nanay、Criollo、SCA Contamana、Nacional、LCT EEN Curaray、Purus）。巴拿马种质表现为Amelonado集群内的混合梯度，介于Amelonado与Criollo之间，并有Nacional、IMC（Iquitos）和Nanay的贡献，整体无显著分支支持（AU=0%）。个体杂合度H_i介于0至75%，均值32%（博卡斯德尔托罗33%，埃雷拉29%），地上叶（34%）与地下根（29%）无显著差异。群体结构显示两区背景迥异：博卡斯德尔托罗以Amelonado、Nacional、Nanay为主，Criollo较少；埃雷拉以Criollo和Amelonado为主，含14份纯（>99%）Criollo个体。根叶组织遗传组贡献无显著差异，但区间差异极显著（F=35.88, p=0.001）。

农场内基因型的维护与交换

遗传图谱显示43%树木为不同来源的砧穗对（即嫁接），57%为单一遗传来源。博卡斯德尔托罗嫁接率（55%）显著高于埃雷拉（16%）（χ²=11.53, p<0.001）。跨160份样本鉴定出10个遗传重复组（GDG），如GDG-3为纯Amelonado，GDG-10为纯Criollo，分布于不同农场甚至区域，代表无性系或近缘同胞。部分GDG对应IDIAP保护的野生优选接穗或CATIE推荐商业品种，反映当地适应种质与引进品种并存。

环境因素与遗传组的关联

基于结构遗传组的PCA中PC1解释约64.8%方差，分离Amelonado与Criollo；PC2解释约25.4%，分离二者与其余组。两区环境差异显著（R²=0.58, p=0.001）。岭回归表明全氮（TN）关联PC1中Criollo减少与Amelonado增加（效应为其他变量2-10倍），PC2中TN关联Amelonado/Criollo减少与其他组增加；pH在叶组织中关联PC1中Criollo减Amelonado增，整体关联PC2中Amelonado/Criollo增与其他组减；有机质（OM）在合并分析中呈类似TN趋势。土壤性质与遗传结构存在显著耦合。

讨论部分总结

研究人员讨论指出，样本杂合度跨度大反映繁殖模式多样，低杂合度纯合群体（如Farm 5的Amelonado、Farm 14的Criollo）遗传变异低易受胁迫，高杂合度混合群体增强适应性。巴拿马种质为Amelonado-Criollo混合并含Nacional等，博卡斯德尔托罗高混合度源于20世纪初UF及CATIE克隆引入，而Guiana等组出现暗示土著社区自然或历史人为传播。埃雷拉纯Criollo罕见，可能源于Las Minas毗邻Veraguas原始林及Ng?be-Buglé文化隔离形成的避难