一种公开的中密度基因分型平台用于美洲山核桃[Carya illinoinensis (Wangenh.) K. Koch]

《The Plant Genome》：A public mid-density genotyping platform for pecan [Carya illinoinensis (Wangenh.) K. Koch]

　　

美洲山核桃[Carya illinoinensis (Wangenh.) K. Koch]是全球栽培第五大树坚果，80%的产量来自美国南部各州。尽管山核桃具有经济和健康效益，但育种者缺乏用于作物改良的基因组工具。山核桃育种群体较小，多数育种项目在采用技术方面存在诸多障碍，尤其是在创建和使用遗传标记面板以进行基于基因组的选种决策所需的成本和专业知识方面。研究人员报告了一个包含3100个位点的DArTag（Diversity Array Technology [DArT]）面板的创建，这些位点分布于二倍体山核桃基因组，用于分子育种和基因组预测。研究人员展示了该面板区分栽培山核桃的关键亲本和奠基种与其他可能用于预育种的Carya物种的能力。研究人员还证明，所得数据可以在双亲群体中创建连锁图谱。该标记面板的创建为山核桃育种项目带来了成本效益高且快速的分型能力，使常规基因分型成为任何山核桃育种者的现实。此外，该平台提供的开放获取功能使得在不同项目、机构和国别之间，基于该标记面板生成的遗传数据集能够进行比较和整合。

这篇发表在《The Plant Genome》上的论文聚焦于美洲山核桃[Carya illinoinensis (Wangenh.) K. Koch]育种中缺乏基因组工具的问题。美洲山核桃是美国和墨西哥本土的重要坚果作物，美国供应全球80%的产量，但其育种进展滞后，主要原因包括幼龄期长（7–10年）、缺乏基因组资源，以及小规模育种群体面临的技术和成本门槛。传统分子标记如简单序列重复（SSR）、扩增片段长度多态性（AFLP）和随机扩增多态性DNA（RAPD）数量有限，而基于基因分型测序（GBS）的单核苷酸多态性（SNP）数据虽多，但缺失率高、重复性差，难以用于持续选择和标记辅助选择（MAS）。为了克服这些障碍，研究人员开发并验证了一个公开的中密度DArTag（Diversity Array Technology [DArT]）基因分型面板，包含3100个位点，旨在为山核桃育种提供低成本、高通量的基因分型解决方案。

研究人员利用一个由48份C. illinoinensis种质组成的多样性面板，通过比对到参考基因组v1.1，发现181,589,925个SNP，并经过严格过滤（质量>20、双等位、次要等位基因频率>5%、缺失率<25%、连锁不平衡剪除r²<0.9、位于低多态性区域、优先选择基因区非重复区域的SNP、Hardy-Weinberg平衡检验p>0.01）获得13,918个高置信SNP。经DArT公司质控后剩余12,817个，进一步通过Breeding Insight精选为3100个位点，确保基因组均匀覆盖，并包含30个与农艺性状潜在关联的标记。测试样本包括376份材料，涵盖16个Carya物种（186份多样性样本和190份F₁双亲群体，亲本为‘Lakota’和‘87MX3-2.11’），叶片组织送DArT进行DNA提取、文库构建和测序。原始数据以缺失等位基因发现计数（MADC）文件提供，包含所有微单倍型（81 bp amplicon）的读数。经BLAST比对和固定等位基因ID分配，最终获得10,199个微单倍型（来自3073个成功扩增的位点）。基于读数的等位基因频率进行主成分分析（PCA）；使用多态性信息含量（PIC）公式评估标记多样性；采用updog软件进行剂量调用（二倍体模式）；利用MAPpoly2构建连锁图谱。

研究结果按原文小标题总结如下：
- **3.1 3K山核桃DArTag分型面板的性能**：对376份二倍体山核桃种质的评估显示，整体样本平均缺失率为9.0%，其中F₁群体平均缺失率仅4.0%，多样性面板为14.2%。90%的样本至少80%的标记位点成功扩增，表明面板检测效率高，重复性好。
- **3.2 主成分分析揭示山核桃群体的遗传模式**：基于读数的PCA显示，前两个主成分（PC）分别解释24.0%和12.4%的遗传方差。F₁群体紧密聚集于两亲本（‘Lakota’和‘87MX3-2.11’）之间，符合双亲遗传预期。多样性样品中，C. illinoinensis种质相对聚集，而不同Carya物种聚类分离明显，种间杂种位于亲本种之间，验证了面板准确区分种内和种间遗传关系的能力。
- **3.3 多样性群体中DArTag位点和SNP的多态性信息含量**：去除单态标记后，对3016个DArTag位点（10,084个微单倍型）的微单倍型PIC值估计范围为0.0001–0.825，平均值0.334；对7442个SNP的SNP-PIC值范围为0.0001–0.375，平均值0.182。超过67%的位点中，微单倍型PIC高于SNP，表明微单倍型能捕获更多遗传变异。
- **3.4 F₁和多样性山核桃群体中序列多态性的分布与比较**：微单倍型分布显示，F₁群体检测到7229个，多样性群体检测到9849个。SNP水平上，剂量调用后多样性群体中有6232个信息性SNP，F₁群体中有1573个。多样性群体中更多标记符合其更广的遗传基础。
- **3.5 山核桃连锁图谱的构建**：利用F₁群体（n=182）和7344个SNP，经过滤（去除单态、异常分离、冗余标记）后剩余1114个标记，最终1053个SNP成功定位到16个连锁群（LG），对应山核桃单倍体染色体数。图谱总长1639.2 cM，LG长度62.3–143.4 cM，平均每cM 0.64个标记，最大间隙21.6 cM（LG 6）。遗传图谱与物理图谱比较显示多数染色体存在重组抑制区域，表明标记分布均匀，为QTL定位和MAS提供了可靠框架。

讨论部分总结了该面板的开放可用性（通过DArT订购），高扩增率（样本水平缺失率<15%），以及适用于遗传图谱构建、MAS、关联研究、基因组选择、重组模式重建、等位基因剂量估计和亲本确认等多种育种应用。PIC值与之前基于RAPD、SSR、indel的研究相比处于合理范围。工作流程周转时间3周，成本约每样本15美元，且面板可通过重新混合寡核苷酸池进行更新（每年一次低费或免费）。未来可通过深度测试建立固定等位基因数据库以改进跨项目数据整合，必要时可补充至7K面板或制作互补3K面板以增加标记密度。

研究结论部分可翻译为：该3K DArTag面板为山核桃育种者提供了一种成本效益高、快速且可重复的基因分型工具，使常规基因分型成为现实，且其开放获取特性促进了不同项目、机构和国别间遗传数据的比较与整合。该面板已显示其在区分亲本与物种、构建连锁图谱以及支持基于基因组的预测育种方面的能力，有望显著加速山核桃品种改良进程。