《Plant Stress》：Deciphering the molecular mechanisms of maize defence against the fall armyworm (Spodoptera frugiperda) through integrated transcriptomics and metabolomics
编辑推荐：
草地贪夜蛾（Fall armyworm, FAW, Spodoptera frugiperda）是一种极具破坏性的多食性害虫，严重威胁全球特别是热带和亚热带农业生态系统中的玉米生产。理解寄主植物抗性的分子基础对于制定持久且可持续的害虫管理策略至关重要。本研究采
草地贪夜蛾（Fall armyworm, FAW, Spodoptera frugiperda）是一种极具破坏性的多食性害虫，严重威胁全球特别是热带和亚热带农业生态系统中的玉米生产。理解寄主植物抗性的分子基础对于制定持久且可持续的害虫管理策略至关重要。本研究采用整合转录组和代谢组学的方法，剖析了在FAW侵害和未侵害条件下，抗性自交系CML338和感性自交系CML451中差异的防御反应。RNA测序（RNA-seq）产生了高质量的测序数据集，比对效率达84–87%，揭示了抗性基因型在侵害后广泛的转录重编程。多变量分析显示了清晰的基因型和处理依赖性分离，主成分解释了超过60%的总方差。差异表达基因（Differentially Expressed Genes, DEGs）显著富集于植物激素信号转导、丝裂原活化蛋白激酶（Mitogen-Activated Protein Kinase, MAPK）级联、转录调控、脂肪酸代谢和次生代谢物生物合成等通路。非靶向液相色谱-高分辨质谱（Untargeted LC-HRMS）代谢组学分析显示，CML338中推定抗性相关代谢物显著积累，包括脂肪酸（亚油酸和花生四烯酸）和萜类化合物（α-松油醇、辣椒玉红素、5-羟基玉米黄质、中黄质和紫黄质）。整合转录组-代谢组相关性分析（Pearson Correlation Coefficient, PCC ≥ 0.85, P ≤ 0.05）表明代谢物积累与关键生物合成基因（包括TPS、DXS、PLD、LOX和FAD）的表达之间存在紧密的协调关系。此外，分子对接（in-silico molecular docking）显示α-松油醇与FAW肠道酶氨基肽酶N1（Aminopeptidase N1）之间存在强结合力（?6.6 kcal mol?1），提示了潜在的相互作用；然而，这仅是计算分析得出的推定预测，需生化实验和生测验证确认。总之，本研究提供了玉米–FAW互作的系统水平理解，并鉴定了用于培育抗虫玉米品种的关键候选基因、通路和代谢物。
论文解读：玉米应对草地贪夜蛾的整合转录代谢防御机制
研究背景
草地贪夜蛾（FAW）作为一种高度多食性的入侵害虫，对全球玉米生产构成严重威胁。目前针对FAW的防治过度依赖化学杀虫剂，导致害虫产生抗药性及杀伤天敌等生态问题。宿主植物抗性（Host Plant Resistance, HPR）作为一种可持续的治理策略日益受到重视，但现代玉米育种因侧重产量提升而导致许多防御性状丢失。尽管已鉴定出部分抗性种质，但其背后的分子机制尚不明确。为此，研究人员在《Plant Stress》发表研究，旨在通过整合转录组学和代谢组学手段，解析玉米对FAW早期防御反应的分子基础，挖掘关键抗性基因与代谢物。
关键技术方法
研究人员选用CIMMYT开发的抗性自交系CML338和感性自交系CML451，在温室条件下设置虫害侵染与未侵染对照，于侵染后48小时采集叶片样品。研究结合Illumina NovaSeq 6000平台的转录组测序（RNA-seq）与非靶向液相色谱-高分辨质谱（LC-HRMS）代谢组学技术，分别获取基因表达谱与代谢物谱。通过加权基因共表达网络分析（Weighted Gene Co-expression Network Analysis, WGCNA）、皮尔逊相关性分析（PCC）整合多组学数据，并利用同源建模与分子对接技术预测关键代谢物与FAW肠道酶的互作潜力。
研究结果
3.1. 玉米基因型的表型再确认
研究人员通过无选择人工接虫实验证实，抗性系CML338的叶片损伤评级（4.8）显著低于感性系CML451（6.2），确立了后续分子分析的遗传材料基础。
3.2. 转录组分析
高质量测序数据的比对效率达84–87%。主成分分析（PCA）显示，PC1（36.8%）区分基因型，PC2（23.7%）区分处理条件，表明遗传背景与虫害胁迫共同驱动了转录组的变异。
3.3. 转录组多变量分析
维恩图分析发现四个比较组间共享29,144个DEGs，同时存在特定条件的独有DEGs。RP vs SP比较组中DEGs数量最多，表明抗性系在受虫害诱导后发生了剧烈的转录重塑。热图与聚类分析进一步揭示了基因表达的动态变化模式。
3.4. 差异表达基因富集分析与通路分类
功能注释显示，抗性相关组成型（Resistance-related constitutive, RRC）基因主要富集于植物激素信号转导、MAPK信号通路和次生代谢物生物合成。抗性相关诱导型（Resistance-related induced, RRI）基因则显著富集于氧化还原酶活性与膜组分相关通路，提示了活性氧（Reactive Oxygen Species, ROS）爆发与细胞膜重构在防御中的重要性。
3.5. 代谢组学分析
代谢组PCA显示PC1（45%）区分基因型，PC2（29.1%）区分处理。鉴定出的RRC代谢物主要涉及脂肪酸类，而RRI代谢物则在虫害诱导后显著增加。关键代谢物包括亚油酸、花生四烯酸以及多种类胡萝卜素衍生物（如5-羟基玉米黄质、辣椒玉红素、紫黄质等）。
3.6. 代谢组学WGCNA
研究人员构建了十个代谢物共表达模块，其中蓝色模块与抗性表型呈最强正相关（r = 0.982）。该模块包含脂肪酸酰基、二萜生物碱、聚酮化合物和萜类化合物，这些化合物均已知具有抗草食作用。
3.7. 转录组与代谢组的整合分析
Sankey图及相关性网络分析揭示了脂肪酸代谢与萜类代谢的高度协调性。研究人员发现TPS、DXS、PLD、LOX和FAD等关键基因的表达与相应代谢物的积累呈强正相关（PCC ≥ 0.85）。
3.8. 分子对接
分子对接结果显示，α-松油醇与FAW中肠氨基肽酶N1（Aminopeptidase N1）的结合能最低（-6.6 kcal/mol），表明其可能通过干扰昆虫消化生理发挥抗性作用，其次是氨基肽酶N和碱性磷酸酶（Alkaline Phosphatase, ALP）。
3.9. qRT-PCR验证
定量逆转录聚合酶链式反应（qRT-PCR）验证了脂肪酸、类黄酮和萜类通路中关键候选基因（如AAMT1、PLD、TPS、DFR）在抗性系中显著上调，与转录组数据高度一致。
讨论与结论
研究人员在讨论中指出，玉米对FAW的防御是一个涉及受体识别、钙信号、MAPK级联、转录因子调控及激素串扰的复杂网络。RLKs、CDPKs及WRKY、MYC等转录因子的激活，驱动了下游脂质衍生信号（如JA）和萜类、黄酮类次生代谢物的生物合成。特别是α-松油醇等单萜类物质，不仅作为直接防御毒素，还可能通过抑制昆虫中肠APN1酶活性干扰营养吸收。该研究从系统生物学层面阐明了玉米–FAW互作的分子机制，鉴定的关键基因（如TPS、LOX、PLD）和代谢物（如α-松油醇、紫黄质）为抗虫玉米的分子标记辅助育种（Marker-Assisted Selection, MAS）提供了重要的靶标资源。
