分根转录组学揭示水稻基因型依赖性根系硼滞留的分子基础

《Environmental and Experimental Botany》：Split-root transcriptomics reveals the molecular basis of genotype-dependent root boron retention in rice

　　

局部硼胁迫（localized boron stress）触发水稻根系快速且局部的系统响应。硼转运通过BOR1、PLT5和DTX16的调控被重塑。基于谷胱甘肽（glutathione, GSH）和过氧化氢酶（catalase, CAT）的防御被快速激活并持续。分根转录组学（split-root transcriptomics）将外部硼响应（external-boron responsive, EBR）与内部硼响应（internal-boron responsive, IBR）模块分离。更广泛的通路抑制是根系硼滞留增强的基础。

**水稻根系硼毒害早期响应及基因型依赖性硼滞留的分子机制解析**

硼是植物生长必需的微量元素，参与细胞壁完整性、膜稳定性和碳水化合物代谢；然而，当土壤中硼浓度超出适宜范围（通常0.5–2.0 mg kg^?1）时，便会对植物产生毒害。硼污染来源包括矿山活动、含硼灌溉水等自然和人为因素，在中国部分矿区土壤硼含量可超过1000 mg kg^?1，严重威胁农业可持续性和粮食安全。过量硼会破坏根系结构、抑制光合作用、阻碍养分吸收并诱发氧化胁迫。尽管土壤改良和培育耐性品种是应对硼毒害的重要策略，但以往研究多集中于长期、均匀硼处理，忽略了两个关键问题：一是胁迫感知后最初数小时内发生的早期系统分子事件，二是根系内部对胁迫感知的空间异质性及其导致的基因型特异性适应策略。水稻（*Oryza sativa* L.）是重要的硼敏感作物，但其硼毒害分子机制尚未充分阐明。为此，研究人员选用两个硼积累能力存在显著差异的水稻品种Bei_Jing_1702（BJ）和Yan_Jing_Ran_Yu_22（YR），利用分根体系（split-root system）结合转录组学，探究局部硼胁迫下12小时内的早期分子响应，区分外部硼感知和内部硼状态相关的转录模块，并揭示两品种根系硼积累差异的分子基础。该研究发表于《Environmental and Experimental Botany》。

研究人员主要采用了以下关键技术方法：首先，对15日龄水稻幼苗进行不同浓度硼处理（0.2–100 mg L^?1），结合形态和生理指标确定40 mg L^?1为后续胁迫浓度；利用电感耦合等离子体质谱（ICP-MS）测定根系和地上部硼浓度，并对分根系统中不同 compartments 的硼分配进行定量。核心方法是转录组测序（RNA-seq），对分根处理后12小时的根系组织构建cDNA文库并进行Illumina测序，利用edgeR鉴定差异表达基因（DEGs），结合基因本体（GO）和KEGG通路富集分析、STRING数据库蛋白互作网络（PPI）以及加权基因共表达网络分析（WGCNA），并采用qRT-PCR验证关键基因表达。样本来源为沈阳农业大学水稻研究所提供的两个水稻品种。

**结果**
**3.1 不同硼浓度下的形态响应**
通过设置0.2–100 mg L^?1硼浓度处理15天，研究发现高硼显著抑制两品种的株高和生物量，且根系比地上部更敏感。从20 mg L^?1开始，总根长、表面积、体积、平均直径和根尖数均显著下降。在40 mg L^?1处理下，YR的根鲜重和根系活跃吸收面积显著高于BJ，表明YR在硼胁迫下能更好地维持根系生理状态。同时，BJ根系硼浓度和积累量始终高于YR，表明两品种存在固有的硼吸收与滞留差异。

**3.2 硼吸收动力学与系统分配**
12小时吸收动力学显示，YR在前6小时具有更高的初始吸收速率，但其最终根系硼浓度反而显著低于BJ。分根实验表明，经硼处理的半边根系中两品种硼浓度相似，但BJ未处理根系中的硼浓度（233 mg kg^?1）几乎是YR（135 mg kg^?1）的两倍；相反，YR向地上部转运的硼更多（301 mg kg^?1 vs. 196 mg kg^?1），且转运系数显著高于BJ（1.1 vs. 0.7）。这说明YR倾向于将硼向地上部转运，而BJ倾向于在根中滞留。

**3.3 硼胁迫下的早期转录组重编程**
RNA-seq共鉴定到1948个DEGs。YR在硼处理下的转录响应更强烈（492个DEGs），而BJ较少。GO和KEGG富集分析显示，最显著富集的GO条目为苯丙烷代谢过程、细胞壁生物发生和谷胱甘肽转移酶活性；KEGG通路主要为谷胱甘肽代谢和苯丙烷生物合成。这表明早期响应涉及细胞壁重塑和抗氧化防御。

**3.4 细胞壁生物合成与组装的抑制**
12小时硼处理下，39个细胞壁相关DEGs中33个（85%）显著下调，包括纤维素合酶（*CESA7*）、半纤维素合酶（*CSLs*）、木质素聚合相关基因（*LACs*, *PERs*）等。PPI网络显示*CESA7*、*GT32*和*LAC12*为中心节点。然而，15天连续胁迫后，纤维素、果胶和木质素含量分别增加16.8%、36.6%和21.0%，半纤维素下降4.54%。这说明早期转录抑制与后期生化积累并非简单对应，细胞壁代谢受时间动态调控。

**3.5 细胞周期阻滞与DNA损伤/修复诱导**
转录组分析发现，G2/有丝分裂特异性细胞周期蛋白（*CCB11*, *CCA12*）、微管去稳定蛋白（*MDP60*）、动粒蛋白（*NUF2*）等细胞周期进程核心基因均下调；DNA复制相关基因如核糖核苷二磷酸还原酶小链（*RIR2*）、起点识别复合物亚基（*ORC3*）和复制蛋白A亚基（*RFA1B*）也显著抑制。相反，9个DNA/RNA损伤修复基因上调。PPI网络识别出以*CCB22*、*CCB11*和*RIR2*为核心模块。结果表明硼毒害在早期即抑制细胞分裂并激活DNA损伤修复。

**3.6 植物激素通路重编程**
硼胁迫显著改变激素稳态：赤霉素（GA）生物合成（*CPS4*, *KOL5*）和信号转导（*CXE17*, *GSR1*）被抑制；生长素（IAA）生物合成（*YUCCA3*）、激活（*ILL9*）和极性运输（*PILS1*, *PIN5c*）上调；脱落酸（ABA）和茉莉酸（JA）合成及信号基因（*NCED5*, *AOS3*, *OPR1*等）激活。水杨酸（SA）、乙烯（ET）和油菜素甾醇（BR）通路也发生复杂变化。PPI网络显示GA合成*CPS4*和SA/JA响应*MPAG1*为关键节点。这表明激素重编程驱动根系从生长导向转向防御导向状态。

**3.7 谷胱甘肽和过氧化氢酶诱导**
12小时硼处理后，多个谷胱甘肽S-转移酶（GST）基因（*GSTU1*, *GST23*）和过氧化氢酶基因*CATA3*显著上调，而编码过氧化物酶（POD）和超氧化物歧化酶（SOD）的基因下调。PPI分析显示*REHYB*为中心枢纽。15天处理后，CAT和SOD活性分别增加37.91%和64.9%，GSH含量升高2.19倍，POD活性下降。这表明CAT和GSH介导的活性氧（ROS）清除是硼毒害下的快速且持续的核心抗氧化系统。

**3.8 转运蛋白重编程**
硼转运系统被重塑：外排转运蛋白基因*BOR1*下调，解毒转运蛋白*DTX16*和多元醇转运蛋白*PLT5*上调。同时，其他养分转运蛋白也发生变化，如铁（*IRT2*）和硫（*SUT12*）转运蛋白下调，磷（*PHO12*）和氮/肽（*PTR53*）转运蛋白上调。这种复杂重校准表明植物在硼胁迫下主动限制硼进入并启动解毒机制。

**3.9 空间转录组比较**
利用分根系统将DEGs分为外部硼响应（EBR）基因集和内部硼响应（IBR）基因集。EBR来自BJ品种硼处理与对照根系的比较（两者内部硼浓度相近），主要反映对外部硼刺激的直接响应；IBR来自两品种未处理根系的比较（内部硼浓度差异显著），反映内源硼稳态差异。EBR模块以快速损伤控制和限制进入为特征，包括细胞增殖抑制和DNA/RNA修复；IBR模块以主动解毒和代谢调整为特征，包括隔离转运蛋白*PLT5*上调以及营养代谢和胁迫蛋白诱导。这揭示了根系对局部和系统信号的功能特化。

**3.10 加权基因共表达网络分析（WGCNA）**
WGCNA显示，硼浓度显著影响共表达网络结构。在YRB vs. YRCK中，正负相关模块分别含2681和1438个基因；在BJ品种相应比较中模块较小。而直接比较BJB与YRB（两者硼浓度相近）时未检出硼相关模块，表明硼浓度是驱动网络差异的主要因素。进一步比较YRB vs. YRCK和BJCK vs. YRCK（两者硼浓度差相似）发现，正相关模块的KEGG富集谱相似，但负相关模块差异显著：BJCK vs. YRCK中负相关模块有71个DEGs分布于24条富集通路，而YRB vs. YRCK中仅3个DEGs。BJ的负相关模块富集到MAPK信号通路、植物激素信号转导、脂肪酸延长、氨基酸代谢和苯丙烷生物合成等，且表达下调幅度更大，表明BJ在硼积累下经历更广泛的通路抑制。

**讨论与结论**
讨论部分总结了硼毒害下水稻的早期硼分配策略、毒害机制（细胞壁抑制、细胞周期阻滞、激素失衡）、适应机制（硼转运系统重塑、GST-CAT抗氧化系统）、空间转录响应以及基因型依赖的根系硼滞留。通过综合转录组和生理数据，提出模型：YR因维持更广泛的应激信号和代谢通路活性，向地上部转运硼的能力较强，根系滞留较少；而BJ因MAPK、激素信号、脂质及氨基酸代谢和苯丙烷通路被广泛抑制，根系生理活性下降，导致硼滞留增多。

本研究结论表明，水稻对硼毒害的适应涉及局部胁迫感知与系统转录协调的整合响应。分根实验证实，局部硼胁迫快速抑制细胞壁生物合成、细胞周期进程和激素平衡，同时激活抗氧化防御、DNA/RNA损伤修复以及转运蛋白重编程。两品种在硼分配策略上显著不同：BJ倾向于根系滞留，YR倾向于向地上部转运。这种差异与转录状态相关，BJ表现出MAPK信号、激素信号、氨基酸和脂质代谢以及苯丙烷通路的广泛抑制，反映其根系生理活性降低；而YR能更有效地维持这些通路，与其较低的根系硼滞留和较高的地上部转运一致。这些发现为理解水稻基因型依赖性硼耐性提供了机制框架，并为在硼污染土壤上改良作物硼转运效率和胁迫适应能力提供了候选通路。