转录组解析揭示小麦对小麦条纹花叶病毒耐受性的分子机制：激素调控与细胞壁强化协同作用

《BMC Genomics》：Transcriptomic profiling provides molecular insights into tolerance mechanisms in wheat to wheat streak mosaic virus (WSMV)

本研究针对小麦卷叶螨（WCM）及其传播的小麦条纹花叶病毒（WSMV）造成的严重产量损失问题，通过田间筛选和转录组分析，发现商用品种Hatcher在无已知抗性基因背景下表现出显著耐受性。研究揭示其通过协调JA、SA、ABA激素通路激活、木质素生物合成途径上调及生长-防御平衡调控实现耐受，为培育持久抗性小麦品种提供了新靶点。

在全球小麦生产中，小麦卷叶螨（Wheat curl mite, WCM）及其传播的病毒病害构成严重威胁。这些微小的螨虫不仅通过直接取食危害小麦，更是小麦条纹花叶病毒（Wheat streak mosaic virus, WSMV）等主要病毒的关键传播媒介。在美国大平原地区，由WCM-WSMV复合体造成的产量损失平均达2-3%，严重田块甚至可高达100%。目前缺乏有效的化学防治手段，主要依赖清除自生小麦、延迟播种以及使用抗性品种等综合管理策略。然而，随着强毒力WCM种群和抗性突破型WSMV分离株的出现，开发更多样化、更持久的抗性来源变得尤为迫切。

传统抗性基因（如Wsm1、Wsm2）的应用常伴随产量损失、连锁累赘及高温下有效性降低等限制。与抗性（resistance）不同，耐受性（tolerance）允许植物在受侵染情况下仍能维持生长和产量，而不降低害虫或病原体水平，这为作物保护提供了新思路。在此背景下，Pingault等人在《BMC Genomics》上发表的研究，对一种名为Hatcher的商业小麦品种展开了深入探究。该品种不含已知的WCM或WSMV抗性基因，却在连续两年的田间筛选中，于高WSMV病害压力下 consistently（持续地）表现出更高的产量，甚至优于携带Wsm1或Wsm2抗性基因的品种。这一反常现象暗示Hatcher可能具有对WCM和WSMV的耐受性。为了揭示其潜在的分子机制，研究人员比较了Hatcher与感病基因型CO15D173R对WCM和WSMV侵染的响应，并进行了转录组学分析。

为开展研究，研究人员主要运用了以下关键技术方法：首先，在田间进行为期两年的随机区组设计试验，对五个小麦品种（包括Hatcher）进行自然病害压力下的筛选，评估病毒发生率（通过DAS-ELISA检测）和产量。其次，在温室控制条件下进行因子设计实验，比较Hatcher和感病系CO15D173R在接种带毒WCM后的农艺性状（分蘖数、株高、穗数、百粒重）。核心的分子生物学技术包括从叶片组织提取总RNA后进行RNA测序（RNA-seq），利用Illumina HiSeq平台进行双端测序。生物信息学分析涉及 reads 的质量控制、比对至中国春小麦参考基因组（v2.1）、转录本重构、差异表达基因（DEGs）鉴定（Cuffdiff软件），以及后续的基因本体（GO）富集分析、KEGG通路分析、加权基因共表达网络分析（WGCNA）和转录因子分析。

Hatcher在高WSMV发病率下保持高产

田间试验结果显示，2015和2016年WSMV发病率在25%至93%之间。尽管Hatcher的病毒发病率在2016年显著高于携带Wsm1和Wsm2抗性基因的品种（Mace和Snowmass 2.0），但其两年平均产量最高，与Snowmass 2.0相当，并显著高于Mace和TAM112（携带抗螨基因Cmc^TAM112）。温室实验进一步证实，在WCM+WSMV侵染下，感病基因型CO15D173R的分蘖数和百粒重显著降低，而Hatcher的分蘖数保持不变，且百粒重反而显著增加，尽管侵染植株的株高有所降低，这表明Hatcher能够将资源从营养生长向生殖生长重新分配，维持关键产量构成因素。

Hatcher和CO15D173R对WCM和WSMV的转录组学响应

RNA测序分析表明，基因型是转录组变异的主要驱动力。与CO15D173R相比，Hatcher对侵染的响应更为强烈，鉴定出5191个差异表达基因（DEGs），且以上调基因为主；而CO15D173R仅鉴定出2814个DEGs，且以下调为主。主成分分析（PCA）显示，侵染处理在Hatcher中引起了显著的转录组偏移。

功能富集分析揭示对WCM和WSMV侵染的不同分子响应

GO富集分析显示，Hatcher中上调的基因显著富集于典型的植物防御反应功能，如苯丙烷类代谢（苯丙氨酸解氨酶PAL活性）、单加氧酶活性、UDP-葡萄糖基转移酶活性、钙离子结合等。而下调基因则富集于叶绿体结构和光合作用相关功能，表明资源从初级代谢向防御的重新分配。相反，CO15D173R中上调基因未富集典型防御类别，而下调基因则显示出光合作用和转录机制核心成分的广泛抑制，暗示其细胞稳态的崩溃。

代谢通路分析揭示对WCM和WSMV侵染响应的差异

MapMan可视化分析表明，Hatcher广泛激活了与结构和细胞防御相关的基因，包括细胞壁生物合成、β-葡聚糖酶、热激蛋白（HSPs）、过氧化物酶、谷胱甘肽-S-转移酶以及JA和ABA等激素通路。而CO15D173R虽然强烈上调了JA、ABA、乙烯（ET）和生长素途径的基因，但许多核心防御通路（如氧化还原调节、次生代谢）却被下调。

木质素生物合成在Hatcher中受WCM和WSMV侵染后转录激活

在分析的221个木质素生物合成相关基因中，有61个在侵染后差异表达，其中50个在Hatcher中特异性上调，覆盖了从p-香豆酰辅酶A到木质素聚合物合成的几乎所有主要酶促步骤。这表明细胞壁强化是Hatcher防御策略的重要组成部分。CO15D173R中这些基因的调控则微乎其微。

KEGG通路分析突出Hatcher中的防御激活和CO15D173R中的代谢紊乱

KEGG通路分析进一步证实了基因型间的差异。在核心代谢通路（碳水化合物、脂质、氨基酸、能量代谢）中，CO15D173R相关基因普遍下调，而Hatcher则基本维持稳定。激素信号通路方面，Hatcher中水杨酸（SA）响应基因（如NPR1）和α-亚麻酸通路（JA生物合成前体）基因被诱导，PR蛋白基因上调，ABA受体PYL基因下调。CO15D173R则表现出SA响应基因抑制、PYL基因上调以及MAPK信号通路中多个防御调节因子（如CNGCs, CDPK, RPM1, HSP90）的下调。这表明Hatcher协调激活了SA和JA介导的信号传导，并维持核心代谢，而CO15D173R则出现代谢崩溃和防御激活不完全的特征。

转录因子的差异调控凸显基因型特异性响应

共鉴定出162个差异表达的转录因子（TFs），涉及36个家族。Hatcher中有70个TFs上调，50个下调，表现出强烈的转录重编程。其中NAC家族TFs（28个基因）响应显著，14个NAC TFs在Hatcher受侵染后特异性高表达。而CO15D173R的TF响应有限。这提示NAC TFs在Hatcher协调防御基因激活中可能起核心作用。

层次聚类揭示与基因型特异性响应相关的不同表达模块

对7394个DEGs进行层次聚类得到9个簇，揭示了基因型特异性和处理响应性的表达模式。例如，簇7中的基因仅在受侵染的Hatcher中高表达，且富集于防御反应、免疫系统过程等功能。簇8中的基因在除受侵染CO15D173R外的所有处理中均表达，其功能涉及生物刺激响应、代谢等，该簇在CO15D173R中的抑制暗示其生物胁迫下代谢和信号功能的受损。

综上所述，本研究通过整合田间表型与转录组数据，系统揭示了小麦品种Hatcher对WCM和WSMV复合胁迫的耐受性分子机制。其主要特征在于：1）精细调控的激素交叉对话，优先表达ABA和SA信号成分，同时约束JA信号的过度激活；2）NAC转录因子和MAPK信号通路的协同激活，驱动下游防御基因（如PR蛋白）的表达；3）木质素生物合成途径的广泛上调，可能通过强化细胞壁物理屏障来阻碍螨虫定殖和病毒传播；4）维持生长相关激素（如赤霉素GA和油菜素内酯BR）通路的活性，实现生长与防御的平衡。这些多层次、协调的防御网络使得Hatcher能够在胁迫下维持代谢稳定和产量形成。该研究不仅为理解植物对虫媒病毒复合胁迫的耐受机制提供了新的见解，而且鉴定出的关键通路和候选基因（如NAC TFs、DnaK分子伴侣、木质素合成相关基因）为通过分子育种培育具有广谱、持久耐受性的小麦新品种奠定了重要的理论基础。