Lentil（Lens culinaris Medik.）抗性淀粉（RS）全基因组关联定位及其营养性状改良的候选基因解析

《Frontiers in Plant Science》：Genome-wide association mapping of resistant starch for improved nutritional traits in lentils

　　

引言：Lentil（Lens culinaris Medik.）是一种营养密集型豆类，在全球粮食安全和人类健康中具有重要作用。除富含蛋白质、膳食纤维和微量营养素外，小扁豆还含有抗性淀粉（resistant starch, RS），其作为膳食纤维的组成成分与改善血糖反应、增强肠道微生物群和降低炎症相关。然而，小扁豆中RS变异的遗传基础仍知之甚少。方法：研究人员在三个环境中对192份多样性小扁豆材料和5个对照品种的RS浓度、总淀粉（total starch, TS）浓度、百粒重（100-seed weight, 100-SW）和小区种子重量（plot seed weight, Plot-SW）进行了定量评价。利用5,583个高质量单核苷酸多态性（single nucleotide polymorphisms, SNPs）开展全基因组关联研究（genome-wide association study, GWAS）。在显著标记-性状关联（marker-trait associations, MTAs）的100-kb窗口内鉴定候选基因，并进行基因本体论（Gene Ontology, GO）和京都基因与基因组百科全书（Kyoto Encyclopedia of Genes and Genomes, KEGG）通路富集分析。结果：研究人员鉴定出9个与RS显著相关的MTAs、3个与TS相关的MTAs、20个与100-SW相关的MTAs和11个与Plot-SW相关的MTAs，可解释2.9%–42.6%的表型方差。候选基因包括参与淀粉生物合成的可溶性淀粉合酶III（soluble starch synthase III, SS3）。功能富集分析显示，生物学过程显著富集于氧化应激响应、解毒、淀粉与支链淀粉生物合成、ATP合成和生长素运输；分子功能显著富集于过氧化物酶、肽基脯氨酰异构酶和UDP-葡萄糖基转移酶活性；细胞组分显著富集于叶绿体、过氧化物酶体、液泡和细胞连接。讨论：这些发现为小扁豆种子中RS积累的遗传基础提供了见解，并鉴定了未来育种改良小扁豆品种营养品质的潜在遗传靶点。

研究背景方面，Lentil（Lens culinaris Medik.）是营养丰富的自花授粉冷季豆类，富含蛋白质、膳食纤维、矿物质和植物化学物，并含有作为膳食纤维功能组分的抗性淀粉（resistant starch, RS）。RS可逃避小肠消化，在结肠中被肠道微生物发酵产生短链脂肪酸（short-chain fatty acids, SCFAs），尤其丁酸，与改善血糖控制、胰岛素敏感性、降低炎症、增强饱腹感和维持肠道微生物平衡相关。尽管RS在水稻、小麦和大麦等作物中已有较多研究，但小扁豆RS积累的遗传基础仍 largely unexplored；小扁豆为二倍体（2n = 14），基因组约4 Gb且结构复杂，复杂数量性状需借助全基因组关联研究（genome-wide association study, GWAS）解析。已有GWAS在豆类中用于抗生物/非生物胁迫、种子营养和农艺性状，但RS相关GWAS仍处于起步阶段。为此，研究人员旨在评价多样性小扁豆面板中RS和总淀粉（total starch, TS）浓度的自然变异，利用GWAS鉴定与RS、TS、百粒重（100-seed weight, 100-SW）和小区种子重量（plot seed weight, Plot-SW）相关的基因组区域，并鉴定参与RS积累的候选基因，为标记辅助选择改良营养品质提供基因组资源。

关键技术方法概括（不超过250字）：研究人员使用来源于USDA-ARS（Pullman, WA）Lentil Single Plant（LSP）种质库的192份小扁豆材料和Avondale、Cedar、Crimson、Merritt、Pardina五个对照品种，在Spillman, WA、Central Ferry, WA及Irwell, Canterbury, NZ三个环境采用增广设计进行田间评价，测定RS、TS、100-SW和Plot-SW；以线性混合模型结合REML估计BLUP和广义遗传力（H²）；基于GBS（genotyping-by-sequencing）的5,583个高质量SNP，经PCA校正后采用BLINK模型开展GWAS，并以Bonferroni校正确定显著MTAs；基于LD和100-kb窗口筛选候选基因，通过BLASTP对Arabidopsis Araport11同源注释，并行GO/KEGG富集和SNP编码区效应注释。

研究结果如下。3.1 表型性状：通过三个环境的多环境表型评价和混合模型分析，研究人员发现RS、TS、100-SW和Plot-SW均存在显著环境效应及环境内区组效应，且各性状环境间分布差异明显；RS在中至高广义遗传力（H² = 0.632），TS遗传力低（H² = 0.182），100-SW遗传力很高（H² = 0.987），Plot-SW中等（H² = 0.440），表明RS等性状受遗传与环境共同影响。3.2 相关性分析：通过Pearson相关分析，研究人员发现RS与橙色种皮百分比（OrangeCP）和粉色种皮百分比（PinkCP）呈负相关，与棕色种皮百分比（BrownCP）和黑色种皮百分比（BlackCP）呈正相关；RS与TS呈正相关（r = 0.37），与Plot-SW呈正相关（r = 0.23），100-SW与RS呈弱负相关（r = -0.34），与TS呈弱负相关（r = -0.19）。3.3 标记-性状关联：通过BLINK模型GWAS和Bonferroni校正，研究人员鉴定出RS相关显著MTAs 9个，分布于1、2、4、5和7号染色体，解释9.1%–38.5%的表型方差，其中多环境BLUP分析中CHR2_24807231解释38.5%表型方差且关联最强；TS相关MTAs 3个，Spillman环境中CHR2_64012776解释42.6%表型方差；100-SW相关MTAs 20个，CHR1_319564320和CHR5_273800854在多个环境重复检测；Plot-SW相关MTAs 11个，表型方差解释率约10.1%–26.8%。3.4 候选基因与功能注释：通过在显著RS相关MTAs的100-kb窗口内筛选并基于Arabidopsis thaliana同源基因注释，研究人员发现候选基因包括淀粉合酶同源基因SS1、SS2、SS3、SS4和颗粒结合淀粉合酶1（granule-bound starch synthase 1, GBSS1）；GO富集显示α-1,4-葡聚糖葡萄糖基转移酶活性和葡萄糖基转移酶活性显著富集，细胞组分中淀粉体显著富集；KEGG富集包括苯丙烷生物合成、次级代谢物生物合成、代谢途径、ABC转运蛋白和玉米素生物合成，GO还富集于过氧化氢代谢、细胞解毒、活性氧代谢、过氧化物酶活性和抗氧化活性。37个独特显著MTAs中，15个位于编码区（10个同义、4个非同义、1个终止获得），14个位于内含子，2个位于3′ UTR，6个位于基因间区；RS相关MTA中CHR4_6631374为预测Val-to-Glu非同义变异，CHR4_389689669和CHR7_509092321为同义变异。最强关联MTA CHR2_24807231距Lcu.2RBY.2g009830约67.6 kb，该基因为Arabidopsis SS3同源基因；CHR5_472805423距两个β-半乳糖苷酶（β-galactosidases, BGALs）基因Lcu.2RBY.5g073830和Lcu.2RBY.5g073840分别约94 kb和96.4 kb；CHR4_389689669与UDP-葡萄糖基转移酶（UDP-glycosyltransferase, UGT）基因串联簇共定位，包括UGT85A23和UGT85A24（Lcu.2RBY.4g060560–580）。

讨论部分总结：讨论指出RS是影响餐后血糖和胰岛素反应的重要碳水化合物，并作为肠道微生物前体底物产生SCFAs，尤其丁酸，从而连接肠道与宿主糖脂代谢。三个环境中RS呈现较宽的地点依赖变异，而TS变异较窄；RS具有中至高广义遗传力，TS遗传力低，提示RS选择潜力及环境互作的重要性。种皮颜色与RS的相关性与深色种皮积累更多酚类物质、限制酶接近淀粉颗粒并减慢水解的已有认识一致；RS与TS正相关、与100-SW负相关，讨论将其与淀粉结构、直链淀粉含量、回生及环境管理敏感性联系起来。GWAS结果支持这些性状的多基因架构；RS多数MTAs为环境特异，两个通过跨环境BLUPs检测，环境特异位点应视为候选关联；CHR5_145398145同时与RS和TS相关，可能反映淀粉生物合成与分配的共享调控。候选基因方面，SS3参与支链淀粉中间链延伸并影响淀粉颗粒结构，与RS3形成相关；BGALs可能通过改变细胞壁孔隙和机械阻力影响淀粉生物可及性；UGTs可能作为上游调控因子影响碳水化合物分配或发育过程。讨论同时强调这些候选基因基于基因组邻近性和功能注释，尚未在小扁豆中功能验证，未来需通过KASP基因分型、独立群体和多环境验证确定稳定性及育种价值。

研究结论翻译：本研究通过整合多环境表型鉴定、混合模型分析、GWAS、候选基因鉴定和功能富集分析，为小扁豆（lentil）RS及相关种子性状的遗传架构提供了新见解。所观察到的RS、种皮颜色和种子重量之间的关系与先前关于多酚积累和淀粉生物合成影响种子组成与功能性的报道一致。混合模型分析显示RS存在显著遗传变异，而TS在不同环境间的总体遗传变异有限。GWAS鉴定出9个与RS显著相关的MTAs，包括位于2号染色体SS3基因附近的一个主要位点；SS3在支链淀粉链延伸和RS形成中具有重要作用。其他MTAs位于涉及细胞壁重塑（BGALs）和糖代谢（UGTs）的基因附近，提示碳水化合物代谢和细胞壁相关过程可能参与RS变异。然而，这些候选基因基于基因组邻近性和功能注释鉴定，尚未在小扁豆中进行功能验证。总体而言，这些发现突出了可能与RS相关的基因组区域和候选基因，并为该性状遗传调控的进一步研究提供了基础。未来研究应聚焦于将这些SNP转化为育种友好的KASP检测，并在独立小扁豆群体中验证其关联；一旦验证，这些标记可促进标记辅助选择，以培育营养强化的小扁豆品种。