编辑推荐：
大豆籽粒蛋白含量在加拿大东西部产区存在显著差异（西低东高1%-5%），其环境调控机制尚不明晰。为解决此问题，研究人员通过多加权基因共表达网络分析（multiWGCNA），解析了不同环境和蛋白质含量下的大豆叶片转录组。研究揭示了由环境和蛋白质含量共同调控的关键共表达模块，这些模块涉及表观遗传调控、翻译调控、光合作用、碳氮代谢及种子蛋白储存液泡填充等生物学过程。这项工作阐明了影响大豆籽粒组成的复杂环境调控网络，为定向改良大豆品质以适应不同产区提供了重要分子线索。
大豆，作为全球最重要的油料和蛋白作物之一，其种子富含的蛋白质是人类和动物营养的重要来源。然而，一个有趣的现象长期困扰着加拿大的大豆种植者和科学家：在加拿大西部产区（如曼尼托巴省和萨斯喀彻温省）种植的大豆，其籽粒蛋白质含量往往比东部产区（如安大略省和魁北克省）种植的大豆低大约1%到5%。这种差异持续了二十多年，但其背后的分子机制却如同笼罩在迷雾之中。是更短的生长季、更少的降水、还是更低的温度在起作用？环境因素是如何“指挥”大豆体内的基因，最终影响种子中蛋白质的积累的？回答这些问题，对于培育适应加拿大不同气候区的高产优质大豆品种，保障国家粮食安全和农业竞争力具有至关重要的意义。
为了拨开这层迷雾，一支由加拿大农业及农业食品部（AAFC）渥太华研发中心领导的研究团队展开了一项雄心勃勃的研究。他们选取了10个具有不同种子蛋白含量背景的早熟大豆品系，在2018年至2021年间，分别种植于一个东部地点（安大略省渥太华）和三个西部地点（曼尼托巴省布兰登、莫登及萨斯喀彻温省萨斯卡通），构建了一个覆盖多年份、多地点、多基因型的复杂研究体系。通过测量种子蛋白含量，并采集大豆鼓粒初期（R5期）的叶片进行RNA测序，研究团队获得了海量的表型和转录组数据。面对如此错综复杂的数据，传统的差异表达分析可能难以捕捉基因间协同工作的整体网络。因此，研究人员祭出了生物信息学分析的“利器”——加权基因共表达网络分析（WGCNA）及其拓展方法多WGCNA（multiWGCNA）。这些方法能够从全转录组中识别出共同表达、功能相关的基因模块，并能将这些模块与特定的样本性状（如产区“东/西”和蛋白含量“低/中/高”）进行关联，从而找出那些与“产区×蛋白含量”这一复合性状显著相关的关键基因共表达网络。
这项系统性的研究最终成功发表在了国际知名期刊《Genome》上。为了深入探索大豆籽粒蛋白环境调控的奥秘，研究人员主要运用了以下几项关键技术方法：首先，他们在加拿大东西部四个地点进行了为期四年（2018-2021）的田间试验，构建了包含10个不同蛋白含量大豆品系的多地点多年份样本队列。其次，他们利用近红外分析仪（Infratec 1241 Grain Analyzer）精确测定了种子蛋白含量。然后，研究团队对R5期叶片样本进行了RNA提取和高通量RNA测序（Illumina平台），并利用Wm82.a2.v1参考基因组进行比对和定量。最后，也是最为核心的分析，是应用加权基因共表达网络分析（WGCNA）和多WGCNA（multiWGCNA）对标准化后的转录组数据进行分析，以识别与“产区”和“蛋白含量”性状显著相关的差异共表达基因模块，并对这些模块进行基因本体（GO）富集分析和京都基因与基因组百科全书（KEGG）通路映射，以阐释其潜在的生物学功能。
研究结果
环境测量与种子蛋白测量
环境数据证实了东西部产区的气候差异：西部生长季平均温度比东部低2.8–6.4°C，平均相对湿度低1.3%–9.2%，总降水量也少于东部。种子蛋白测量结果与长期观察一致：在大多数年份和品系中，西部种植的大豆籽粒蛋白含量显著低于东部。直观展示了各品系在不同地点和年份的平均蛋白含量，西部（蓝色）普遍低于东部（黄色）。
PCA、尺度自由拓扑和平均连通性
对转录组数据的主成分分析（PCA）显示，东部和西部样本的基因表达谱存在明显分离，表明环境对大豆整体转录状态有深远影响。通过尺度自由拓扑分析筛选，2019年和2021年的数据集符合网络构建标准，被用于后续的WGCNA/multiWGCNA分析。
模块识别与差异模块表达
通过multiWGCNA分析，研究人员从2019年和2021年的网络中分别识别出44个和63个共表达基因模块。经过严格筛选，他们重点关注了两个与“蛋白×产区”性状显著相关的模块：2019-浅青色（lightcyan）模块和2021-浅黄色（lightyellow）模块。这两个模块的基因在不同产区及不同蛋白含量样本中的共表达模式存在显著差异。
2019年与蛋白×产区相关的模块
2019-浅青色模块包含74个基因，其表达在西部产区的中、低蛋白样本中上调，而在东部产区的同类样本中下调。展示了该模块的差异共表达热图、样本特征基因表达条形图及按性状分组表达箱线图。该模块的枢纽基因（hub gene）编码一个网状蛋白（reticulon），其他高连通性基因涉及蔗糖合酶、NADH–细胞色素b5还原酶和HVA22样蛋白等。功能分析表明，该模块富含应激反应（如抗坏血酸过氧化物酶、热激蛋白）、碳代谢（蔗糖合酶、苹果酸脱氢酶）、氮同化（天冬酰胺-tRNA连接酶）以及病原体反应（类奇甜蛋白、葡聚糖酶）相关基因。展示了该模块内基因的共表达关系网络。
2021年与蛋白×产区相关的模块
2021-浅黄色模块包含70个基因，其表达在东部产区（尤其是高蛋白样本）中上调，在西部产区则普遍下调。展示了该模块的分析结果。该模块的枢纽基因编码核糖体小亚基蛋白S12，模块内富含核糖体亚基基因（19个）、光合系统蛋白基因、染色质调节蛋白（如DEK结构域蛋白）以及组蛋白去甲基酶基因等。展示了其共表达网络。功能富集分析显示，该模块与染色质重塑、组蛋白结合、翻译以及光合作用等过程密切相关。
研究结论与讨论
综合以上结果，本研究通过WGCNA/multiWGCNA这一强大的生物信息学工具，成功解析了环境调控大豆籽粒蛋白质积累的复杂分子网络。研究发现的两个关键模块揭示了截然不同但又相辅相成的调控机制。
2019-浅青色模块指向了一种应激响应驱动的代谢重编程机制。西部产区相对严苛的环境（低温、少雨）可能激发了更强的生物与非生物胁迫响应。模块中HVA22样蛋白的上调可能通过抑制囊泡运输来延迟蛋白质储存液泡的充实，而应激相关的碳氮代谢基因（如蔗糖合酶、天冬酰胺-tRNA连接酶）的差异共表达，则暗示资源可能被重新分配以应对胁迫，而非全力进行种子蛋白质合成，这直接导致了西部大豆籽粒蛋白含量的降低。
2021-浅黄色模块则揭示了表观遗传与翻译水平上的精密调控。该模块在东部产区高蛋白样本中的高表达，关联到染色质重塑活性（如DEK蛋白、组蛋白修饰酶）和核糖体生物合成的上调。这表明东部适宜的环境可能通过表观遗传修饰“打开”有利于蛋白质合成和光合作用的基因网络，并通过增加翻译机器（核糖体）的丰度来高效生产蛋白质，从而支持更高的籽粒蛋白积累。此外，模块中光合作用相关基因的共表达也提示，东部环境下可能具有更强的光合能力，为蛋白质合成提供了更充足的碳骨架和能量。
重要意义 本研究的突破性在于，它超越了寻找单个差异表达基因的传统思路，从基因网络协同工作的视角，系统性地揭示了环境通过调控多个互相关联的生物学过程（胁迫响应、碳氮代谢、表观遗传、翻译效率等）来影响大豆籽粒蛋白质含量的整合机制。这不仅为理解“西低东高”这一长期农学现象提供了坚实的分子生物学解释，更为大豆育种家提供了宝贵的靶点信息。未来，通过分子标记辅助选择或基因编辑技术，针对这些关键模块中的枢纽基因（如网状蛋白、HVA22样蛋白、核糖体蛋白、DEK蛋白等）进行调控，有望培育出在不同环境条件下都能稳定保持高蛋白含量的大豆新品种，从而提升加拿大乃至全球大豆生产的品质与韧性。这项研究也展示了multiWGCNA在解析复杂农艺性状与环境互作方面的强大能力，为其他作物的类似研究提供了典范。
