The Asian Honey Bee (Apis cerana) is Significantly in Decline

Conservation genomics analysis reveals recent population decline and possible causes in bumblebee Bombus opulentus

Bumblebees are a genus of pollinators (Bombus) that play important roles in natural ecosystem and agricultural production. Several bumblebee species have been recorded as under population decline, and the proportion of species experiencing population decline within subgenus Thoracobombus is higher than average. Bombus opulentus is 1 species in Thoracobombus, but little is known about its recent population dynamics. Here, we employed conservation genomics methods to investigate the population dynamics of B. opulentus during the recent past and identify the likely environmental factors that may cause population decline. Firstly, we placed the scaffold-level of B. opulentus reference genome sequence onto chromosome-level using Hi-C technique. Then, based on this reference genome and whole-genome resequencing data for 51 B. opulentus samples, we reconstructed the population structure and effective population size (Ne) trajectories of B. opulentus and identified genes that were under positive selection. Our results revealed that the collected B. opulentus samples could be divided into 2 populations, and 1 of them experienced a recent population decline; the declining population also exhibited lower genetic diversity and higher inbreeding levels. Genes related to high-temperature tolerance, immune response, and detoxication showed signals of positive selection in the declining population, suggesting that climate warming and pathogen/pesticide exposures may contribute to the decline of this B. opulentus population. Taken together, our study provided insights into the demography of B. opulentus populations and highlighted that populations of the same bumblebee species could have contrasting Ne trajectories and population decline could be caused by a combination of various stressors.

Combined Toxic Effects of Lead and Glyphosate on Apis cerana cerana

Glyphosate (GY) is the most widely used herbicide in agriculture worldwide. Lead is a common heavy metal in the natural environment. Honeybees, as pollinators, are exposed to these pollutants. So far, few reports have evaluated the toxic effects of GY mixed with heavy metals on honeybees (Apis cerana cerana). This study found that the acute toxicity of lead (LC50 = 1083 mg/L) is much greater than that of GY (LC50 = 4764 mg/L) at 96 h. The acute toxicities of the mixed substances were as follows: LC50 = 621 mg/L of lead and LC50 = 946 mg/L of GY. The combination of lead and GY was more toxic than either of the individual substances alone. Compared to the individual toxicity, combined treatment significantly affected the bees' learning and cognitive abilities and changed the relative expression of genes related to immune defense and detoxification metabolism in A. c. cerana. The combination of lead and GY seriously affected the behavior and physiology of the studied honeybees. This study provides basic data for further research on the combined effects of GY and heavy metals on bee health. It also serves as a reference for effective colony protection.

Enhancing Honey Bee Health: Evaluating Pollen Substitute Diets in Field and Cage Experiments

Honey bees (Apis mellifera L.) play vital roles as agricultural pollinators and honey producers. However, global colony losses are increasing due to multiple stressors, including malnutrition. Our study evaluated the effects of four pollen substitute diets (Diet 1, Diet 2, Diet 3, and Control) through field and cage experiments, analyzing 11 parameters and 21 amino acids. Notably, Diet 1 demonstrated significantly superior performance in the field experiment, including the number of honey bees, brood area, consumption, preference, colony weight, and honey production. In the cage experiment, Diet 1 also showed superior performance in dried head and thorax weight and vitellogenin (vg) gene expression levels. Canonical discriminant and principle component analyses highlighted Diet 1's distinctiveness, with histidine, diet digestibility, consumption, vg gene expression levels, and isoleucine identified as key factors. Arginine showed significant correlations with a wide range of parameters, including the number of honey bees, brood area, and consumption, with Diet 1 exhibiting higher levels. Diet 1, containing apple juice, soytide, and Chlorella as additive components, outperformed the other diets, suggesting an enhanced formulation for pollen substitute diets. These findings hold promise for the development of more effective diets, potentially contributing to honey bee health.

Biotic and abiotic stresses on honeybee health

Honeybees are the most critical pollinators providing key ecosystem services that underpin crop production and sustainable agriculture. Amidst a backdrop of rapid global change, this eusocial insect encounters a succession of stressors during nesting, foraging, and pollination. Ectoparasitic mites, together with vectored viruses, have been recognized as central biotic threats to honeybee health, while the spread of invasive giant hornets and small hive beetles also increasingly threatens colonies worldwide. Cocktails of agrochemicals, including acaricides used for mite treatment, and other pollutants of the environment have been widely documented to affect bee health in various ways. Additionally, expanding urbanization, climate change, and agricultural intensification often result in the destruction or fragmentation of flower-rich bee habitats. The anthropogenic pressures exerted by beekeeping management practices affect the natural selection and evolution of honeybees, and colony translocations facilitate alien species invasion and disease transmission. In this review, the multiple biotic and abiotic threats and their interactions that potentially undermine bee colony health are discussed, while taking into consideration the sensitivity, large foraging area, dense network among related nestmates, and social behaviors of honeybees.

Decline in Honeybees and Its Consequences for Beekeepers and Crop Pollination in Western Nepal

In understudied regions of the world, beekeeper records can provide valuable insights into changes in pollinator population trends. We conducted a questionnaire survey of 116 beekeepers in a mountainous area of Western Nepal, where the native honeybee Apis cerana cerana is kept as a managed bee. We complemented the survey with field data on insect-crop visitation, a household income survey, and an interview with a local lead beekeeper. In total, 76% of beekeepers reported declines in honeybees, while 86% and 78% reported declines in honey yield and number of beehives, respectively. Honey yield per hive fell by 50% between 2012 and 2022, whilst the number of occupied hives decreased by 44%. Beekeepers ranked climate change and declining flower abundance as the most important drivers of the decline. This raises concern for the future food and economic security of this region, where honey sales contribute to 16% of total household income, and where Apis cerana cerana plays a major role in crop pollination, contributing more than 50% of all flower visits to apple, cucumber, and pumpkin. To mitigate further declines, we promote native habitat and wildflower preservation, and using well-insulated log hives to buffer bees against the increasingly extreme temperature fluctuations.

A Bee Trp-Arg Dense Peptide with Antiproliferation Efficacy against the Prostate Cancer Cell Line DU145

Prostate cancer accounts for 14% of male cancer-related fatalities in the UK. Given the challenges associated with hormone-based therapies in the context of androgen-independent prostate cancer, there is an imperative need for research into anticancer drugs. N0821, a peptide belonging to the Trp-Arg dense region and derived from the homologous region of various bee species, shows substantial potential for an anticancer effect. Both MTT assays and 3D spheroid assays were conducted to substantiate its antiproliferation potential and strongly indicated the antiproliferation effect of N0820 (WWWWRWWRKI) and N0821 (YWWWWRWWRKI). Notably, the mechanism underlying this effect is related to the downregulation of CCNA2 and the upregulation of CCNE1. Cell cycle arrest results from the reduction of CCNA2 in the S/G2 phase, leading to the accumulation of CCNE1. Our peptides were predicted to make an α-helix structure. This can act as an ion channel in the cell membrane. Therefore, we analyzed genes implicated in the influx of calcium ions into the mitochondria. Trp-Arg dense-region peptides are known for their antibacterial properties in targeting cell membranes, making the development of resistance less likely. Hence, further research in this area is essential and promising.

Pan-genome analysis highlights the role of structural variation in the evolution and environmental adaptation of Asian honeybees

The Asian honeybee, Apis cerana, is an ecologically and economically important pollinator. Mapping its genetic variation is key to understanding population-level health, histories and potential capacities to respond to environmental changes. However, most efforts to date were focused on single nucleotide polymorphisms (SNPs) based on a single reference genome, thereby ignoring larger scale genomic variation. We employed long-read sequencing technologies to generate a chromosome-scale reference genome for the ancestral group of A. cerana. Integrating this with 525 resequencing data sets, we constructed the first pan-genome of A. cerana, encompassing almost the entire gene content. We found that 31.32% of genes in the pan-genome were variably present across populations, providing a broad gene pool for environmental adaptation. We identified and characterized structural variations (SVs) and found that they were not closely linked with SNP distributions; however, the formation of SVs was closely associated with transposable elements. Furthermore, phylogenetic analysis using SVs revealed a novel A. cerana ecological group not recoverable from the SNP data. Performing environmental association analysis identified a total of 44 SVs likely to be associated with environmental adaptation. Verification and analysis of one of these, a 330 bp deletion in the Atpalpha gene, indicated that this SV may promote the cold adaptation of A. cerana by altering gene expression. Taken together, our study demonstrates the feasibility and utility of applying pan-genome approaches to map and explore genetic feature variations of honeybee populations, and in particular to examine the role of SVs in the evolution and environmental adaptation of A. cerana.

Effects of the insecticide imidacloprid on the post-embryonic development of the honey bee Apis mellifera (Hymenoptera: Apidae)

The widespread use of pesticides in agriculture has been linked to declines in bee populations worldwide. Imidacloprid is a widely used systemic insecticide that can be found in the pollen and nectar of plants and has the potential to negatively impact the development of bee larvae. We investigated the effects of oral exposure to a realistic field concentration (20.5 ng g-1) of imidacloprid on the midgut and fat body of Apis mellifera worker larvae. Our results showed that larvae exposed to imidacloprid exhibited changes in the midgut epithelium, including disorganization of the brush border, nuclear chromatin condensation, cytoplasm vacuolization, and release of cell fragments indication cell death. Additionally, histochemical analysis revealed that the midgut brush border glycocalyx was disorganized in exposed larvae. The fat body cells of imidacloprid-exposed larvae had a decrease in the size of lipid droplets from 50 to 8 μm and increase of 100 % of protein content, suggesting possible responses to the stress caused by the insecticide. However, the expression of de cdc20 gene, which plays a role in cell proliferation, was not affected in the midgut and fat body of treated larvae. These results suggest that imidacloprid negatively affects non-target organs during the larval development of A. mellifera potentially impacting this important pollinator species.

Genetic differentiation and local adaptation of the Japanese honeybee, Apis cerana japonica

We examine the population genetic structure and divergence among the regional populations of the Japanese honeybee, Apis cerana japonica, by re-sequencing the genomes of 105 individuals from the three main Japanese islands with diverse climates. The genetic structure results indicated that these individuals are distinct from the mainland Chinese A. cerana samples. Furthermore, population structure analyses have identified three genetically distinct geographic regions in Japan: Northern (Tohoku-Kanto-Chubu districts), Central (Chugoku district), and Southern (Kyushu district). In some districts, "possible non-native" individuals, likely introduced from other regions in recent years, were discovered. Then, genome-wide scans were conducted to detect candidate genes for adaptation by two different approaches. We performed a population branch statistics (PBS) analysis to identify candidate genes for population-specific divergence. A latent factor mixed model (LFMM) was used to identify genes associated with climatic variables along a geographic gradient. The PBSmax analysis identified 25 candidate genes for population-specific divergence whereas the LFMM analysis identified 73 candidate genes for adaptation to climatic variables along a geographic gradient. However, no common genes were identified by both methods.

Acute and chronic effects of sublethal neonicotinoid thiacloprid to Asian honey bee (Apis cerana cerana)

Pesticide pollution is one of the most important factors for global bee declines. Despite many studies have revealed that the most important Chinese indigenous species，Apis cerana, is presenting a high risk on exposure to neonicotinoids, the toxicology information on Apis cerana remain limited. This study was aimed to determine the acute and chronic toxic effects of thiacloprid (IUPAC name: {(2Z)-3-[(6-Chloro-3-pyridinyl)methyl]-1,3-thiazolidin-2-ylidene}cyanamide) on behavioral and physiological performance as well as genome-wide transcriptome in A. cerana. We found the 1/5 LC50 of thiacloprid significantly impaired learning and memory abilities after both acute and chronic exposure, nevertheless, has no effects on the sucrose responsiveness and phototaxis climbing ability of A. cerana. Moreover, activities of detoxification enzyme P450 monooxygenases and CarE were increased by short-term exposure to thiacloprid, while prolonged exposure caused suppression of CarE activity. Neither acute nor chronic exposure to thiacloprid altered honey bee AChE activities. To further study the potential defense molecular mechanisms in Asian honey bee under pesticide stress, we analyzed the transcriptomes of honeybees in response to thiacloprid stress. The transcriptomic profiles revealed consistent upregulation of immune- and stress-related genes by both acute or chronic treatments. Our results suggest that the chronic exposure to thiacloprid produced greater toxic effects than a single administration to A. cerana. Altogether, our study deepens the understanding of the toxicological characteristic of A. cerana against thiacloprid, and could be used to further investigate the complex molecular mechanisms in Asian honey bee under pesticide stress.

Survey Results of Honey Bee Colony Losses in Winter in China (2009-2021)

There is growing concern that massive loss of honey bees can cause serious negative effects on biodiversity and ecosystems. Surveys of colony losses have been performed worldwide to monitor the dynamic changes and health status of honey bee colonies. Here, we present the results of surveys regarding winter colony losses from 21 provinces in China from 2009 to 2021, with a total of 1,744,324 colonies managed by 13,704 beekeepers. The total colony losses were low (9.84%; 95% Confidence Interval (CI): 9.60-10.08%) but varied among years, provinces, and scales of apiaries. As little is known about the overwintering mortality of Apis cerana, in this study, we surveyed and compared the loss rates between Apis mellifera and A. cerana in China. We found colonies of A. mellifera suffered significantly lower losses than A. cerana in China. Larger apiaries resulted in higher losses in A. mellifera, whereas the opposite was observed in A. cerana. Furthermore, we used generalized linear mixed-effects models (GLMMs) to evaluate the effects of potential risk factors on winter colony losses and found that the operation size, species, migration, migration×species interaction, and queen problems were significantly related to the loss rates. New queens can increase their colony overwintering survival. Migratory beekeepers and large operations reported lower loss rates.

Preliminary Survey of Pathogens in the Asian Honey Bee (Apis cerana) in Thailand

Widespread parasites, along with emerging threats, globalization, and climate change, have greatly affected honey bees' health, leading to colony losses worldwide. In this study, we investigated the detection of biotic stressors (i.e., viruses, microsporidian, bacteria, and fungi) in Apis cerana by surveying the colonies across different regions of Thailand (Chiang Mai in the north, Nong Khai and Khon Kaen in the northeast, and Chumphon and Surat Thani in the south, in addition to the Samui and Pha-ngan islands). In this study, we detected ABPV, BQCV, LSV, and Nosema ceranae in A. cerana samples through RT-PCR. ABPV was only detected from the samples of Chiang Mai, whereas we found BQCV only in those from Chumphon. LSV was detected only in the samples from the Samui and Pha-ngan islands, where historically no managed bees are known. Nosema ceranae was found in all of the regions except for Nong Khai and Khon Kaen in northeastern Thailand. Paenibacillus larvae and Ascosphaera apis were not detected in any of the A. cerana samples in this survey. The phylogenetic tree analysis of the pathogens provided insights into the pathogens' movements and their distribution ranges across different landscapes, indicating the flow of pathogens among the honey bees. Here, we describe the presence of emerging pathogens in the Asian honey bee as a valuable step in our understanding of these pathogens in terms of the decline in eastern honey bee populations.

Characterisation of the heat shock protein Tid and its involvement in stress response regulation in Apis cerana

Objective: The impact of various environmental stresses on native Apis cerana cerana fitness has attracted intense attention in China. However, the defence responses of A. cerana cerana to different stressors are poorly understood. Here, we aimed to elucidate the regulatory mechanism mediated by the tumorous imaginal discs (Tid) protein of A. cerana cerana (AccTid) in response to stressors. Methods: We used some bioinformatics softwares to analyse the characterisation of Tid. Then, qRT-PCR, RNA interference and heat resistance detection assays were used to explore the function of Tid in stress response in A. cerana cerana. Results: AccTid is a homologous gene of human Tid1 and Drosophila Tid56, contains a conserved J domain and belongs to the heat shock protein DnaJA subfamily. The level of AccTid induced expression was increased under temperature increases from 40°C to 43°C and 46°C, and AccTid knockdown decreased the heat resistance of A. cerana cerana, indicating that the upregulation of AccTid plays an important role when A. cerana cerana is exposed to heat stress. Interestingly, contrary to the results of heat stress treatment, the transcriptional level of AccTid was inhibited by cold, H₂O₂ and some agrochemical stresses and showed no significant change under ultraviolet ray and sodium arsenite stress. These results suggested that the requirement of A. cerana cerana for Tid differs markedly under different stress conditions. In addition, knockdown of AccTid increased the mRNA levels of some Hsps and antioxidant genes. The upregulation of these Hsps and antioxidant genes may be a functional complement of AccTid knockdown. Conclusion: AccTid plays a crucial role in A. cerana cerana stress responses and may mediate oxidative damage caused by various stresses. Our findings will offer fundamental knowledge for further investigations of the defence mechanism of A. cerana cerana against environmental stresses.

Conservation Genomic Analysis of the Asian Honeybee in China Reveals Climate Factors Underlying Its Population Decline

The Asian honeybee, Apis cerana, is one of the most important native pollinators in Asia. Asian honeybees were believed to be under significant decline in China based on a report in 2005. On the contrary, a recent survey revealed that Asian honeybee populations in China are stable and even slightly increased in some regions. Therefore, the declining status of A. cerana populations in China is still unclear. Taking advantage of the abundant, publicly available genomic data for Asian honeybees in China, we employed conservation genomics methods to understand if Asian honeybee populations in China are declining and what the underlying climate factors are. We reconstructed the changes of effective population size (Ne) within the recent past for 6 population groups of Asian honeybees and found out that only one of them (population in Bomi, Tibet) showed a consistently declining Ne from the last 100 generations to 25 generations. Selective sweep analysis suggests that genes related to the tolerance of low temperatures and strong ultraviolet radiation are under selection in the declining population, indicating that these two climate factors most likely underlie the decline of BM populations during the recent past. Our study provides insights into the dynamic changes of Asian honeybee populations in China and identifies climate factors that underlie its population decline, which is valuable for the conservation of this important pollinator.

Honey DNA metabarcoding revealed foraging resource partitioning between Korean native and introduced honey bees (Hymenoptera: Apidae)

Honey DNA metabarcoding provides information of floral sources of honey and foraging plant preferences of honey bees. We evaluated the floral composition of honey from two different species of honey bees, Apis cerana honey (ACH) and A. mellifera honey (AMH) in a mixed apiary located in a semi-forest environment to understand the floral preference and level of interspecific competition on floral resource. Three honey samples were collected from different hives of each species in mid-August. In total, 56 plant taxa were identified across the honey samples and among them, 38 taxonomic units were found in ACH compared with a total of 33 in AMH. The number of major plants (> 1% of reads) in honey samples was 9 and 11 in ACH and AMH respectively indicating the higher diversity of plant taxa in AMH. 23 taxonomic units were found exclusively in ACH, 18 taxonomic units were found only in AMH and 15 taxonomic units were shared between ACH and AMH indicating that 73% of the taxonomic units were present only in honey originated from one of the honeybee species. Qualitative and quantitative analyses of the shared major plants revealed the division of floral resource between these co-existing honey bee species pointing to a low level of interspecific competition between these two important pollinators.

Periodically taken photographs reveal the effect of pollinator insects on seed set in lotus flowers

Understanding of pollination systems is an important topic for evolutionary ecology, food production, and biodiversity conservation. However, it is difficult to grasp the whole picture of an individual system, because the activity of pollinators fluctuates depending on the flowering period and time of day. In order to reveal effective pollinator taxa and timing of visitation to the reproductive success of plants under the complex biological interactions and fluctuating abiotic factors, we developed an automatic system to take photographs at 5-s intervals to get near-complete flower visitation by pollinators during the entire flowering period of selected flowers of Nelumbo nucifera and track the reproductive success of the same flowers until fruiting. Bee visits during the early morning hours of 05:00-07:59 on the second day of flowering under optimal temperatures with no rainfall or strong winds contributed strongly to seed set, with possible indirect negative effects by predators of the pollinators. Our results indicate the availability of periodic and consecutive photography system in clarifying the plant-pollinator interaction and its consequence to reproductive success of the plant. Further development is required to build a monitoring system to collect higher-resolution time-lapse images and automatically identify visiting insect species in the natural environment.

Epidemiology and risk factors of self-reported systemic allergic reactions to a Hymenoptera venom in beekeepers worldwide: a protocol for a systematic review of observational studies

INTRODUCTION

Systemic allergic reaction (SAR) to a Hymenoptera venom is a potentially life-threatening disorder. The rate of SAR between beekeepers in comparison with a healthy individual is different. The risk for an SAR is particularly high in beekeepers due to their persistent or seasonal exposure to the stinging Hymenoptera. We aim to provide a critical appraisal and a synthesis of evidence-based data from epidemiological observational studies, focusing on SARs to a Hymenoptera venom and the associated risk factors for SARs in beekeepers worldwide.

METHODS AND ANALYSIS

Searching will include seven electronic databases for published studies without language restrictions, from inception up to 3 August 2021, and it will be rerun for all electronic databases prior publication. Only epidemiological observational studies in beekeepers will be included. The risk of bias in the included studies will be appraised by using the Joanna Briggs Institute Critical Appraisal Checklist for Analytical Cross-Sectional Studies and the Newcastle-Ottawa Scale, adapted for cross-sectional studies. For the certainty of evidence, the Grading of Recommendations Assessment, Development and Evaluation approach will be used. Qualitative synthesis will be presented in a tabulated format with the selected characteristics across primary studies and the main outcome of interest. A meta-analysis is planned to be performed if there will be a sufficient number of homogeneous studies with complete data. The Preferred Reporting Items for Systematic Review and Meta-Analysis Protocols 2015 statement will guide the reporting of this systematic literature review.

ETHICS AND DISSEMINATION

No ethics approval is needed to conduct the systematic literature review since it will be solely based on the published literature. Findings will be disseminated through the relevant conferences, peer-review and open-access journals.PROSPERO registration numberCRD42021260922.

Full length transcriptomes analysis of cold-resistance of Apis cerana in Changbai Mountain during overwintering period

Apis cerana in Changbai Mountain is an ecological type of Apis cerana, which is an excellent breeding material with cold-resistant developed by long-term natural selection under the ecological conditions. However, the physiological and molecular mechanisms of Changbai Mountain population under cold stress are still unclear. In this study, the Nanopore sequencing was carried out for the transcriptome of Apis cerana in Changbai Mountain in the coldest period of overwintering, which will provide a reference to the cold-resistant mechanism. We determined 5,941 complete ORF sequences, 1,193 lncRNAs, 619 TFs, 10,866 SSRs and functional annotations of 11,599 new transcripts. Our results showed that the myosin family and the C2H2 zinc finger protein transcription factor family possibly have significant impacts on the response mechanism of cold stress during overwintering. In addition, the cold environment alters genes expression profiles in honeybees via different AS and APA mechanisms. These altered genes in Hippo, Foxo, and MARK pathways help them counter the stress of cold in overwinter period. Our results might provide clues about the response of eastern honeybees to extreme cold, and reflect the possible genetic basis of physiological changes.

Differential Expression Characterisation of the Heat Shock Proteins DnaJB6, DnaJshv, DnaJB13, and DnaJB14 in Apis cerana cerana Under Various Stress Conditions

As key pollinators, bees are frequently exposed to multiple environmental stresses and have developed crucial mechanisms by which they adapt to these stressors. However, the molecular bases mediated at the gene level remain to be discovered. Here, we found four heat shock protein DnaJB subfamily genes, DnaJB6, DnaJshv, DnaJB13, and DnaJB14, from Apis cerana cerana, that all have J domains in their protein sequences. The expression levels of DnaJB6 and DnaJshv were upregulated by different degrees of heat stress, and the transcript level of DnaJB14 was gradually upregulated as the degree of heat stress increased, while the mRNA level of DnaJB13 was downregulated at multiple time points during heat stress treatment. The mRNA levels of all four DnaJBs were upregulated by cold and UV stress. In addition, the expression levels of DnaJB6, DnaJshv and DnaJB13 were reduced under abamectin, imidacloprid, cypermethrin, bifenthrin, spirodiclofen, and methomyl stresses. The transcript level of DnaJB14 was decreased by imidacloprid, cypermethrin, spirodiclofen, and methomyl exposure but increased by abamectin and bifenthrin exposure. These results indicate that the demand of A. cerana cerana for these four DnaJBs differs under various stress conditions. To further explore the role of DnaJBs in the stress response, we successfully silenced DnaJshv and DnaJB14. The content of protein carbonyl was increased, while the content of VC, the enzymatic activities of CAT, GST, and SOD, the mRNA levels of many antioxidant-related genes, and the total antioxidant capacity were reduced after knockdown of DnaJshv and DnaJB14 in A. cerana cerana. These results indicate that silencing DnaJshv and DnaJB14 increases oxidative damage and decreases the antioxidant ability of A. cerana cerana. Taken together, our results demonstrate that DnaJB6, DnaJshv, DnaJB13, and DnaJB14 are differentially expressed under stress conditions and play crucial roles in response to various stressors, possibly through the antioxidant signalling pathway. These findings will be conducive to understanding the molecular basis of bee responses to environmental stresses and are beneficial for improving bee protection.

Geometric morphology and population genomics provide insights into the adaptive evolution of Apis cerana in Changbai Mountain

Background

Exploration of adaptive evolution of organisms in response to environmental change can help to understand the evolutionary history of species and the underlying mechanisms of adaptation to local environments, thus guiding future conservation programmes. Before the introduction of Apis mellifera in China, eastern honey bees (Apis cerana) were the only species used for beekeeping in this region. In the mountains of Changbai, populations of A. cerana are considered a distinct ecotype of the species which formed through the distinct selective pressures in this area over time.

Result

We performed a measure of 300 wing specimens of eastern honey bees and obtained the geometric morphological variation in the wing of A. cerana in Changbai Mountain. A total of 3,859,573 high-quality SNP loci were yielded via the whole-genome resequencing of 130 individuals in 5 geographic regions.

Conclusion

Corresponding geometric morphology and population genomics confirmed the particularity of the A. cerana in Changbai Mountain. Genetic differentiation at the subspecies level exists between populations in Changbai Mountain and remaining geographic regions, and a significant reduction in the effective population size and an excessive degree of inbreeding may be responsible for a substantial loss of population genetic diversity. Candidate genes potentially associated with cold environmental adaptations in populations under natural selection were identified, which may represent local adaptations in populations. Our study provided insights into the evolutionary history and adaptation of A. cerana in Changbai Mountain, as well as its conservation.

Development of a Kit for Rapid Immunochromatographic Detection of Sacbrood Virus Infecting Apis cerana (AcSBV) Based on Polyclonal and Monoclonal Antibodies Raised against Recombinant VP1 and VP2 Expressed in Escherichia coli

A Korean isolate of the sacbrood virus infecting Apis cerana (AcSBV-Kor) is the most destructive honeybee virus, causing serious economic damage losses in Korean apiculture. To address this, here, we attempted to develop an assay for the rapid detection of AcSBV-Kor based on immunochromatographic detection of constituent viral proteins. Genes encoding VP1 and VP2 proteins of AcSBV-Kor were cloned into an expression vector (pET-28a) and expressed in Escherichia coli BL21(DE3). During purification, recombinant VP1 (rVP1) and VP2 (rVP2) proteins were found in the insoluble fraction, with a molecular size of 26.7 and 24.9 kDa, respectively. BALB/c mice immunized with the purified rVP1 and rVP2 produced polyclonal antibodies (pAbs) such as pAb-rVP1 and pAb-rVP2. Western blot analysis showed that pAb-rVP1 strongly reacted with the homologous rVP1 but weakly reacted with heterologous rVP2. However, pAb-rVP2 strongly reacted not only with the homologous rVP2 but also with the heterologous rVP1. Spleen cells of the immunized mice fused with SP2/0-Ag14 myeloma cells produced monoclonal antibodies (mAbs) such as mAb-rVP1-1 and mAb-rVP2-13. Western blot analysis indicated that pAb-rVP1, pAb-rVP2, mAb-rVP1-1, and mAb-rVP2-13 reacted with AcSBV-infected honeybees and larvae as well as the corresponding recombinant proteins. These antibodies were then used in the development of a rapid immunochromatography (IC) strip assay kit with colloidal gold coupled to pAb-rVP1 and pAb-rVP2 at the conjugate pad and mAb-rVP1-1 and mAb-rVP2-13 at the test line. One antibody pair, pAb-rVP1/mAb-VP1-1, showed positive reactivity as low as 1.38 × 10³ copies, while the other pair, pAb-rVP2/mAb-VP2-13, showed positive reactivity as low as 1.38 × 10⁴ copies. Therefore, the antibody pair pAb-rVP1/mAb-VP1-1 was selected as a final candidate for validation. To validate the detection of AcSBV, the IC strip tests were conducted with 50 positive and 50 negative samples and compared with real-time PCR tests. The results confirm that the developed IC assay is a sufficiently sensitive and specific detection method for user-friendly and rapid detection of AcSBV.

The Current State of the Protected Apis mellifera mellifera Population in Russia: Hybridization and Nosematosis

The Southern Urals of Russia are the habitat of one of the surviving populations of the dark forest bee-the Burzyan population of Apis mellifera mellifera. In this study, we present the results of the subspecies identification of bee colonies in the Altyn-Solok Nature Reserve in the Southern Ural Mountains using the intergenic mtDNA COI-COII locus and the assessment of the prevalence of nosematosis. Analysis of the mtDNA COI-COII intergenic locus in the studied sample showed that 30.4% of the colonies belong to the lineage C. The PCR diagnostics of nosematosis in 92 colonies selected from different sectors of the Altyn-Solok Nature Reserve showed that about half of the analyzed colonies were infected with Nosema apis. Nosema ceranae was found in eight colonies. Both of these factors can lead to the extinction of this population of the dark forest bee.

Genetic divergence and functional convergence of gut bacteria between the Eastern honey bee Apis cerana and the Western honey bee Apis mellifera

•

The inter-species diversity of A. cerana and A. mellifera core gut bacteria was revealed.

•

Core bacterial species of A. cerana and A. mellifera are distinctive in function.

•

Functional profile of overall gut community of A. cerana and A. mellifera are similar.

•

Metabolome showed that A. cerana and A. mellifera gut bacteria have similar metabolic capability.

•

A. cerana and A. mellifera core gut bacteria have no strict host specificity.

Introduction

The functional relevance of intra-species diversity in natural microbial communities remains largely unexplored. The guts of two closely related honey bee species, Apis cerana and A. mellifera, are colonised by a similar set of core bacterial species composed of host-specific strains, thereby providing a good model for an intra-species diversity study.

Objectives

We aim to assess the functional relevance of intra-species diversity of A. cerana and A. mellifera gut microbiota.

Methods

Honey bee workers were collected from four regions of China. Their gut microbiomes were investigated by shotgun metagenomic sequencing, and the bacterial compositions were compared at the species level. A cross-species colonisation assay was conducted, with the gut metabolomes being characterised by LC-MS/MS.

Results

Comparative analysis showed that the strain composition of the core bacterial species was host-specific. These core bacterial species presented distinctive functional profiles between the hosts. However, the overall functional profiles of the A. cerana and A. mellifera gut microbiomes were similar; this was further supported by the consistency of the honey bees’ gut metabolome, as the gut microbiota of different honey bee species showed rather similar metabolic profiles in the cross-species colonisation assay. Moreover, this experiment also demonstrated that the gut microbiota of A. cerana and A. mellifera could cross colonise between the two honey bee species.

Conclusion

Our findings revealed functional differences in most core gut bacteria between the guts of A. cerana and A. mellifera, which may be associated with their inter-species diversity. However, the functional profiles of the overall gut microbiomes between the two honey bee species converge, probably as a result of the overlapping ecological niches of the two species. Our findings provide critical insights into the evolution and functional roles of the mutualistic microbiota of honey bees and reveal that functional redundancy could stabilise the gene content diversity at the strain-level within the gut community.

Using Citizen Science to Scout Honey Bee Colonies That Naturally Survive Varroa destructor Infestations

Simple Summary

Citizen Science is a valuable resource that can substantially contribute to the conservation of biodiversity. However, its use in honey bee research has remained minimal. The Survivors Task Force of the COLOSS association created and promoted an online surveying tool with the aim of identifying potential cases of Western honey bee, Apis mellifera, populations that are surviving infestations with ectoparasitic mites Varroa destructor without control measures by beekeepers. The reports suggest that there could be twice as many naturally surviving colonies worldwide than are currently known. The survey also shows that citizens can be readily engaged through social media, personal networks, and promotional campaigns to gather valuable and previously inaccessible data. These reports of surviving honey bee colonies will now be validated through the new initiative Honey Bee Watch, a global and multi-year Citizen Science project founded to connect citizens, beekeepers, and scientists. This will enable to increase scientific knowledge, mitigate honey bee colony losses, and develop education and conservation campaigns.

Abstract

Citizen Science contributes significantly to the conservation of biodiversity, but its application to honey bee research has remained minimal. Even though certain European honey bee (Apis mellifera) populations are known to naturally survive Varroa destructor infestations, it is unclear how widespread or common such populations are. Such colonies are highly valuable for investigating the mechanisms enabling colony survival, as well as for tracking the conservation status of free-living honey bees. Here, we use targeted Citizen Science to identify potentially new cases of managed or free-living A. mellifera populations that survive V. destructor without mite control strategies. In 2018, a survey containing 20 questions was developed, translated into 13 languages, and promoted at beekeeping conferences and online. After three years, 305 reports were collected from 28 countries: 241 from managed colonies and 64 from free-living colonies. The collected data suggest that there could be twice as many naturally surviving colonies worldwide than are currently known. Further, online and personal promotion seem to be key for successful recruitment of participants. Although the survivor status of these colonies still needs to be confirmed, the volume of reports and responses already illustrate how effectively Citizen Science can contribute to bee research by massively increasing generated data, broadening opportunities for comparative research, and fostering collaboration between scientists, beekeepers, and citizens. The success of this survey spurred the development of a more advanced Citizen Science platform, Honey Bee Watch, that will enable a more accurate reporting, confirmation, and monitoring of surviving colonies, and strengthen the ties between science, stakeholders, and citizens to foster the protection of both free-living and managed honey bees.

Viral impacts on honey bee populations: A review

Honey bee is vital for pollination and ecological services, boosting crops productivity in terms of quality and quantity and production of colony products: wax, royal jelly, bee venom, honey, pollen and propolis. Honey bees are most important plant pollinators and almost one third of diet depends on bee's pollination, worth billions of dollars. Hence the role that honey bees have in environment and their economic importance in food production, their health is of dominant significance. Honey bees can be infected by various pathogens like: viruses, bacteria, fungi, or infested by parasitic mites. At least more than 20 viruses have been identified to infect honey bees worldwide, generally from Dicistroviridae as well as Iflaviridae families, like ABPV (Acute Bee Paralysis Virus), BQCV (Black Queen Cell Virus), KBV (Kashmir Bee Virus), SBV (Sacbrood Virus), CBPV (Chronic bee paralysis virus), SBPV (Slow Bee Paralysis Virus) along with IAPV (Israeli acute paralysis virus), and DWV (Deformed Wing Virus) are prominent and cause infections harmful for honey bee colonies health. This issue about honey bee viruses demonstrates remarkably how diverse this field is, and considerable work has to be done to get a comprehensive interpretation of the bee virology.

Gene reuse facilitates rapid radiation and independent adaptation to diverse habitats in the Asian honeybee

Animals with recent shared ancestry frequently adapt in parallel to new but similar habitats, a process often underlined by repeated selection of the same genes. Yet, in contrast, few examples have demonstrated the significance of gene reuse in colonization of multiple disparate habitats. By analyzing 343 genomes of the widespread Asian honeybee, Apis cerana, we showed that multiple peripheral subspecies radiated from a central ancestral population and adapted independently to diverse habitats. We found strong evidence of gene reuse in the Leucokinin receptor (Lkr), which was repeatedly selected in almost all peripheral subspecies. Differential expression and RNA interference knockdown revealed the role of Lkr in influencing foraging labor division, suggesting that Lkr facilitates collective tendency for pollen/nectar collection as an adaptation to floral changes. Our results suggest that honeybees may accommodate diverse floral shifts during rapid radiation through fine-tuning individual foraging tendency, a seemingly complex process accomplished by gene reuse.

Functional and transcriptomic analyses of the NF-Y family provide insights into the defense mechanisms of honeybees under adverse circumstances

As predominant pollinators, honeybees are important for crop production and terrestrial ecosystems. Recently, various environmental stresses have led to large declines in honeybee populations in many regions. The ability of honeybees to respond to these stresses is critical for their survival. However, the details of the stress defense mechanisms of honeybees have remained elusive. Here, we found that the Nuclear Factor Y (NF-Y) family (containing NF-YA, NF-YB, and NF-YC) is a novel stress mediator family that regulates honeybee environmental stress resistance. NF-YA localized in the nucleus, NF-YB accumulated in the cytoplasm, and NF-YC presented in both the nucleus and cytoplasm. NF-YC interacted with NF-YA and NF-YB in vitro and in vivo, and the nuclear import of NF-YB relied on its interaction with NF-YC. We further found that the expression of NF-Y was induced under multiple stress conditions. In addition, NF-Y regulated many stress responses and antioxidant genes at the transcriptome-wide level, and knockdown of NF-Y repressed the expression of stress-inducible genes, particularly LOC108003540 and LOC107994062, under adverse circumstances. Silencing NF-Y lowered honeybee stress resistance by reducing total antioxidant capacity and enhancing oxidative impairment. Collectively, these results indicate that NF-Y plays important roles in stress responses. Our study sheds light on the underlying defense mechanisms of honeybees under environmental stress.

Transcriptomic and metabolomic landscape of the molecular effects of glyphosate commercial formulation on Apis mellifera ligustica and Apis cerana cerana

Understanding the causes of the decline in bee population has attracted intensive attention worldwide. The indiscriminate use of agrochemicals is a persistent problem due to their physiological and behavioural damage to bees. Glyphosate and its commercial formulation stand out due to their wide use in agricultural areas and non-crop areas, such as parks, railroads, roadsides, industrial sites, and recreational and residential areas, but the mode of action of glyphosate on bees at the molecular level remains largely unelucidated. Here, we found that the numbers of differentially expressed genes and metabolites under glyphosate commercial formulation (GCF) stress were significantly higher in Apis cerana cerana than in Apis mellifera ligustica. Despite these differences, the number of differentially expressed transcripts increased following an increase in the GCF treatment time in both A. cerana cerana and A. mellifera ligustica. GCF exerted adverse impacts on the immune system, digestive system, nervous system, amino acid metabolism, carbohydrate metabolism, growth and development of both bee species by influencing their key genes and metabolites to some extent. The expression of many genes involved in immunity, agrochemical detoxification and resistance, such as antimicrobial peptides, cuticle proteins and cytochrome P450 families, was upregulated by GCF in both bee species. Collectively, our results indicate that both A. cerana cerana and A. mellifera ligustica strive to mitigate the pernicious effects caused by GCF by regulating detoxification and immune systems. Moreover, A. cerana cerana might be better able to withstand the toxic effects of GCF with lower fitness costs than A. mellifera ligustica. Our work will contribute to elucidating the deleterious physiological and behavioural impacts of GCF on bees.

Substantial Genetic Progress in the International Apis mellifera carnica Population Since the Implementation of Genetic Evaluation

Simple Summary

The Apis mellifera carnica subspecies of the honeybee is known for its gentleness and good honey yield. In the early 20th century, systematic breeding efforts began. Breeding progress was slow before the introduction of modern techniques of genetic evaluation in the mid 1990s. Here, the results of the official breeding value estimation in BeeBreed.eu are analyzed. From about 2000 onward, breeding progress accelerated. The result is a considerable gain in honey yield and desirable properties without increased inbreeding coefficients. The success of A. m. carnica breeding shows the potential of genetic evaluation.

Abstract

The Apis mellifera carnica subspecies of the honeybee has long been praised for its gentleness and good honey yield before systematic breeding efforts began in the early 20th century. However, before the introduction of modern techniques of genetic evaluation (best linear unbiased prediction, BLUP) and a computerized data management in the mid 1990s, genetic progress was slow. Here, the results of the official breeding value estimation in BeeBreed.eu are analyzed to characterize breeding progress and inbreeding. From about the year 2000 onward, the genetic progression accelerated and resulted in a considerable gain in honey yield and desirable properties without increased inbreeding coefficients. The prognostic quality of breeding values is demonstrated by a retrospective analysis. The success of A. m. carnica breeding shows the potential of BLUP-based breeding values and serves as an example for a large-scale breeding program.

The effects of diploid male production on honey bee colony evolution and survival

The order Hymenoptera includes most of the eusocial species on the planet. Correlated is the fact that many of the social species within the order are haplodiploid and use complementary sex determination (CSD) to determine the sex of offspring. CSD is the mechanism by why single sex alleles within an organism result in male development (haploid) and mismatched sex alleles develop into females (diploids). Related to this is the production of diploid males: fertilized eggs with matched sex alleles which develop as male instead of female. Honey bees are no exception to this, and as their numbers continue to suffer globally and their genetic diversity lowers, the effects of diploid male production (DMP) may pose an increased risk to the survival of bee colonies. In the present study, we develop a model for diploid male production in a honey bee colony and show that with ample resources, this phenomena has little effect on a colony's health, but there is a limit to the sustainability of a colony suffering from diploid male production. We use our model to show that there were likely no great evolutionary pressures against CSD and DMP in wild honey bees as its effects on colony health in the wild would have been negligible but increased environmental hazards such as pesticides and monoculture crops increase the effects of DMP on colony health.

Long-Term Evaluation of Breeding Scheme Alternatives for Endangered Honeybee Subspecies

Modern breeding structures are emerging for European honeybee populations. However, while genetic evaluations of honeybees are becoming increasingly well understood, little is known about how selection decisions shape the populations' genetic structures. We performed simulations evaluating 100 different selection schemes, defined by selection rates for dams and sires, in populations of 200, 500, or 1000 colonies per year and considering four different quantitative traits, reflecting different genetic parameters and numbers of influential loci. Focusing on sustainability, we evaluated genetic progress over 100 years and related it to inbreeding developments. While all populations allowed for sustainable breeding with generational inbreeding rates below 1% per generation, optimal selection rates differed and sustainable selection was harder to achieve in smaller populations and for stronger negative correlations of maternal and direct effects in the selection trait. In small populations, a third or a fourth of all candidate queens should be selected as dams, whereas this number declined to a sixth for larger population sizes. Furthermore, our simulations indicated that, particularly in small populations, as many sires as possible should be provided. We conclude that carefully applied breeding provides good prospects for currently endangered honeybee subspecies, since sustainable genetic progress improves their attractiveness to beekeepers.

Analyses of the function of DnaJ family proteins reveal an underlying regulatory mechanism of heat tolerance in honeybee

There is clear evidence of severe honeybee declines in recent years, and parallel declines of plant community and crop productivity that rely on them. Different stresses, including heat stress, are among the primary drivers of this decline. However, the mechanisms by which honeybees respond to heat stress are elusive. Though heat shock proteins (Hsps) play important roles in heat stress response, the function of DnaJs (a subfamily of Hsps) is unclear. Here, we aimed to determine the underlying regulatory mechanism of honeybees to heat stress mediated by DnaJs. We found that several DnaJ genes, including DnaJA1, DnaJB12 and DnaJC8, are key for honeybee heat tolerance. DnaJA1 and DnaJB12 are cytoplasmic proteins, and DnaJC8 is a nuclear protein. The expression of DnaJA1, DnaJB12 and DnaJC8 was induced at different levels under short-term and long-term heat stress. Phenotypic analysis indicated that DnaJA1, DnaJB12 and DnaJC8 knockdown attenuated honeybee heat resistance. In addition, DnaJA1 participated in the heat stress response by upregulating many heat-inducible genes at the transcriptome-wide level, especially LOC108002668 and LOC107995148. Importantly, the upregulation of LOC108002668 and LOC107995148 was significantly repressed under heat stress when DnaJA1 was knocked down. We also found that knockdown of DnaJA1, DnaJB12 and DnaJC8 decreased antioxidant defense ability and increased the degree of oxidative damage in the honeybee. Taken together, our results indicate that DnaJ genes play important roles under heat stress in the honeybee. Overexpression of DnaJ genes may protect honeybees from heat stress-induced injuries and increase their survival rate.

Population Genomics Provide Insights into the Evolution and Adaptation of the Eastern Honey Bee (Apis cerana)

The mechanisms by which organisms adapt to variable environments are a fundamental question in evolutionary biology and are important to protect important species in response to a changing climate. An interesting candidate to study this question is the honey bee Apis cerana, a keystone pollinator with a wide distribution throughout a large variety of climates, that exhibits rapid dispersal. Here, we resequenced the genome of 180 A. cerana individuals from 18 populations throughout China. Using a population genomics approach, we observed considerable genetic variation in A. cerana. Patterns of genetic differentiation indicate high divergence at the subspecies level, and physical barriers rather than distance are the driving force for population divergence. Estimations of divergence time suggested that the main branches diverged between 300 and 500 Ka. Analyses of the population history revealed a substantial influence of the Earth's climate on the effective population size of A. cerana, as increased population sizes were observed during warmer periods. Further analyses identified candidate genes under natural selection that are potentially related to honey bee cognition, temperature adaptation, and olfactory. Based on our results, A. cerana may have great potential in response to climate change. Our study provides fundamental knowledge of the evolution and adaptation of A. cerana.

Population genomics and morphological features underlying the adaptive evolution of the eastern honey bee (Apis cerana)

Background

The adaptation of organisms to changing environments is self-evident, with the adaptive evolution of organisms to environmental changes being a fundamental problem in evolutionary biology. Bees can pollinate in various environments and climates and play important roles in maintaining the ecological balance of the earth.

Results

We performed an analysis of 462 Apis cerana (A. cerana) specimens from 31 populations in 11 regions and obtained 39 representative morphological features. We selected 8 A. cerana samples from each population and performed 2b-RAD simplified genome sequencing. A total of 11,506 high-quality single nucleotide polymorphism (SNP) loci were obtained. For these SNPs, the minor allele frequency (MAF) was > 1%, the average number of unique labels for each sample was 49,055, and the average depth was 72.61x. The ratios of the unique labels of all samples were 64.27–86.33%.

Conclusions

Using 39 morphological characteristics as the data set, we proposed a method for the rapid classification of A. cerana. Using genomics to assess population structure and genetic diversity, we found that A. cerana has a large genetic difference at the ecotype level. A comparison of A. cerana in North China revealed that some physical obstacles, especially the overurbanization of the plains, have isolated the populations of this species. We identified several migration events in North China and Central China. By comparing the differences in the environmental changes in different regions, we found that A. cerana has strong potential for climate change and provides a theoretical basis for investigating and protecting A. cerana.

The genomic basis of adaptation to high-altitude habitats in the eastern honey bee (Apis cerana)

The eastern honey bee (Apis cerana) is of central importance for agriculture in Asia. It has adapted to a wide variety of environmental conditions across its native range in southern and eastern Asia, which includes high-altitude regions. eastern honey bees inhabiting mountains differ morphologically from neighbouring lowland populations and may also exhibit differences in physiology and behaviour. We compared the genomes of 60 eastern honey bees collected from high and low altitudes in Yunnan and Gansu provinces, China, to infer their evolutionary history and to identify candidate genes that may underlie adaptation to high altitude. Using a combination of F_ST -based statistics, long-range haplotype tests and population branch statistics, we identified several regions of the genome that appear to have been under positive selection. These candidate regions were strongly enriched for coding sequences and had high haplotype homozygosity and increased divergence specifically in highland bee populations, suggesting they have been subjected to recent selection in high-altitude habitats. Candidate loci in these genomic regions included genes related to reproduction and feeding behaviour in honey bees. Functional investigation of these candidate loci is necessary to fully understand the mechanisms of adaptation to high-altitude habitats in the eastern honey bee.