昆虫学报 ›› 2022, Vol. 65 ›› Issue (2): 187-196.doi: 10.16380/j.kcxb.2022.02.007
蔡宗兵1,#, 张文德1,#, 隆琦1, 余岢骏1, 孙明会1, 吴鹰1, 许雅静1, 刘佳美1, 郭意龙1, 徐细建3, 陈大福1,2,*, 郭睿1,2,*
CAI Zong-Bing1,#, ZHANG Wen-De1,#, LONG Qi1, YU Ke-Jun1, SUN Ming-Hui1, WU Ying1, XU Ya-Jing1, LIU Jia-Mei1, GUO Yi-Long1, XU Xi-Jian3, CHEN Da-Fu1,2,*, GUO Rui1,2,*
摘要: 【目的】本研究拟利用已获得的蜜蜂球囊菌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位点与关键生物学过程的潜在关联。