Wen Feng Bai1,2†
Yuanzhen Liu
Jay D. Evans
Qiang Huang- 1Honeybee Research Institute, Jiangxi Agricultural University, Nanchang, China
- 2Jiangxi Province Key Laboratory of Honeybee Biology and Beekeeping, Jiangxi Agricultural University, Nanchang, China
- 3Periodicals Agency, Jiangxi Agricultural University, Nanchang, China
- 4Department of Medical Biochemistry and Microbiology, Science for Life Laboratory, Uppsala University, Uppsala, Sweden
- 5Environment and Plant Protection Institute, Chinese Academy of Tropical Agricultural Sciences, Haikou, China
- 6USDA-ARS Bee Research Laboratory, Beltsville, MD, United States
The small hive beetle (SHB), a social parasite of beehives, is native to sub-Saharan Africa and has spread to America, Europe, and Australia. Recently, these beetles invaded China, causing widespread colony collapses in the honeybee, Apis cerana. In this study, single nucleotide polymorphisms (SNPs) were identified in the beetle genome from its native range (Africa), a region that was invaded by SHBs nearly 30 years ago (America), and more recent invasions (Asia). The beetles in the United States formed the earliest branch and show signs of two decades of gene flow and local adaptation to differentiate this population from the native ones. The beetles in China were deep branched and showed the highest fixation index when compared to the US populations. The number of SNPs in overexpressed genes was significantly higher than the transcriptome. Gene-expression profiles presented here distinguish the characters between adult and larvae SHBs.
Introduction
The small hive beetle (Aethina tumida Murray, hereafter SHB) is a parasite of bee colonies. This beetle belongs to the family Nitidulidae; members of this family mainly feed on decaying vegetables, ripened fruit, and sap (Mckenna et al., 2015; Neumann et al., 2016). SHBs can survive on a wide range of food but thrive in the beehive that provides shelter and protein-rich food (Cuthbertson et al., 2013; Neumann et al., 2016). The adult beetles are attracted by beehive volatiles and invade the hive around dusk (Torto et al., 2005; Torto et al., 2007; Graham et al., 2011). Once inside the beehive, adult beetles employ a “sit and wait” strategy until an opportunity for reproduction arises (Neumann et al., 2015). The adult beetles lay eggs in cracks, which hatch in approximately 72 h (Neumann et al., 2016). Larvae are the most damaging stage for the beehive as they tunnel through the combs and ferment the honey, which attracts other SHBs (Benda et al., 2008). The larvae can feed on pollen, honey, and bee brood, which is also a vector of bee viruses (Eyer et al., 2009; Huang et al., 2019; Huwiler et al., 2020). At the later stage, the larvae crawl out of the beehive and pupate in soil, which may take up to months depending on the temperature and humidity (Neumann et al., 2016; Cornelissen et al., 2020). Once adult SHBs emerge from the soil, they search for and fly to beehives individually or in swarms (Neumann and Elzen, 2004).
The SHB is endemic in sub-Saharan Africa, and the damage to local honeybees is minor, as the native bees can better guard the entrance and the comb (Neumann et al., 2016). Outside the native range, SHBs were first reported in the United States in 1996 and caused considerable damage to apiculture (Hood, 2000; van Engelsdorp et al., 2007). SHBs were then reported in Canada and Australia in 2002 (Gillespie et al., 2003; Clay, 2006), Mexico in 2007 (Del Valle Molina, 2007), Italy and the Philippines in 2014 (Palmeri et al., 2014; Cervancia et al., 2016), and South Korea in 2019 (Mohamadzade Namin et al., 2019). In China, SHBs were first found to infest Eastern honeybee colonies in 2017 (Zhao et al., 2020). Using a fragment of mitochondrial DNA, the SHBs in China and the Philippines formed a cluster that was distant from beetles collected in other areas, leaving an unresolved invasion source (Liu et al., 2021a). In this study, we used single nucleotide polymorphisms (SNPs) at the genome level to reconstruct the invasion routes of SHBs in China based on a phylogenetic analysis. In addition, the characteristically expressed genes were distinguished between larvae and adults. The number of SNPs in the highly expressed genes was further quantified and compared with the number of SNPs in the remaining genes to infer gene selection.
Materials and Methods
Sample Collection
The small hive beetle was first observed to parasitize honeybee colonies in south China in 2017. In July 2019, 23 SHBs, including 8 larvae and 15 adults, were randomly collected from five collapsed honeybee (Apis cerana) colonies in the coastal city of Haikou. The SHBs were stored in a −80°C freezer.
RNA Extraction and Gene Expression Analysis
Out of 23 SHBs, 15 beetles were used for RNA extraction, including 8 larvae and 7 adults. Total RNA was extracted from individual beetle using TRIzol. The RNAseq libraries were prepared for each beetle individually and sequenced on Illumina HiSeq 2000. In total, 15 libraries were sequenced. The quality of RNA-seq reads was controlled using Fastp with default parameters (Chen et al., 2018). The reads were aligned to the SBH genome (GCA_001937115.1) using Hisat2 with default parameters (Kim et al., 2015; Evans et al., 2018). The output files were compressed, sorted, and indexed using Samtools (Li et al., 2009). The variance of replicates was used to calculate significantly regulated genes with the edgeR package and adjusted for multiple comparisons with FDR (Robinson et al., 2010; R Core Team, 2013). The protein sequences of significantly regulated genes were used to query KEGG databases to retrieve the putative functions and involved pathways (Kanehisa and Goto, 2000).
DNA Extraction, Variant Calling, and Phylogenetic Analysis
DNA was extracted from 8 adult beetles using DNAzol individually. DNA sequencing libraries were prepared for each beetle and sequenced on Illumina HiSeq 2000. In total, 8 libraries were sequenced. In addition, DNA sequencing reads of 48 SHBs were obtained from NCBI for phylogenetic analysis, including samples from South Africa (N = 12), United States (N = 9), Tanzania (N = 6), Burkina Faso (N = 12), and Liberia (N = 9). Low-quality reads were filtered using Fastp with default parameters (Chen et al., 2018). DNA reads were then mapped to the SBH genome (GCA_001937115.1) using BWA with default parameters (Li and Durbin, 2009). The variants were called using the Picard-GATK-SNPEFF pipeline (Van der Auwera et al., 2013). Then the SNPs of 57 samples were integrated to generate a single gVCF file by CombineGVCFs. The high-quality variants were extracted using GenotypeGVCF and SelectVariants functions in GATK (McKenna et al., 2010; Christmas et al., 2021). The output VCF files were converted to PHYLIP files for genetic phylogenetic analysis using VCF2PHYLIP (Ortiz, 2019). A maximum-likelihood tree was constructed using IQTree with 1000 bootstrap replicates (Nguyen et al., 2015). As the host, the honeybee (Apis mellifera) was used as an outgroup to root the tree. The fixation index FST was pairwise calculated among the beetles collected in Africa, the United States, and China using PoPoolation2 (Kofler et al., 2011). The distribution of significantly regulated genes in adults and larvae was compared using Pearson’s Chi-squared test, R (R Core Team, 2013). The number of SNPs per gene was compared using a T-test, R. To infer the population structure of the studied beetles, principal component analysis (PCA) was performed using Plink, and the ancestry was estimated using Admixture (Purcell et al., 2007; Zhou et al., 2011).
Results and Discussion
Phylogenetic Analysis Suggests Burkina Faso as the Original Source
On average, 54,358,906 paired reads (151 bp per read) were aligned to the SHB genome, and 5,400,117 SNPs were called in an individual beetle collected in China. After normalizing SHBs collected from Africa and America, 4,541,776 SNPs were identified in all samples, which were used for the phylogenetic analysis. The beetles collected from the United States were clustered into two groups, which formed the earliest branch of studied beetles (Figure 1). The results supported the previous studies that two haplotypes were formed due to multiple intrusions in the United States (Evans et al., 2003; Lounsberry et al., 2010). In addition, the early branch also indicates that the multiple intrusion routes and transporting pollination activities facilitated novel genotypes in the United States, leading to differentiated haplotypes. Indeed, we found that the beetles in the United States showed the highest fixation index with the ones in China compared with those in Africa (Table 1). In our study, the beetles were primarily clustered per country, which reflects that geographic proximity and substantial dispersal of beetles were not detected within their native region. However, occasional intrusions from Burkina Faso to the neighboring countries were indicated, which might be due to human-mediated transportation (Idrissou et al., 2019; Liu et al., 2021b). The PCA recapitulated the occasional dispersal of beetles from Burkina Faso to other regions (Figure 2). In our previous study, the SHBs in China were closely related to the ones in the Philippines based on a fragment of mitochondrial DNA, suggesting that the beetles that invaded the two countries were from the same source (Liu et al., 2021a). However, those beetles cannot be clustered with the ones in their native region. In the PCA plot, the beetles in China were distant from those of other regions. In the phylogenetic tree, the beetles in China and the ones in Burkina Faso formed a cluster, suggesting that the SHBs that invaded China shared a common ancestor with those in Burkina Faso. However, it does not necessarily indicate that SHBs in China were directly imported from Burkina Faso. The ancestry inference analysis also supported the close relationship between the beetles in China and those in Burkina Faso (Figure 2). In Asia, the beetles were also reported to invade South Korea. Using mitochondrial DNA, it was found that the beetles in Korea likely originated from the United States (Mohamadzade Namin et al., 2019). It seems clear that SHBs invaded Asia through multiple invasion paths.
FIGURE 1. Phylogenetic tree of 56 small hive beetles crossing native and novel regions. The tree was built on 4,541,776 SNPs along the genome. The beetles collected from the same region are highlighted with color. The beetles in the United States were clustered into two clusters. The beetles collected in China formed the deepest branch, which was closely related to those in Burkina Faso. The value on the branch indicates the bootstrap confidence. Honeybee was used as an outgroup to root the tree. The scale bar for branch length is shown at the bottom.
TABLE 1. Pairwise fixation index FST of the three beetle sources along the genome (mean ± SE). The beetles collected from China and the United States showed the highest divergence.
FIGURE 2. Structure analysis of the beetles. (A) PCA plot of the studied beetles. The beetles in China were distant from the ones in native and established regions. (B) Ancestral population inference of the studied beetles. The beetles in China may share an ancestor with the ones in Burkina Faso.
Gene Expression of Larvae and Adults
On average, 29,249,729 paired reads (151 bp per read) were aligned to SHBs. Out of 13,656 annotated protein-coding genes, 12,068 genes were detected in adults, and 11,042 genes were detected in larvae (Figure 3). Among the expressed genes, 203 genes were significantly upregulated in larvae, and 460 genes were significantly upregulated in adults, which significantly deviated from random (Pearson’s Chi-squared test, df = 1, p < 0.001). By quantifying the number of SNPs in each gene, we found that the 663 overexpressed genes in adults and larvae showed a significantly higher number of SNPs than the remaining genes (t-test, p < 0.05). Using the Kyoto Encyclopedia of Genes analysis, the distribution of genes in six main pathways was significantly different between adults and larvae (Pearson’s Chi-squared test, df = 5, p < 0.001). Among the six main pathways, metabolism was associated with the highest number of genes in both adults (31.0%) and larvae (40.1%). The larvae showed a higher number of upregulated genes associated with translation, replication, and repair than the adults. In the white-striped longhorn beetle, more genes were upregulated in adults than in larvae, and the antennae showed tissue-specific expressed olfactory genes (Yang et al., 2018). Similarly, in the coffee berry borer beetle, more genes were significantly upregulated in adults than in larvae. As expected, higher expression levels of larval-specific cuticle-binding proteins and chitinases were found in larvae (Noriega et al., 2019). In another pest beetle, Tribolium castaneum, a higher number of upregulated genes were found in larvae than in adults. Large numbers of ribosomal proteins were included, indicative of protein production variation between the life stages (Perkin and Oppert, 2019). In our data, the relative abundance of genes involved in genetic information processing was two-fold higher in larvae than in adults, which indicates the fast growth at this stage as supported by highly expressed metabolic genes.
FIGURE 3. Gene expression of small hive beetles. (A) Venn diagram of expressed genes in SHB adults and larvae. In the beetle genome, 13,656 protein-coding genes were annotated. Overall, 78% genes were expressed in both adults and larvae. (B) Significantly upregulated genes in adults. In total, 460 genes were significantly upregulated in adults compared with larvae SHBs, and the genes were enriched in environmental information processing. (C) Significantly upregulated genes in larvae. Overall, 203 genes were significantly upregulated in larvae compared with adults and the genes associated with translation, replication, and repair showed a high relative abundance compared with adults.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) are as follows: NCBI PRJNA782806 and PRJNA776042.
Author Contributions
WB conducted the experiment. QH designed the experiment. YL annotated the phylogenetic tree. JL and WH collected the beetles. WB, JL, YL, WH, JE, and QH organized the manuscript.
Funding
The project is supported by the initiation package of Jiangxi Agricultural University 050014/923230722 (QH) and the Science and Technology Program Project of Hainan Province ZDYF2021XDNY282 (WH).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
The authors appreciate the technical support of Dr. Lizhen Zhang.
References
Benda, N. D., Boucias, D., Torto, B., and Teal, P. (2008). Detection and Characterization ofKodamaea Ohmeriassociated with Small Hive beetleAethina Tumidainfesting Honey Bee Hives. J. Apicultural Res. 47, 194–201. doi:10.1080/00218839.2008.11101459
Cervancia, C. R., de Guzman, L. I., Polintan, E. A., Dupo, A. L. B., and Locsin, A. A. (2016). Current Status of Small Hive Beetle Infestation in the Philippines. J. Apicultural Res. 55, 74–77. doi:10.1080/00218839.2016.1194053
Chen, S., Zhou, Y., Chen, Y., and Gu, J. (2018). Fastp: an Ultra-fast All-In-One FASTQ Preprocessor. Bioinformatics 34, i884–i890. doi:10.1093/bioinformatics/bty560
Christmas, M. J., Jones, J. C., Olsson, A., Wallerman, O., Bunikis, I., Kierczak, M., et al. (2021). Genetic Barriers to Historical Gene Flow between Cryptic Species of Alpine Bumblebees Revealed by Comparative Population Genomics. Mol. Biol. Evol. 38, 3126–3143. doi:10.1093/molbev/msab086
Cornelissen, B., Neumann, P., and Ellis, J. D. (2020). Successful Pupation of Small Hive Beetle, Aethina Tumida (Coleoptera: Nitidulidae), in Greenhouse Substrates. J. Econ. Entomol. 113, 3032–3034. doi:10.1093/jee/toaa224
Cuthbertson, A. G. S., Wakefield, M. E., Powell, M. E., Marris, G., Anderson, H., Budge, G. E., et al. (2013). The Small Hive Beetle Aethina Tumida: A Review of its Biology and Control Measures. Curr. Zool. 59, 644–653. doi:10.1093/czoolo/59.5.644
Del Valle Molina, J. A. (2007). Small Hive Beetle Infestation (Aethina Tumida) in Mexico: Immediate Notification Report. Ref OIE, 6397.
Evans, J. D., McKenna, D., Scully, E., Cook, S. C., Dainat, B., Egekwu, N., et al. (2018). Genome of the Small Hive Beetle (Aethina Tumida, Coleoptera: Nitidulidae), a Worldwide Parasite of Social Bee Colonies, Provides Insights into Detoxification and Herbivory. Gigascience 7, giy138. doi:10.1093/gigascience/giy138
Evans, J. D., Pettis, J. S., Hood, W. M., and Shimanuki, H. (2003). Tracking an Invasive Honey Bee Pest: Mitochondrial DNA Variation in North American Small Hive Beetles. Apidologie 34, 103–109. doi:10.1051/apido:2003004
Eyer, M., Chen, Y. P., Schäfer, M. O., Pettis, J., and Neumann, P. (2009). Small Hive beetle,Aethina Tumida, as a Potential Biological Vector of Honeybee Viruses. Apidologie 40, 419–428. doi:10.1051/apido:2008051
Gillespie, P., Staples, J., King, C., Fletcher, M. J., and Dominiak, B. (2003). Small Hive Beetle Aethina Tumida (Murray) (Coleoptera: Nitidulidae) in New South Wales. Gen. Appl. Entomol. 32, 5–7.
Graham, J. R., Ellis, J. D., Carroll, M. J., and Teal, P. E. A. (2011). Aethina Tumida (Coleoptera: Nitidulidae) Attraction to Volatiles Produced by Apis mellifera (Hymenoptera: Apidae) and Bombus Impatiens (Hymenoptera: Apidae) Colonies. Apidologie 42, 326–336. doi:10.1007/s13592-011-0017-9
Hood, W. M. (2000). Overview of the Small Hive beetle,Aethina Tumida, in North America. Bee World 81, 129–137. doi:10.1080/0005772x.2000.11099483
Huang, Q., Lopez, D., and Evans, J. D. (2019). Shared and Unique Microbes between Small Hive Beetles (Aethina Tumida) and Their Honey Bee Hosts. Microbiologyopen 8, e899. doi:10.1002/mbo3.899
Huwiler, M., Papach, A., Cristina, E., Yañez, O., Williams, G. R., and Neumann, P. (2020). Deformed Wings of Small Hive Beetle Independent of Virus Infections and Mites. J. Invertebr. Pathol. 172, 107365. doi:10.1016/j.jip.2020.107365
Idrissou, F. O., Huang, Q., Yañez, O., and Neumann, P. (2019). International Beeswax Trade Facilitates Small Hive Beetle Invasions. Sci. Rep. 9, 10665. doi:10.1038/s41598-019-47107-6
Kanehisa, M., and Goto, S. (2000). KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 28, 27–30. doi:10.1093/nar/28.1.27
Kim, D., Langmead, B., and Salzberg, S. L. (2015). HISAT: a Fast Spliced Aligner with Low Memory Requirements. Nat. Methods 12, 357–360. doi:10.1038/nmeth.3317
Kofler, R., Pandey, R. V., and Schlötterer, C. (2011). PoPoolation2: Identifying Differentiation between Populations Using Sequencing of Pooled DNA Samples (Pool-Seq). Bioinformatics 27, 3435–3436. doi:10.1093/bioinformatics/btr589
Li, H., and Durbin, R. (2009). Fast and Accurate Short Read Alignment with Burrows-Wheeler Transform. Bioinformatics 25, 1754–1760. doi:10.1093/bioinformatics/btp324
Li, H., Handsaker, B., Wysoker, A., Fennell, T., Ruan, J., Homer, N., et al. (2009). The Sequence Alignment/Map Format and SAMtools. Bioinformatics 25, 2078–2079. doi:10.1093/bioinformatics/btp352
Liu, Y., Han, W., Gao, J., Su, S., Beaurepaire, A., Yañez, O., et al. (2021a). Out of Africa: Novel Source of Small Hive Beetles Infesting Eastern and Western Honey Bee Colonies in China. J. Apicultural Res. 60, 108–110. doi:10.1080/00218839.2020.1816686
Liu, Y., Henkel, J., Beaurepaire, A., Evans, J. D., Neumann, P., and Huang, Q. (2021b). Comparative Genomics Suggests Local Adaptations in the Invasive Small Hive Beetle. Ecol. Evol. 11, 15780–15791. doi:10.1002/ece3.8242
Lounsberry, Z., Spiewok, S., Pernal, S. F., Sonstegard, T. S., Hood, W. M., Pettis, J., et al. (2010). Worldwide Diaspora of Aethina Tumida (Coleoptera: Nitidulidae), a Nest Parasite of Honey Bees. Ann. Entomol. Soc. Am. 103, 671–677. doi:10.1603/an10027
McKenna, A., Hanna, M., Banks, E., Sivachenko, A., Cibulskis, K., Kernytsky, A., et al. (2010). The Genome Analysis Toolkit: a MapReduce Framework for Analyzing Next-Generation DNA Sequencing Data. Genome Res. 20, 1297–1303. doi:10.1101/gr.107524.110
Mckenna, D. D., Wild, A. L., Kanda, K., Bellamy, C. L., Beutel, R. G., Caterino, M. S., et al. (2015). The Beetle Tree of Life Reveals that C Oleoptera Survived End‐ P Ermian Mass Extinction to Diversify during the C Retaceous Terrestrial Revolution. Syst. Entomol. 40, 835–880. doi:10.1111/syen.12132
Mohamadzade Namin, S., Koh, Y., Osabutey, A. F., and Jung, C. (2019). Invasion Pathway of the Honeybee Pest, Small Hive Beetle, Aethina Tumida (Coleoptera: Nitidulidae) in the Republic of Korea Inferred by Mitochondrial DNA Sequence Analysis. J. Asia-Pacific Entomol. 22, 963–968. doi:10.1016/j.aspen.2019.07.008
Neumann, P., and Elzen, P. J. (2004). The Biology of the Small Hive Beetle (Aethina Tumida, Coleoptera: Nitidulidae): Gaps in Our Knowledge of an Invasive Species. Apidologie 35, 229–247. doi:10.1051/apido:2004010
Neumann, P., Naef, J., Crailsheim, K., Crewe, R. M., and Pirk, C. W. W. (2015). Hit‐and‐run Trophallaxis of Small Hive Beetles. Ecol. Evol. 5, 5478–5486. doi:10.1002/ece3.1806
Neumann, P., Pettis, J. S., and Schäfer, M. O. (2016). Quo Vadis Aethina Tumida? Biology and Control of Small Hive Beetles. Apidologie 47, 427–466. doi:10.1007/s13592-016-0426-x
Nguyen, L.-T., Schmidt, H. A., von Haeseler, A., and Minh, B. Q. (2015). IQ-TREE: A Fast and Effective Stochastic Algorithm for Estimating Maximum-Likelihood Phylogenies. Mol. Biol. Evol. 32, 268–274. doi:10.1093/molbev/msu300
Noriega, D. D., Arias, P. L., Barbosa, H. R., Arraes, F. B. M., Ossa, G. A., Villegas, B., et al. (2019). Transcriptome and Gene Expression Analysis of Three Developmental Stages of the Coffee berry Borer, Hypothenemus Hampei. Sci. Rep. 9, 12804. doi:10.1038/s41598-019-49178-x
Ortiz, E. M. (2019). vcf2phylip v2.0: Convert a VCF Matrix into Several Matrix Formats for Phylogenetic Analysis (v2.0). Zenodo. doi:10.5281/zenodo.2540861
Palmeri, V., Scirtò, G., Malacrinò, A., Laudani, F., and Campolo, O. (2014). A Scientific Note on a New Pest for European Honeybees: First Report of Small Hive Beetle Aethina Tumida, (Coleoptera: Nitidulidae) in Italy. Apidologie 46, 527–529. doi:10.1007/s13592-014-0343-9
Perkin, L. C., and Oppert, B. (2019). Gene Expression in Tribolium castaneum Life Stages: Identifying a Species-specific Target for Pest Control Applications. PeerJ 7, e6946. doi:10.7717/peerj.6946
Purcell, S., Neale, B., Todd-Brown, K., Thomas, L., Ferreira, M. A. R., Bender, D., et al. (2007). PLINK: a Tool Set for Whole-Genome Association and Population-Based Linkage Analyses. Am. J. Hum. Genet. 81, 559–575. doi:10.1086/519795
R Core Team (2013). R: A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation Statistical Computing. Available at: http://www.R-project.org/.
Robinson, M. D., McCarthy, D. J., and Smyth, G. K. (2010). edgeR: a Bioconductor Package for Differential Expression Analysis of Digital Gene Expression Data. Bioinformatics 26, 139–140. doi:10.1093/bioinformatics/btp616
Torto, B., Boucias, D. G., Arbogast, R. T., Tumlinson, J. H., and Teal, P. E. A. (2007). Multitrophic Interaction Facilitates Parasite-Host Relationship between an Invasive Beetle and the Honey Bee. Proc. Natl. Acad. Sci. U.S.A. 104, 8374–8378. doi:10.1073/pnas.0702813104
Torto, B., Suazo, A., Alborn, H., Tumlinson, J. H., and Teal, P. E. A. (2005). Response of the Small Hive Beetle (Aethina Tumida) to a Blend of Chemicals Identified from Honeybee (Apis mellifera) Volatiles. Apidologie 36, 523–532. doi:10.1051/apido:2005038
Van der Auwera, G. A., Carneiro, M. O., Hartl, C., Poplin, R., Del Angel, G., Levy-Moonshine, A., et al. (2013). From FastQ Data to High Confidence Variant Calls: the Genome Analysis Toolkit Best Practices Pipeline. Curr. Protoc. Bioinformatics 43, 11–33. doi:10.1002/0471250953.bi1110s43
van Engelsdorp, D., Underwood, R., Caron, D., and Hayes, J. J. (2007). Estimate of Managed colony Losses in the winter of 2006-2007: A Report Commissioned by the Apiary Inspectors of America. Am. Bee J. 147, 599–603.
Yang, H., Cai, Y., Zhuo, Z., Yang, W., Yang, C., Zhang, J., et al. (2018). Transcriptome Analysis in Different Developmental Stages of Batocera Horsfieldi (Coleoptera: Cerambycidae) and Comparison of Candidate Olfactory Genes. PLoS One 13, e0192730. doi:10.1371/journal.pone.0192730
Zhao, H., Yang, S., Liu, J., Huang, W., Ji, C., Ren, Q., et al. (2020). First Detection of Small Hive Beetle Aethina Tumida Murray (Coleoptera: Nitidulidae) Infesting Eastern Honeybee, Apis cerana Fabricius (Hymenoptera: Apidae), in China. Sociobiology 67, 126–128. doi:10.13102/sociobiology.v67i1.4381
Keywords: small hive beetle, phylogeny, SNP, evolution, honeybee
Citation: Bai WF, Liu J, Liu Y, Han W, Evans JD and Huang Q (2022) Phylogenetic Analysis of Small Hive Beetles From Native to Introduced Populations. Front. Genet. 13:900795. doi: 10.3389/fgene.2022.900795
Received: 21 March 2022; Accepted: 15 April 2022;
Published: 19 May 2022.
Edited by:
Rui Guo, Fujian Agriculture and Forestry University, ChinaCopyright © 2022 Bai, Liu, Liu, Han, Evans and Huang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Qiang Huang, cWlhbmctaHVhbmdAbGl2ZS5jb20=
†These authors have contributed equally to this work