本研究旨在基于已获得的第三代纳米孔全长转录组数据对东方蜜蜂微孢子虫Nosema ceranae基因的可变剪接(alternative splicing,AS)和可变多聚腺苷酸化(alternative polyadenylation,APA)进行分析.通过Astalavista软件鉴定东方蜜蜂微孢子虫基因的AS事件类型.采用TAPIS pipeline对东方蜜蜂微孢子虫基因的APA位点进行鉴定.利用MEME软件分析所有全长转录本的poly(A)剪接位点上游50bp的序列特征并鉴定motif.共鉴定出发生于5个基因的5次AS事件,包括1次可变供体位点(alternative donor site,ADS)和4次内含子保留(intron retention,IR).鉴定出233个基因发生APA,其中含有1个poly(A)剪接位点的基因数量最多(143,61.37％).序列特征分析结果显示,东方蜜蜂微孢子虫全长转录本的3'UTR的上下游表现出明显的碱基倾向性,U在3'UTR的上游富集,而A在3'UTR的下游富集.此外,在东方蜜蜂微孢子虫全长转录本的poly(A)剪接位点上游鉴定到3个motif(AAUAAA、UGAUGC和GCGACG).研究结果揭示了东方蜜蜂微孢子虫转录组的复杂性,完善了现有的参考基因组和转录组注释,也为深入探究不同剪接异构体的分子功能提供了基础.此外,本研究为其他真菌的AS和APA研究提供了有益的思路和方法借鉴.

作者：陈华枝；范小雪；范元婵；王杰；祝智威；蒋海宾；张文德；隆琦；熊翠玲；郑燕珍；付中民；徐国钧；陈大福；郭睿

来源：菌物学报 2021 年 40卷 1期