引言：花青素是彩色玉米籽粒中重要的天然色素和生物活性化合物。尽管彩色鲜食玉米因其营养和功能价值日益受到关注，但不同颜色玉米籽粒中花青素积累的分子基础仍不完全清楚。
方法：本研究以四种籽粒颜色对比鲜明的鲜食玉米品种为材料，包括白色（WN 2000）、黄色（HN）、多色（CN1）和黑色（HTN188），旨在探究花青素积累的机制。研究人员测量了不同发育阶段的总花青素含量，并进行了综合代谢组和转录组分析，以鉴定与籽粒色素沉着相关的花青素相关代谢物、差异表达基因和候选调控因子。
结果：代谢组学分析鉴定出49种花青素相关代谢物/特征，其中黑玉米品种HTN188中矢车菊素、芍药素和天竺葵素相关化合物显著富集。转录组学分析在四个品种间鉴定出12,557个差异表达基因，这些基因主要富集在苯丙烷生物合成、类黄酮生物合成、谷胱甘肽代谢、果糖和甘露糖代谢以及乙醛酸和二羧酸代谢途径中。若干花青素生物合成基因和调控转录因子候选者在富含花青素的籽粒中表达上调，包括DFR（Zm00001d011438）、MYB（Zm00001d003052）、bHLH（Zm00001d015990）和NAC（Zm00001d012508）。特别是，DFR（Zm00001d011438）的表达与总花青素含量呈强正相关。加权基因共表达网络分析进一步鉴定了与花青素积累相关的基因模块和候选中枢调控因子。
讨论：这些发现为不同颜色鲜食玉米籽粒中花青素积累提供了比较多组学视角。结果表明，黑玉米中高花青素积累与花青素相关代谢物的协同富集以及关键结构基因和转录因子的上调有关。本研究为未来的功能研究和花青素导向的鲜食玉米育种提供了候选基因和调控信息。
不同颜色玉米籽粒花青素积累机制的多组学研究解读
研究背景与意义
玉米是全球重要的主食作物，而鲜食玉米因其高直链淀粉含量及丰富的营养成分备受青睐。近年来，随着健康意识的提升，富含花青素等生物活性物质的彩色鲜食玉米逐渐成为研究热点。花青素是一类广泛存在于植物中的水溶性类黄酮多酚，不仅赋予植物组织红、紫、蓝等色泽，还具有抗氧化、抗炎、心血管保护等多种生理功效。然而，尽管已有研究表明花青素的积累受转录水平的严格调控，涉及苯丙烷/类黄酮途径的结构基因及MYB-bHLH-WD40（MBW）等转录因子复合体，但关于不同颜色玉米籽粒中花青素积累的深层分子机制，尤其是针对多种颜色表型的比较多组学研究仍相对匮乏。为了填补这一空白，研究人员对不同颜色鲜食玉米品种进行了综合分析，以期揭示其色素积累的分子基础，为彩色玉米的分子育种提供理论依据。该研究成果发表在了《Frontiers in Genetics》期刊上。
关键技术方法
研究人员选取了四个具有代表性的鲜食玉米品种：白色（WN 2000）、黄色（HN）、多色（CN1）和黑色（HTN188）。田间试验采用随机完全区组设计，在吐丝后15、21、27和33天（DAS）收获籽粒。代谢组学方面，利用液相色谱-串联质谱（LC-MS/MS）技术对花青素相关代谢物进行定性和定量分析。转录组学方面，提取籽粒总RNA并进行测序（RNA-seq），通过生物信息学手段筛选差异表达基因（DEGs）。此外，还结合了加权基因共表达网络分析（WGCNA）、转录因子预测以及实时荧光定量PCR（qRT-PCR）验证等技术手段，系统地解析了花青素代谢通路及相关调控网络。
研究结果
表型特征与总花青素含量
对不同颜色玉米品种的表型和总花青素含量测定显示，彩色甜玉米的花青素水平显著高于白色和黄色品种，其中黑玉米HTN188的花青素积累量最高。随着籽粒发育，着色品种（CN1和HTN188）的总花青素含量显著增加，而白色（WN 2000）和黄色（HN）籽粒则维持在较低水平且变化甚微。值得注意的是，虽然WN2000和HN在花青素含量上无显著差异，但两者在表型上存在明显区别，这暗示了除花青素外，其他色素（如类胡萝卜素）可能也参与了籽粒颜色的塑造。
样本全局表达模式分析
基于代谢组数据的PCA分析表明，所有生物学重复样本聚类紧密，重现性好。主成分1（PC1）解释了90.40%的方差，能清晰地将HTN188与其他品种区分开，说明HTN188具有最独特的代谢谱。而WN2000和HN聚在一起，形成一个“低色素/低花青素”簇，CN1则位于两者之间，表现出独特的代谢差异。
不同颜色玉米品种的代谢物谱分析
通过LC-MS/MS分析，共检测到49种花青素相关代谢物/特征，主要包括矢车菊素（16种）、芍药素（8种）、天竺葵素（7种）、飞燕草素（6种）、矮牵牛素（6种）和锦葵素（4种）。维恩图分析显示，HTN188拥有最多的独特差异积累代谢物（29种），且在与其他品种的比较中表现出最多的上调代谢物。这表明HTN188在花青素相关代谢物的富集程度上最为显著，其黑色籽粒表型可能与锦葵素、飞燕草素和矢车菊素的广泛积累有关。
转录组测序、组装与统计
高质量的RNA-seq数据成功获得了总计97.26 Gb的Clean Data，各样本的GC含量平均达49.59%，Q30碱基占比超过95.27%。比对到玉米参考基因组的平均比对率为88.69%，其中唯一比对率为83.45%。这些数据证明了转录组测序的高质量，适用于后续分析。
不同颜色玉米品种的转录组分析
在六个品种对的比较中，共鉴定出12,557个DEGs。其中，CN1与WN2000之间的DEGs数量最多（7,642个），HTN188与HN之间最少（3,423个）。韦恩图分析揭示了107个共有的DEGs和3,195个独特的DEGs。K-means聚类分析将DEGs分为9个簇，其中第5簇的基因表达量从白色到深色品种逐渐增加，这与花青素的积累趋势一致，可作为区分色素相关表型的潜在标记。
差异表达基因的功能注释
GO分析显示，DEGs主要富集在“生物过程”、“代谢过程”和“细胞过程”等条目。KEGG富集分析则指出，DEGs显著富集在“苯丙烷生物合成”、“谷胱甘肽代谢”、“果糖和甘露糖代谢”、“乙醛酸和二羧酸代谢”以及“类黄酮生物合成”等途径。特别是，“类黄酮生物合成”途径在深色品种（如HTN188与CN1）中更为活跃，而“植物激素信号转导”途径则在浅色品种（如WN2000）中表现突出。
花青素生物合成途径关键DEGs的鉴定
在花青素生物合成途径中，共有35个DEGs被鉴定出来，涵盖了苯丙氨酸解氨酶（PAL）、查尔酮合酶（CHS）、二氢黄酮醇4-还原酶（DFR）等关键酶编码基因。皮尔逊相关性分析表明，22个基因与总花青素含量呈正相关，13个呈负相关。其中，基因Zm00001d011438（编码DFR）的表达与总花青素含量呈极强正相关，是关联花青素积累的关键候选基因。
WGCNA模块与主要性状的相关性分析
通过WGCNA构建了共表达网络，识别出17个表达模块。模块与性状的相关性分析发现，MEblack模块与天竺葵素相关性状呈极显著负相关，而MEmagenta模块则与矢车菊素相关性状、花青素合酶（ANS）和UDP-黄酮苷葡萄糖转移酶（UFGT）的表达呈显著正相关。这表明MEmagenta模块包含了与花青素生物合成和调控相关的候选基因。
花青素表达的转录因子调控预测
研究人员在转录本中鉴定出6,928个转录因子（TFs），主要属于MYB、C2C2、AP2/ERF等家族。通过与35个花青素途径DEGs的互作网络分析，识别出一个中枢结构基因DFR（Zm00001d011438）和三个候选中枢转录因子：MYB（Zm00001d003052）、bHLH（Zm00001d015990）和NAC（Zm00001d012508）。这些中枢基因的表达模式与总花青素含量的变化趋势平行。
代谢组与转录组数据联合分析
以WN2000为对照的代谢物-基因相关性分析揭示了不同品种中特定代谢物与基因的关联。例如，在HTN188中，多个基因（如ANS、DFR、CHI等）与矢车菊素、芍药素和矮牵牛素的生物合成相关。ANS、CHI和4-香豆酸辅酶A连接酶（4CL）基因在所有比较组中均存在，突显了它们在花青素积累中的核心作用。
RNA-seq数据的qRT-PCR验证
为了验证转录组数据的可靠性，研究人员选取了10个基因进行qRT-PCR检测。结果显示，这些候选基因的相对表达趋势与RNA-seq数据完全一致，线性回归分析也证实了两者之间存在强相关性，从而验证了转录组数据集的可信度。
讨论与结论总结
讨论部分指出，不同颜色玉米籽粒的颜色差异不仅源于花青素的积累，还可能涉及类胡萝卜素等其他色素的贡献。代谢组学结果证实，HTN188含有最丰富多样的花青素化合物，尤其是矢车菊素、芍药素和天竺葵素的显著富集决定了其黑色表型。在分子机制上，深色玉米更倾向于激活苯丙烷/类黄酮途径，而浅色玉米则偏向谷胱甘肽代谢途径，这种次级代谢网络的资源分配权衡反映了植物在生长与防御间的平衡。研究特别强调了DFR（Zm00001d011438）作为关键结构基因的重要性，其表达量与花青素总量高度相关。同时，MYB、bHLH和NAC等转录因子可能作为上游调控因子，通过MBW复合体等机制调控花青素的生物合成。
结论部分表明，黑玉米品种HTN188的花青素含量显著高于其他品种。综合多组学分析揭示，HTN188的黑色籽粒表型与其矢车菊素、芍药素和天竺葵素相关化合物的增加积累密切相关。花青素生物合成途径中结构基因（如F3H、DFR和UFGT）的高表达，特别是DFR基因Zm00001d011438的表达，与花青素积累紧密相关。WGCNA分析进一步锁定了与花青素积累相关的共表达模块及MYB、bHLH、NAC等候选转录因子。这项研究为理解鲜食玉米花青素积累提供了多组学视角，并为未来的功能验证和遗传改良提供了宝贵的候选基因资源。
打赏
