意大利蜜蜂工蜂中肠piRNA的鉴定与分析

【目的】Piwi蛋白互作RNA(PIWI-interacting RNA, piRNA)在昆虫的发育和免疫等重要生物学过程发挥重要的调控作用。本研究旨在丰富西方蜜蜂Apis mellifera的piRNA信息,并为进一步探究差异表达piRNA(differentially expressed piRNA, DEpiRNA)调控意大利蜜蜂Apis mellifera ligustica工蜂中肠发育的分子机理提供基础。【方法】基于已获得的意大利蜜蜂7日龄(Am7)和10日龄(Am10)工蜂中肠的small RNA(sRNA)组学数据,对piRNA进行预测和分析。将质控后的sRNA数据比对西方蜜蜂参考基因组,再将比对上的序列标签(tags)进一步比对数据库以滤除rRNA和tRNA等小非编码RNA(non-coding RNA,ncRNA),进而根据piRNA的长度特征鉴定piRNA。采用TPM(tags per million)算法对piRNA表达量进行计算和归一化处理。根据|log₂ fold change|≥1且P≤0.05的标准筛选Am7 vs Am10 比较组的DEpiRNA。通过相关软件预测DEpiRNA的靶mRNA并进行GO和KEGG数据库注释。利用Cytoscape软件对DEpiRNA-mRNA调控网络进行可视化。通过Stem-loop RT-PCR对随机选取的Am7与Am10两组共有的6个piRNA的表达进行验证。利用RT-qPCR对随机挑选的6个DEpiRNA的表达趋势进行验证。【结果】从意大利蜜蜂工蜂中肠中共鉴定到596个piRNA,长度介于24~33 nt,且Am7和Am10组中不同长度分布范围的piRNA的首位碱基偏向性具有明显差异。在Am7 vs Am10比较组共筛选出41个DEpiRNA,其中piR-ame-11093, piR-ame-1111451,piR-ame-190949和piR-ame-932156可分别靶向1 195, 1 018, 4 040和1 063条mRNA。上述靶mRNA可分别注释到45个功能条目和45条通路。Stem-loop RT-PCR验证结果显示6个piRNA(piR-ame-1084826, piR-ame-11093, piR-ame-14476, piR-ame-24995, piR-ame-39500和piR-ame-774987)均真实表达。RT-qPCR结果显示共有6个DEpiRNA (piR-ame-1084826,piR-ame-11093, piR-ame-14476, piR-ame-24995, piR-ame-39500和piR-ame-774987)的表达趋势与测序数据中的表达趋势一致,证实了本研究中sRNA-seq数据的可靠性和piRNA差异表达趋势的真实性。【结论】本研究从意大利蜜蜂工蜂中肠中鉴定到596个piRN,长度介于24~33 nt。不同长度分布范围的piRNA的首位碱基偏向性具有明显差异;piR-ame-11093, piR-ame-1111451,piR-ame-190949和piR-ame-932156具有靶向调控基因表达参与工蜂中肠发育过程的潜力。