蜜蜂球囊菌菌丝和孢子中毒力因子相关全长转录本的鉴定及分析

目的]蜜蜂球囊菌Ascosphaera apis是一种专性侵染蜜蜂幼虫而导致白垩病的致死性真菌病原.基于前期已获得的高质量纳米孔(Nanopore)长读段测序数据,对蜜蜂球囊菌菌丝(AaM)和孢子(AaS)中的毒力因子相关全长转录本进行鉴定和分析,为毒力因子相关剪接异构体的功能研究提供参考信息和基础.方法]利用BLAST工具,将AaM和AaS中的所有全长转录本比对Nr数据库,以鉴定蜜蜂球囊菌的毒力因子(几丁质酶、脂肪酶、水解酶和蛋白酶)相关的全长转录本.使用minimap2软件,将AaM和AaS中的全长转录本序列与蜜蜂球囊菌参考基因组注释的已知转录本序列进行比对,将比对到参考基因组的全长转录本进行归一化处理,再通过每百万计数法(Counts per million,CPM)计算毒力因子相关全长转录本的表达量.利用百迈克云平台的相关工具绘制转录本的表达量聚类热图.通过IGV浏览器对部分毒力因子相关全长转录本结构进行可视化.结果]在AaM鉴定到毒力因子相关的367个基因及407个全长转录本,包括12条几丁质酶相关全长转录本,48条脂肪酶相关全长转录本,289条水解酶相关全长转录本,58条蛋白酶相关全长转录本.在AaS鉴定到毒力因子相关367个基因及400个全长转录本,包括14条几丁质酶相关全长转录本,63条脂肪酶相关全长转录本,267条水解酶相关全长转录本,56条蛋白酶相关全长转录本.另外,AaM和AaS特有的毒力因子(几丁质酶、脂肪酶)相关全长转录本分别有0条和17条,共有的毒力因子相关全长转录本有60条.进一步分析发现,蜜蜂球囊菌的部分毒力因子基因可通过可变剪接形成多条剪接异构体.结论]共鉴定到蜜蜂球囊菌毒力因子(几丁质酶、脂肪酶、水解酶和蛋白酶)相关的367个基因和486条全长转录本;相比于蜜蜂球囊菌参考基因组注释的转录本,绝大多数毒力因子基因对应的全长转录本数量更多且结构更为复杂.研究结果丰富了蜜蜂球囊菌毒力因子相关基因和转录本的注释信息,为毒力因子相关剪接异构体的功能研究提供了基础,也为白垩病防控提供了潜在靶点.