整合生理学与转录组学揭示谷子（Setaria italica）响应LiCl胁迫的适应机制

《Food and Energy Security》：Integrative Physiology and Transcriptomics Reveal the Adaptive Mechanisms of Foxtail Millet (Setaria italica) in Response to LiCl Stress

　　

随着锂离子电池的普及，锂已成为一种重要的环境污染物，然而其在植物中的毒性和耐受机制仍知之甚少。本研究探讨了C4模式作物谷子（Setaria italica）对LiCl胁迫的生理和分子响应。生理学分析揭示了一种浓度依赖性效应：低浓度Li⁺可激活抗氧化系统，而暴露于50 mg/L LiCl则触发剧烈的氧化爆发，导致抗氧化酶（SOD、POD、CAT）活性受抑、脂质过氧化和显著的生长抑制。对品种“晋谷21”的转录组分析鉴定了1562个共同调控的差异表达基因（DEGs），表明Li⁺胁迫扰乱Na⁺/K⁺稳态并使代谢途径重编程，包括支链氨基酸降解和植物-病原互作途径的上调。为验证这些发现，研究人员分析了5个具有差异耐受性的额外品种。比较分析表明，“冀谷22”观察到的强耐受性与关键基因的强烈诱导相关——特别是转录因子SiBHLH148、液泡转运蛋白SiNHX1和脂质转运蛋白SiDIR1——表明它们在维持ROS稳态和离子区室化中发挥关键作用。这些结果阐明了谷子适应LiCl的分子基础，并为培育耐受锂污染的作物提供了关键遗传靶点。

谷子响应LiCl胁迫的整合生理学与转录组学研究解读

随着锂离子电池在新能源领域的广泛应用，锂在生产、使用和废弃过程中大量释放到环境中，已成为一种新兴的土壤污染物。锂主要以Li⁺离子形式存在，具有高迁移性和生物有效性，易被植物吸收并沿食物链传递，威胁生态系统和人体健康。然而，植物对锂的毒性响应和耐受机制尚不清楚。谷子（Setaria italica）是一种基因组紧凑、生育期短、对非生物胁迫具有较强抗性的C4模式作物，适合用于研究锂胁迫的生理和分子机制。为此，研究人员以谷子为材料，系统研究了LiCl胁迫下的适应机制。相关论文发表在《Food and Energy Security》。

研究人员以“晋谷21”为材料，设置0、0.1、5、50 mg/L四个LiCl浓度处理4天，测定生长形态、叶绿素含量、生物量、H₂O₂含量、超氧化物歧化酶（SOD）、过氧化物酶（POD）、过氧化氢酶（CAT）活性及丙二醛（MDA）含量；同时进行转录组测序（RNA-seq），分析差异表达基因（DEGs）及其富集通路；并进一步在5个谷子品种（“豫谷1”、“冀谷39”、“冀谷22”、“大同40”、“长谷16”）中验证候选基因的表达。品种种子由山西农业大学张海英博士提供。主要技术方法包括：利用1/2 MS液体培养基进行LiCl处理；使用商业试剂盒测定生理生化参数；基于Illumina NovaSeq 6000平台进行RNA-seq，用fastp质控、HISAT2比对参考基因组、HTSeq-count定量，以p<0.05且|log₂(fold change)|>1筛选DEGs；用定量实时PCR（qRT-PCR）验证关键基因表达。所有实验设3个生物学重复和3个技术重复。

结果部分如下：

3.1 不同锂浓度对“晋谷21”生长形态、生理生化特性的影响。通过表型观察和生理指标测定，发现低浓度（0.1或5 mg/L）LiCl处理未引起显著表型变化，但50 mg/L处理导致叶片斑驳、卷曲和萎蔫，叶绿素含量和总生物量显著下降。生化分析显示，0.1 mg/L处理显著升高SOD、POD、CAT活性；5 mg/L时H₂O₂和MDA积累，但酶活性仍上调；50 mg/L时H₂O₂过量积累，抗氧化酶活性骤降，MDA剧增，表明抗氧化系统崩溃和严重膜脂过氧化。

3.2 不同浓度LiCl处理下谷子转录组测序及DEGs分析。转录组测序共获得78.84 Gb clean data，Q30比例96.38%~98.13%，比对率96.72%~97.50%。与对照相比，0.1、5、50 mg/L处理分别有3427、3210、3833个DEGs，其中1562个DEGs在所有处理中共同变化，构成核心调控网络。

3.3 GO富集分析和KEGG通路分析。基因本体论（GO）分析显示，低浓度下富集“对几丁质的响应”、“茉莉酸介导信号通路调控”等生物学过程；中浓度富集“亮氨酸分解代谢”、“昼夜节律”等；高浓度富集“蛋白磷酸化”、“对化学物质的响应”、“蛋白丝氨酸/苏氨酸激酶活性”等。京都基因与基因组百科全书（KEGG）分析显示，“缬氨酸、亮氨酸和异亮氨酸降解”在所有处理中均显著富集，而“谷胱甘肽代谢”、“植物激素信号转导”、“植物-病原互作”等通路在不同浓度下被激活。

3.4 LiCl胁迫损害植物-病原互作和Na⁺/K⁺稳态通路。通过筛选共同DEGs，研究人员绘制了效应子触发免疫（ETI）和系统获得抗性（SAR）相关基因表达谱，发现50 mg/L处理诱导SiRGA3介导的ETI和SiNPR1相关的SA信号，但SiDIR1、SiTGA、SiPR1s未同步上调，表明下游SAR可能不完整。此外，LiCl胁迫下SiKEA3上调，而SiHKT11、SiCHX15和SiNHX1下调，提示Na⁺/K⁺相关转运系统受到干扰。

3.5 五个谷子品种在LiCl胁迫下的形态表型。对“豫谷1”、“冀谷39”、“冀谷22”、“大同40”、“长谷16”进行LiCl处理，发现“长谷16”最敏感，“冀谷22”和“大同40”耐受性最强，其中“冀谷22”在高浓度下叶绿素含量显著增加。综合“晋谷21”表型，六个品种耐受性排序为：“冀谷22”>“大同40”>“冀谷39”>“晋谷21”>“豫谷1”>“长谷16”。

3.6 不同谷子品种中LiCl胁迫相关基因的表达。qRT-PCR验证了5个候选基因（SiDIR1、SiRGA3、SiKEA3、SiBHLH148、SiNHX1）在“晋谷21”中的表达模式与转录组一致。耐性品种“冀谷22”在50 mg/L下SiDIR1和SiBHLH148强诱导，SiNHX1在“大同40”中上调；而敏感品种“豫谷1”和“长谷16”中SiDIR1和SiNHX1下调更明显。结果表明，SiDIR1、SiNHX1和SiBHLH148的强诱导与LiCl耐受性正相关。

总结讨论部分，研究人员认为LiCl胁迫通过ROS积累抑制谷子生长。高浓度Li⁺导致抗氧化系统崩溃和膜脂过氧化，这是生长抑制的主要驱动因素。转录组分析显示LiCl胁迫重编程大量基因表达，并激活植物-病原互作信号通路，但SAR反应不完全，可能是一种防御与能量消耗之间的权衡。与典型盐胁迫不同，Li⁺没有激活常规的Na⁺/K⁺稳态响应，而是抑制SiHKT11、SiCHX15和SiNHX1等转运蛋白表达，干扰离子区室化。基于品种差异，耐性品种中SiBHLH148、SiNHX1和SiDIR1的诱导程度更高，这些基因可能分别通过促进ROS清除、液泡区室化和脂质信号传导增强耐锂性。

研究结论：对“晋谷21”的生理和转录组分析表明，高浓度LiCl（50 mg/L）通过损害ROS清除能力抑制生长。LiCl暴露改变了参与支链氨基酸分解、昼夜节律、甘油酯代谢、植物-病原互作和离子稳态的基因表达。不同耐受性品种的比较表明，SiDIR1、SiNHX1和SiBHLH148可能通过调控ROS稳态和锂转运而贡献于LiCl耐受性。这些发现为理解谷子锂胁迫响应和培育耐锂品种提供了基础。