[目的]蜜蜂球囊菌(Ascosphaera apis)是一种专性侵染蜜蜂幼虫的致死性真菌病原.本研究旨在利用PacBio单分子实时(single molecule real-time,SMRT)测序技术对蜜蜂球囊菌孢子(AaS)中基因的可变剪切(alternative splicing,AS)和可变多聚腺苷酸化(alternative polyadenylation,APA)以及长链非编码RNA(long non-coding RNA,lncRNA)进行鉴定和分析,进而揭示蜜蜂球囊菌孢子中转录组的复杂性.[方法]采用Suppa软件对蜜蜂球囊菌孢子中基因的AS事件进行鉴定.通过RT-PCR对不同类型的AS事件进行验证.采用TAPIS pipeline对蜜蜂球囊菌孢子基因的APA位点进行鉴定.利用MEME软件分析孢子全长转录本的poly(A)剪接位点上游50 bp的序列特征并鉴定motif.联用CPC和CNCI软件和比对Swiss-prot数据库的方法预测lncRNA,取三者的交集作为高可信度的lncRNA集合.进一步比较lncRNA和mRNA的转录本长度,外显子数量与长度,内含子长度,GC含量,AS事件数量.[结果]在蜜蜂球囊菌孢子中共鉴定到2 609次AS事件,包括 1 227 次(47.03％)内含子保留(retained intron,RI),842 次(32.27％)可变 3'剪切(alternative 3'splice sites,A3),415 次(15.91％)可变 5'剪切(alternative 5'splice sites,A5),85 次(3.26％)可变起始外显子(alternative first exons,AF),35 次(1.34％)外显子跳跃(skipping exon,SE),4 次(0.15％)可变末端外显子(alternative last exons,AL)和 1 次(0.04％)互斥外显子(mutually exclusive exons,MX).通过RT-PCR证实了不同类型AS事件(RI、A5、SE和A3)的可靠性.共鉴定出5 552个基因含APA位点,其中含有＞5个APA位点的基因数量最多,达到2197个,其次为含有1个APA位点的基因,数量为1 149个;此外,含有2、3、4和5个APA位点的基因数量分别为782、596、477和351个.蜜蜂球囊菌孢子中全长转录本的上下游序列具有明显的碱基倾向性,U和A分别在转录本的下游和上游富集.共鉴定到953条lncRNA.与mRNA相比,上述lncRNA的外显子数量更少且长度更短,内含子长度更短,转录本长度更短,GC含量更低,AS事件数量更少.[结论]本研究全面解析了 AaS中的AS、APA和lncRNA,研究结果丰富了蜜蜂球囊菌的基础生物学信息,深入揭示了蜜蜂球囊菌孢子中转录组的复杂性,并为进一步探究不同剪接异构体在病原孢子和病原侵染中的分子功能打下了基础.

蜜蜂球囊菌;PacBio测序;可变剪接;可变多聚腺苷酸化;长链非编码RNA;剪接异构体

孙明会;刘佳美;王思懿;朱乐冉;王紫馨;叶亚萍;钱加珺;顾小雨;徐细建;陈大福;郭睿;徐国钧

福建农林大学动物科学学院(蜂学学院),福建福州350002;江西省养蜂研究所,江西南昌330000;福建农林大学蜂疗研究所,福建福州350002

微生物学报

[1]孙明会;刘佳美;王思懿;朱乐冉;王紫馨;叶亚萍;钱加珺;顾小雨;徐细建;陈大福;郭睿;徐国钧-.结合三代测序与二代测序技术揭示蜜蜂球囊菌孢子转录组的复杂性)[J].微生物学报,2022(08):2981-2994

蜜蜂球囊菌,AaS,Suppa,TAPIS

三代测序,二代测序技术,孢子,转录组,Ascosphaera,apis,侵染,蜜蜂幼虫,死性,真菌病,PacBio,单分子,子实,single,molecule,real,SMRT,可变剪切,alternative,splicing,AS,可变多聚腺苷酸化,polyadenylation,APA,长链非编码,long,coding,lncRNA,pipeline,MEME,全长,转录本,bp,序列特征,motif,CPC,CNCI,Swiss,prot,交集,可信度,子长,件数,内含子保留,retained,intron,RI,splice,sites,A3,A5,变起,first,exons,AF,外显子跳跃,skipping,SE,last,AL,互斥,mutually,exclusive,MX,上下游,碱基,倾向性,基础生物学,剪接异构体,可变剪接

0.295963