整合QTL与GWAS揭示玉米气孔密度和大小多基因座遗传基础及ZmbHLH140与ZmCESA6功能验证

《The Plant Journal》：Multi-locus genetic basis of stomatal density and size in maize revealed by integrated QTL and GWAS

　　

气孔模式是植物光合性能与水分利用效率的基本方面，但玉米等C₄禾本科作物的遗传基础仍认识不足。研究人员整合了三个重组自交系（recombinant inbred line, RIL）群体的连锁分析与一个多样化玉米面板的全基因组关联研究（genome-wide association study, GWAS），以解析六个气孔性状的遗传结构。该方法鉴定出33个数量性状位点（quantitative trait loci, QTL）和51个关联位点，并揭示气孔密度（stomatal density, SD）与大小之间存在显著负遗传相关，多个基因组区间显示两类性状效应共定位。通过整合已报道的细胞类型特异性转录组数据，研究人员优先筛选出高置信候选基因，包括保守气孔调控因子以及先前未表征的、参与细胞壁生物合成等过程的候选基因。对候选基因ZmbHLH140和ZmCESA6的乙基甲磺酸（ethyl methanesulfonate, EMS）诱导突变体进行功能分析表明，ZmbHLH140功能缺失降低SD，而ZmCESA6突变具有相反效应。此外，ZmbHLH140位点的选择特征与气孔调控在玉米驯化中的作用一致。总之，研究结果建立了玉米气孔发育的综合遗传框架，突出了塑造气孔模式的保守与禾本科特异性机制，并为未来禾谷类作物气孔发育功能研究提供了基础。

气孔密度（stomatal density, SD）与气孔大小（stomatal size, SS）是植物光合、水分利用和逆境适应相关的重要性状，但其在玉米等C₄禾本科作物中的遗传基础尚不清楚。已有气孔发育框架主要来自模式植物拟南芥，涉及EPF类信号肽、TMM、ER类受体激酶及SPCH、MUTE、FAMA等bHLH转录因子；禾本科还存在副卫细胞形成、气孔成行排列及肌动蛋白动态等特有机制。玉米中SD与SS常呈负相关，二者独立优化的遗传约束尚待解析；以往研究多依赖单一作图群体或缺少功能基因组整合。为此，研究人员整合多个RIL群体的连锁作图、多样化关联面板的GWAS、细胞类型特异性单核转录组数据优先排序及EMS突变体功能验证，构建玉米气孔性状变异的综合遗传框架，并鉴定ZmbHLH140和ZmCESA6等候选基因，为禾谷类气孔发育机制解析和未来改良提供基础。

关键技术方法：研究人员在三个RIL群体BYK（K22×By815，163系）、BYD（DE3×By815，142系）、SZ（Z58×SK，105系；Pan et al.提供）中开展连锁作图，在262份玉米自交系（Yang et al.提供；含中国、CIMMYT和美国种质）中开展GWAS（混合线性模型），2023年海南三亚种植。表型采用穗位叶中部气孔印迹法。候选优先整合已发表单核转录组数据；EMS突变体ZmbHLH140（EMS4-008220）和ZmCESA6（EMS5-056202，B73背景）经回交、自交后以同胞株系比较。群体选择分析使用81份墨西哥类玉米、71份小颖类玉米和507份玉米自交系，计算π、F_ST、Tajima's D和XP-CLR。

研究结果：

遗传变异：研究人员在三个RIL群体和关联面板中量化SD、气孔长度（stomatal length, SL）、气孔宽度（stomatal width, SW）、气孔面积（stomatal area, SA）、气孔行数（row number, RN）和气孔排列比（stomatal arrangement ratio, SAR）六个性状，发现各性状呈连续变异并近似正态分布；SD与RN强正相关、与SL、SW和SA负相关，SAR与其他性状无显著相关，关联面板保持相同相关结构。该结果支持玉米SD与SS之间存在稳健的密度—大小权衡。

气孔性状QTL鉴定：通过三群体连锁分析，研究人员初始鉴定43个非冗余QTL区间、55个性状—QTL关联，最终保留33个稳健QTL区间，单个QTL解释5.08%—16.87%的表型方差，其中14个PVE大于10%；GWAS鉴定128个显著非冗余单核苷酸多态性（single-nucleotide polymorphism, SNP），定义51个非冗余关联位点，48个与既往气孔或干旱响应研究区间重叠。连锁与关联分析共同支持6个QTL区间，表明存在主要效应位点、跨群体重复检测和两类方法共定位的结构化遗传架构。

QTL区间内候选基因优先排序：研究人员结合已知气孔调控因子同源性、功能注释和细胞类型特异性表达，分别保留40个关联分析和67个连锁分析核心候选基因，归纳为保守气孔调控因子、细胞壁生物合成与修饰、信号与代谢基因三类；前者包括SPCH同源bHLH、EPFL和SCR相关GRAS成员，细胞壁类包括CESA6、IRX15A及多糖修饰酶。基于此，研究人员选择ZmbHLH140和ZmCESA6进行功能验证。

候选基因调控气孔性状的功能验证：EMS突变体分析显示，Zmbhlh140纯合突变体SD、RN和气孔指数（stomatal index, SI）显著降低，SL增加而SW和SA不变；候选基因关联分析鉴定两个非同义SNP与SD显著关联，单倍型间SD存在差异，群体选择统计显示ZmbHLH140位点多态性降低、Tajima's D为负及F_ST/XP-CLR中度升高，但未超过全基因组经验阈值，提示与驯化选择相关但非定论。zmcesa6纯合突变体SD和RN显著增加，SL、SW和SA降低，SI不变，且未见明确选择信号。结果支持SD受细胞身份调控网络与细胞壁介导生长约束共同塑造，且不同基因呈现完全或部分权衡效应。

讨论部分总结：研究人员指出，SD与SS的负相关具有遗传相关证据，多个QTL可同时影响密度和大小；但仅凭QTL共定位不能区分真多效性与紧密连锁。研究表型来自单一田间环境，6个共享区间缺少区间重叠的正式置换检验，且262份面板受玉米快速连锁不平衡（linkage disequilibrium, LD）衰减限制，因此共享区间仅为提示性；与Xie et al.（2021）独立作图研究比较，33个稳健QTL中5个、51个关联位点中16个重叠。SAR表现为遗传上相对独立的性状。保守调控因子和细胞壁生物合成基因的反复检出支持所揭示架构；ZmbHLH140和ZmCESA6功能验证显示基因依赖性权衡，ZmCESA6表现完全权衡，ZmbHLH140表现部分效应。局限包括单发育阶段验证、成熟叶未验证、QTL群体SI无法量化、气孔性状与抗旱或产量之间尚无直接表型与遗传相关分析，ZmbHLH140选择证据为提示性而非确定性。

研究结论翻译：总之，该研究通过整合连锁作图、关联分析、转录组优先排序和功能验证，建立了玉米气孔性状变异的综合遗传框架。研究结果提示，SD与SS之间的经典权衡可能受多效性遗传位点、保守发育调控因子和细胞壁介导的机械控制影响，但多效性程度具有基因依赖性，某些情况下表现为完全权衡（ZmCESA6），另一些情况下表现为部分效应（ZmbHLH140）。通过揭示这些过程如何互作以塑造气孔模式，该工作推进了对禾本科植物气孔发育的理解，并为后续对研究中鉴定到的遗传调控因子开展功能验证提供了基础。