中华蜜蜂酚氧化酶和丝氨酸蛋白酶相关基因 及其全长转录本发掘与分析

为深入开展中华蜜蜂相关研究提供参考信息和基础，利用已获得的纳米孔(Nanopore)测序数据对中华蜜蜂Apis cerana cerana的酚氧化酶(phenoloxidase, PO)和丝氨酸蛋白酶(serine protease, SP)相关基因和全长转录本进行发掘和分析。通过Blast工具将中华蜜蜂所有全长转录本比对KEGG和Nr数据库以筛选出PO和SP相关基因和全长转录本。采用gffcompare软件将PO和SP相关全长转录本与东方蜜蜂参考基因组(ACSNU-2. 0)已注释基因进行比较以优化结构。使用Astalavista软件鉴定PO和SP相关基因的可变剪接(alternative splicing, AS)事件，再利用IGV浏览器进行结构可视化。通过PCR验证AS事件的真实性。利用TAPIS pipeline预测和分析可变多聚腺苷酸化(alternative polyadenylation, APA)位点，并通过TBtool软件鉴定APA位点上游的基序(motif)。鉴定到中华蜜蜂PO和SP相关的42个基因与146条转录本。优化了16个参考基因组已注释PO和SP相关基因的结构，其中5′端延长和3′端延长的基因均有7个，5′端和3′端同时延长的基因有2个。鉴定到PO和SP相关的10个基因的389次AS事件，其中最丰富的AS类型是5′端可变剪接。PCR结果证实了2次AS事件的真实性。共鉴定到34个PO和SP相关基因含有1个及以上的APA位点，其中含有3个APA位点的基因数量最多；在APA位点上游鉴定到多个motif，一致性序列为：GGHKSYWSHHTRATWTCNBHDMRRYWYRTNYTVACNGCKGCDCAYTGYR。鉴定到中华蜜蜂PO和SP相关的42个基因和146条全长转录本，优化了东方蜜蜂参考基因组已注释的PO和SP相关基因结构，并发掘出PO和SP相关基因的389次AS事件和237个APA位点。