《Plants》：Multi-Omics Analysis Reveals Hormonal, miRNA-Associated, and NLR-Linked Aphid Defense in Lycium barbarum
编辑推荐：
蚜虫侵扰制约了枸杞生产，但其防御架构仍不清晰。研究人员整合了比较转录组学、植物激素谱分析、小RNA测序以及全基因组核苷酸结合富亮氨酸重复（NLR）分析，以解析抗性枸杞（宁杞-10号）与易感枸杞之间的蚜虫防御机制。蚜虫生物测定证实，抗性枸杞对蚜虫施加了显著的适合
蚜虫侵扰制约了枸杞生产，但其防御架构仍不清晰。研究人员整合了比较转录组学、植物激素谱分析、小RNA测序以及全基因组核苷酸结合富亮氨酸重复（NLR）分析，以解析抗性枸杞（宁杞-10号）与易感枸杞之间的蚜虫防御机制。蚜虫生物测定证实，抗性枸杞对蚜虫施加了显著的适合度代价，降低了其存活率、寿命和繁殖力。转录组比较揭示，抗性种质中差异表达基因数量较少但富含更多防御相关功能。加权基因共表达网络分析（WGCNA）鉴定出一个应激相关的青绿色模块（288个基因，r=0.981），包含类受体激酶和丝裂原活化蛋白激酶（MAPK）组分，桥接了胞外信号感知与激素信号传导。植物激素谱分析显示，水杨酸（SA）持续积累至72小时，而茉莉酸（JA）在24小时短暂达到峰值，导致SA/JA比值升高，表明信号转向以SA为主导。MiR398a-3p和miR319-3p的动态变化与这一SA/JA转变相吻合，提示其可能与激素协调调控存在关联，并为验证提供了候选靶标。NLR分析鉴定出151个候选基因，其中18个在SA主导信号下上调的抗性枸杞偏向性NLR，暗示其与SA信号存在联系。综上所述，这些发现支持了一个整合了SA信号、转录重编程以及偶联的miRNA和NLR应答的多层次蚜虫防御框架。
研究背景与研究意义
枸杞属植物是全球重要的药食同源经济树种，其中宁夏枸杞作为中国西北地区的主导栽培种，年产值已达数十亿元。然而，蚜虫危害严重制约枸杞产业发展，尤其在宁夏和青海主产区，高峰期可造成20%-40%的产量损失。当前主要依赖化学杀虫剂进行防控，但长期使用已导致抗药性蚜虫种群进化、天敌群落破坏以及果实农药残留等问题。因此，迫切需要基于天然防御机制的分子理解来发展可持续管理策略。
植物对蚜虫的防御由多个调控层面协同完成。植物激素，特别是水杨酸（salicylic acid, SA）和茉莉酸（jasmonic acid, JA），充当核心信号枢纽，SA常主导应对韧皮部取食昆虫的反应，而JA驱动伤口和草食动物诱导的防御程序。此外，SA还能直接抑制植食性昆虫的免疫反应。在转录后水平，微小RNA（microRNA, miRNA）增加了额外的调控层，特定miRNA如miR393、miR319、miR398和miR156可调节激素信号和防御基因表达。同时，激素信号也可能与核苷酸结合富亮氨酸重复（nucleotide-binding leucine-rich repeat, NLR）免疫受体的应答相关联。
尽管这些SA/JA和miRNA模块在被子植物中广泛保守，但它们在枸杞这类多年生木本植物中的部署机制仍知之甚少。与一年生模式植物不同，枸杞作为多年生灌木必须在防御与多年生长和果实生产之间取得平衡，这可能对防御信号施加了独特的选择压力。此外，两种亲缘关系密切但蚜虫抗性截然不同的枸杞种质——抗性宁夏枸杞（宁杞-10号）和易感枸杞——为解析该经济重要属内保守防御模块的部署方式提供了天然的比较系统。然而，关于枸杞蚜虫防御的分子研究仍处于起步阶段，缺乏涵盖转录组重编程、植物激素动态、miRNA相关转录后调控以及NLR免疫受体应答的综合视角。
为此，本研究整合了比较转录组学、小RNA测序、靶向植物激素谱分析和全基因组NLR分析，系统表征了抗性宁夏枸杞相对于易感枸杞的蚜虫防御机制，旨在解析抗性种质中激素和转录后调控的时间动态，鉴定种质偏向性的NLR候选基因，为枸杞抗蚜育种提供分子框架和优先候选靶标。
主要关键技术方法
研究人员采用了以下关键技术方法：（1）年龄-阶段两性生命表分析，用于量化蚜虫在两种枸杞上的适合度参数；（2）比较转录组学分析，对宁夏枸杞和枸杞在蚜虫侵染72小时后的叶片样本进行RNA测序，采用DESeq2进行差异表达基因（differentially expressed gene, DEG）鉴定，并使用加权基因共表达网络分析（weighted gene co-expression network analysis, WGCNA）构建共表达模块；（3）靶向液相色谱-串联质谱（LC-MS/MS）植物激素定量分析，对抗性宁夏枸杞在0、24和72小时的激素浓度进行时间序列检测；（4）小RNA测序，对抗性宁夏枸杞在0、8、24、36、48和72小时六个时间点的18个小RNA文库进行测序，鉴定已知和新颖miRNA，并使用TargetFinder预测靶基因；（5）基于隐马尔可夫模型（hidden Markov model, HMM）的全基因组NLR基因鉴定与分类，从宁夏枸杞参考基因组中系统识别候选NLR基因。
研究结果
2.1 蚜虫侵染诱导抗性和易感枸杞产生对比鲜明的转录程序
通过年龄-阶段两性生命表分析，研究人员确认抗性宁夏枸杞显著降低了蚜虫的存活率、寿命和繁殖力，净生殖率（R0）、内禀增长率（rm）和平均世代时间（T）均显著低于易感枸杞。在侵染40天后，抗性枸杞仅表现轻微褪绿和有限卷叶，而易感枸杞则出现明显叶片变形、广泛褪绿和高密度蚜虫定殖。转录组分析发现，抗性枸杞在蚜虫胁迫下有3441个DEG（1980个上调，1461个下调），而易感枸杞有8397个DEG（5083个上调，3314个下调）。功能富集分析显示，抗性枸杞中苯丙烷代谢过程和植物-病原互作途径显著富集，且DEG更集中于防御相关通路。热图分析显示，防御基因PR1（LOC132605191，FC=8.38）、WRKY70（LOC132605276，FC=5.37）和CHS（LOC132640283，FC=-1.30）表现出种质偏向性应答。
2.2 WGCNA将青绿色共表达模块与胁迫处理和防御应答相关联
WGCNA分析鉴定了10个共表达模块，其中青绿色模块（288个基因）与胁迫处理的相关性最强（r=0.981，p<0.0001）。该模块的枢纽基因包括类CRK25受体样激酶（receptor-like kinase, RLK）候选基因（LOC132607248）、二萜生物合成相关基因以及多个WRKY相关转录因子。基因本体（Gene Ontology, GO）过代表分析显示，防御应答和应激应答在该模块枢纽基因中显著富集，京都基因与基因组百科全书（Kyoto Encyclopedia of Genes and Genomes, KEGG）分析也确认了植物-病原互作的显著富集。RT-qPCR验证了PR4、LOX、GST、WRKY70和PTI5五个基因的表达趋势与RNA-seq结果一致。
2.3 NLR库显示出种质偏向性的转录激活
基于NB-ARC结构域搜索，研究人员在宁夏枸杞参考基因组中鉴定出151个候选NLR基因，主要包括其他结构域架构（113个）、CC-NBS-LRR/NB-ARC（22个）、TIR-NBS-LRR（12个）和NB-ARC-only（4个）。比较分析显示，蚜虫侵染后，47个NLR在抗性枸杞中显著上调，76个在易感枸杞中显著上调，其中29个为共享，18个为抗性枸杞偏向性，47个为易感枸杞偏向性。这18个抗性枸杞偏向性候选NLR在抗性枸杞中显著上调但在易感枸杞中未达到相同显著性阈值，编码TIR-NBS-LRR和CC-NBS-LRR/NB-ARC结构域架构的蛋白质。
2.4 持续的SA积累和瞬时的JA爆发定义了激素动态
靶向植物激素谱分析显示，在抗性宁夏枸杞中，游离SA从0小时的33.24 ng/g上升至24小时的1215.40 ng/g（36.6倍），并在72小时继续增加至1456.17 ng/g（43.8倍）。SA糖苷（SAG）平行积累（72小时达9.9倍），表明通过糖基化进行活性SA稳态调节。相比之下，JA和JA-Ile在24小时呈现瞬时爆发（分别较基线增加3.6倍和2.8倍），随后迅速下降。描述性SA/JA比值随时间推移持续升高。脱落酸（abscisic acid, ABA）水平呈现“对勾形”波动，在24小时下降后反弹；L-苯丙氨酸丰度在24小时下降后逐渐恢复，两者动态趋势大致与瞬时JA/JA-Ile爆发相反。
2.5 协调的miRNA表达动态提示转录后调控
小RNA测序从18个文库中鉴定出253个miRNA，包括84个已知和169个新颖miRNA。差异表达分析显示，stu-miR398a-3p在24小时从2224 TPM增至23673 TPM（10.64倍），stu-miR408b-3p约增加2.01倍，两者峰值时间与JA爆发相吻合。stu-miR156e在36小时达到最高平均丰度（相对0小时为2.09倍），而stu-miR319-3p在24小时增加1.60倍，但到72小时降至基线以下。TargetFinder预测了这四个miRNA的候选靶转录本，包括stu-miR156e预测靶向SPL17-like和SPL12-like转录本，stu-miR319-3p预测靶向TCP4和PCF5。
讨论与结论
讨论部分指出，本研究通过整合比较转录组学、时间序列植物激素谱分析、小RNA测序和全基因组NLR分析，在统一生物学模型中解析了枸杞-蚜虫系统的多组学框架。关键发现是抗性宁夏枸杞激活了更少的DEG（3441个）但更富集于防御通路，这与有效防御可能依赖于通路激活的精准度和聚焦度而非响应基因数量的概念一致。WGCNA鉴定的青绿色模块包含RLK和MAPK信号组分，表明RLK介导的胞外信号感知与WRKY驱动的转录重编程在抗性种质中紧密协调。
植物激素时间序列分析揭示了SA的持续积累（72小时达43.8倍）和JA/JA-Ile的瞬时爆发（24小时达3.6倍和2.8倍），SA/JA比值持续升高，符合经典SA-JA拮抗框架。ABA和L-苯丙氨酸的相反动态轨迹提示早期防御中存在可能的代谢重构：ABA在24小时的瞬时下降可能解除其对JA信号的拮抗作用，而L-苯丙氨酸的消耗可能反映其被用于苯丙烷衍生防御代谢物的合成。
小RNA时间序列显示，stu-miR398a-3p和stu-miR319-3p在24小时达到峰值，与JA爆发和早期SA积累时间耦合，提示这些miRNA可能作为激素信号的下游或平行调节因子参与防御微调。NLR分析鉴定的18个抗性枸杞偏向性候选NLR在SA主导信号下优先上调，暗示这些NLR候选基因可能与SA通路协同发挥作用。
研究结论部分总结如下：本研究通过整合比较转录组学、植物激素谱分析、小RNA测序和全基因组NLR分析，在自然比较系统中对枸杞蚜虫防御进行了多组学解析。结果定义了一个以植物激素动力学为中心、miRNA和NLR作为协作调控层的时序协调防御框架。蚜虫攻击后，抗性宁夏枸杞激活了比易感枸杞更少的DEG（3441对8397），但这些基因在防御通路中富集程度更高，符合更聚焦的防御应答模式。抗性种质的植物激素谱分析揭示了清晰的时间分裂：JA和JA-Ile在24小时瞬时爆发，而SA逐步积累至72小时，SA/JA比值升高，提示SA主导信号的趋势。在转录后水平，stu-miR398a-3p和stu-miR319-3p在24小时达到峰值，与JA爆发和早期SA积累时间耦合。在NLR层面，鉴定出151个候选NLR基因，其中18个抗性枸杞偏向性NLR在SA主导信号下优先上调。综合这些发现，研究支持了一个多层模型，其中RLK/MAPK介导的感知启动了以两个时间上不同的阶段展开的防御应答：24小时的早期阶段，以瞬时JA爆发伴随ABA/L-苯丙氨酸代谢重构和miRNA调控为特征；72小时的持续阶段，以SA逐步积累驱动SA主导信号并与NLR候选基因上调相吻合。该框架为理解枸杞蚜虫防御提供了分子基础，并为抗蚜育种确定了优先候选靶标。
