《Journal of Advanced Research》：Multi-omics dissection of panicle development reveals a GSE9-centered regulatory network controlling rice grain shape
编辑推荐：
摘要：水稻穗发育是决定产量与品质的关键农艺性状，但其潜在的分子网络，特别是由de novo（从头起源）基因调控的网络仍缺乏充分表征。为解决此空白，研究人员对de novo基因GSE9（GRAIN SHAPE ON CHROMOSOME 9）在穗发育期间的转录组
摘要：水稻穗发育是决定产量与品质的关键农艺性状，但其潜在的分子网络，特别是由de novo（从头起源）基因调控的网络仍缺乏充分表征。为解决此空白，研究人员对de novo基因GSE9（GRAIN SHAPE ON CHROMOSOME 9）在穗发育期间的转录组、代谢组和蛋白质组进行了整合分析。结果表明，GSE9缺失在关键阶段激活了广泛的转录重编程。功能分析揭示GSE9可能指导一个层级调控网络，可能通过调节关键转录因子（如ERF、WRKY和bZIP）来整合激素信号通路，包括赤霉素（GA）和脱落酸（ABA）。同时，GSE9缺陷导致广泛的代谢失调，特别是涉及苯丙烷和黄酮类生物合成的次生代谢，影响颖壳特性。应用于代谢组数据的机器学习（KANMB，Kolmogorov-Arnold Network for identifying Metabolic Biomarkers）方法鉴定了一组核心代谢物和共表达基因，表明GSE9协同调节淀粉代谢与次生代谢产物生物合成。蛋白质组分析进一步证实了这些通路的改变。关键的是，GSE9敲除诱导了全基因组转录-翻译解偶联，特异性地损害了同步的mRNA-蛋白协调性，而非简单地引入转录与翻译间的时间延迟。进一步分析揭示此效应源于转录主导（transcription-dominant）的调控模式，可能伴随翻译补偿（translational compensation），尽管缺乏直接证据。结论：本研究描绘了GSE9决定粒形的多层次调控机制，强调了基因编辑广泛的网络级效应，为基础研究和精准作物育种提供了关键见解。
研究背景与意义
亚洲栽培稻（Oryza sativa L.）是全球过半人口的主食，其粒形不仅是籼粳亚种分化的形态标志，也是影响产量、外观品质和市场价值的重要农艺性状，受植物激素信号通路复杂网络调控。De novo（从头起源）基因虽被视为进化创新驱动力并在真核生物生理过程中起作用，但其在主要作物形态分化中的功能尚不清楚。近期研究发现de novo基因GSE9（GRAIN SHAPE ON CHROMOSOME 9）起源于野生稻非编码区，通过单核苷酸变异产生起始密码子翻译为107氨基酸肽，调控水稻粒形及颖壳细胞扩增与增殖，且在粳稻中保守而在籼稻中为祖先非编码等位基因。然而GSE9调控粒形的分子机制不明。鉴于其在穗发育早期高表达，研究人员假设其通过调节颖壳细胞影响粒形，下游通路未明，故采用多组学解析其效应。该研究发表于《Journal of Advanced Research》，揭示了de novo基因在作物形态建成的网络整合机制，为高产育种提供资源。
主要关键技术方法
研究人员以野生型ZH11、GSE9敲除突变体（Cas9-1）及过表达株系（OE-1）为材料，在7个关键发育阶段（从孕穗至乳熟期）采集穗组织进行多组学测序。关键技术包括：转录组学（RNA-seq）采用Illumina NovaSeq 6000测序与DESeq2鉴定差异表达基因（DEG），加权基因共表达网络分析（WGCNA）解析模块；代谢组学采用LC-MS非靶向分析鉴定差异积累代谢物（DAM）并经KANMB机器学习框架筛选核心代谢物；蛋白质组学采用数据非依赖采集（DIA）质谱与DIA-NN定量分析差异表达蛋白（DEP）；整合分析通过计算Spearman秩相关评估mRNA-蛋白转录-翻译协调性，并结合生化测定验证苯丙烷途径酶活性及淀粉、黄酮、木质素含量。
研究结果
GSE9 disruption triggers global transcriptomic reprogramming during panicle development
研究人员对野生型ZH11、过表达OE-1及敲除Cas9-1在7个阶段的穗进行RNA-seq，PCA显示样本按发育阶段及基因型分离。DEG分析显示发育进程是基因表达变化主要驱动因子，而GSE9操作在阶段c（pre-heading）、b和g诱导最广泛的转录重编程，阶段c尤为关键（Cas9-1 vs ZH11达8320个DEG），表明其为调控网络活跃窗口。
Core co-expression modules reveal disrupted hormone signaling pathways
研究人员对13986个DEG进行WGCNA鉴定21个模块，Green模块富集赤霉素（GA）代谢与生物合成过程，Turquoise模块富集激素信号转导；鉴定906个差异表达转录因子（TF），聚类分析显示Cas9-1中簇III、V、VII（富含ERF、WRKY、bZIP家族）在早中期表达严重紊乱，GO富集于GA信号转导及ABA响应，表明GSE9可能通过调制核心TF网络整合激素信号，qRT-PCR验证相关系数0.87–0.95。
Metabolomic profiling reveals secondary metabolic dysregulation in gse9 mutants
代谢组检出2816种代谢物，鉴定2593种差异积累代谢物（DAM）。Cas9-1在大多数阶段DAM数量多于野生型，阶段a显著下调为主；K-means聚为8簇，基因-代谢物共响应分析鉴定13836个DEG与代谢物共调控。TF-代谢物共表达网络（|PCC|>0.8）显示228个TF与749个DAM相关，核心TF Os03t0180800-01与细胞壁加固代谢物（如tricin衍生物、Momilactone C）负相关，暗示GSE9抑制壁硬化代谢物维持塑性。
KANMB?based integration of metabolome and transcriptome identifies a GSE9?regulated metabolic network central to panicle development
研究人员应用KANMB模型从代谢组识别439个关键代谢物，关联964个高相关DEG（|PCC|>0.8），GO/KEGG富集于淀粉代谢、次生代谢物生物合成（苯丙烷、黄酮、二萜）。生化验证显示Cas9-1在d–e阶段黄酮（1.48倍）、木质素（1.42倍）及PAL/C4H/4CL活性（1.55–1.62倍）显著升高，淀粉降低（0.95倍）；OE-1反之。证实GSE9抑制苯丙烷代谢保壁塑性并促进淀粉合成，且影响激素信号受体激酶及糖转运蛋白（SWEET）表达。
Proteomic validation identifies key pathway alterations
DIA蛋白质组定量16048个蛋白，鉴定2287个差异表达蛋白（DEP），K-means聚为6簇，簇1/4/5/6富集植物激素信号转导，簇2/3富集黄酮生物合成。Cas9-1与ZH11在早期（a,c）及晚期（g）DEP差异显著，证明遗传修饰阶段性重塑蛋白子集。
GSE9 orchestrates dynamic transcriptional-translational coordination and its disruption leads to genome-wide decoupling
研究人员整合转录组与蛋白质组计算Spearman相关，野生型ZH11平均mRNA-蛋白同步相关0.379（高相关基因35.2%）；Cas9-1降至0.305（高相关仅20.6%），同步相关（Corr_sync）显著下降（0.3487 vs 0.4286, P<0.0001），而时滞相关（Corr_lag）无差异，表明GSE9缺失损害实时协调性而非时间延迟。调控分类显示47.8%基因为转录主导型（transcription-dominant），401基因为补偿/拮抗型（compensatory/antagonistic），提示GSE9以转录控制为主，翻译缓冲维持蛋白稳态。
讨论与结论翻译
讨论：GSE9作为穗发育关键期的激素信号与转录重编程中枢，WGCNA及TF分析表明其通过ERF/WRKY/bZIP间接整合GA与ABA信号（GSE9定位于质膜非核，不直接互作TF），阶段c为决定颖壳尺寸的窗口。代谢层面GSE9作为代谢抑制因子，缺失导致苯丙烷/黄酮异常激活增加壁刚性限制横向扩张致粒细长；过表达维持壁塑性形成圆粒。相较GS3、GW2等经典调控因子，GSE9为de novo起源无序肽，诱导全基因组转录-翻译解偶联，属高阶整合者。ML整合指出其双管齐下协调淀粉与次生代谢，反映de novo基因快速嵌入现有网络的特征。转录-翻译解偶联源于同步性丧失而非时序偏移，转录主导模式下翻译补偿缓冲扰动，囊泡运输等物流过程去耦合促成形态缺陷；需注意单株系潜在脱靶影响。
结论：本研究描绘了GSE9决定粒形的多层次调控机制，强调了基因编辑广泛的网络级效应，为基础研究和精准作物育种提供了关键见解。
