Next Article in Journal
Investigating Virus–Host Interactions in Cultured Primary Honey Bee Cells
Next Article in Special Issue
Workflows for Rapid Functional Annotation of Diverse Arthropod Genomes
Previous Article in Journal
Cytochrome P450s Are Essential for Insecticide Tolerance in the Endoparasitoid Wasp Meteorus pulchricornis (Hymenoptera: Braconidae)
Previous Article in Special Issue
The USDA-ARS Ag100Pest Initiative: High-Quality Genome Assemblies for Agricultural Pest Arthropod Research
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Brief Report

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

1
Laboratory of Evolutionary and Functional Genomics, School of Life Sciences, Chongqing University, Chongqing 401331, China
2
Chongqing General Station of Animal Husbandry Technology Extension, Chongqing 401121, China
3
Department of Zoology, Faculty of Engineering and Applied Sciences, Faisalabad Campus, Riphah International University, Faisalabad 38000, Pakistan
4
Entomology Section, Agriculture Research Institute D. I. Khan, Dera Ismail Khan 29120, Pakistan
5
Department of Entomology, The University of Agriculture Peshawar, Peshawar 25120, Pakistan
6
Faculty of Animal Science and Technology, Yunnan Agricultural University, Kunming 650201, China
7
Institute of Apicultural Research, Chinese Academy of Agricultural Sciences, Beijing 100093, China
*
Authors to whom correspondence should be addressed.
These authors contribute equally.
Submission received: 29 May 2021 / Revised: 6 July 2021 / Accepted: 15 July 2021 / Published: 16 July 2021
(This article belongs to the Special Issue Insect Genomics)

Abstract

:

Simple Summary

The Asian honeybee, Apis cerana, is used for honey production and pollination services in Pakistan. However, its genome sequence is still unknown. We collected A. cerana samples from its main rearing region in Pakistan and performed whole genome sequencing. We obtained a remarkably complete genome sequence for A. cerana in Pakistan, from which we identified a total of 11,864 protein-coding genes. Phylogeny analysis indicated an unexpectedly close relationship between A. cerana 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 A. cerana in Pakistan.

Abstract

In Pakistan, Apis cerana, 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 A. cerana 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 A. cerana 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 A. cerana 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 A. cerana in Pakistan.

1. Introduction

Among all honeybee species, the three species, Apis cerana, A. florea, and A. dorsata, are endemic to Pakistan. The traditional beekeepers have been using A. cerana for honey production for the last two centuries because of its ability to cope with environmental challenges [1,2]. Currently, the average honey production of A. cerana has ranged from 7 to 9 kg/colony, and on an annual basis, 60 tons of honey are produced from A. cerana in Pakistan [2,3]. A. cerana could also be used for pollination services to improve crop production [4]. The first time, PARC (Pakistan Agricultural Research Council) made an effort to use A. cerana as pollinators, which has also promoted increasing awareness among the farmers to use honeybee colonies as pollinators in Pakistan [5].
Currently, introduced honeybee species Apis mellifera are used for honey production and pollination services in Pakistan. However, such an introduction presents a severe threat to local biodiversity [6]. For example, in the 1980s, the tracheal mite attacking A. cerana might be due to the mite transfer from A. mellifera [2]. Hence, priority should be given to local species for honey production and pollination services, to avoid such incidents. A. cerana beekeeping is again being encouraged, especially in the mountainous regions where it is better adapted than A. mellifera and has no problems with Varroa destructor [1]. However, the genomic makeup of the commonly used A. cerana in Pakistan is still unknown, which will impede its improvement and conservation. In addition, the natural distribution of A. cerana occurs across southern and southeastern Asia, and extends to the Afghanistan region in the west [7]. There are different climatic conditions across the distribution range of this species, which will promote the genetic differentiation of A. cerana. However, little is known about the phylogenetic relationship between A. cerana in Pakistan and those in other countries.
In this study, we collected A. cerana samples from the main cerana-keeping region in Pakistan. We then generated and characterized its genome sequence, and based on this we inferred the phylogenetic relationship between A. cerana in Pakistan and those in other countries.

2. Materials and Methods

Drones of A. cerana were collected from one colony reared at an apiary located in Khyber Pakhtunkhwa, a main cerana-keeping region in Pakistan (34.0851029° E & 71.6149948° N), in 2018. Genomic DNA was purified from each drone with the use of a Qiagen Gentra Puregene Tissue Kit. To avoid contamination, the abdomen of each drone was removed before DNA extraction.
We took a genome sequencing and assembly method as described before [8]. In brief, one single drone of A. cerana was used to extract genomic DNA, based on which one fragment library (insert size: 450 bp) was generated by using NEBNext® Ultra™ DNA Library Prep Kit for Illumina® (NEB, Ipswich, MA, USA). The obtained fragment library was sequenced on the Illumina HiSeq 2500 platform (read length: 250 bp), from which overlapping paired-end shotgun reads (2 × 250 bp) were produced. Genomic DNA purified from five drones of the same A. cerana colony was used to produce large-insert jumping libraries with four different sizes (4 kb, 6 kb, 8 kb, and 10 kb, respectively), following a protocol as described before [9]. The obtained jumping libraries were sequenced on the Illumina HiSeq X Ten platform, with a read length of 150 bp, from which paired-end shotgun reads (2 × 150 bp) were generated. To produce genome assembly, at first, the 250 bp overlapping paired-end reads from the fragment library were assembled using the software DISCOVAR de novo v52488, an assembler performing well at assembling insect genomes [10], to produce continuous sequences (contigs). Then, shotgun reads from the four jumping libraries were used for scaffolding the contigs using the software BESST v2.2.6 [11]. Software BUSCO v3 [12] was employed to evaluate the completeness of the genome assembly by using 4415 Benchmarking Universal Single-Copy Orthologs (BUSCOs) (dataset: hymenoptera_odb9).
Protein-coding genes of the A. cerana genome were annotated by the MAKER computational pipeline, which was based on ab initio gene predictions, transcript evidence, and homologous protein evidence [13]. Ab initio gene prediction software includes SNAP [14], GENEMARK [15] and AUGUSTUS [16]. Transcript evidence came from predicted transcripts of A. cerana (GenBank: GCF_001442555.1_ACSNU-2.0_rna.fna) and A. mellifera (GenBank: GCF_003254395.2_Amel_HAv3.1_rna.fna). Protein evidence came from annotated Apoidea proteins that are publicly available in GenBank (last accessed on 29 January 2018).
To provide functional clues for the predicted protein-coding genes, the protein sequences that they encode were used as queries to do BLASTp (in BLAST suite v2.28 [17]) against the Swiss-Prot database (last accessed on 28 January 2018). In addition, program InterproScan-5 [18] was used to identify protein domains and GO terms that are associated with the predicted protein-coding genes. The quality of gene annotation was evaluated by comparing predicted protein sequences of the A. cerana genome to 4415 BUSCOs, using the software BUSCO v3 (lineage dataset: hymenoptera_odb9) [12].
Transposable elements (TEs) were de novo identified by RepeatModeller2 [19]. Using the repeat library obtained from RepeatModeller2, A. cerana genome assembly was analyzed with RepeatMasker (http://www.repeatmasker.org (accessed on 28 January 2018) to generate a summary of its TE landscape. Tandem Repeat Finder [20] was employed to identify tandem repeats in the A. cerana genome.
The predicted protein sequences for A. mellifera (GenBank assembly: GCA_003254395.2), A. cerana from Korea (GenBank assembly: GCA_001442555.1), A. cerana from China (GenBank assembly: GCA_002290385.1), A. florea (GenBank assembly: GCA_000184785.2), A. dorsata (GenBank assembly: GCA_000469605.1), and Bombus terrestris (GenBank assembly: GCA_000214255.1) were downloaded from GenBank. Proteins for A. cerana in Pakistan are from this study. Proteins for A. laboriosa were from the National Genomics Data Center, Beijing Institute of Genomics, Chinese Academy of Sciences (https://bigd.big.ac.cn/ (accessed on 21 October 2020); assembly accession: GWHAOTM00000000). Only the longest isoform of each protein-coding gene was kept for downstream analysis. Protein sequences of the eight sequenced bee genomes were used to infer orthologous groups (gene families) by using Orthofinder v2.4.0 [21].
To construct the molecular phylogeny for these seven honey bees that have a sequenced genome (A. mellifera, A. cerana from Korea, A. cerana from China, A. florea, A. dorsata, A. laboriosa and A. cerana from Pakistan), with one bumblebee (B. terrestris) serving as the outgroup, universal single-copy orthologs inferred by Orthofinder [21] were extracted. Protein sequences from those universal single-copy orthologs were multiply aligned with program MAFFT [22], respectively, followed by alignment trimming by software BMGE [23]. Alignments shorter than 50 amino acids or with more than 20% gaps were removed. The remaining alignments were concatenated for each of the eight bee species, resulting in eight super-sequences. Software ModelFinder [24] was used to select the best-fitting amino acid substitution model. IQ-TREE version 1.6.1 [25] and RAxML-NG version 1.0.2 [26] were employed to construct maximum likelihood trees for these eight bees, with 1000 bootstrap replicates, respectively.

3. Results

3.1. Genome Sequencing, Assembly and Annotation

Genomic DNA extracted from six males of the same A. cerana colony were used to construct one fragment library and four jumping libraries for whole genome sequencing (see Methods). We obtained a total of 22.67 Gb of overlapping paired-end shotgun reads (2 × 250 bp) from the fragment library (Table 1), which were used to produce contigs. A further 5.95 Gb of shotgun reads were obtained from the four jump libraries (Table 1), which were used to scaffold the contigs.
The obtained genome assembly had a total length of 214.44 Mb, with a GC content of 32.77%. The assembly had a scaffold N50 of 2.85 Mb and a contig N50 of over 311 Kb (Table 1). To evaluate the completeness of the resulting genome assembly, we compared genes present in the assembly to 4415 hymenopteran BUSCOs. We found that 99.5% of BUSCOs (99% of them are complete) were present in the genome assembly (Table 1), suggesting that we have obtained a remarkably complete genome sequence for reared A. cerana in Pakistan.
A total of 11,864 protein-coding genes were predicted using the MAKER pipeline (Table 1). Of them, 6750 genes were assigned at least one Gene Ontology (GO) term, and 8813 genes were annotated with at least one protein domain (Table 1). BUSCO analysis suggested a high-quality of genome annotation (Table 1). The assembled A. cerana genome contains 6.08% of repetitive sequences, and 0.7% of the genome are recognized as transposable elements (TEs). Mariner DNA transposons represent the most abundant TEs, consistent with previous reports for A. cerana in other countries [27,28,29].

3.2. Phylogeny Analysis of A. cerana in Pakistan

To understand the phylogenetic position of reared A. cerana in Pakistan, we constructed the genome-scale molecular phylogeny for seven honeybee genomes, including three Asian bees from different geographical regions (A. mellifera, A. florea, A. dorsata, A. laboriosa, A. cerana from Korea, A. cerana from China and A. cerana from Pakistan), with one bumblebee species (B. terrestris) serving as the outgroup. Protein sequences of these eight bees were downloaded and clustered, with single copy orthologs being delineated (see Methods). A total of 6101 single copy genes were identified across the eight bees and were multiply aligned. After trimming the alignments and getting rid of low-quality alignments, a total of 5234 single-copy genes remained, which were concatenated for each species, resulting in eight super-sequences. The super-alignment contained 2,980,871 amino-acid sites, with 44,493 distinct site patterns.
IQ-TREE [25] and RAxML-NG [26] were used to construct maximum likelihood trees for the eight bees, with the best-fitting amino acid substitution model JTT+F+R3. Both softwares yielded exactly the same topology, with high supporting values (Figure 1). From the phylogeny, we could see that A. cerana in Pakistan consistently clustered with A. cerana from China with a short branch, which formed a sister group to A. cerana from Korea.

4. Discussion

At present, there are nine species that have been commonly recognized in the genus Apis, which including A. andreniformis, A. cerana, A. dorsata, A. florea, A. koschevnikovi, A. laboriosa, A. mellifera, A. nigrocincta and A. nuluensis [30]. Phylogenetic analyses using mitochondrial genes strongly supported the basic topology that groups the nine honeybee species into three major clusters: giant bees (A. dorsata and A. laboriosa), dwarf bees (A. andreniformis and A. florea) and cavity-nesting bees (A. mellifera, A. cerana, A. koschevnikovi, A. nuluensis and A. nigrocincta) [30,31] A. mellifera and A. cerana are both widely used for pollination services and honey production. Based on a recent analysis, Apis species originated from tropical regions, and out of the nine Apis species, only A. mellifera and A. cerana have expanded northward into the temperate zone, which may explain their large-scale application in agriculture [31].
In this study, we have generated a high-quality reference genome for A. cerana in Pakistan (Table 1), which will form a foundation for future research, including resequencing and population genomic studies for functional gene positioning and cloning, and thus will facilitate the genetic improvement and conservation of A. cerana in Pakistan. Besides the A. cerana genome in Pakistan that was generated in this study, there are two other A. cerana genomes with publicly available genome assemblies and annotations: one is for A. cerana samples collected in Korea [27], and the other is for A. cerana samples collected in China [28]. Phylogenetic analysis based on whole genomes will yield more accurate results than analysis based solely on mitochondrial genes, so these A. cerana genomes will be valuable for understanding the origin and evolution of bees.
Unexpectedly, we found a very close relationship between A. cerana in Pakistan and A. cerana in China (Figure 1). Pakistan represents the west end of the A. cerana natural distribution range [32], so it is not impossible that A. cerana in the two counties are similar. However, the Himalayan mountains lie between Pakistan and China, which could serve as a natural barrier to prevent the mixing of A. cerana between the two countries. In contrast, no such natural barriers exist between China and Korea. Therefore, it is quite un-expected that the phylogenetic relationship between the A. cerana of China and Pakistan is closer than that of China and Korea (Figure 1). Although phylogeny based solely on bee species with genome sequences could be incomplete, it might be also possible that A. cerana in Pakistan was introduced via trade or exchange from China, especially considering the fact that during 1980–1983, almost all of the A. cerana colonies kept in modern and traditional hives were killed by an epidemic of tracheal mites in Pakistan [2].

Author Contributions

Conceptualization, Z.Z. and C.S.; methodology, H.T. and C.S.; resources, M.N., H.A., M.S., and H.K.; data curation, C.S.; writing—original draft preparation, H.T., C.S. and M.N.; writing—review and editing, Z.Z.; supervision, Z.Z. and C.S.; funding acquisition, C.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Fundamental Research Funds of the Chinese Academy of Agricultural Sciences (grant numbers: Y2019XK13 and Y2021XK16).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The raw Illumina shotgun reads and assembled genome sequence for Apis cerana in Pakistan have been deposited in National Center for Biotechnology Information (NCBI) under the BioProject number PRJNA669617.

Acknowledgments

The authors would like to thank Zachary Huang for his critical reading and improvement on the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Khan, K.A.; Ansari, M.J.; Al-Ghamdi, A.; Sharma, D.; Ali, H. Biodiversity and relative abundance of different honeybee species (Hymenoptera: Apidae) in Murree-Punjab, Pakistan. J. Entomol. Zool. Stud. 2014, 2, 324–327. [Google Scholar]
  2. Khan, K. Beekeeping in Pakistan (History, Potential, and Current Status). Preprints 2020, 2020070503. [Google Scholar] [CrossRef]
  3. Ahmad, R.; Camphor, E.; Ahmed, M. Factors affecting honey yield of the oriental bee, Apis cerana in Pakistan. Pak. J. Agric. Res. 1983, 4, 190–197. [Google Scholar]
  4. Koetz, A.H. Ecology, Behaviour and Control of Apis cerana with a Focus on Relevance to the Australian Incursion. Insects 2013, 4, 558–592. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Jasra, A.W.; Sabira, A.; Kasi, M.A. Apple pollination problems in Balochistan, Pakistan. Int. J. Agric. Biol. 2001, 3, 210–213. [Google Scholar]
  6. Chantawannakul, P.; De Guzman, L.I.; Li, J.; Williams, G.R. Parasites, pathogens, and pests of honeybees in Asia. Apidologie 2016, 47, 301–324. [Google Scholar] [CrossRef]
  7. Ilyasov, R.; Youn, H.G.; Lee, M.-L.; Kim, K.W.; Proshchalykin, M.; Lelej, A.; Takahashi, J.-I.; Kwon, H.W. Phylogenetic Relationships of Russian Far-East Apis cerana with Other North Asian Populations. J. Apic. Sci. 2019, 63, 289–314. [Google Scholar] [CrossRef] [Green Version]
  8. Sun, C.; Huang, J.; Wang, Y.; Zhao, X.; Su, L.; Thomas, G.W.C.; Zhao, M.; Zhang, X.; Jungreis, I.; Kellis, M.; et al. Genus-Wide Characterization of Bumblebee Genomes Provides Insights into Their Evolution and Variation in Ecological and Behavioral Traits. Mol. Biol. Evol. 2021, 38, 486–501. [Google Scholar] [CrossRef]
  9. Heavens, D.; Accinelli, G.G.; Clavijo, B.; Clark, M.D. A method to simultaneously construct up to 12 differently sized Illumina Nextera long mate pair libraries with reduced DNA input, time, and cost. Biotechniques 2015, 59, 42–45. [Google Scholar] [CrossRef] [Green Version]
  10. Love, R.R.; Weisenfeld, N.I.; Jaffe, D.B.; Besansky, N.J.; Neafsey, D.E. Evaluation of DISCOVAR de novo using a mosquito sample for cost-effective short-read genome assembly. BMC Genom. 2016, 17, 187. [Google Scholar] [CrossRef] [Green Version]
  11. Sahlin, K.; Vezzi, F.; Nystedt, B.; Lundeberg, J.; Arvestad, L. BESST-Efficient scaffolding of large fragmented assemblies. BMC Bioinform. 2014, 15, 1–11. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Waterhouse, R.M.; Seppey, M.; Simão, F.A.; Manni, M.; Ioannidis, P.; Klioutchnikov, G.; Kriventseva, E.V.; Zdobnov, E.M. BUSCO Applications from Quality Assessments to Gene Prediction and Phylogenomics. Mol. Biol. Evol. 2018, 35, 543–548. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Cantarel, B.L.; Korf, I.; Robb, S.M.; Parra, G.; Ross, E.; Moore, B.; Holt, C.; Alvarado, A.S.; Yandell, M. MAKER: An easy-to-use annotation pipeline designed for emerging model organism genomes. Genome Res. 2008, 18, 188–196. [Google Scholar] [CrossRef] [Green Version]
  14. Korf, I. Gene finding in novel genomes. BMC Bioinform. 2004, 5, 59. [Google Scholar] [CrossRef] [Green Version]
  15. Lomsadze, A.; Ter-Hovhannisyan, V.; Chernoff, Y.; Borodovsky, M. Gene identification in novel eukaryotic genomes by self-training algorithm. Nucleic Acids Res. 2005, 33, 6494–6506. [Google Scholar] [CrossRef] [PubMed]
  16. Stanke, M.; Schöffmann, O.; Morgenstern, B.; Waack, S. Gene prediction in eukaryotes with a generalized hidden Markov model that uses hints from external sources. BMC Bioinform. 2006, 7, 62. [Google Scholar] [CrossRef] [Green Version]
  17. Altschul, S.F.; Madden, T.L.; Schäffer, A.A.; Zhang, J.; Zhang, Z.; Miller, W.; Lipman, D.J. Gapped BLAST and PSI-BLAST: A new generation of protein database search programs. Nucleic Acids Res. 1997, 25, 3389–3402. [Google Scholar] [CrossRef] [Green Version]
  18. Jones, P.; Binns, D.; Chang, H.-Y.; Fraser, M.; Li, W.; McAnulla, C.; McWilliam, H.; Maslen, J.; Mitchell, A.; Nuka, G.; et al. InterProScan 5: Genome-scale protein function classification. Bioinformatics 2014, 30, 1236–1240. [Google Scholar] [CrossRef] [Green Version]
  19. Flynn, J.M.; Hubley, R.; Goubert, C.; Rosen, J.; Clark, A.G.; Feschotte, C.; Smit, A.F. RepeatModeler2 for automated genomic discovery of transposable element families. Proc. Natl. Acad. Sci. USA 2020, 117, 9451–9457. [Google Scholar] [CrossRef]
  20. Benson, G. Tandem repeats finder: A program to analyze DNA sequences. Nucleic Acids Res. 1999, 27, 573–580. [Google Scholar] [CrossRef] [Green Version]
  21. Emms, D.; Kelly, S. OrthoFinder: Phylogenetic orthology inference for comparative genomics. Genome Biol. 2019, 20, 1–14. [Google Scholar] [CrossRef] [Green Version]
  22. Katoh, K.; Misawa, K.; Kuma, K.; Miyata, T. MAFFT: A novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res. 2002, 30, 3059–3066. [Google Scholar] [CrossRef] [Green Version]
  23. Criscuolo, A.; Gribaldo, S. BMGE (Block Mapping and Gathering with Entropy): A new software for selection of phylogenetic informative regions from multiple sequence alignments. BMC Evol. Biol. 2010, 10, 210. [Google Scholar] [CrossRef] [Green Version]
  24. Kalyaanamoorthy, S.; Minh, B.Q.; Wong, T.K.F.; Von Haeseler, A.; Jermiin, L.S. ModelFinder: Fast model selection for accurate phylogenetic estimates. Nat. Methods 2017, 14, 587–589. [Google Scholar] [CrossRef] [Green Version]
  25. Lam-Tung, N.; Schmidt, H.A.; von Haeseler, A.; Minh, B.Q. IQ-TREE: A Fast and Effective Stochastic Algorithm for Estimating Maximum-Likelihood Phylogenies. Mol. Biol. Evol. 2015, 32, 268–274. [Google Scholar]
  26. Kozlov, A.M.; Darriba, D.; Flouri, T.; Morel, B.; Stamatakis, A. RAxML-NG: A fast, scalable and user-friendly tool for maximum likelihood phylogenetic inference. Bioinformatics 2019, 35, 4453–4455. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Park, D.; Jung, J.W.; Choi, B.-S.; Jayakodi, M.; Lee, J.; Lim, J.; Yu, Y.; Choi, Y.-S.; Lee, M.-L.; Park, Y.; et al. Uncovering the novel characteristics of Asian honey bee, Apis cerana, by whole genome sequencing. BMC Genom. 2015, 16, 1. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Diao, Q.; Sun, L.; Zheng, H.; Zeng, Z.; Wang, S.; Xu, S.; Zheng, H.; Chen, Y.; Shi, Y.; Wang, Y.; et al. Genomic and transcriptomic analysis of the Asian honeybee Apis cerana provides novel insights into honeybee biology. Sci. Rep. 2018, 8, 1–14. [Google Scholar] [CrossRef]
  29. Wang, Z.-L.; Zhu, Y.; Yan, Q.; Yan, W.; Zheng, H.; Zeng, Z. A Chromosome-Scale Assembly of the Asian Honeybee Apis cerana Genome. Front. Genet. 2020, 11, 279. [Google Scholar] [CrossRef] [PubMed]
  30. Koeniger, N.; Koeniger, G.; Smith, D. Phylogeny of the Genus Apis; Springer: Berlin/Heidelberg, Germany, 2011; pp. 23–50. [Google Scholar]
  31. Ji, Y. The geographical origin, refugia, and diversification of honey bees (Apis spp.) based on biogeography and niche modeling. Apidologie 2021, 52, 367–377. [Google Scholar] [CrossRef]
  32. Radloff, S.E.; Hepburn, C.; Hepburn, H.R.; Fuchs, S.; Hadisoesilo, S.; Tan, K.; Engel, M.; Kuznetsov, V. Population structure and classification ofApis cerana. Apidologie 2010, 41, 589–601. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Maximum likelihood trees constructed by IQ-TREE (A) and RAxML-NG (B). Numbers on the nodes indicate bootstrap values.
Figure 1. Maximum likelihood trees constructed by IQ-TREE (A) and RAxML-NG (B). Numbers on the nodes indicate bootstrap values.
Insects 12 00652 g001
Table 1. The results of genome sequencing, assembly and annotation for A. cerana in Pakistan.
Table 1. The results of genome sequencing, assembly and annotation for A. cerana in Pakistan.
Genome Sequencing
Read number (Million)Read length (bp)Total read length (Gb)
Fragment library90.66 25022.67
Jump libraries39.651505.95
Genome Assembly
Genome assembly sizeScaffold N50 (Mb)Contig N50 (Kb)BUSCO
214.442.85311.1399.50%
Genome Annotation
Protein-coding gene numberGenes with a GO termGenes with a protein domainBUSCO
11,8646750881397.90%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tan, H.; Naeem, M.; Ali, H.; Shakeel, M.; Kuang, H.; Zhang, Z.; Sun, C. Genome Sequence of the Asian Honeybee in Pakistan Sheds Light on Its Phylogenetic Relationship with Other Honeybees. Insects 2021, 12, 652. https://0-doi-org.brum.beds.ac.uk/10.3390/insects12070652

AMA Style

Tan H, Naeem M, Ali H, Shakeel M, Kuang H, Zhang Z, Sun C. Genome Sequence of the Asian Honeybee in Pakistan Sheds Light on Its Phylogenetic Relationship with Other Honeybees. Insects. 2021; 12(7):652. https://0-doi-org.brum.beds.ac.uk/10.3390/insects12070652

Chicago/Turabian Style

Tan, Hongwei, Muhammad Naeem, Hussain Ali, Muhammad Shakeel, Haiou Kuang, Ze Zhang, and Cheng Sun. 2021. "Genome Sequence of the Asian Honeybee in Pakistan Sheds Light on Its Phylogenetic Relationship with Other Honeybees" Insects 12, no. 7: 652. https://0-doi-org.brum.beds.ac.uk/10.3390/insects12070652

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop