《Stresses》：Investigating the Effects of Conventional and No-Tillage Cultivation Methods on Plant Physiological Processes Using Genome-Wide Transcriptomic Analysis
编辑推荐：
尽管对可持续农业系统的需求日益增长，但耕作实践的长期影响仍存在争议。免耕（NT）系统具有若干潜在优势，包括改善土壤结构、增强土壤生物活性和减少环境胁迫，但其应用也可能带来杂草控制困难或作物产量不稳定等挑战。尽管先前的研究已广泛调查了NT系统对土壤和作物的影响，
尽管对可持续农业系统的需求日益增长，但耕作实践的长期影响仍存在争议。免耕（NT）系统具有若干潜在优势，包括改善土壤结构、增强土壤生物活性和减少环境胁迫，但其应用也可能带来杂草控制困难或作物产量不稳定等挑战。尽管先前的研究已广泛调查了NT系统对土壤和作物的影响，但对植物细胞适应机制，特别是基因表达和生化响应的了解仍然有限。本研究旨在利用整合转录组学和生化方法，比较常规耕作（CT）和NT系统对向日葵植株的影响。研究人员对来自三个不同田间地点的叶片样本进行了基于下一代测序的全基因组转录组分析，并辅以与选定代谢过程相关的生化参数测量。探索性转录组分析表明，与CT系统相比，在NT系统中生长的植物发生了与初级代谢过程相关的若干基因表达变化，包括光合作用、细胞呼吸、碳水化合物代谢以及某些氨基酸的生物合成。同时，研究人员观察到转录模式表明几个次级代谢途径活性增加，这可能与环境胁迫的适应机制有关。对总可溶性糖、粗蛋白、总酚和总黄酮的测定为转录组结果所指示的变化提供了独立的生化支持。研究结果表明，耕作系统影响向日葵的细胞调控过程。在适应免耕环境的过程中，植物能够同时维持基本代谢过程并激活防御机制，这可能有助于适应改变的生长条件。该研究有助于更好地理解耕作系统对植物分子和生理影响的后果。
**论文解读：免耕系统对向日葵转录组与生化响应的多地点整合分析**
**研究背景与问题**
全球气候变化、粮食需求增长及农业资源压力加剧，推动可持续作物生产系统的开发。传统集约化耕作长期导致土壤结构退化、有机质损失及生物活性下降，因此土壤保护技术（如免耕系统）受到广泛关注。免耕（NT）系统可改善土壤结构、提高保水能力、促进微生物活性并有助于碳封存，但其农艺效益并非普遍适用，受气候、土壤类型、作物种类及栽培技术显著影响。既往研究多聚焦于土壤物理、化学和微生物属性及作物产量，但对植物细胞调控机制——尤其是基因表达模式与生化响应——如何因土壤环境改变而直接受影响，仍知之甚少。植物通过激活胁迫响应基因、重组代谢途径及修饰抗氧化防御机制等复杂分子过程适应环境变化，但长期田间条件下NT系统影响植物适应性的分子机制尚不明确。
**研究目的与意义**
本研究假设长期耕作制度变化对植物施加环境选择压力，改变与胁迫响应和代谢相关的基因表达过程，并反映在可测量的生化差异上。为此，研究人员在匈牙利三个相邻田间地点（Ecseg、Kiszombor、Vityapuszta，均属黑钙土，温带大陆性气候，年均温约10°C，年降水600–650 mm），对长期（10年）免耕与常规耕作（CT）下的向日葵（品种SY Excellio）进行整合转录组学和生化分析。研究揭示了免耕系统如何通过基因表达重编程影响植物初级与次级代谢，并阐明位点特异性响应规律，为理解可持续农业中作物的分子适应机制提供新见解。论文发表在《Stresses》。
**主要技术方法**
样本来源：匈牙利三个田间地点（Ecseg、Kiszombor、Vityapuszta）的向日葵叶片，每个地点分别设常规耕作和免耕处理，各采集4个生物学重复（共6个混合RNA样本）。主要方法：RNA提取及mRNA富集后，采用NEXTFLEX试剂盒构建RNA-seq文库，在Illumina NovaSeq 6000平台进行低深度测序；使用Trinity软件进行de novo转录组组装，通过OmicsBox平台（CloudBlast、Blast2GO、eggNOG）进行基因鉴定与功能注释；利用NOISeq算法（概率>0.9，且|M|≥1）进行差异表达分析；通过KEGG和Plant Reactome数据库进行通路富集。生化参数测定包括总可溶性糖（TSS）、粗蛋白（CPC）、总酚（TPC）和总黄酮（TFC），并采用单因素和双因素方差分析（ANOVA）进行统计检验。
**研究结果**
**2.1 全基因组转录组分析**
通过de novo组装获得66,372个转录本（43,749个基因），平均长度995 bp。三个地点的差异表达分析显示：Ecseg地点的免耕样本中，1098个基因上调、3295个下调；Kiszombor地点314个上调、508个下调；Vityapuszta地点1299个上调、2629个下调。基因富集分析发现，所有三个地点免耕处理下，与代谢和解剖结构发育相关的生物过程（BP）均被激活。联合通路分析表明，在Ecseg，上调最显著的通路包括精氨酸和脯氨酸代谢、苯丙氨酸代谢、甘油酯代谢、β-丙氨酸代谢及苯丙氨酸、酪氨酸和色氨酸生物合成；下调最显著的通路包括半胱氨酸和甲硫氨酸代谢、卡尔文循环碳固定、甘氨酸/丝氨酸/苏氨酸代谢、半乳糖代谢、戊糖和葡萄糖醛酸相互转化、甘油磷脂代谢、嘌呤代谢、卟啉代谢、丙酮酸代谢、淀粉和糖代谢、萜类骨架构建及糖酵解-糖异生。Kiszombor地点上调过程包括氧化磷酸化、嘌呤代谢、萜类骨架生物合成、核黄素代谢、嘧啶代谢及多种植物次生代谢物生物合成；下调过程包括淀粉和糖代谢、甘油磷脂代谢、戊糖和葡萄糖醛酸相互转化及半乳糖代谢。Vityapuszta地点上调过程包括糖酵解/糖异生、嘌呤代谢、丙酮酸代谢、半胱氨酸和甲硫氨酸代谢、淀粉和糖代谢、甘油磷脂代谢、萜类骨架生物合成、甘氨酸/丝氨酸/苏氨酸代谢、卡尔文循环碳固定及苯丙烷生物合成；下调过程包括氨基糖和核苷酸糖代谢、卟啉代谢及鞘脂代谢。文氏图显示，多数基因表达变化由耕作方式引起，少部分由环境因素引起。
**2.2 实验室生化参数测定**
在Ecseg，免耕导致可溶性糖含量、粗蛋白浓度、酚类化合物和类黄酮含量均增加，表明初级代谢增强且次生代谢启动，与转录组结果一致。Kiszombor的生化参数变化趋势类似但差异较小，显示代谢重排及适应机制启动。Vityapuszta的生化参数变化趋势与Ecseg和Kiszombor相似，但处理间差异不显著，且产量最高。双因素方差分析表明，地点和耕作方式对TSS、CPC、TPC、TFC均有显著影响，且交互作用显著（除TPC的交互作用不显著），说明免耕效应部分依赖于环境条件。
**2.3 收获后产量研究**
产量结果显示，Ecseg免耕区产量显著低于常规耕作（P
1%显著），Kiszombor免耕区产量略低（P
10%显著），Vityapuszta免耕区产量与常规耕作无显著差异。表明免耕对产量的影响高度位点特异性。
**2.4 耕作与地点的综合效应**
双因素方差分析进一步确认，粗蛋白含量（CPC）受地点、耕作方式及其交互作用强烈影响（P
2,12=202.68, 631.62, 158.68，均P<0.001），免耕在Ecseg和Vityapuszta显著增加蛋白含量，在Kiszombor变化较小。总酚含量（TPC）受地点和耕作方式显著影响，但交互作用不显著，表明免耕在所有地点均一致提高酚含量。总黄酮含量（TFC）受地点、耕作方式及交互作用显著影响，免耕在所有地点提高黄酮含量，但程度因地点而异（Vityapuszta最突出）。总体而言，免耕系统并非引起统一响应，而是与地点环境相互作用，植物通过动态调控生化反应适应耕作变化。
**讨论与结论**
三个地点的转录组模式呈现不同生理状态，但共同指向免耕触发的资源重新分配策略。Ecseg呈现代谢极化：上调氨基酸、脂质和芳香代谢（如脯氨酸、精氨酸、苯丙烷类）指示防御和氮储存策略，同时碳固定、能量代谢、核苷酸和卟啉通路强烈下调，表明光合和生长功能停滞，植物转向适应而非生长。Kiszombor表现为高能量需求但资源受限的适应：上调氧化磷酸化、辅因子产生及次生代谢，下调碳水化合物和细胞壁构建，表明植物依赖内部储备维持能量和防御。Vityapuszta呈现相对平衡的代谢：光合碳固定、糖酵解、嘌呤、脂质、氨基酸及次生代谢均上调，仅氨基糖、卟啉、鞘脂代谢轻微下调，说明植物在维持生长的同时增强防御。
**结论翻译**：综合Ecseg、Kiszombor和Vityapuszta地点的基因表达模式，向日葵采用高度可塑的代谢策略，对由耕作方式差异导致的不同环境条件作出强烈响应。三个地点并非仅在代谢活性强度上不同，而是代表不同的生理状态，但汇合于共同的适应性生存策略。Ecseg表现为代谢极化，免耕诱导强烈的防御、胁迫耐受和氮储存策略（脯氨酸、精氨酸、苯丙烷类），同时碳固定和能量代谢显著下调，植物放弃生长而转向适应。Kiszombor呈现高能量需求但资源受限的适应状态，上调能量和辅因子产生过程，下调碳水化合物和细胞壁构建，植物依赖内部储备维持防御。Vityapuszta呈现协调的活跃代谢网络，光合作用、糖酵解、嘌呤、脂质、氨基酸及次生代谢均上调，仅少数过程微调，表明植物在优化资源的同时保持生长。总之，尽管存在显著的位点特异性代谢调控差异，三个田间地点揭示了免耕响应的共同适应性生存策略：免耕持续促进资源从初级生长过程向次生代谢和防御相关途径的重新分配，指示一种保守的、地点无关的生存响应。这些发现表明向日葵并非被动承受环境影响，而是根据局部条件在转录水平动态重组代谢，同时维持共同的适应框架。这种代谢灵活性可能有助于向日葵在可持续农业中的恢复力和作物稳定性。
