基于PacBio测序数据的蜜蜂球囊菌SNP与InDel位点发掘及分析

[目的]本研究拟利用已获得的蜜蜂球囊菌Ascosphaera apis菌丝的PacBio单分子实时(single molecule real-time,SMRT)测序数据对蜜蜂球囊菌的单核苷酸多态性(single nucleotide polymorphism,SNP)和插入缺失(insertion-deletion,InDel)突变位点进行发掘和分析,旨在丰富蜜蜂球囊菌的SNP和InDel位点信息,并为新型分子标记的开发和应用提供基础.[方法]采用SAMtools对蜜蜂球囊菌菌丝PacBio SMRT测序数据进行全长转录本识别,再利用ANNOVAR软件将识别的全长转录本与蜜蜂球囊菌参考基因组(assembly AAP 1.0)进行比对以检测和分析SNP位点和InDel位点.通过相关生物信息学软件将上述SNP位点和InDel位点基因分别比对GO和KEGG数据库进行功能和通路注释.[结果]在蜜蜂球囊菌菌丝中共鉴定到6 743个SNP位点,主要分布于外显子,包括6 091个纯合位点和652个杂合位点;其中发生转换和颠换的SNP位点分别有4 887和1 856个;SNP位点密码子突变类型中数量最多的是同义单核苷酸突变;SNP位点基因可注释到34个GO条目和76个KEGG通路.共鉴定到597个InDel位点,包括349个纯合位点和248个杂合位点,这些InDel位点在基因下游区分布的数量最多;InDel位点密码子突变类型中非移码插入最为丰富;InDel位点基因可在39个GO条目和87条KEGG通路.[结论]鉴定到蜜蜂球囊菌的大量SNP和InDel位点,明确了 SNP和InDel位点的突变类型、基因组功能元件分布和密码子突变类型,并揭示了 SNP和InDel位点与关键生物学过程的潜在关联.