The authors' names are listed in alphabetical order with the exception of the corresponding author. All authors' contributions are equal. SummaryThe term "quality" in relation to queens and drones refers to certain quantitative physical and / or behavioural characters. It is generally believed that a high quality queen should have the following physical characteristics: high live weight; high number of ovarioles; large size of spermatheca; high number of spermatozoa in spermatheca; and be free from diseases and pests. It is, however, also known that the performance of a honey bee colony is the result of its queen's function as well as of that of the drones that mated with her. These two approaches are often considered together and give a general picture of the queen production technique and selection. Here we describe the most common and well known anatomical, physiological, behavioural and performance characters related to the queens, as measured in different European countries: the live weight of the virgin queen (Bulgaria); the live weight of the laying queen (Bulgaria, Italy); the diameter and volume of spermatheca (Bulgaria, Greece, Slovenia); the number of ovarioles (Greece, Italy, Slovenia); the weight of ovaries (Slovenia); the number of spermatozoa in spermatheca (Italy, Poland, Slovenia); the brood pattern (Bulgaria, Greece); the egg laying ability / fecundity (Bulgaria); the brood production (Croatia, Serbia); the colony population development (Croatia, Serbia, Slovakia); the honey production (Croatia, Denmark, Serbia, Slovakia); the hygienic behaviour (Croatia, Denmark, Serbia, Slovakia); the defence behaviour (Croatia); the calmness / sitting on the comb (Croatia, Denmark); and swarming (Croatia, Denmark). The data presented fit well with the findings of the same characters in the literature, and in general they support the argument for the term "quality characters". Especially for the weight of the queen, the number of ovarioles, the volume of the spermatheca and the number of spermatozoa, data per country proved its own accuracy by repetition through the years. We also report that when instrumentally inseminated queens are kept under mass production conditions (in 338 Hatjina et al.

Background The honey bee (Apis mellifera) is an ecologically and economically important species that provides pollination services to natural and agricultural systems. The biodiversity of the honey bee in parts of its native range is endangered by migratory beekeeping and commercial breeding. In consequence, some honey bee populations that are well adapted to the local environment are threatened with extinction. A crucial step for the protection of honey bee biodiversity is reliable differentiation between native and nonnative bees. One of the methods that can be used for this is the geometric morphometrics of wings. This method is fast, is low cost, and does not require expensive equipment. Therefore, it can be easily used by both scientists and beekeepers. However, wing geometric morphometrics is challenging due to the lack of reference data that can be reliably used for comparisons between different geographic regions. Findings Here, we provide an unprecedented collection of 26,481 honey bee wing images representing 1,725 samples from 13 European countries. The wing images are accompanied by the coordinates of 19 landmarks and the geographic coordinates of the sampling locations. We present an R script that describes the workflow for analyzing the data and identifying an unknown sample. We compared the data with available reference samples for lineage and found general agreement with them. Conclusions The extensive collection of wing images available on the Zenodo website can be used to identify the geographic origin of unknown samples and therefore assist in the monitoring and conservation of honey bee biodiversity in Europe.

A diverse supply of pollen is an important factor for honey bee health, but information about the pollen diversity available to colonies at the landscape scale is largely missing. In this COLOSS study, beekeeper citizen scientists sampled and analyzed the diversity of pollen collected by honey bee colonies. As a simple measure of diversity, beekeepers determined the number of colors found in pollen samples that were collected in a coordinated and standardized way. Altogether, 750 beekeepers from 28 different regions from 24 countries participated in the two-year study and collected and analyzed almost 18,000 pollen samples. Pollen samples contained approximately six different colors in total throughout the sampling period, of which four colors were abundant. We ran generalized linear mixed models to test for possible effects of diverse factors such as collection, i.e., whether a minimum amount of pollen was collected or not, and habitat type on the number of colors found in pollen samples. To identify habitat effects on pollen diversity, beekeepers’ descriptions of the surrounding landscape and CORINE land cover classes were investigated in two different models, which both showed that both the total number and the rare number of colors in pollen samples were positively affected by ‘urban’ habitats or ‘artificial surfaces’, respectively. This citizen science study underlines the importance of the habitat for pollen diversity for bees and suggests higher diversity in urban areas.

Apis cerana and Apis mellifera are important honey bee species in Asia. A. cerana populations are distributed from a cold, sharply continental climate in the north to a hot, subtropical climate in the south. Due to the Sacbrood virus, almost all A. cerana populations in Asia have declined significantly in recent decades and have recovered over the past five years. This could lead to a shift in the gene pool of local A. cerana populations that could affect their sustainability and adaptation. It was assumed that adaptation of honey bees could be observed by comparative analysis of the sequences of genes involved in development, labor division, and caste differentiation, such as the gene Vitellogenin VG. The VG gene nucleotide sequences were used to assess the genetic structure and signatures of adaptation of local populations of A. cerana from Korea, Russia, Japan, Nepal, and China. A. mellifera samples from India and Poland were used as the outgroup. The signatures of adaptive selection were found in the local population of A. cerana using VG gene sequence analysis based on Jukes–Cantor genetic distances, cluster analysis, dN/dS ratio evaluation, and Tajima’s D neutrality test. Based on analysis of the VG gene sequences, Apis cerana koreana subspecies in the Korean Peninsula were subdivided into three groups in accordance with their geographic localization from north to south. The VG gene sequences are acceptable tools to study the sustainability and adaptation of A. cerana populations.

Apis cerana and Apis mellifera, are very important honey species for agriculture in Asian countries. In recent decades, A. cerana populations have sharply declined in all Asian countries as a result of Sacbrood Virus infection and have now recovered to their original size. It can change the genetic structure of local populations of A. cerana. We used the nuclear gene Vitellogenin VG to assess the genetic structure of local populations of A. cerana and the signature of adaptive selection. We performed a population genetic analysis of the honey bees A. cerana from South Korea in comparison with A. cerana samples from Russia, Japan, Nepal, and China. The sequences of the gene VG of a closely related honey bee species, A. mellifera, from India and Poland were used as outgroup samples. A comparative analysis of northern and southern A. cerana populations was performed. The signatures of positive adaptive selection were found in the local population of A. cerana. We performed the Tajima's neutrality D test for A. cerana populations from different local populations based on the gene VG exon sequences. All A. cerana populations showed signs of population size expansion following the possible recent decline in population sizes. The local populations of A. c. koreana were subdivided according to their geographical distribution into southern, northern, and central Korean clusters. The gene VG exon sequences can be used as informative markers for monitoring the changes in genetic structure and adaptation to the environment processes in A. cerana populations.