《Frontiers in Plant Science》：Crosstalk between glutathione and melatonin in chromium detoxification in sweet potato revealed by high-throughput sequencing and physio-biochemical profiling
编辑推荐：
本研究综述了外源谷胱甘肽（GSH）与褪黑素（MT）缓解甘薯铬（Cr）毒性的分子机制。通过转录组学结合生理生化分析，揭示了两种分子在转录重编程、通路调控（如MAPK信号、谷胱甘肽代谢、植物激素信号转导）及关键靶点（如GA20OX1、HEMA1、TIP2-1）激活方面的差异，阐明了GSH如何通过强化以氧化还原为核心的防御网络，更有效地增强作物对铬污染的耐受性，为作物育种和基因工程提供了重要的分子靶点。
引言：铬污染与作物胁迫的挑战
六价铬（Cr(VI)）是一种广泛存在的环境污染物，对植物生产力、光合作用和氧化还原平衡造成严重破坏。甘薯作为重要的全球粮食作物，其根茎易从土壤中吸收重金属，面临铬胁迫导致的生长抑制、光合效率下降和氧化损伤。为寻找经济可持续的生物缓解策略，外源应用具有生物活性的小分子，如谷胱甘肽（GSH）和褪黑素（MT），已成为增强植物重金属耐受性的有效途径。先前生理生化研究已表明，GSH在改善甘薯生长、降低铬积累和氧化损伤方面优于MT，但其背后差异化的分子机制尚不明确。本研究旨在通过整合转录组学与生物信息学分析，系统比较并阐明GSH和MT在甘薯铬胁迫应答中的分子对话。
材料与方法：从实验设计到数据分析的系统框架
本研究采用甘薯品种‘海大HD7791’，在严格的水培条件下进行。实验设置四个处理组：对照（CK）、铬胁迫（Cr）、褪黑素缓解（MC，1 μM MT预处理后加Cr）、谷胱甘肽缓解（GC，2 mM GSH预处理后加Cr）。在相同发育阶段采集叶片样本，进行RNA提取、文库构建，并在Illumina HiSeq X Ten平台上进行高通量测序。生信分析流程包括：使用FastQC和Trimmomatic进行原始数据质量控制，通过HISAT2将高质量reads比对到甘薯参考基因组，利用DESeq2进行差异表达基因（DEGs）分析（阈值：|log2FC| ≥ 1, FDR ≤ 0.05）。对DEGs进行GO和KEGG功能富集分析，通过visCluster进行表达谱聚类，并利用WGCNA构建与关键生理生化性状（如生物量、光合参数、氧化损伤标志物、铬积累量）相关的基因共表达网络，从中挖掘核心调控模块与枢纽基因。最后，通过RT-qPCR对八个候选基因的表达模式进行验证，确认转录组数据的可靠性。
结果解析：揭示GSH与MT差异化的转录调控蓝图
3.1 转录组数据概览与假设验证
样本的主成分分析（PCA）与层次聚类显示，GC处理组在转录谱上明显与Cr和MC组分离，而MC组与Cr组聚类更近，表明GSH引发了更显著且独特的转录重编程。样本间表达相关性分析进一步支持了这一点：GC与Cr的相关系数（r = 0.92–0.94）低于MC与Cr的系数（r = 0.97），表明GSH对铬胁迫响应谱的修饰作用更强。
3.2 差异表达基因的识别与模式
在六组对比中，共鉴定出7,734个DEGs。铬胁迫（Cr vs CK）诱导了5,601个DEGs。直接比较显示，GSH（GC vs Cr）重编程了824个基因，其中82.5%为下调；而MT（MC vs Cr）仅影响319个基因。GC与MC对比（GC vs MC）有533个DEGs，其中64.5%在GC中表达更高。这些数据定量证实了GSH是比MT更强的转录调控因子。
3.3 差异表达基因的生物学功能
GO富集分析揭示了共同与特异的应答模式。所有处理（Cr、GC、MC vs CK）均导致光合作用、碳固定、细胞壁组织等生长相关功能的基因被显著抑制。在应激应答方面，Cr和MC主要激活了外源物跨膜运输、毒素代谢等通路。而GSH的独特之处在于，其特异性富集了与活性氧（ROS）应答、乙醛代谢过程、解毒作用等相关的GO条目，突出了其对氧化还原和解毒网络的强化。
3.4 差异表达基因参与的代谢通路
KEGG通路分析显示，铬胁迫上调了MAPK信号、谷胱甘肽代谢、α-亚麻酸代谢等防御通路，同时下调了光合作用、碳固定等初级代谢通路。GSH处理（GC vs CK）在共享这些变化的基础上，额外显著富集了半胱氨酸和蛋氨酸代谢、糖酵解/糖异生、丙酮酸代谢以及内质网中的蛋白质加工等通路，这与其直接增强硫同化和中心碳代谢以支持解毒的功能一致。相比之下，MT处理（MC vs Cr）则更偏向于富集植物激素信号转导和膜/囊泡运输等通路。
3.5 聚类分析鉴定谷胱甘肽相关基因模块
将所有DEGs进行表达谱聚类，得到四个主要模块（C1-C4）。其中，C1包含的铬诱导基因能被GSH有效抑制回基础水平；C4则代表了一组被GSH进一步强化的铬诱导基因，与ROS和H2O2应答、丙酮酸代谢、内质网蛋白质加工等功能相关，构成了GSH特异的保护性模块。而C3模块中的基因在MT处理下表达进一步增强，定义了一个MT倾向的解毒模块。
3.6 与生理生化性状关联模块的鉴定
WGCNA分析构建了基因共表达网络，并鉴定出与关键性状显著相关的模块。其中，蓝色模块的表达特征与植株生物量、水分状态、光合性能呈极显著正相关，与膜损伤、氧化标志物（如MDA、H2O2）及铬的积累转运呈极显著负相关。该模块在GSH处理下上调最明显，代表了一个核心的铬耐受性模块，囊括了GSH驱动下生长与光合恢复、氧化损伤降低的分子网络。
3.7 与性状改善相关的候选枢纽基因挖掘
从蓝色模块中，基于节点连接度筛选出前30个枢纽基因。这些基因功能多样，包括：
- •
光合与叶绿素合成：HEMA1（谷氨酰-tRNA还原酶1）、CAB4、PSAN等，直接关联叶绿素含量和光系统效率的恢复。
- •
氧化还原与解毒：过氧化物酶基因（响应氧化应激、H2O2分解），与ASA-GSH循环协同清除ROS。
- •
运输与渗透调节：TIP2-1（液泡膜内在蛋白），可能参与液泡区隔化；PHT1-4（磷酸盐转运蛋白）；TPPA（海藻糖磷酸磷酸酶）。
- •
生长、信号与结构调节：GA20OX1/GA20OX2（赤霉素生物合成）、HAT4（HD-ZIP转录因子）、XTH1（木葡聚糖内转葡糖基酶/水解酶）等，整合了激素信号、细胞壁重塑以支持生物量恢复。
3.8 RT-qPCR验证RNA-seq结果
对八个候选基因（如金属转运蛋白NRAMP5、谷胱甘肽S-转移酶HSP26-A、转录因子WRKY70和ERF113）进行RT-qPCR验证，其表达变化趋势与RNA-seq数据高度一致（Pearson相关系数 r = 0.915），证实了转录组分析结果的可靠性。
讨论与结论：机制阐释与展望
本研究表明，铬胁迫在甘薯中触发了一个核心应激特征，即抑制光合与生长相关初级代谢，同时激活基础解毒与氧化防御机制。GSH和MT虽然都能提供保护，但其分子作用模式存在显著差异。MT主要通过调控植物激素信号转导和细胞膜运输通路来微调应激响应。而GSH则表现出更强的转录调控能力，它独特地强化了以硫同化（半胱氨酸和蛋氨酸代谢）和谷胱甘肽代谢为核心的解毒网络，直接提升了植物合成植物螯合肽（PCs）的能力和ROS清除效率，从而构建了一个更高效、以氧化还原为中心的防御性代谢网络。通过WGCNA鉴定的30个枢纽基因，是连接GSH处理、关键生理生化性状改善与转录重编程的核心分子节点，为未来通过分子育种或基因工程培育铬耐受型甘薯品种提供了明确且可操作的高价值靶点。综上，这项研究从系统水平揭示了GSH在甘薯铬解毒中的优越性及其分子基础，为作物重金属胁迫的绿色防控策略奠定了重要的理论依据。
