Tetracycline-induced gut community dysbiosis and Israeli Acute Paralysis Virus infection synergistically negatively affect honeybees.

Antibiotics are frequently employed to control bacterial diseases in honeybees, but their broad-spectrum action can disrupt the delicate balance of the gut microbiome, leading to dysbiosis. This imbalance in the gut microbiota of honeybees adversely affects their physiological health and weakens their resistance to pathogens, including viruses that significantly threaten honeybee health. In this study, we investigated whether tetracycline-induced gut microbiome dysbiosis promotes the replication of Israeli acute paralysis virus (IAPV), a key virus associated with colony losses and whether IAPV infection exacerbates gut microbiome dysbiosis. Our results demonstrated that tetracycline-induced gut microbiome dysbiosis increases the susceptibility of honeybees to IAPV infection. The viral titer in worker bees with antibiotic-induced gut microbiome dysbiosis prior to IAPV inoculation was significantly higher than in those merely inoculated with IAPV. Furthermore, we observed a synergistic effect between tetracycline and IAPV on the disruption of the honeybee gut microbiome balance. The progression of IAPV replication could, in turn, exacerbate antibiotic-induced gut microbiome dysbiosis in honeybees. Our research provides novel insights into the role of the gut microbiota in host-virus interactions, emphasizing the complex interplay between antibiotic use, gut microbiome health, and viral susceptibility in honeybees. We highlight the crucial role of a balanced gut microbiota in honey bees for their immune response against pathogens and emphasize the importance of careful, safe antibiotic use in beekeeping to protect these beneficial microbes.

Chemical Analyses and Antimicrobial Activity of Nine Kinds of Unifloral Chinese Honeys Compared to Manuka Honey (12+ and 20+).

Honey has good antimicrobial properties and can be used for medical treatment. The antimicrobial properties of unifloral honey varieties are different. In this study, we evaluated the antimicrobial and antioxidant activities of nine kinds of Chinese monofloral honeys. In addition, headspace gas chromatography-ion mobility spectrometry (HS-GC-IMS) technology was used to detect their volatile components. The relevant results are as follows: 1. The agar diffusion test showed that the diameter of inhibition zone against Staphylococcus aureus of Fennel honey (21.50 ± 0.41 mm), Agastache honey (20.74 ± 0.37 mm), and Pomegranate honey (18.16 ± 0.11 mm) was larger than that of Manuka 12+ honey (14.27 ± 0.10 mm) and Manuka 20+ honey (16.52 ± 0.12 mm). The antimicrobial activity of Chinese honey depends on hydrogen peroxide. 2. The total antioxidant capacity of Fennel honey, Agastache honey, and Pomegranate honey was higher than that of other Chinese honeys. There was a significant positive correlation between the total antioxidant capacity and the total phenol content of Chinese honey (r = 0.958). The correlation coefficient between the chroma value of Chinese honey and the total antioxidant and the diameter of inhibition zone was 0.940 and 0.746, respectively. The analyzed dark honeys had better antimicrobial and antioxidant activities. 3. There were significant differences in volatile components among Fennel honey, Agastache honey, Pomegranate honey, and Manuka honey. Hexanal-D and Heptanol were the characteristic components of Fennel honey and Pomegranate honey, respectively. Ethyl 2-methylbutyrate and 3-methylpentanoic acids were the unique compounds of Agastache honey. The flavor fingerprints of the honey samples from different plants can be successfully built using HS-GC-IMS and principal component analysis (PCA) based on their volatile compounds. Fennel honey, Agastache honey, and Pomegranate honey are Chinese honey varieties with excellent antimicrobial properties, and have the potential to be developed into medical grade honey.

Plant microRNAs in larval food regulate honeybee caste development.

The major environmental determinants of honeybee caste development come from larval nutrients: royal jelly stimulates the differentiation of larvae into queens, whereas beebread leads to worker bee fate. However, these determinants are not fully characterized. Here we report that plant RNAs, particularly miRNAs, which are more enriched in beebread than in royal jelly, delay development and decrease body and ovary size in honeybees, thereby preventing larval differentiation into queens and inducing development into worker bees. Mechanistic studies reveal that amTOR, a stimulatory gene in caste differentiation, is the direct target of miR162a. Interestingly, the same effect also exists in non-social Drosophila. When such plant RNAs and miRNAs are fed to Drosophila larvae, they cause extended developmental times and reductions in body weight and length, ovary size and fecundity. This study identifies an uncharacterized function of plant miRNAs that fine-tunes honeybee caste development, offering hints for understanding cross-kingdom interaction and co-evolution.

Occurrence of multiple honeybee viruses in the ectoparasitic mites Varroa spp. in Apis cerana colonies.

In this study, we investigated the prevalence of honeybee viruses in Varroa destructor and Varroa underwoodi infesting Apis cerana colonies in China. Deformed wing virus (DWV) was the most prevalent virus in these two mite species, followed by Israeli acute paralysis virus (IAPV), Black queen cell virus (BQCV), Kashmir bee virus (KBV), Chronic bee paralysis virus (CBPV), Apis mellifera filamentous virus (AmFV) and Sacbrood virus (SBV) in V. destructor, while in V. underwoodi, it was followed by CBPV, AmFV, BQCV, IAPV and KBV. In addition, multiple viruses were commonly detectable in both mite species.

A New Propolis Type from Changbai Mountains in North-east China: Chemical Composition, Botanical Origin and Biological Activity.

Propolis is a bee product with a wide range of biological activities and its chemical compounds depend highly on the type of plant accessible to the bees. The Changbai Mountains are a major mountain range in Northeast China and are one of the major bee product-producing areas in China. In this study, we evaluated the total phenolic acids and flavonoid contents as well as the antioxidant activity of propolis sampled from the Changbai Mountains area (CBM). We identified the major compounds and qualified their contents by HPLC-ESI/MS and HPLC-UV, and found that the content of p-coumaric acid and an unknown peak (CBE) in CBM propolis was higher than in propolis from other parts of China. The unknown compound CBE was isolated, purified, and identified as benzyl p-coumarate by MS and NMR. Possible plant sources of CBM propolis are Populus davidiana dode and Populus simonii Carr, which widely distributed in the Changbai Mountains area. CBM propolis is a new propolis type, that could be an excellent raw material for health foods and pharmaceuticals.

Authentication of Apis cerana Honey and Apis mellifera Honey Based on Major Royal Jelly Protein 2 Gene.

In Asia, honey is mainly produced by Apis mellifera and Apis cerana. However, the price of A. cerana honey is usually much higher than A. mellifera honey. Seeing considerable profits, some dishonest companies and beekeepers mislabel A. mellifera honey as A. cerana honey or incorporate A. mellifera honey into A. cerana honey. In the present study, we developed methods to discriminate A. cerana honey from A. mellifera honey based on the MRJP2 (major royal jelly protein 2) gene. Two pairs of species-specific primers were designed. The amplification products of A. cerana and A. mellifera were 212 and 560 bp, respectively. As little as one percent incorporation of A. mellifera honey in the mixture can be detected by duplex PCR. Additionally, another method based on the melt curve analysis using the same primers was also developed, allowing a rapid discrimination of real-time PCR product of different species. Our study shows that the entomological authentication of honey samples can be identified by nuclear genes other than mitochondrial genes and this extends the possibility of gene selection in identification. The authentication system we proposed could be a useful tool for discriminating A. cerana honey from A. mellifera honey.

Sublethal doses of neonicotinoid imidacloprid can interact with honey bee chemosensory protein 1 (CSP1) and inhibit its function.

As a frequently used neonicotinoid insecticide, imidacloprid can impair the chemoreceptive behavior of honey bees even at sublethal doses, while the physiochemical mechanism has not been further revealed. Here, multiple fluorescence spectra, thermodynamic method, and molecular docking were used to study the interaction and the functional inhibition of imidacloprid to the recombinant CSP1 protein in Asian honey bee, Apis cerana. The results showed that the fluorescence intensity (λem = 332 nm) of CSP1 could be significantly quenched by imidacloprid in a dynamic mode. During the quenching process, ΔH > 0, ΔS > 0, indicating that the acting forces of imidacloprid with CSP1 are mainly hydrophobic interactions. Synchronous fluorescence showed that the fluorescence of CSP1 was mainly derived from tryptophan, and the hydrophobicity of tryptophan decreased with the increase of imidacloprid concentration. Molecular docking predicted the optimal pose and the amino acid composition of the binding process. Circular dichroism (CD) spectra showed that imidacloprid reduced the α-helix of CSP1 and caused the extension of the CSP1 peptide chain. In addition, the binding of CSP1 to floral scent ß-ionone was inhibited by nearly 50% of the apparent association constant (KA) in the presence of 0.28-2.53 ng/bee of imidacloprid, and the inhibition rate of nearly 95% at 3.75 ng/bee of imidacloprid at sublethal dose level. This study initially revealed the molecular physiochemical mechanism that sublethal doses of neonicotinoid still interact and inhibit the physiological function of the honey bees' chemoreceptive system.

Evidence of Apis cerana Sacbrood virus Infection in Apis mellifera.

Sacbrood virus(SBV) is one of the most destructive viruses in the Asian honeybee Apis cerana but is much less destructive in Apis mellifera In previous studies, SBV isolates infecting A. cerana(AcSBV) and SBV isolates infecting A. mellifera(AmSBV) were identified as different serotypes, suggesting a species barrier in SBV infection. In order to investigate this species isolation, we examined the presence of SBV infection in 318A. mellifera colonies and 64A. cerana colonies, and we identified the genotypes of SBV isolates. We also performed artificial infection experiments under both laboratory and field conditions. The results showed that 38A. mellifera colonies and 37A. cerana colonies were positive for SBV infection. Phylogenetic analysis based on RNA-dependent RNA polymerase (RdRp) gene sequences indicated that A. cerana isolates and most A. mellifera isolates formed two distinct clades but two strains isolated fromA. mellifera were clustered with theA. cerana isolates. In the artificial-infection experiments, AcSBV negative-strand RNA could be detected in both adult bees and larvae ofA. mellifera, although there were no obvious signs of the disease, demonstrating the replication of AcSBV inA. mellifera Our results suggest that AcSBV is able to infectA. melliferacolonies with low prevalence (0.63% in this study) and pathogenicity. This work will help explain the different susceptibilities ofA. cerana and A. melliferato sacbrood disease and is potentially useful for guiding beekeeping practices.

High Royal Jelly-Producing Honeybees (Apis mellifera ligustica) (Hymenoptera: Apidae) in China.

China is the largest producer and exporter of royal jelly (RJ) in the world, supplying >90% of the global market. The high production of RJ in China is principally owing to the high RJ-producing lineage of honeybees (Apis mellifera ligustica Spinola, 1806) established by beekeepers in the 1980s. We describe the development of high royal jelly-producing honeybees and the management of this lineage today. Previous research and recent advances in the genetic characterization of this lineage, and the molecular markers and mechanisms associated with high RJ production are summarized. The gaps in our knowledge and prospects for future research are also highlighted.

Antioxidant activities and molecular mechanisms of the ethanol extracts of Baccharis propolis and Eucalyptus propolis in RAW64.7 cells.

Context Numerous studies have reported that propolis possesses strong antioxidant activities. However, their antioxidant molecular mechanisms are unclear. Objective We utilize ethanol extracts of Chinese propolis (EECP) as a reference to compare ethanol extracts of Eucalyptus propolis (EEEP) with ethanol extracts of Baccharis propolis (EEBGP) based on their antioxidant capacities and underlying molecular mechanisms. Materials and methods HPLC and chemical analysis are utilized to evaluate compositions and antioxidant activities. ROS-eliminating effects of EEBGP (20-75 µg/mL), EEEP (1.25-3.75 µg/mL) and EECP (1.25-5 µg/mL) are also determined by flow cytometry analysis. Moreover, we compared antioxidant capacities by determining their effects on expressions of antioxidant genes in RAW264.7 cells with qRT-PCR, western blot and confocal microscopy analysis. Results EEBGP mainly contains chlorogenic acid (8.98 ± 0.86 mg/g), kaempferide (11.18 ± 8.31 mg/g) and artepillin C (107.70 ± 10.86 mg/g), but EEEP contains 10 compositions, whereas EECP contains 17 compositions. Meantime, although EEEP shows DPPH (IC50 19.55 ± 1.28), ABTS (IC50 20.0 ± 0.31) and reducing power (2.70 ± 0.08 mmol TE/g) better than EEBGP's DPPH (IC50 43.85 ± 0.54), ABTS (IC50 38.2 ± 0.33) and reducing power (1.53 ± 0.05 mmol TE/g), EEBGP exerts much higher ROS inhibition rate (40%) than EEEP (under 20%). Moreover, EEBGP strengthen antioxidant system by activating p38/p-p38 and Erk/p-Erk kinase via accelerating nucleus translocation of Nrf2. EEEP and EECP improve antioxidant gene expression only via Erk/p-Erk kinase-Nrf2 signalling pathway. Discussion and conclusion EEBGP and EEEP exert antioxidant activities via different molecular mechanisms, which may depend on chemical compositions.

Detection of Spiroplasma melliferum in honey bee colonies in the US.

Spiroplasma infections in honey bees have been reported in Europe and Asia quite recently, due to intensive studies on the epidemiology of honey bee diseases. The situation in the US is less well analyzed. Here, we examined the honey bee colonies in Beltsville, MD, where Spiroplasmamelliferum was originally reported and found S. melliferum infection in honey bees. Our data showed high variation of S. melliferum infection in honey bees with a peak prevalence in May during the course of one-year study period. The colony prevalence increased from 5% in February to 68% in May and then decreased to 25% in June and 22% in July. Despite that pathogenicity of spiroplasmas in honey bee colonies remains to be determined, our results indicated that spiroplasma infections need to be included for the consideration of the impacts on honey bee health.

Spore Loads May Not be Used Alone as a Direct Indicator of the Severity of Nosema ceranae Infection in Honey Bees Apis mellifera (Hymenoptera:Apidae).

Nosema ceranae Fries et al., 1996, a microsporidian parasite recently transferred from Asian honey bees Apis cerana F., 1793, to European honey bees Apis mellifera L., 1758, has been suspected as one of the major culprits of the worldwide honey bee colony losses. Spore load is a commonly used criterion to describe the intensity of Nosema infection. In this study, by providing Nosema-infected bees with sterilized pollen, we confirmed that pollen feeding increased the spore loads of honey bees by several times either in the presence or absence of a queen. By changing the amount of pollen consumed by bees in cages, we showed that spore loads increased with an increase in pollen consumption. Nosema infections decrease honey bee longevity and transcription of vitellogenin, either with or without pollen feeding. However, the reduction of pollen consumption had a greater impact on honey bee longevity and vitellogenin level than the increase of spore counts caused by pollen feeding. These results indicate that spore loads may not be used alone as a direct indicator of the severity of N. ceranae infection in honey bees.

Evaluation of cage designs and feeding regimes for honey bee (Hymenoptera: Apidae) laboratory experiments.

The aim of this study was to improve cage systems for maintaining adult honey bee (Apis mellifera L.) workers under in vitro laboratory conditions. To achieve this goal, we experimentally evaluated the impact of different cages, developed by scientists of the international research network COLOSS (Prevention of honey bee COlony LOSSes), on the physiology and survival of honey bees. We identified three cages that promoted good survival of honey bees. The bees from cages that exhibited greater survival had relatively lower titers of deformed wing virus, suggesting that deformed wing virus is a significant marker reflecting stress level and health status of the host. We also determined that a leak- and drip-proof feeder was an integral part of a cage system and a feeder modified from a 20-ml plastic syringe displayed the best result in providing steady food supply to bees. Finally, we also demonstrated that the addition of protein to the bees' diet could significantly increase the level ofvitellogenin gene expression and improve bees' survival. This international collaborative study represents a critical step toward improvement of cage designs and feeding regimes for honey bee laboratory experiments.

The Braconid Syntretomorpha szaboi Papp Is Posing a Great Threat to the Eastern Honeybee, Apis cerana.

The expansion of pathogen distribution may result in a new threat to the host. The braconid Syntretomorpha szaboi Papp is an obligate parasite that targets Apis cerana, the Eastern honeybee, engaging in endoparasitism by ovipositing eggs inside the host bee. Although S. szaboi has been documented in India and in various regions across China, its epidemiological data are notably lacking. In this study, we summarized the distribution of S. szaboi based on the available literature and described the symptoms of infested honeybee workers. We also investigated the infestation rate in 36 apiaries in Zhejiang Province, China, after a new occurrence of the parasite was reported in these regions in 2020. A rapid increase in infestation rate was found from the year 2021 to 2022, reaching 53.88% at the colony level of the sampled colonies in the Jinhua and Wenzhou apiaries. The infestation rate at an individual level in positive colonies reached an average of 26%. A monthly survey showed high seasonal variation in S. szaboi infestation, with the peak occurring from May to August. These results suggest that S. szaboi poses a great threat to A. cerana. Further research is needed to elucidate its epidemiology and pathology and to develop disease prevention and control strategies.

TALEN-mediated homologous-recombination-based fibroin light chain in-fusion expression system in Bombyx mori.

Silkworm was the first domesticated insect and has important economic value. It has also become an ideal model organism with applications in genetic and expression studies. In recent years, the use of transgenic strategies has made the silkworm silk gland an attractive bioreactor for the production of recombinant proteins, in particular, piggyBac-mediated transgenes. However, owing to differences in regulatory elements such as promoters, the expression levels of exogenous proteins have not reached expectations. Here, we used targeted gene editing to achieve site-specific integration of exogenous genes on genomic DNA and established the fibroin light chain (FibL) in-fusion expression system by TALEN-mediated homology-directed recombination. First, the histidine-rich cuticular protein (CP) was successfully site-directed inserted into the native FibL, and the FibL-CP fusion gene was correctly transcribed and expressed in the posterior silk gland under the control of the endogenous FibL promoter, with a protein expression level comparable with that of the native FibL protein. Moreover, we showed based on molecular docking that the fusion of FibL with cuticular protein may have a negative effect on disulfide bond formation between the C-terminal domain of fibroin heavy chain (FibH) and FibL-CP, resulting in abnormal spinning and cocoon in homozygotes, indicating a significant role of FibL in silk protein formation and secretion. Our results demonstrate the feasibility of using the FibL fusion system to express exogenous proteins in silkworm. We expect that this bioreactor system will be used to produce more proteins of interest, expanding the application value of the silk gland bioreactor.

An innovative molecular approach towards the cost-effective entomological authentication of honey.

Honey authentication and traceability are crucial not only for economic purposes but also for ensuring safety. However, the widespread adoption of cutting-edge technologies in practical applications has been hampered by complex, time-consuming sample pre-treatment processes, the need for skilled personnel, and substantial associated expenses. This study aimed to develop a simple and cost-effective molecular technique to verify the entomological source of honey. By utilizing newly designed primers, we successfully amplified the mitochondrial 16S ribosomal RNA gene of honey bees from honey, confirming the high quality of the extracted DNA. Employing RFLP analysis with AseI endonuclease, species-specific restriction patterns were generated for honey derived from six closely related honey bees of the Apis genus. Remarkably, this method was proven equally effective in identifying heat-treated and aged honey by presenting the same RFLP profiles as raw honey. As far as we know, this is the initial research of the simultaneous differentiation of honey from closely related honey bee species using the restriction endonuclease AseI and mitochondrial 16S rRNA gene fragments. As a result, it holds tremendous potential as a standardized guideline for regulatory agencies to ascertain the insect origins of honey and achieve comprehensive traceability.

Genome sequencing and comparative genomics of honey bee microsporidia, Nosema apis reveal novel insights into host-parasite interactions.

BACKGROUND: The microsporidia parasite Nosema contributes to the steep global decline of honey bees that are critical pollinators of food crops. There are two species of Nosema that have been found to infect honey bees, Nosema apis and N. ceranae. Genome sequencing of N. apis and comparative genome analysis with N. ceranae, a fully sequenced microsporidia species, reveal novel insights into host-parasite interactions underlying the parasite infections. RESULTS: We applied the whole-genome shotgun sequencing approach to sequence and assemble the genome of N. apis which has an estimated size of 8.5 Mbp. We predicted 2,771 protein- coding genes and predicted the function of each putative protein using the Gene Ontology. The comparative genomic analysis led to identification of 1,356 orthologs that are conserved between the two Nosema species and genes that are unique characteristics of the individual species, thereby providing a list of virulence factors and new genetic tools for studying host-parasite interactions. We also identified a highly abundant motif in the upstream promoter regions of N. apis genes. This motif is also conserved in N. ceranae and other microsporidia species and likely plays a role in gene regulation across the microsporidia. CONCLUSIONS: The availability of the N. apis genome sequence is a significant addition to the rapidly expanding body of microsprodian genomic data which has been improving our understanding of eukaryotic genome diversity and evolution in a broad sense. The predicted virulent genes and transcriptional regulatory elements are potential targets for innovative therapeutics to break down the life cycle of the parasite.

Geographical influences on content of 10-hydroxy-trans-2-decenoic acid in royal jelly in China.

The content of 10-hydroxy-trans-2-decenoic acid (10-HDA), a marker compound in royal jelly (RJ), is the most important criterion in grading RJ for commercial trade and varies with its origin. To identify the effect of geographical origin on 10-HDA content in RJ, 138 samples were collected from 19 provinces of China (divided into three groups) produced by either Apis mellifera ligustica Spinola, 1806 or a hybrid of A. m. ligustica and Apis mellifera carnica Pollman, 1879 and analyzed for moisture, sugar, crude protein, ash, acid, and 10-HDA concentration. The results show that RJ from western China has a significantly higher 10-HDA level (2.01 +/- 0.05%) than those from northeastern (1.87 +/- 0.05%) and eastern (1.75 +/- 0.03%) China. RJ secreted by hybrid bees contained more 10-HDA (1.89 +/- 0.03%) than that secreted by A. m. ligustica (1.78 +/- 0.03%). The 10-HDA content of RJ produced during flowering of rape (Brassica campestris L.), lime (Tilia amurensis Ruprecht), and vitex (Vitex negundo L. variety heterophylla (Franch.) Rehder) was 1.92, 1.80, and 1.68%, respectively. The results would be helpful during the process of price determination of RJ by providing some basis of geographical, bee strain, and botanical information for commercial trade.

Genomic signatures underlying the oogenesis of the ectoparasitic mite Varroa destructor on its new host Apis mellifera.

INTRODUCTION: Host shift of parasites may have devastating effects on the novel hosts. One remarkable example is that of the ectoparasitic mite Varroa destructor, which has shifted its host from Eastern honey bees (Apis cerana) to Western honey bees (Apis mellifera) and posed a global threat to apiculture. OBJECTIVES: To identify the genetic factors underlying the reproduction of host-shifted V. destructor on the new host. METHODS: Genome sequencing was conducted to construct the phylogeny of the host-shifted and non-shifted mites and to screen for genomic signatures that differentiated them. Artificial infestation experiment was conducted to compare the reproductive difference between the mites, and transcriptome sequencing was conducted to find differentially expressed genes (DEGs) during the reproduction process. RESULTS: The host-shifted and non-shifted V. destructor mites constituted two genetically distinct lineages, with 15,362 high-FST SNPs identified between them. Oogenesis was upregulated in host-shifted mites on the new host A. mellifera relative to non-shifted mites. The transcriptomes of the host-shifted and non-shifted mites differed significantly as early as 1h post-infestation. The DEGs were associated with nine genes carrying nonsynonymous high-FST SNPs, including mGluR2-like, Lamb2-like and Vitellogenin 6-like, which were also differentially expressed, and eIF4G, CG5800, Dap160 and Sas10, which were located in the center of the networks regulating the DEGs based on protein-protein interaction analysis. CONCLUSIONS: The annotated functions of these genes were all associated with oogenesis. These genes appear to be the key genetic determinants of the oogenesis of host-shifted mites on the new host. Further study of these candidate genes will help elucidate the key mechanism underlying the success of host shifts of V. destructor.

Population Structure, Demographic History, and Adaptation of Giant Honeybees in China Revealed by Population Genomic Data.

There have been many population-based genomic studies on human-managed honeybees (Apis mellifera and Apis cerana), but there has been a notable lack of analysis with regard to wild honeybees, particularly in relation to their evolutionary history. Nevertheless, giant honeybees have been found to occupy distinct habitats and display remarkable characteristics, which are attracting an increased amount of attention. In this study, we de novo sequenced and then assembled the draft genome sequence of the Himalayan giant honeybee, Apis laboriosa. Phylogenetic analysis based on genomic information indicated that A. laboriosa and its tropical sister species Apis dorsata diverged ∼2.61 Ma, which supports the speciation hypothesis that links A. laboriosa to geological changes throughout history. Furthermore, we re-sequenced A. laboriosa and A. dorsata samples from five and six regions, respectively, across their population ranges in China. These analyses highlighted major genetic differences for Tibetan A. laboriosa as well as the Hainan Island A. dorsata. The demographic history of most giant honeybee populations has mirrored glacial cycles. More importantly, contrary to what has occurred among human-managed honeybees, the demographic history of these two wild honeybee species indicates a rapid decline in effective population size in the recent past, reflecting their differences in evolutionary histories. Several genes were found to be subject to selection, which may help giant honeybees to adapt to specific local conditions. In summary, our study sheds light on the evolutionary and adaptational characteristics of two wild giant honeybee species, which was useful for giant honeybee conservation.