《Genome》：Integration of consensus WGCNA and differential expression analysis reveals stable heat shock protein networks in soybean adaptation to Western Canada
编辑推荐：
大豆（Glycine max (L.) Merr.）向加拿大西部种植区的扩展受到独特非生物胁迫因素的制约，包括较短生长季节和干旱。热休克蛋白（heat shock proteins，HSPs）是胁迫适应的关键调控因子；然而，在大田条件下表征其表达常因显著的年度
大豆（Glycine max (L.) Merr.）向加拿大西部种植区的扩展受到独特非生物胁迫因素的制约，包括较短生长季节和干旱。热休克蛋白（heat shock proteins，HSPs）是胁迫适应的关键调控因子；然而，在大田条件下表征其表达常因显著的年度环境变异性而受到干扰。在本研究中，研究人员整合了差异表达分析、一致性加权基因共表达网络分析（consensus Weighted Gene Co-expression Network Analysis，WGCNA）以及严格的蛋白质-蛋白质相互作用组建模，利用了2018–2021年共4年的田间数据，涉及在加拿大东部和西部种植的10个基因型。通过过滤时间噪声，研究人员鉴定出一个稳定的、区域特异性模块，经计算预测该模块可解析为一个高度协调的物理伴侣复合体。对该网络的亚细胞表征揭示了一个靶向线粒体和叶绿体的稳健伴侣核心。关键的是，高置信度预测鉴定出线粒体HSP60枢纽蛋白与一个核结构域重排甲基转移酶之间的物理相互作用，提示该伴侣复合体可能将局部细胞器防御与系统性表观遗传调控联系起来，以保障代谢能力。这些发现刻画了一个经计算建模的、依赖于地理位置的胁迫应答机器，并为旨在提高大豆在高纬度地区适应性的育种项目提供了理想的遗传靶标。
大豆（Glycine max (L.) Merr.）在加拿大农业中的地位日益重要，其高蛋白和高油含量的籽粒广泛应用于饲料、生物燃料和人类消费。将大豆种植向加拿大西部扩展是加拿大育种者的重要经济目标，但西部地区更长的夏季日照、较短的生长季节以及较少的降水量，使得现有品种在该区域常表现为产量和籽粒品质下降。先前研究已证实，西部大豆的籽粒蛋白含量和油含量以及产量均显著低于东部地区，而由地理差异引起的非生物胁迫是造成这一品质差距的主要因素。因此，解析这些环境下调控胁迫响应的遗传机制，对培育适应西部地区的品种至关重要。
热休克蛋白（heat shock protein，HSP）是一类负责协调其他蛋白质折叠的多样蛋白家族，在非生物胁迫下尤为关键。HSP家族按分子量从HSP10到HSP100分为多种类型，并广泛参与大豆的胁迫应答。例如，HSP20参与胁迫级联反应中的蛋白聚集，也涉及大豆胞囊线虫防御；HSP40与干旱、盐害和重金属胁迫下的稳态维持相关，并与开花和荚发育有关；HSP70通过与HSP40直接互作，在蛋白质折叠网络中发挥多重作用。然而，由于非生物胁迫和植物发育过程的遗传基础复杂，加之田间环境中降水、温度等变量在年份和地点间显著波动，稳定遗传信号常被掩盖，识别驱动西部大豆区域适应性的特定HSP网络十分困难。标准差异表达（differential expression，DE）分析捕捉到的往往是年度暂时性波动，难以分离多年间持续存在的潜在生物学网络。为克服这种年份效应，需要采用系统生物学方法过滤时间噪声，识别在不同环境中稳定保留的共表达模式。基于此，研究人员整合差异表达分析与一致性加权基因共表达网络分析（consensus Weighted Gene Co-expression Network Analysis，WGCNA），利用4年田间数据，旨在过滤混淆性年际变异，鉴定稳定且区域特异性的HSP共表达模块，为适应加拿大西部市场的大豆育种提供遗传资源。
该研究使用了2018–2021年连续4年田间数据，涉及10个大豆基因型，种植于加拿大东部（渥太华）和西部（莫登、布兰登、萨斯卡通）共4个地点。叶片样品于第5生殖期（R5）采集，每品系3次重复。RNA经提取和质检后，采用Illumina HiSeq 4000平台进行双端测序（2×150 bp）；原始读段经GenPipes流程处理，比对到Glycine max Wm82.a2.v1参考基因组。差异表达分析采用DESeq2中的似然比检验（likelihood ratio test，LRT），以年份为控制变量；一致性网络利用WGCNA包基于方差稳定化计数矩阵构建，软阈值设为14，以各年份为独立输入集；蛋白质-蛋白质相互作用（protein-protein interaction，PPI）预测采用PIPE4引擎，并以拟南芥数据集为训练集，后验概率阈值设为99%；亚细胞定位预测综合使用WoLF PSORT和DeepLoc-2.0工具。
**RNA sequencing quality control（RNA测序质量控制）**：测序和质控结果显示，平均每个文库获得29,313,314条读段，绝大多数碱基的Phred质量分数不低于36，质控后平均保留98.0%的读段，最低单样本保留率为90.9%，确认数据集具有高质量。
**Regional agronomic phenotypes（区域农艺表型）**：通过线性混合效应模型分析4年表型数据，发现西部大豆产量（东部均值2653.36 kg/ha vs西部2160.94 kg/ha，p=5.60E-8）、籽粒蛋白含量（43.63% vs40.88%，p=3.37E-50）和油含量（20.73% vs19.97%，p=7.84E-14）均显著低于东部，证实西部种植条件诱导了可量化的胁迫效应。
**HSP family analysis（HSP家族分析）**：基于多源注释在基因组中鉴定出539个HSP基因，其中329个在东部和西部之间显著差异表达。各HSP家族在差异表达群体中的相对比例与全基因组中比例总体相似。
**Differential expression analysis（差异表达分析）**：全转录组水平上，地理效应共产生27,184个差异表达基因，而年度时间效应产生42,245个差异表达基因，表明年际环境波动是巨大的变异来源。在329个区域差异表达HSP中，42个显著下调，17个显著上调；主成分分析显示第一主成分（39%）对应年份的连续变化，说明即便经统计控制，年效应仍较大，需通过共识网络进一步过滤。
**Weighted gene correlation network analysis（加权基因共表达网络分析）**：基于4年数据构建的共识网络共包含36,313个基因，鉴定出161个共表达模块。模块特征基因与农艺性状的相关性分析未发现显著关联，反映出最终收获性状在年际间波动较大。
**HSP regional co-expression networks（HSP区域共表达网络）**：对161个模块进行富集分析，仅模块“indianred4”显著富集区域差异表达HSP基因（校正p<0.01）。该模块含92个基因，其中23个为区域差异表达HSP，其余69个为共伴侣和其他蛋白折叠相关基因。基于拓扑重叠矩阵的共表达网络可视化显示该模块内部高度协调。
**Module analyses（模块分析）**：PPI预测显示模块“indianred4”内部存在67条高置信度预测物理互作，并解析为三个互不重叠的亚复合体。亚细胞定位预测发现该模块包括靶向线粒体的两个HSP70和一个HSP60，以及靶向叶绿体的HSP20；其中HSP20（Glyma.02G052200）同时具有线粒体（0.78）和质体（0.94）双靶向信号。此外，高置信度预测识别出线粒体HSP60（Glyma.12G013100）与核结构域重排甲基转移酶（DRM1同源物，Glyma.02G035700）之间的物理互作，提示该伴侣复合体可能将局部细胞器防御与系统性表观遗传调控相连接。
讨论部分指出，共识WGCNA通过仅保留在所有年份输入网络中方向一致的共表达关系，有效滤除了年际环境噪声，是处理多年田间转录组数据的关键步骤。模块“indianred4”不仅富集差异表达HSP基因，还包含多种辅助伴侣，如含有四肽重复（tetratricopeptide repeat，TPR）结构域的HSP70-HSP90组织蛋白（HOP3）同源物、热诱导肽基脯氨酰异构酶（ROF1）同源物，以及线粒体GrpE共伴侣基因。该模块的top10核心基因包括HSP20、HSP70、HSP90、HSP100以及硫醇肽酶等共伴侣，表明多种HSP协同参与西部环境胁迫应答。亚细胞定位结果支持模块优先保护能量转导细胞器（线粒体和叶绿体），以维持R5期源器官的代谢功能。模块特征基因与最终收获性状缺乏显著相关，可能是由于共识网络统计极为保守，且R5期取样与最终成熟（R8）之间存在后期环境波动。特别值得注意的是，线粒体HSP60与DRM甲基转移酶之间的预测互作暗示一个跨区室机制：线粒体伴侣可能保护核表观遗传修饰酶，从而维持能量细胞器相关基因的正常表达。该假说仍需酵母双杂交、双分子荧光互补或免疫共沉淀等体内空间验证，但为理解胁迫应答中局部保护与全局调控的整合提供了新视角。
研究结论为：该研究表明，当标准差异表达分析应用于多年田间数据时，易受年际环境变异的混淆。通过将似然比检验与共识WGCNA整合，并利用严格的物理互作组模型进行筛选，研究人员成功过滤时间噪声，鉴定出模块“indianred4”——一个稳定的、由区域依赖性热休克蛋白组成的物理互作复合体。该模块经计算建模为一个独特的伴侣机器，并特异性地保护能量转导细胞器。关键在于，线粒体HSP60枢纽与DRM甲基转移酶之间的互作预测提示，该复合体可能桥接局部细胞器防御与系统性表观遗传调控，从而保障代谢能力。因此，这一依赖地理位置的胁迫应答谱系突出了伴侣复合体在环境胁迫下维持植物发育的结构性作用，为大豆在加拿大西部扩展背景下的育种提供了理想的候选基因。
