Acta Entomologica Sinica ›› 2022, Vol. 65 ›› Issue (2): 187-196.doi: 10.16380/j.kcxb.2022.02.007

• RESEARCH PAPERS • Previous Articles     Next Articles

Exploration and analysis of SNP and InDel sites in Ascosphaera apis based on PacBio sequencing data

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,*   

  1.  (1. College of Animal Sciences (College of Bee Science), Fujian Agriculture and Forestry University, Fuzhou 350002, China; 2. Apitherapy Research Institute, Fujian Agriculture and Forestry University, Fuzhou 350002, China; 3. Apicultural Research Institute of Jiangxi Province, Nanchang 330000, China)
  • Online:2022-02-20 Published:2022-03-10

Abstract: 【Aim】 In this study, single nucleotide polymorphism (SNP) and insertion-deletion (InDel) mutation sites were explored and analyzed using previously obtained PacBio single molecule realtime (SMRT) sequencing data of Ascosphaera apis mycelia, aiming to enrich information about SNP and InDel sites in A. apis and to provide a foundation for the development and application of new molecular markers. 【Methods】 SAMtools was employed to survey the full-length transcripts of the PacBio SMRT sequencing data of A. apis mycelia, followed by detection and analysis of SNP and InDel sites by alignment of the full-length transcripts to the reference genome of A. apis (assembly AAP 1.0) with ANNOVAR software. Genes containing SNP and InDel sites were aligned to GO database and KEGG database using related bioinformatic software to gain corresponding function and pathway annotations. 【Results】 A total of 6 743 SNP sites including 6 091 homozygous sites and 652 heterozygous sites were identified in A. apis mycelium and mainly distributed in exons, and the numbers of SNP sites with transition and transversion were 4 887 and 1 856, respectively. Additionally, the most abundant mutation type of the codon of SNP was synonymous single nucleotide mutation. Moreover, genes containing SNP sites could be annotated to 34 GO terms and 76 KEGG pathways. In total, 597 InDel sites including 349 homozygous sites and 248 heterozygous sites were identified and mainly distributed in genic downstream region, and non-frame shift insertion was the most abundant codon mutation type. Furthermore, genes containing InDel sites could be annotated to 39 GO terms and 87 KEGG pathways. 【Conclusion】 We identified a large quantity of SNP and InDel sites of A. apis, pinpointed their mutation type, distribution in genomic functional element and codon mutation type, and uncovered the potential relationship between SNP and InDel sites and the key biological processes.

Key words:  Ascosphaera apis, third-generation sequencing technology, single molecule real-time sequencing, single nucleotide polymorphism, insertion and deletion