结合三代测序与二代测序技术揭示蜜蜂球囊菌孢子转录组的复杂性

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

来源：微生物学报, 2022, 62(08): 2981-2994.

DOI：10.13343/j.cnki.wsxb.20210711

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

【目的】蜜蜂球囊菌(Ascosphaera apis)是一种专性侵染蜜蜂幼虫的致死性真菌病原。本研究旨在利用PacBio单分子实时(singlemoleculereal-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)剪接位点上游50bp的序列特征并鉴定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'剪切(alternative3'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%)可变末端外显子(alternativelastexons,AL)和1次(0.04%)互斥外显子(mutuallyexclusiveexons,MX)。通过RT-PCR证实了不同类型AS事件(RI、A5、SE和A3)的可靠性。共鉴定出5 552个基因含APA位点，其中含有>5个APA位点的基因数量最多，达到2 197个，其次为含有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，研究结果丰富了蜜蜂球囊菌的基础生物学信息，深入揭示了蜜蜂球囊菌孢子中转录组的复杂性，并为进一步探究不同剪接异构体在病原孢子和病原侵染中的分子功能打下了基础。

• 单位
  动物科学学院; 江西省养蜂研究所; 福建农林大学