基于RNA-seq的木霉长链非编码RNA的生物信息学预测及其重寄生相关性分析

目的]本研究旨在全基因组范围内识别木霉(Trichoderma guizhouense) NJAU 4742(NJAU 4742)的长链非编码RNA(long noncoding RNA,lncRNA),探究lncRNA在重寄生过程中可能参与的调控作用。方法]用链特异性RNA-seq技术,对重寄生过程中与病原菌互作接触前、后以及独立培养的木霉菌进行转录组测序;构建生物信息学流程,识别lncRNA并用RSEM和DESeq2软件分析lncRNA在重寄生过程中的表达情况;对lncRNA的靶标预测和表达情况进行分析,探索lncRNA在重寄生过程中可能参与的调控作用。结果]在木霉NJAU 4742中识别了1 676个lncRNA,其中包含1 049个基因间型lncRNA,590个反义型lncRNA,32个正义型lncRNA以及5个内含子型lncRNA。与编码基因相比,lncRNA的外显子数量偏少,序列长度偏短,表达量偏低,在基因组上的跨度偏短。靶标预测结果显示:1 496个lncRNA能够靶向2 269个蛋白编码基因,其中1 492个lncRNA以顺式作用形式靶向2 262个编码基因,4个lncRNA以反式作用形式靶向7个编码基因。GO功能分类结果显示:代谢过程(metabolic process)、催化活性(catalytic activity)和细胞过程(cellular process)是lncRNA靶标分布数量较多的3个类别。KEGG通路分析结果显示:信号转导(signal transduction)、转运与分解代谢(transport and catabolism)和碳水化合物代谢(carbohydrate metabolism)是lncRNA靶标分布数量较多的3类通路。进一步分析发现:147个lncRNA靶向编码碳水化合物活性酶和蛋白酶的基因以及次生代谢物合成相关的基因,其中30个lncRNA在重寄生过程中表达水平发生显著变化,有10个lncRNA的表达和靶标基因显著相关。结论]木霉在NJAU 4742中存在长链非编码RNA,部分成员参与对病原菌重寄生过程的调控。