基于GBS全基因组关联分析挖掘柑橘果实品质性状的关键候选基因

时间：2025年11月22日

来源：BMC Genomics

编辑推荐：

本研究针对柑橘育种中果实品质性状遗传基础复杂、传统育种效率低的问题，通过基因分型测序(GBS)和全基因组关联分析(GWAS)对82份柑橘种质资源的23个果实品质性状进行遗传解析。研究鉴定出44,229个高质量SNP，发现354个显著数量性状核苷酸(QTN)，其中36个为多效性QTN。在5Mb连锁不平衡(LD)区域内鉴定出12个与中果皮宽度(MS)、果汁含量(JCE)等关键性状相关的候选基因（如具有氧化还原酶活性的Cs2g_pb011080和Cs7g_pb024190），为柑橘品质育种提供了分子标记和基因资源。

柑橘作为全球重要的经济作物，其果实品质直接影响市场竞争力。然而，柑橘育种面临诸多挑战：多年生作物幼年期长导致育种周期漫长；高度杂合性、无融合生殖和自交不亲和等生物学特性使得传统育种效率低下；果实品质性状多为数量性状，受多基因控制，遗传机制复杂。这些因素严重制约了柑橘品种改良进程。为了突破这些瓶颈，研究人员将目光投向现代分子育种技术。近日发表在《BMC Genomics》的研究通过整合基因分型测序(GBS)和全基因组关联分析(GWAS)技术，系统解析了柑橘果实品质性状的遗传基础。

研究人员首先收集了82份代表性柑橘种质资源，包括甜橙、柚、柠檬、葡萄柚等不同类型，在2020-2022年间系统评估了23个果实品质性状，包括14个形态性状、3个理化性状和6个营养性状。利用GBS技术生成高质量SNP标记，采用FarmCPU和BLINK多模型进行GWAS分析，显著提高了QTN检测效率。

研究团队通过GBS技术获得了215 million reads，覆盖38.78Gb序列数据，平均每个样本产生2.8 million reads。经过严格质控，最终获得44,229个高质量SNP，均匀分布在柑橘9条染色体上。群体结构分析将82份种质分为6个亚群，与地理起源高度一致，其中柚、 mandarin和 lemon/lime类群遗传背景相对纯净，而杂交种类群显示出明显的基因混合现象。

表型变异分析揭示显著遗传多样性

所有23个性状均符合正态分布，适合进行GWAS分析。果重(FW)、单果种子数(ANSF)等性状变异系数(CV)超过78%，表明群体具有丰富的遗传多样性。遗传力分析显示各性状遗传力在72%-99.1%之间，说明性状主要受遗传因素控制。相关性分析发现果实直径(FD)与果实长度(FL)呈极显著正相关(r=0.92)，而雌雄同体花数(FTH)与叶宽(LLW)呈显著负相关(r=-0.48)。

多效性QTN挖掘与候选基因鉴定

GWAS共检测到354个显著QTN，其中36个为多效性QTN，同时控制多个性状。这些多效性QTN主要与中果皮宽度(MS)、种子性状和叶片形态性状相关。特别值得注意的是，在第5号染色体上发现了一个多效性QTN热点区域，同时控制MS、JCE和ANSF等多个重要品质性状。

研究人员以显著SNP位点为中心，在5Mb LD区域内进行候选基因筛选，共鉴定出12个候选基因。其中，Cs2g_pb011080和Cs7g_pb024190两个基因编码氧化还原酶，可能通过调控柑橘苦味物质代谢直接影响果实风味品质。Cs4g_pb016990编码锰离子结合蛋白，参与细胞壁代谢，可能与果实质地形成相关。Cs5g_pb003860编码GDP-岩藻糖O-岩藻糖基转移酶，参与糖代谢途径，可能影响果实糖分积累。

性状特异性遗传机制解析

针对不同类别性状，GWAS揭示了其特异性遗传架构：果实特征性状检测到172个显著SNP，其中果重(FW)相关QTN解释的表型变异最大；种子性状关联分析发现92个显著位点，其中chr1_5129339位点与种子重量(SW)显著相关；营养性状鉴定出52个显著关联位点，叶片长宽比(RLW)相关位点主要分布在第8号染色体。

研究还发现一些有趣现象：虽然酸度(ACD)是重要的品质指标，但GWAS未检测到显著关联信号，可能与基因型-环境互作效应有关。此外，嫁接砧木的使用可能通过调节接穗基因表达影响表型测定，这为后续研究提供了重要启示。

群体结构与系统进化分析

ADMIXTURE分析显示K=6时交叉验证误差最小(0.24)，将群体分为6个亚群：葡萄柚、mandarin、柠檬/ lime、杂交种、柚和混合群。AMOVA分析表明72.2%的遗传变异存在于亚群间，27.8%存在于亚群内，说明群体具有明显的遗传分化。系统进化树将性状分为两大簇，第一簇包含生化性状和营养性状，第二簇包含果实特征性状，这种聚类模式与性状的生物学功能高度一致。