小麦条锈病抗性相关基因间长非编码RNA（lincRNA）的鉴定与表达分析

《Frontiers in Plant Science》：Identification and expression analysis of intergenic long non-coding RNAs (lincRNA) associated with stripe rust resistance in wheat

研究背景：小麦（Triticum aestivum）作为重要作物，其生产长期受条锈菌（Puccinia striiformis f. sp. tritici, Pst）严重威胁。尽管已有研究关注蛋白编码基因在免疫中的作用，但基因间长非编码RNA（lincRNA）在该病理系统中的功能景观仍不明确。本研究旨在通过全基因组水平鉴定小麦叶片在条锈菌侵染过程中的lincRNA，并探讨其潜在的调控机制。研究方法：研究人员利用链特异性RNA-seq数据，对小麦品种铭贤169接种条锈菌后0、24、48、120小时共四个时间点的叶片转录组进行分析。通过严格的生物信息学流程（包括转录本组装、编码潜能评估、结构域筛选等）鉴定高可信度lincRNA。进一步采用差异表达分析、功能富集、共表达分析、lincRNA-miRNA-mRNA网络预测以及基于AlphaFold的蛋白-RNA结构预测等手段，解析候选lincRNA的调控特征。研究结果：共鉴定出2,736个表达的lincRNA，其中515个在侵染期间差异表达。差异表达lincRNA邻近100 kb范围内的基因显著富集于防御相关通路，如MAPK信号通路和植物-病原互作通路。共表达分析揭示了lincRNA与邻近防御相关mRNA的协调表达模式，提示可能存在顺式调控关联。lincRNA-miRNA-mRNA网络预测将65个差异表达lincRNA鉴定为21个保守miRNA的潜在内源靶标模拟子（endogenous target mimics），展示了广泛的竞争性内源RNA（ceRNA）调控网络。尤为重要的是，基于AlphaFold的蛋白-RNA结构预测表明，共表达的lincRNA MSTRG.56980能与RGA1型抗病蛋白TraesCS6B03G0179100形成稳定复合体。研究结论：该研究首次提供了小麦lincRNA与免疫受体之间可能存在物理相互作用的计算预测。这些整合性发现显著拓展了对小麦免疫中非编码调控网络的理解，并为优先筛选可提高条锈病抗性的候选lincRNA提供了基础资源。

论文解读：小麦条锈病抗性相关基因间长非编码RNA（lincRNA）的鉴定与表达分析

一、研究背景与问题

植物作为固着生物，必须不断适应复杂多变的环境以保障生长和繁殖。非编码RNA，尤其是长非编码RNA（lncRNA），在植物环境响应调控中扮演关键角色。根据基因组位置，lncRNA可分为基因间lincRNA、天然反义转录本（NAT）、内含子非编码RNA等类型。其中，lincRNA转录自蛋白编码基因之间的非编码区域，近年被证实是植物免疫应答的重要调控因子。例如，拟南芥中的ELENA1通过与中介体复合物亚基互作激活PR1转录从而增强对丁香假单胞菌的抗性；水稻中的ALEX1可激活茉莉酸介导的免疫。然而，在小麦与条锈菌（Puccinia striiformis f. sp. tritici, Pst）互作体系中，lincRNA的系统鉴定和功能机制研究仍十分有限。现有报道多集中于NLR受体和病程相关蛋白等蛋白编码组分，对lincRNA的分类、调控模式及其与mRNA网络的关联缺乏深入解析。由于条锈菌毒力小种进化迅速，揭示非编码调控层对于作物抗病育种具有迫切意义。此外，ceRNA假说和AlphaFold等工具的出现，为研究lincRNA-miRNA-mRNA互作网络及lincRNA与蛋白的空间互作提供了新途径，但这些方法在小麦-Pst体系中的系统应用尚属空白。

二、研究内容与结论

研究人员利用来自NCBI SRA数据库（登录号PRJNA975982）的小麦叶组织链特异性RNA-seq数据，聚焦感病品种铭贤169在接种条锈菌后0、24、48、120小时的转录组，开展了全基因组lincRNA鉴定。通过严格的过滤流程，共获得2,736个表达的lincRNA，其中515个在侵染期间差异表达。结构特征分析显示，小麦lincRNA比mRNA更短、外显子数更少、表达水平更低，符合植物lincRNA的典型特征。差异表达动态分析显示早期有适度诱导，后期呈现显著下调，表明可能存在免疫激活与负反馈调节的时序协调。功能富集分析表明，差异表达lincRNA邻近100 kb范围内的基因富集于MAPK信号通路、光合作用、植物-病原互作等防御相关通路。共表达分析发现lincRNA与邻近防御相关mRNA（如PR蛋白、转运蛋白）具有协调表达模式，并通过qRT-PCR验证了多个lincRNA-mRNA对的表达动态。ceRNA网络预测鉴定出65个差异表达lincRNA可作为21个保守miRNA的潜在内源靶标模拟子（eTM），呈现多对多的复杂调控架构。关键的是，AlphaFold蛋白-RNA结构预测显示，lincRNA MSTRG.56980可与RGA1型抗病蛋白TraesCS6B03G0179100形成稳定复合体，预测界面包含碱基配对和堆积互作。该研究首次提供了小麦lincRNA与免疫受体可能发生物理互作的计算证据，为理解小麦免疫中的非编码调控网络提供了重要资源。

三、主要关键技术方法

本研究采用以下核心技术手段：1）基于HISAT2和StringTie的基因组引导转录组组装，结合Gffcompare筛选“u”类基因间转录本；2）联合CPC2、CNCI、PLEK三种工具评估编码潜能，并通过HMMER、BLASTX及开放阅读框（ORF）长度过滤排除蛋白编码转录本；3）DESeq2差异表达分析（校正p值≤0.05，|log2倍数变化|≥2）；4）基于psRNATarget平台预测lincRNA-miRNA-mRNA互作网络，使用Cytoscape可视化；5）AlphaFold Server进行蛋白-RNA复合物结构预测，PyMOL分析结合界面；6）qRT-PCR验证候选lincRNA-mRNA对表达（内参为Ta2291）。样本队列来源为NCBI SRA PRJNA975982，包含铭贤169在0、24、48、120小时三个生物学重复的转录组数据；qRT-PCR验证使用扬麦158品种。

四、研究结果

lincRNA的系统鉴定与特征：通过对197,328个组装转录本进行分级筛选，鉴定出30,681个基因间转录本（占15.55%），经结构过滤、编码潜能筛选、蛋白结构域和同源比对过滤后，最终获得2,736个表达的lincRNA，其中515个在条锈菌侵染期间差异表达。lincRNA结构特征分析显示其外显子数较少、转录本较短且表达丰度较低，符合植物lncRNA的典型特征。

lincRNA差异表达动态：侵染早期（6小时）有12个上调、32个下调；12小时上调增至56个、下调18个；24小时出现显著下调（221个下调、38个上调）；6-24小时窗口内共检测到19个上调和11个下调，12-24小时区间变化最剧烈（317个下调、36个上调），表明lincRNA在侵染不同阶段呈现动态调控模式。

lincRNA邻近基因功能富集：834个邻近mRNA的GO分析显示，生物学过程主要富集于细胞过程和代谢过程，分子功能主要涉及结合和催化活性；KEGG分析显示核糖体、光合作用和D-氨基酸代谢是显著富集通路，同时包含MAPK信号通路和植物-病原互作等防御相关通路。

lincRNA-mRNA共表达验证：qRT-PCR结果表明MSTRG.18458与其邻近基因Ta2B1、Ta2B2、Ta2B3在侵染全程下调；MSTRG.34447邻近基因Ta3A1和Ta1D2在24小时降低后48小时回升；MSTRG.56980邻近基因Ta6B1和Ta3D1随处理时间表达增加，其中Ta6B1编码RGA1型抗病蛋白，含Rx N端结构域（PF18052）。

ceRNA网络预测：psRNATarget筛选发现108个miRNA与515个差异表达lincRNA存在潜在互作，最终65个差异表达lincRNA被预测为21个保守miRNA的潜在内源靶标模拟子，形成多对多的lincRNA-miRNA-mRNA调控网络。

蛋白-RNA复合物结构预测：AlphaFold预测MSTRG.56980与TraesCS6B03G0179100可形成稳定复合体，结合界面包含局部碱基配对和堆积互作，蛋白和RNA分别以蓝色和绿色显示，界面残基以红色突出，为lincRNA与NLR蛋白物理互作提供了结构证据。

五、讨论与结论总结

讨论部分指出，小麦lincRNA的结构与表达特征与植物lncRNA的普遍规律一致，其侵染后期的显著下调可能代表一种防止过度免疫激活的反馈机制。邻近基因富集于MAPK信号、光合作用及植物-病原互作通路，提示lincRNA可能参与协调免疫信号传导与代谢重编程。ceRNA网络揭示lincRNA可能作为miRNA海绵整合多条防御信号通路。AlphaFold预测的MSTRG.56980与抗病蛋白的稳定复合体，为该lincRNA参与NLR蛋白调控提供了结构层面的计算线索。值得注意的是，Ta6B1的KEGG自动注释K15078（SLX1核酸内切酶）与其含有的PF18052（Rx_N）结构域不一致，被判定为假阳性注释。研究结论强调：该研究系统鉴定了2,736个表达lincRNA，其中515个差异表达，揭示了lincRNA-mRNA顺式调控模块和ceRNA调控网络，并首次计算预测了小麦lincRNA与免疫受体的物理互作。这些发现为小麦条锈病抗性育种中候选lincRNA的功能验证奠定了资源基础。未来需要通过实验验证lincRNA与蛋白的直接结合、解析其调控NLR信号的分子机制，并将RNA策略应用于作物真菌病害防控。