整合转录组学与生理学分析揭示籼稻和粳稻中汞耐受性亚种偏好的调控机制

《Environmental and Experimental Botany》：Integrative Transcriptomic and Physiological Analysis Uncovers Subspecies-preferred Regulatory Mechanisms of Mercury Tolerance in Indica and Japonica Rice

　　

研究人员通过整合表型分析、生理测量和根转录组分析，筛选了具有对比汞（Hg）耐受性的籼稻和粳稻品种，以研究其对汞胁迫的亚种偏好性响应。本研究的目的是表征籼稻和粳稻幼苗在汞暴露下的生理响应和汞积累模式，在全基因组水平上鉴定其差异响应背后的关键代谢途径和转录调控网络，并确定与汞胁迫响应相关的核心候选基因和转录因子。这些发现为水稻对汞适应性的遗传多样性提供了新见解，并为低汞积累水稻品种的分子育种奠定了基础。

**论文解读：整合转录组与生理学分析揭示籼稻和粳稻汞耐受性亚种偏好调控机制**

**1. 研究背景与问题**

汞（Hg）是一种全球性重金属污染物，具有强神经毒性和生物放大效应。在水稻（*Oryza sativa* L.）种植中，长期淹水导致根际厌氧环境，促进无机汞向甲基汞（MeHg）的微生物转化，而稻米摄入是人群甲基汞暴露的主要途径，尤其在贵州等汞矿区，食用受污染稻米是当地居民汞暴露的主要风险，对孕妇及新生儿神经发育构成威胁。因此，阐明水稻对汞的吸收、转运及解毒机制对保障食品安全和公共健康至关重要。

植物对重金属胁迫已进化出多种防御策略，包括限制吸收、液泡区隔、螯合解毒及抗氧化清除。水稻的籼稻（*indica*）和粳稻（*japonica*）两个主要栽培亚种在形态、生理和胁迫耐受性上存在显著分化。已有研究报道了籼粳稻对镉（Cd）和铅（Pb）的差异积累与耐受性，但两者对汞胁迫响应的系统比较，尤其是全基因组转录水平上的比较，尚属空白。因此，开展此项研究旨在解析籼粳稻在汞胁迫下防御策略的分子差异。

**2. 主要技术方法**

研究人员从实验室保存的164份水稻种质中，通过20 μmol L^-1 HgCl₂胁迫7天的初步筛选，结合隶属函数分析和主成分分析（PCA），选取了耐汞品种（籼稻NH701、粳稻NHH）和敏感品种（籼稻N1311R、粳稻WN）。随后对幼苗进行20 μmol L^-1 HgCl₂处理7天，采用以下关键方法：冷原子荧光光谱法（CV-AFS）测定总汞（THg）含量；分光光度法测定丙二醛（MDA）含量、超氧化物歧化酶（SOD）、过氧化物酶（POD）、过氧化氢酶（CAT）活性、还原型谷胱甘肽（GSH）含量及谷胱甘肽S-转移酶（GST）活性；利用Illumina NovaSeq 6000平台对根样本进行转录组测序，通过DESeq2鉴定差异表达基因（DEGs），并利用加权基因共表达网络分析（WGCNA）识别关键模块，所有样本均来自上述处理的三个生物学重复。

**3. 研究结果**

**3.1. 汞胁迫对水稻幼苗生长的影响**
汞胁迫显著抑制所有品种的生长，但耐性品种（NH701和NHH）的株高、根长和生物量减少幅度小于敏感品种（N1311R和WN），其中耐性籼稻NH701的抑制程度最小。

**3.2. 汞积累与分配的亚种偏好模式**
耐汞品种的根汞含量最高，其中NH701的根汞浓度比敏感籼稻N1311R高30.57%，而其地上部汞浓度仅为N1311R的54.89%。NH701的转运系数（TF）最低，显著低于其他品种，表明其具有最强的根汞滞留和限制向上转运的能力。

**3.3. 品种特异性抗氧化响应**
汞胁迫诱导所有品种根系MDA含量上升，但敏感品种增幅更大。POD活性在粳稻中激活更显著，CAT和SOD活性在敏感品种中补偿性上调更明显，显示粳稻偏好过氧化物酶介导的抗氧化防御。

**3.4. 汞胁迫下的品种特异性转录重编程**
每个品种均鉴定出大量DEGs，其中NH701最多（6,234个）。共有1,495个DEGs为四个品种共享，同时存在大量品种特异性DEGs，表明转录响应具有品种和亚种偏好性。

**3.5. 籼粳稻在汞胁迫下保守与分化的功能途径**
所有品种均富集“氧化应激响应”、“谷胱甘肽代谢”等途径，但耐性品种更富集抗氧化防御，敏感品种更富集初级代谢。亚种偏好方面，籼稻显著富集“谷胱甘肽代谢过程”、“茉莉酸代谢过程”；粳稻显著富集“热响应”、“海藻糖代谢过程”和“木质素代谢过程”。KEGG分析显示，“ABC转运蛋白”途径在籼稻耐性品种NH701中更富集，“黄酮类生物合成”在粳稻耐性品种NHH中更富集。

**3.6. 籼稻中硫代谢与谷胱甘肽合成的协同激活增强汞螯合解毒**
所有品种中“谷胱甘肽代谢”途径均被激活，但籼稻中硫酸盐吸收和同化基因（如*OsSULTR3;6*、*OsATPS*、*OsOAS-TL*）及谷胱甘肽合成相关基因（*OsGSTU4*、*OsGR2*）协同上调更显著，GSH含量和GST活性增幅也更大。这为籼稻构建了从硫酸盐到半胱氨酸再到谷胱甘肽的完整解毒代谢流。

**3.7. 粳稻激活苯丙烷生物合成途径与木质素合成以构建物理屏障**
苯丙烷生物合成途径在粳稻中更强烈激活，尤其是木质素合成分支基因（*OsHCT4*、*OsCCR10*、*OsCCR18*）和过氧化物酶（PRX）家族基因在粳稻中上调幅度显著高于籼稻。木质素含量测定证实粳稻在汞胁迫下木质化程度更高，表明其偏好通过细胞壁强化和过氧化物酶介导的活性氧（ROS）清除来构建物理屏障。

**3.8. ABC转运蛋白家族与ZIP金属转运蛋白的差异调控**
ABC转运蛋白基因（尤其是*OsABCC1*、*OsABCG36*）在籼稻中诱导更强烈，与籼稻较低的TF和较高的根汞浓度一致。而ZIP家族基因（如*OsZIP5*、*OsZIP9*）在籼稻中抑制更广泛，可能限制汞的流入，但敏感品种N1311R的更强抑制反映了其更严重的胞内汞毒性。

**3.9. 通过WGCNA识别潜在的关键转录因子**
WGCNA分析鉴定出MEturquoise模块与汞胁迫高度正相关，且富集于谷胱甘肽代谢和苯丙烷生物合成途径。该模块内*OsWRKY24*、*OsNAC4*、*OsHSFA4D*等转录因子在耐汞品种中强烈诱导，其中*OsWRKY24*与谷胱甘肽代谢和苯丙烷途径关键基因高度共表达，可能作为调控枢纽。

**3.10. RT-qPCR验证**
通过RT-qPCR对RNA-seq结果中关键基因的表达模式进行了验证，两者趋势一致。

**4. 讨论与结论**

**讨论**：研究人员指出，耐汞籼稻NH701通过硫代谢-谷胱甘肽合成-ABC转运蛋白介导的液泡区隔这一“螯合-区隔”协同模块，实现了高效的根汞滞留和低转运系数。而粳稻则偏好通过细胞壁木质化（木质素沉积）和过氧化物酶（PRX）介导的抗氧化防御，形成“物理屏障”策略，但该策略在限制汞向上转运方面效率较低。*OsWRKY24*作为潜在调控枢纽，可能通过双机制（激活谷胱甘肽合成酶和苯丙烷途径关键酶）协调防御响应。两亚种在汞分配策略上的差异，为进一步培育低汞积累水稻品种提供了靶基因（如*OsGSTU4*、*OsABCC1*）。

**结论翻译**：在本研究中，将籼稻（NH701、N1311R）和粳稻（NHH、WN）幼苗暴露于20 μmol L^-1 HgCl₂处理7天。汞胁迫显著抑制所有品种的生长，但耐性基因型比敏感基因型表现出更少的生物量减少和更强的根汞滞留能力，尤其是耐性籼稻品种。转录组分析揭示，两个亚种的差异表达基因主要富集在保守途径中，包括“谷胱甘肽代谢”、“苯丙烷生物合成”和“氨基酸生物合成”。尽管存在这种共享富集，但两个亚种表现出了不同防御模块的优先激活。在粳稻基因型（NHH和WN）中，上调最显著的基因主要与抗氧化系统和木质素合成相关；这种转录偏好通过生化测量得到证实，即与籼稻基因型（NH701和N1311R）相比，粳稻中汞诱导的木质素积累显著更高。在籼稻中，硫代谢和谷胱甘肽合成相关基因被优先激活，同时伴随ABC转运蛋白的强诱导和ZIP家族基因的广泛抑制。生理测量进一步证实了这些转录模式：粳稻品种表现出过氧化物酶（POD）和过氧化氢酶（CAT）活性的更大增加，而籼稻品种则表现出谷胱甘肽（GSH）含量和谷胱甘肽S-转移酶（GST）活性的更高提升，与各亚种观察到的转录偏好一致。相应地，在测试的四个品种中，耐汞籼稻品种NH701表现出最高的根汞浓度和最低的转运系数，而粳稻品种表现出中等至较高的TF值。总体上，这些发现表明，在汞胁迫下，粳稻主要通过激活抗氧化酶和加强细胞壁木质化来减轻氧化损伤。相比之下，籼稻通过协调的硫同化、谷胱甘肽螯合、ABC转运蛋白介导的液泡区隔以及限制ZIP介导的汞流入，实现了高效的根汞滞留。防御策略的早期分化构成了转录和生理水平上的亚种偏好响应模式，而整体耐受水平主要由基因型特异性螯合-区隔网络的协调效率决定，如耐汞籼稻品种NH701所体现的。