Genome Sequence of the Asian Honeybee in Pakistan Sheds Light on Its Phylogenetic Relationship with Other Honeybees.

Hongwei Tan, Muhammad Naeem, Hussain Ali, Muhammad Shakeel, Haiou Kuang, Ze Zhang, Cheng Sun
Author Information
  1. Hongwei Tan: Laboratory of Evolutionary and Functional Genomics, School of Life Sciences, Chongqing University, Chongqing 401331, China.
  2. Muhammad Naeem: Department of Zoology, Faculty of Engineering and Applied Sciences, Faisalabad Campus, Riphah International University, Faisalabad 38000, Pakistan.
  3. Hussain Ali: Entomology Section, Agriculture Research Institute D. I. Khan, Dera Ismail Khan 29120, Pakistan.
  4. Muhammad Shakeel: Department of Entomology, The University of Agriculture Peshawar, Peshawar 25120, Pakistan.
  5. Haiou Kuang: Faculty of Animal Science and Technology, Yunnan Agricultural University, Kunming 650201, China.
  6. Ze Zhang: Laboratory of Evolutionary and Functional Genomics, School of Life Sciences, Chongqing University, Chongqing 401331, China. ORCID
  7. Cheng Sun: Institute of Apicultural Research, Chinese Academy of Agricultural Sciences, Beijing 100093, China. ORCID

Abstract

In Pakistan, the Asian Honeybee, has been used for honey production and pollination services. However, its genomic makeup and phylogenetic relationship with those in other countries are still unknown. We collected samples from the main cerana-keeping region in Pakistan and performed whole genome sequencing. A total of 28 Gb of Illumina shotgun reads were generated, which were used to assemble the genome. The obtained genome assembly had a total length of 214 Mb, with a GC content of 32.77%. The assembly had a scaffold N50 of 2.85 Mb and a BUSCO completeness score of 99%, suggesting a remarkably complete genome sequence for in Pakistan. A MAKER pipeline was employed to annotate the genome sequence, and a total of 11,864 protein-coding genes were identified. Of them, 6750 genes were assigned at least one GO term, and 8813 genes were annotated with at least one protein domain. Genome-scale phylogeny analysis indicated an unexpectedly close relationship between in Pakistan and those in China, suggesting a potential human introduction of the species between the two countries. Our results will facilitate the genetic improvement and conservation of in Pakistan.

Keywords

References

  1. BMC Genomics. 2016 Mar 05;17:187 [PMID: 26944054]
  2. BMC Bioinformatics. 2014 Aug 15;15:281 [PMID: 25128196]
  3. Insects. 2013 Oct 21;4(4):558-92 [PMID: 26462524]
  4. Biotechniques. 2015 Jul 01;59(1):42-5 [PMID: 26156783]
  5. Front Genet. 2020 Mar 27;11:279 [PMID: 32292419]
  6. Genome Res. 2008 Jan;18(1):188-96 [PMID: 18025269]
  7. BMC Genomics. 2015 Jan 02;16:1 [PMID: 25553907]
  8. Sci Rep. 2018 Jan 16;8(1):822 [PMID: 29339745]
  9. Mol Biol Evol. 2015 Jan;32(1):268-74 [PMID: 25371430]
  10. Nucleic Acids Res. 1997 Sep 1;25(17):3389-402 [PMID: 9254694]
  11. Nat Methods. 2017 Jun;14(6):587-589 [PMID: 28481363]
  12. Mol Biol Evol. 2018 Mar 1;35(3):543-548 [PMID: 29220515]
  13. BMC Bioinformatics. 2004 May 14;5:59 [PMID: 15144565]
  14. Nucleic Acids Res. 2002 Jul 15;30(14):3059-66 [PMID: 12136088]
  15. Mol Biol Evol. 2021 Jan 23;38(2):486-501 [PMID: 32946576]
  16. Bioinformatics. 2014 May 1;30(9):1236-40 [PMID: 24451626]
  17. BMC Evol Biol. 2010 Jul 13;10:210 [PMID: 20626897]
  18. Nucleic Acids Res. 2005 Nov 28;33(20):6494-506 [PMID: 16314312]
  19. Nucleic Acids Res. 1999 Jan 15;27(2):573-80 [PMID: 9862982]
  20. Proc Natl Acad Sci U S A. 2020 Apr 28;117(17):9451-9457 [PMID: 32300014]
  21. Bioinformatics. 2019 Nov 1;35(21):4453-4455 [PMID: 31070718]
  22. BMC Bioinformatics. 2006 Feb 09;7:62 [PMID: 16469098]
  23. Genome Biol. 2019 Nov 14;20(1):238 [PMID: 31727128]

Grants

  1. Y2019XK13 and Y2021XK16/the Fundamental Research Funds of Chinese Academy of Agricultural Sciences

Word Cloud

Created with Highcharts 10.0.0PakistangenometotalgenesAsianusedrelationshipcountriessequencingassemblyMbsuggestingsequenceleastonephylogenyhoneybeehoneyproductionpollinationservicesHowevergenomicmakeupphylogeneticstillunknowncollectedsamplesmaincerana-keepingregionperformedwhole28GbIlluminashotgunreadsgeneratedassembleobtainedlength214GCcontent3277%scaffoldN50285BUSCOcompletenessscore99%remarkablycompleteMAKERpipelineemployedannotate11864protein-codingidentified6750assignedGOterm8813annotatedproteindomainGenome-scaleanalysisindicatedunexpectedlycloseChinapotentialhumanintroductionspeciestworesultswillfacilitategeneticimprovementconservationGenomeSequenceHoneybeeShedsLightPhylogeneticRelationshipHoneybeesApiscerana

Similar Articles

Cited By (1)