粳稻保守盐胁迫共表达网络揭示枢纽基因连接性与转录诱导的部分解偶联

《Frontiers in Plant Science》：Aconserved salt-stress co-expression network in japonica rice reveals partially dissociated hub gene connectivity and transcriptional induction

　　

摘要

引言：土壤盐分通过破坏离子稳态和渗透稳态严重损害水稻（Oryza sativa L.）的生长和籽粒产量。尽管转录组分析已推进了对水稻盐胁迫响应的理解，但共表达网络在不同遗传背景间的保守性以及网络连接性与转录诱导之间的关系仍缺乏充分表征。

方法：研究人员对三个粳稻品种在对照（0 mM）、中度（85.6 mM）和重度（136.9 mM）NaCl胁迫下进行了整合的RNA测序和加权基因共表达网络分析（WGCNA）。

结果：所有三个品种均表现出基础盐耐受性。差异表达分析在中度和重度胁迫下分别鉴定出2,067和4,900个非冗余差异表达基因（DEGs），其中分别有205和1,093个DEGs在三个品种间共享。WGCNA鉴定出32个共表达模块。棕色模块表现出对盐胁迫最强的剂量依赖性响应，并富集于苯丙烷生物合成、植物激素信号转导和MAPK信号通路；其特征基因（eigengene）在各品种间高度一致（平均r = 0.987），与共享的叶片水平转录程序一致。在棕色模块基因中，模块内连接性与盐诱导的倍数变化呈正相关但仅中等相关（Pearson r = 0.377）：OsDREB1A是连接性最高的基因，但并非诱导最强的基因，而LOC_Os10g35080是诱导最强的基因，但其连接性排名第14位。在15个节点的高渗盐度响应子网络中，没有单个枢纽基因占主导地位，且OsDREB1A处于外围位置。

讨论：连接性与诱导性因此反映了盐胁迫响应的互补方面，用于分子育种的候选枢纽基因应通过这两个标准联合优先排序，并在不同遗传背景中进行验证。

粳稻保守盐胁迫共表达网络：枢纽基因连接性与转录诱导的部分解偶联

一、研究背景与问题

土壤盐渍化是限制水稻（Oryza sativa L.）生产的重要非生物胁迫因素，盐胁迫通过破坏离子稳态和渗透稳态抑制植株生长并降低籽粒产量。尽管RNA测序（RNA-Seq）技术已广泛应用于水稻盐胁迫转录组研究，但以往工作多集中于单个基因或通路，对共表达网络在不同遗传背景中的保守性，以及网络连接性与转录诱导之间的关系仍缺乏系统认识。因此，研究人员设计了包含三个粳稻品种、三种NaCl浓度的比较转录组实验，以期鉴定跨品种保守的盐响应共表达模块，解析其功能通路，并评估枢纽基因连接性与转录诱导的关联。

二、研究内容与意义

研究人员以三个具有基础盐耐受性的粳稻品种WYJ 31、LJ 11和XN 8333为材料，在三叶期分别进行0、85.6和136.9 mM NaCl处理10天，采集叶片进行转录组测序。通过差异表达分析和加权基因共表达网络分析（WGCNA），研究人员鉴定出32个共表达模块，其中棕色模块对盐胁迫表现出最强的剂量依赖性响应，且其模块特征基因（eigengene，即模块表达矩阵的第一主成分）在三个品种间高度一致（平均r = 0.987）。进一步分析显示，模块内连接性与盐诱导倍数变化仅呈中等程度正相关（Pearson r = 0.377），提示二者反映盐胁迫响应的互补方面。该研究为通过联合连接性和诱导性两个标准筛选分子育种候选枢纽基因提供了依据。

三、主要技术方法

本研究利用三个粳稻品种（WYJ 31、LJ 11、XN 8333）的水培幼苗，在三叶期进行0、85.6、136.9 mM NaCl处理10天，每处理3个生物学重复，共27个叶片样本，由BGI Genomics Co., Ltd.完成Illumina NovaSeq 6000双端测序。数据经fastp质控、HISAT2比对到水稻MSU7参考基因组、StringTie定量。使用DESeq2鉴定差异表达基因（DEGs），clusterProfiler进行GO（基因本体）和KEGG（京都基因与基因组百科全书）通路富集分析。WGCNA以软阈值β=6构建共表达网络，动态树切割鉴定模块，通过模块特征基因相关性评估跨品种一致性。

四、研究结果

3.1 盐胁迫下的表型响应

通过株高和叶片数等表型测量，发现三个品种在S5和S8处理下均出现生长抑制，但未死亡，表明它们具备基础盐耐受性。WYJ 31在胁迫下保留较多叶片和较高株高，LJ 11在S8下株高降幅最大，但品种间差异有限。

3.2 RNA-Seq数据质量评估

27个样本共获得182.49 Gb clean data，平均每样本6.76 Gb，Q30比例均高于86.84%，比对率介于91.56%–95.44%。主成分分析（PCA）显示，处理梯度沿PC1排列，品种间存在转录组差异，但所有品种的盐处理样本在PC2上共同低于对照，提示存在共享的胁迫响应。

3.3 盐胁迫下差异表达基因的转录组分析

各品种DEGs数量在S5下为1,046–3,523个，S8下为3,116–7,283个；三个品种共享的DEGs在S5和S8下分别为205和1,093个。GO富集显示共享DEGs主要涉及非生物刺激响应、水分剥夺响应等；S8特有富集涉及离子转运、光合碳固定、氨基酸和谷胱甘肽代谢等，表明重度胁迫激活额外机制。

3.4 共表达网络构建与模块鉴定

基于27个样本的WGCNA共鉴定出32个模块（87–1,458个基因）。棕色、绿色和蓝色模块的特征基因（eigengene）值随盐浓度升高而增加，其中棕色模块响应最强。棕色模块的特征基因值在品种间具有高度一致性（平均r=0.987），表明该模块代表保守的叶片水平转录程序。

3.5 棕色模块功能富集

GO富集显示棕色模块最显著富集“响应水分剥夺”（GO:0009414）。KEGG富集到16条通路，其中植物激素信号转导（ko04075）、MAPK信号通路（ko04016）、苯丙烷生物合成（ko00940）等达到高度显著，涉及ABA、生长素、乙烯等激素分支及Ca²⁺信号组分，提示该模块整合渗透信号、结构防御和能量代谢。

3.6 枢纽基因鉴定与子网络分析

按模块内连接性（kIM）排名前20的枢纽基因均受S8上调，但连接性与诱导倍数变化仅中等相关（r=0.377）。OsDREB1A（LOC_Os06g45184）是连接性最高的基因，但诱导倍数较低；LOC_Os10g35080（功能未知）诱导倍数最高（203.2倍），连接性仅排第14。在由15个节点组成的GO:0042538高渗盐度响应子网络中，没有单一枢纽主导，OsDREB1A处于外围位置，说明该功能子网络由多基因协同调控。

五、讨论总结

讨论部分指出，棕色模块富集的通路支持其作为渗透信号、结构防御和能量代谢的整合节点。连接性与诱导性仅部分偶联，因此候选枢纽基因应联合两个标准筛选。棕色模块的跨品种保守性提示存在共享的叶片水平盐胁迫程序，但品种间沿PC1的分离表明遗传背景影响转录组景观。研究局限性包括仅分析了单一发育阶段、叶片组织和一个时间点，且表型指标有限。研究人员提出工作模型：盐胁迫信号经MAPK和ABA通路传递，激活转录因子和下游效应子，促进渗透调节、苯丙烷类细胞壁加固及糖酵解供能。

六、研究结论

本研究利用RNA-Seq和WGCNA表征了三个粳稻品种的盐胁迫响应。尽管品种间存在明显的转录组差异，棕色模块在所有品种中保守，并对逐渐增强的盐胁迫呈剂量依赖性响应。该模块可能通过ABA介导的渗透信号、苯丙烷类细胞壁加固和糖酵解ATP供给来协调适应。包括OsDREB1A、OsAP37、OsSRO1c和LOC_Os10g35080在内的枢纽基因是功能验证的