LncRNA驱动的免疫网络桥接PTI（Pattern-Triggered Immunity，模式触发免疫）与ETI（Effector-Triggered Immunity，效应子触发免疫）应答以增强Vigna mungo（黑绿豆）对MYMIV（Mungbean Yellow Mosaic India Virus，印度豇豆黄花叶病毒）的抗性

《Polar Science》：LncRNA-Driven Immune Networks Bridge PTI and ETI Responses to Reinforce MYMIV Resistance in Vigna mungo

　　

Mungbean Yellow Mosaic India Virus（MYMIV，印度豇豆黄花叶病毒）诱导的黄化花叶病是一种极具破坏性的病害，严重制约Vigna mungo（黑绿豆）生产。尽管宿主防御涉及大规模转录重编程，但lncRNA（long non-coding RNA，长链非编码RNA）在塑造这种多层抗性中的调节作用仍鲜有报道。本研究首次报道了MYMIV侵染的抗性（VM-R）和感性（VM-S）基因型中lncRNA的综合图谱，揭示了与相反病害耐受表型相关的不同分子特征。比较分析发现VM-R和VM-S中分别有637和574个差异表达lncRNA（DEL，Differentially Expressed LncRNA），主要以cis方式分别调控858和932个基因，少数以trans方式分别作用于360和313个靶标。功能富集显示不同模式：VM-R优先激活SA（Salicylic Acid，水杨酸）信号、ROS（Reactive Oxygen Species，活性氧）稳态及过敏性反应（HR，Hypersensitive Response），桥接PTI（Pattern-Triggered Immunity，模式触发免疫）与ETI（Effector-Triggered Immunity，效应子触发免疫）以强化抗病毒防御；而VM-S表现出广泛的代谢重编程。加权基因共表达网络分析（WGCNA，Weighted Gene Co-expression Network Analysis）将DEL划分为27个共表达模块，其中MEgreen和MEpurple模块分别与VM-R和VM-S显著关联，突显基因型特异性共表达框架。竞争性内源RNA（ceRNA，Competing Endogenous RNA）网络揭示独特调控架构：VM-R中7个lncRNA作为靶标模拟物吸附7个miRNA以调控11个靶基因，VM-S中调控7个靶mRNA。研究人员提出，VM-R中ceRNA对VmLnc.2876.15–miR5675激活TIR-NB-LRR介导防御，而VM-S中VmLnc.3702.1表达降低无法吸附miR482，导致RING/U-Box水平下降并损害免疫。机器学习模型识别出具有区分抗感互作强预测潜力的DEL及ceRNA模块。qPCR验证了关键DEL、靶基因及相关ceRNA元件，强调其在lncRNA介导免疫调节中的作用。综上，这些发现阐明lncRNA在协调MYMIV抗性中的机制作用，为提升作物抗逆性与产量提供策略。

本研究聚焦黑绿豆（Vigna mungo）感染MYMIV（Mungbean Yellow Mosaic India Virus，印度豇豆黄花叶病毒）后lncRNA（long non-coding RNA，长链非编码RNA）介导的免疫调控机制，首次在豆科作物中系统解析了lncRNA通过桥接PTI（Pattern-Triggered Immunity，模式触发免疫）与ETI（Effector-Triggered Immunity，效应子触发免疫）协同增强抗病毒抗性的分子网络。目前植物病毒抗性研究中lncRNA的功能多限于模式物种，豆科作物中lncRNA如何调控MYMIV抗性及PTI–ETI协同机制尚属空白。研究人员整合抗性品种VM84（VM-R）与感性品种T9（VM-S）的MYMIV侵染转录组，鉴定MYMIV响应lncRNA并构建其cis/trans靶标、WGCNA（Weighted Gene Co-expression Network Analysis，加权基因共表达网络分析）模块及ceRNA（Competing Endogenous RNA，竞争性内源RNA）网络，结合机器学习筛选抗性标志lncRNA并用qPCR验证。结论表明抗性基因型中特异诱导的lncRNA通过顺式调控免疫相关转录因子（TF，Transcription Factor）、NLR（Nucleotide-Binding Leucine-Rich Repeat，核苷酸结合亮氨酸富集重复受体）类抗病基因及ceRNA机制吸附miRNA解除对TIR-NB-LRR等的抑制，从而协同激活PTI与ETI通路并维持ROS稳态与HR（Hypersensitive Response，过敏性反应），为豆科抗病毒分子育种提供新靶点。本文解读依据原文结果浓缩总结。

主要关键技术方法

研究人员选用已发表的MYMIV侵染黑绿豆抗性（VM84，VM-R）与感性（T9，VM-S）基因型RNA-seq数据集（含NCBI BioProject: PRJNA283940、PRJNA288413、PRJNA789694），经FastQC与fastp质控后用HISAT2比对至Vigna mungo参考基因组（ASM1342719v1），StringTie进行参考引导de novo组装并由Cuffmerge整合注释。用FEELnc联合CPC2、CNCI及Rfam过滤得到高置信lncRNA并分类（lincRNA、lncNAT、genic、intronic）。DESeq2筛选|log₂FC|≥2且padj<0.05的差异表达lncRNA（DEL）。cis靶基因取lncRNA上下游100kb内蛋白编码基因（PCG，Protein-Coding Gene），trans靶用LncTar预测并经PCC（Pearson Correlation Coefficient）|PCC|≥0.7过滤。GO与KEGG富集用AgriGO v2.0与KAAS。WGCNA包构建共表达模块并关联表型。psMimic与psRobot预测lncRNA–miRNA–mRNA的ceRNA关系。TensorFlow/Keras构建ANN（Artificial Neural Network，人工神经网络）模型并用SHAP解析特征重要性。温室接种MYMIV后提取总RNA，qPCR以VmACT为内参验证DEL及靶基因表达。

3.1. Systematic profiling uncovers MYMIV-responsive lncRNA landscape in Vigna mungo

经多步过滤从93,198条转录本中获得4,746条高置信lncRNA（2,976条lincRNA、1,009条lncNAT、651条genic lncRNA、110条intronic lncRNA），均匀分布于11条染色体；其中1,309条（27.5%）与CANTATAdb/GreeNC中已知lncRNA同源，Vigna属内保守度高于其他豆科。lncRNA表达量（TPM）低于mRNA、GC含量（37.4%）较低、长度多<300nt且多含2个外显子，符合植物lncRNA典型特征。

3.2. Structural and genomic features define the Vigna mungo lncRNA repertoire

（此部分与上段合并描述结构特征，原文无独立新结论小节）

3.3. Differential expression of lncRNAs underpins genotype-specific responses to MYMIV

MYMIV侵染vs对照鉴定VM-R中637个DEL（336上调、301下调），VM-S中574个DEL（213上调、361下调）；仅少量重叠（36个共同上调/下调），40个DEL在VM-R上调而在VM-S下调，表明lncRNA重编程具强基因型特异性。DEL主要源自基因间区，其次为反义链。

3.4. Differentially expressed lncRNAs (DELs) modulate target genes via Cis-Trans mechanisms

预测VM-R中569个DEL以cis调控858个mRNA，68个DEL以trans调控360个mRNA，形成1,238个lncRNA–mRNA对；VM-S中544个DEL cis调控932个mRNA，30个DEL trans调控313个mRNA，形成1,206对；少数PCG受cis+trans双重调控，说明lncRNA具多层次基因调控能力。

3.5. Functional enrichment of DEL-targets reveals lncRNA-mediated immune networks

VM-R中DEL靶基因显著富集于"response to biotic stimulus"（GO:0009607）、"immune response"（GO:0006955）、SA信号、ROS稳态、MAPK信号（KO:04016）、植物–病原互作（KO:04626）及苯丙烷/类黄酮合成；VM-S中靶基因富集于光合作用（GO:0009765）、碳代谢（KO:01200）、淀粉与蔗糖代谢（KO:00500），表明感性品种发生代谢扰动而非免疫激活。

3.6. WGCNA uncovers regulatory lncRNA modules driving MYMIV-response

27个共表达模块中MEgreen（VM-R关联，PCC>0.5）含105个DEL与709个DEG（Differentially Expressed Gene），富集PTI/ETI及激素信号；MEpurple（VM-S关联）含29个DEL与235个DEG，富集基础代谢。MEgreen中VmLnc.22558.3预测cis调控CC-NB-LRR，VmLnc.18958.1上调RLK（Receptor-Like Kinase，类受体激酶），多个DEL共调控WRKY、NAC、bHLH、TGA7等免疫TF及PR5、TIR-NB-LRR等抗病组分。

3.7. Artificial neural network (ANN) based modelling identifies lncRNA signatures of MYMIV resistance

基于MEgreen/MEpurple模块DEL表达的ANN分类模型AUC–ROC达0.88；SHAP分析显示VmLnc.33854.1（靶向WD40）、VmLnc.423.2（靶向GDSL lipase）、VmLnc.42513.6（靶向WRKY）等为抗感分型关键标志lncRNA。

3.8. qPCR validation confirms lncRNA–target expression dynamics

选取9个DEL进行qPCR，表达趋势与RNA-seq一致：VM-R中VmLnc.42248.9/42513.6/42403.2/28498.3上调并伴随NAC、WRKY、bHLH、TGA7升高；VmLnc.37199.1/839.4/18958.1诱导CC-NB-LRR、STK（Serine/Threonine Kinase，丝氨酸/苏氨酸激酶）、RLK；VmLnc.41325.1与PR5、VmLnc.423.2与GDSL lipase共上调；VmLnc.33854.1在VM-S中显著下调。

3.9. Regulatory ceRNA interactions modulate host immunity against MYMIV

VM-R中MEgreen模块7个lncRNA作为ceRNA吸附miR166/miR426/miR866/miR4397/miR5035/miR5658/miR5675，其中VmLnc.2867.15吸附miR5675解除对TIR-NB-LRR（TNL，TIR-NB-LRR class NLR）的抑制；VM-S中VmLnc.3702.1低表达无法吸附miR482致使RING/U-Box E3 ligase被沉默。ceRNA模块训练的ANN模型AUC–ROC=0.722，qPCR证实VmLnc.2867.15诱导伴随TNL上调，VmLnc.3702.1下调伴随RING/U-Box下调。

讨论部分总结（翻译研究结论）

研究人员提出lncRNA是Vigna mungo抗病毒免疫的核心整合因子，在抗性基因型中MYMIV侵染协调诱导cis与trans作用lncRNA，通过调控PRR（Pattern-Recognition Receptor，模式识别受体）/NLR及下游防御组分（WRKY、NAC、bHLH、STK、GDSL lipase、PR蛋白）强化PTI与ETI，促进SA信号、ROS爆发及HR；同时ceRNA网络缓冲miRNA对TIR-NB-LRR及RING/U-Box等的抑制以维持免疫基因表达。感性基因型缺失该lncRNA-ceRNA调控致免疫相关转录受抑及代谢重编程偏向碳代谢，削弱防御。本研究首次在豆科植物中证明lncRNA参与植物–病毒抗性，所鉴定的DEL与ceRNA模块为黑绿豆抗MYMIV分子标记辅助育种及基因组改良提供理论与候选靶点。需注意ceRNA互作为计算预测假设模型，尚需进一步功能验证。