全基因组关联研究解析花生籽粒抗黄曲霉菌的基因组区域与候选基因

《Theoretical and Applied Genetics》：Genome-wide association studies unravel genomic regions and candidate genes associated with Aspergillus flavus resistance in peanut kernels

为应对花生易受黄曲霉菌（Aspergillus flavus）侵染并产生强致癌性黄曲霉毒素（AFB1、AFB2等）的食品安全挑战，本研究通过多环境表型评价结合高密度SNP标记，对353份花生种质进行全基因组关联分析（GWAS）。研究鉴定出10个显著SNP关联、3个稳定数量性状位点（QTL）以及3个关键候选基因（如编码TIR-NBS-LRR类抗病蛋白的基因），为花生抗黄曲霉育种提供了新的遗传资源和分子标记，对推动抗性品种选育、保障花生产业健康发展具有重要意义。

　　

花生，这种富含油脂和蛋白质的作物，不仅是重要的食用油来源，也是广受欢迎的零食。然而，在田间和储藏过程中，它常常受到一种名为黄曲霉菌（Aspergillus flavus）的真菌侵袭。这种真菌的可怕之处在于，它能在花生籽粒中产生一类被称为黄曲霉毒素（Aflatoxin）的强致癌物质，其中以黄曲霉毒素B₁（AFB₁）的毒性和产量最高。这类毒素污染严重威胁食品安全和人类健康，尤其是增加了肝癌等疾病的发病风险。长期以来，人们主要依靠物理化学方法去毒，但这些方法效果有限、成本高且可能带来化学残留。因此，从根源上培育自身就能抵抗黄曲霉侵染或抑制毒素合成的花生品种，被认为是经济、有效且可持续的解决方案。然而，大多数现有花生品种仍缺乏足够的抗性，挖掘新的抗性遗传资源和关键基因，成为花生育种中亟待突破的瓶颈。

为了破解这一难题，一项发表在《Theoretical and Applied Genetics》上的研究，对353份具有丰富遗传多样性的花生种质资源发起了“攻坚战”。研究人员的目标很明确：第一，精准评估这些种质对黄曲霉菌胁迫的抗性水平，寻找“双抗”（既抗侵染又抗产毒）的优异材料；第二，利用先进的基因组学手段，定位控制抗性的关键遗传位点，并筛选出候选基因，为分子育种提供“导航图”。

为达成研究目标，研究人员整合运用了多项关键技术。首先，他们对一个包含353份花生种质（涵盖5个植物学类型）的自然群体进行了全基因组重测序，获得了935,231个高质量单核苷酸多态性（SNP）标记，构建了高密度遗传图谱。其次，在三个不同环境（河南原阳2020、商丘2020和2021）下，对所有这些种质进行了人工接种黄曲霉菌的表型鉴定。他们采用体外种子定殖（IVSC） assay测定侵染指数，并利用高效液相色谱-质谱联用（HPLC-MS）技术精准定量了籽粒中的AFB₁和AFB₂含量。最后，通过整合多环境表型数据与全基因组SNP数据，运用混合线性模型（MLM）进行了全基因组关联分析（GWAS），以寻找与抗性显著相关的遗传标记。

研究结果

表型变异

研究发现，花生籽粒对黄曲霉菌侵染的抗性（以侵染指数衡量）和产毒抗性（以AFB₁/AFB₂含量衡量）在种质间存在广泛变异，且表型呈连续正态分布，表明这些性状由多基因控制。三个环境下的侵染指数以及AFB₁/AFB₂含量之间均存在极显著正相关，且遗传力（H²）均高于0.90，说明遗传因素是表型变异的主因，且抗性表现较为稳定。研究共鉴定出13份在三个环境中均表现稳定的“双抗”种质，其中C203和C206表现尤为突出。侵染指数与AFB₁、AFB₂含量高度正相关，且AFB₁与AFB₂含量也紧密相关，暗示它们可能受共同遗传机制调控。

比较分析

对不同植物学类型和生长习性的比较发现，属于var. hypogaea类型和匍匐/半匍匐生长习性的种质材料，其侵染指数和黄曲霉毒素含量的中位数显著较低，且在鉴定出的13份抗性种质中占比很高。这表明在未来的抗性筛选中，可以优先关注这些类型的种质资源。

全基因组关联研究

利用MLM模型进行的GWAS分析，共检测到10个显著的标记-性状关联（MTAs）。其中，1个SNP在至少两个环境中稳定与侵染抗性相关，2个SNP同时与AFB₁和AFB₂含量相关。值得注意的是，与侵染抗性相关的SNP位点，和与产毒抗性相关的SNP位点在物理位置上相互独立，这从基因组层面支持了“侵染抗性”和“产毒抗性”是相对独立遗传的性状这一观点。

稳定QTL分类与候选基因鉴定

根据连锁不平衡（LD）衰减估计，以显著SNP上下游各70 kb的区间定义了候选数量性状位点（QTL）。最终确定了3个稳定的QTL，分别命名为qII_A05（与侵染抗性相关）、qAF_A03和qAF_A05（均与产毒抗性相关）。在这些QTL区间内，结合基因功能注释和已有研究，初步筛选出13个候选基因。进一步的单倍型（Haplotype）分析验证了其中3个关键候选基因：Arahy.7046BI.1（编码一个含有PH/START结构域的蛋白）、Arahy.TX2FLU.1（编码一个类受体激酶）和Arahy.ASR4JM.1（编码一个TIR-NBS-LRR类抗病蛋白）。它们的单倍型与抗性表型存在显著关联。其中，与侵染抗性相关的两个基因的有利单倍型（Hap2）能使侵染指数平均降低27%以上，且在种质库中频率高达85.43%和90.6%；与产毒抗性相关的基因的有利单倍型（Hap1）能降低约26-29%的黄曲霉毒素含量，频率高达94.12%。这些高频率的有利单倍型意味着它们易于在育种中被直接利用。

结论与讨论

本研究系统地揭示了花生抗黄曲霉菌胁迫的遗传基础。通过大规模种质鉴定，获得了包括C203、C206在内的多份优异“双抗”资源，为育种提供了直接亲本。通过GWAS定位了3个新的稳定QTL，拓展了该性状的已知遗传图谱。更重要的是，通过单倍型分析验证了3个功能明确的候选基因，它们分别涉及脂质信号转导（PH/START蛋白）、病原模式识别（类受体激酶，RLK）和效应子触发免疫（TIR-NBS-LRR蛋白），这为了解花生应对黄曲霉侵染的复杂免疫机制提供了重要线索。这些基因的有利单倍型在现有种质中频率很高，便于通过分子标记辅助选择（MAS）快速导入优良品种中。

该研究的发现具有多重重要意义。在实践层面，它提供了可直接用于育种的抗性种质、稳定的QTL区间和分子标记，将显著加速抗黄曲霉花生品种的选育进程。在策略层面，研究指出优先筛查var. hypogaea类型和匍匐型种质，能更高效地富集抗性资源，为种质资源工作者提供了明确方向。在科学层面，研究鉴定的候选基因关联了植物免疫的核心通路，为后续深入解析花生与黄曲霉菌互作的分子机制奠定了坚实基础。总之，这项工作为从源头上控制花生黄曲霉毒素污染、保障食品安全和产业健康发展，提供了坚实的遗传学基础和宝贵的育种资源。