《Scientific Reports》：Cell wall remodeling and inositol metabolism coexpression modules associated with nut size in Carya illinoinensis cvs. ‘Mahan’ and ‘Tiny tim’
编辑推荐：
本研究针对山核桃坚果大小这一重要农艺性状的分子机制尚不明确的问题，通过开展'Mahan'和'Tiny Tim'两个品种的时序转录组分析，利用WGCNA技术鉴定出与细胞壁生物合成、肌醇代谢相关的关键共表达模块。研究发现'Mahan'品种中细胞壁组织模块和'Tiny Tim'中蛋白酶解模块分别与坚果大小显著相关，为山核桃品种改良提供了新的分子靶点和理论依据。
山核桃作为一种原产于美国和墨西哥的珍贵坚果树种，其全球市场规模已超过20亿美元。在极为有限的育种周期内，坚果大小始终是品种改良的首要目标性状。然而，关于山核桃坚果个体发育的分子基础及其大小调控机制的研究仍相对匮乏。除了坚果大小问题外，山核桃果实在整个生长季节还面临诸多生理紊乱的挑战，这使得明确各生长阶段的关键基因成为培育新品种和制定管理措施的必要前提。
为深入解析山核桃果实发育的分子机制并挖掘大果表型的候选基因，研究人员开展了一项针对'Mahan'和'Tiny Tim'两个基因型的时序转录组研究。通过加权基因共表达网络分析(WGCNA)技术，将转录本按功能聚类，并利用模块相关性分析鉴定枢纽转录本，从而筛选出可能驱动这些功能簇的关键因子。
研究团队采用的主要技术方法包括：对两个山核桃品种('Mahan'和'Tiny Tim')进行系统采样，利用高通量测序技术获取不同发育阶段的转录组数据，运用WGCNA算法构建基因共表达网络，通过模块-性状相关性分析筛选关键功能模块，并采用拓扑特征识别方法确定各模块中的枢纽基因。
共表达模块鉴定与特征分析
通过WGCNA共识别出18个基因共表达模块，其中蓝色模块(包含1，247个转录本)与'Mahan'品种的坚果大小呈显著正相关(r=0.89， p<0.001)。该模块富集于细胞壁生物合成和重构相关通路，包含扩展蛋白(expansins)和木葡聚糖内转葡糖基酶(xyloglucan endotransglucosylase/hydrolases)等关键基因。
品种特异性模块功能解析
在'Mahan'品种中，黄色模块(包含893个转录本)显著富集于肌醇代谢通路(r=0.85， p<0.01)，包含肌醇加氧酶(myo-inositol oxygenase)和磷酸肌醇激酶(phosphoinositide kinases)等基因。而在'Tiny Tim'品种中，褐色模块(包含1，105个转录本)与坚果大小负相关(r=-0.82， p<0.01)，该模块主要涉及蛋白酶解过程和脱落酸(abscisic acid， ABA)信号响应。
枢纽基因网络拓扑分析
通过计算模块内连通性(intramodular connectivity)和模块成员值(module membership)，鉴定出各关键模块中的核心枢纽基因。在'Mahan'的蓝色模块中，CcXTH1(木葡聚糖内转葡糖基酶)表现出最高的连接度(k=28)，而在黄色模块中，CcMIOX2(肌醇加氧酶同源物)的模块成员值达到0.98。
代谢通路交叉验证
进一步分析发现，细胞壁重塑模块与肌醇代谢模块存在显著的功能互作。肌醇衍生的信号分子如磷脂酰肌醇(phosphatidylinositol)可能通过调控细胞壁松弛酶活性影响果实膨大过程，这为理解山核桃坚果大小形成的多通路协同调控提供了新视角。
本研究系统揭示了山核桃坚果大小形成的分子调控网络，发现细胞壁动态重塑与肌醇代谢信号在决定最终坚果大小中的核心作用。研究不仅为山核桃分子育种提供了有价值的候选基因标记，更重要的是建立了多组学数据整合分析果树性状的研究范式。这些发现对解决山核桃产业中果实大小不均和生理障碍等问题具有重要实践意义，为未来开展精准育种和栽培管理创新奠定了理论基础。论文中鉴定的关键模块和枢纽基因为进一步解析木本植物果实发育的进化保守性和物种特异性提供了重要线索。
