Evolutionary history of the honey bee Apis mellifera inferred from mitochondrial DNA analysis

Evolutionary Origin and Genetic Diversity of the Pannonian Ecotype of Apis mellifera carnica Colonies in Hungary Based on Mitochondrial DNA and Microsatellite Markers

Honey bees are crucial to both the ecosystem and the economy. However, they are subject to different influences that can lead to a loss of genetic diversity. In this study, we used mitochondrial DNA information and nuclear microsatellite markers to compare worker individuals that strictly meet the morphological breed standard of the Pannonian bee in Hungary to those with morphological disorders (yellow color of the abdomen). Additionally, this study involves Carniolan colonies from two European countries and other bee subspecies as a reference group that might have crossed into the Pannonian bee. As for the mitochondrial DNA, the combined assessment of COI and 16S genes identified six haplotypes. Based on the tRNAleu-cox2 intergenic region (E2/H2), our samples belonged to the C evolutionary lineage. According to the microsatellite data, the level of inbreeding was low in all groups investigated, and only the genotypes of the Pannonian bee showed significant deviation from the Hardy-Weinberg equilibrium state. Cluster analysis and the Discriminant Analysis of Principal Components showed that bees that failed the morphological breed identification had started to diverge genetically from those meeting the breed standards, becoming more similar to the Carniolan bee. Our findings suggest that the genetic status of the Pannonian bees investigated in this study is satisfactory. However, in order to maintain an adequate level of diversity, periodic genetic monitoring of the colonies is necessary.

Genetic diversity and novel haplotypes of Apis mellifera jemenitica on the Arabian Peninsula: insights from mtDNA markers

The genetic diversity of Apis mellifera jemenitica populations collected from the Arabian Peninsula (Saudi Arabia, Yemen, and Oman), Jordan, and Ethiopia, was examined using three mtDNA markers: 1- Cytochrome b (Cyt b), 2- Cytochrome c oxidase I (COI) and 3- The intergenic region located between the cytochrome c oxidase I & II (COI-COII). DNA was extracted from 44 samples, amplified for each region using classic PCR, and the resulting amplicons were sequenced using Sanger technology at both ends. Sequences were verified and aligned, and Maximum-Likelihood phylogenetic analyses were conducted with reference sequences from other subspecies. The in silico DraI mtDNA COI-COII (DmCC) test was applied to the COI-COII sequences to identify evolutionary lineages and haplotypes. Moreover, COI-COII haplotype network analyses were conducted to assess the intra- and inter-genetic relationships between samples and references. Based on the Cyt b marker, most samples cluster within the African lineage (A) near lamarckii and syriaca (Sub-lineage Z) subspecies. Few samples from Ethiopia and Yemen were closely related to simensis and scutellata clades. The COI gene separated jemenitica samples (Bootstrap = 97) from subspecies of other lineages (C and O). The DmCC test revealed a P0Q2 structure in the intergenic region for all samples, with a distinct 18 bp deletion in the P0 element observed in two Ethiopian and one Yemeni samples, suggesting litorea or simensis origin. A total of 13 COI-COII haplotypes were identified, among which 8 haplotypes were novel: Saudi Arabia (1), Yemen (3), Oman (1), and Ethiopia (3), with a haplotype diversity (H) of 0.980. Furthermore, molecular-variance parsimony in COI-COII confirmed a distant genetic relationship between Ethiopian samples versus samples of the Arabian Peninsula. The haplotype network analysis suggests a higher intra-jemenitica diversity than previously understood with a syriaca ancestry to this clade. These findings offer crucial insights into the conservation of A. m. jemenitica and its role in preserving biodiversity in arid ecosystems. Additionally, the data enhance our understanding of the genetic diversity of A. m. jemenitica and its evolutionary connections with other neighboring African subspecies.

COI Metabarcoding as a Novel Approach for Assessing the Honey Bee Source of European Honey

Honey is a widely consumed food product frequently subjected to adulteration, with the mislabeling of its botanical or geographical origin being one of the most common practices. Determining the entomological origin of honey is particularly challenging but of high relevance for ensuring its authenticity, especially for products with protected designation of origin (PDO) labels. This study presents a novel DNA metabarcoding approach targeting a highly informative 406 bp fragment of the cytochrome c oxidase I (COI) gene to differentiate among the three major mitochondrial lineages (A, M, and C) of honey bees (Apis mellifera L.) native to Europe. The target region was selected based on the calculated fixation index (FST), which is frequently used in Population Genetics as a measure of differentiation between populations. The approach was validated with 11 honey samples of known entomological origin and applied to 44 commercial honeys from 13 countries. The approach demonstrated high sensitivity, accurately identifying the entomological origin of honey, including samples produced by honey bees of varying ancestries, which could not be resolved by previous methods based on real-time PCR coupled with high-resolution melting (PCR-HRM) analysis. The results demonstrate the effectiveness of COI metabarcoding in verifying honey authenticity and highlight the predominance of C-lineage honey bees in the production of commercial honeys from northwestern Europe. This finding suggests a limited presence of the native M-lineage ancestry, underscoring the need for conservation efforts.

Inferring Long-Term and Short-Term Determinants of Genetic Diversity in Honey Bees: Beekeeping Impact and Conservation Strategies

Bees are vital pollinators in natural and agricultural landscapes around the globe, playing a key role in maintaining flowering plant biodiversity and ensuring food security. Among the honey bee species, the Western honey bee (Apis mellifera) is particularly significant, not only for its extensive crop pollination services but also for producing economically valuable products such as honey. Here, we analyzed whole-genome sequence data from four Apis species to explore how honey bee evolution has shaped current diversity patterns. Using Approximate Bayesian Computation, we first reconstructed the demographic history of A. mellifera in Europe, finding support for postglacial secondary contacts, therefore predating human-mediated transfers linked to modern beekeeping. However, our analysis of recent demographic changes reveals significant bottlenecks due to beekeeping practices, which have notably affected genetic diversity. Black honey bee populations from conservatories, particularly those on islands, exhibit considerable genetic loss, highlighting the need to evaluate the long-term effectiveness of current conservation strategies. Additionally, we observed a high degree of conservation in the genomic landscapes of nucleotide diversity across the four species, despite a divergence gradient spanning over 15 million years, consistent with a long-term conservation of the recombination landscapes. Taken together, our results provide the most comprehensive assessment of diversity patterns in honey bees to date and offer insights into the optimal management of resources to ensure the long-term persistence of honey bees and their invaluable pollination services.

Deciphering the variation in cuticular hydrocarbon profiles of six European honey bee subspecies

The Western honey bee (Apis mellifera) subspecies exhibit local adaptive traits that evolved in response to the different environments that characterize their native distribution ranges. An important trait is the cuticular hydrocarbon (CHC) profile, which helps to prevent desiccation and mediate communication. We compared the CHC profiles of six European subspecies (A. m. mellifera, A. m. carnica, A. m. ligustica, A. m. macedonica, A. m. iberiensis, and A. m. ruttneri) and investigated potential factors shaping their composition. We did not find evidence of adaptation of the CHC profiles of the subspecies to the climatic conditions in their distribution range. Subspecies-specific differences in CHC composition might be explained by phylogenetic constraints or genetic drift. The CHC profiles of foragers were more subspecies-specific than those of nurse bees, while the latter showed more variation in their CHC profiles, likely due to the lower desiccation stress exerted by the controlled environment inside the hive. The strongest profile differences appeared between nurse bees and foragers among all subspecies, suggesting an adaptation to social task and a role in communication. Foragers also showed an increase in the relative amount of alkanes in their profiles compared to nurses, indicating adaptation to climatic conditions.

Distribution of honey bee mitochondrial DNA haplotypes in an Italian region where a legislative act is protecting the Apis mellifera ligustica subspecies

The conservation of the genetic integrity of Apis mellifera subspecies has emerged as an important objective. In 2019, the Emilia-Romagna region became the first Italian regional authority to issue a law specifically addressing the protection of the native Apis mellifera ligustica subspecies. In this study we analysed a highly informative portion of the mitochondrial DNA (mtDNA), widely used for assessing genetic diversity of honey bee populations. By analysing 1143 honey bees sampled after the introduction of this law, we provided a snapshot of the distribution of mtDNA haplotypes in this region. The two most frequent mtDNA haplotypes were C1 (characteristic of A. m. ligustica) and C2 (characteristic of A. m. carnica), reported in 86.5% and 11.0% of the analysed bees, respectively. About 1.3% and 1.1% of the analysed bees carried mtDNA haplotypes of the A and M lineages (haplotypes A1a, A1e, A4, A26, A65 and two novel ones, A2w and A6a; M3, M3a, M4 and M79). Continued genetic monitoring will be important to assess the impact of this regional law over the coming years. Based on the obtained results, we recommend a more stringent policy to prevent the erosion of the genetic integrity of the native subspecies A. m. ligustica.

New COI-COII mtDNA Region Haplotypes in the Endemic Honey Bees Apis mellifera intermissa and Apis mellifera sahariensis (Hymenoptera: Apidae) in Algeria

The practice of beekeeping in Algeria is of great cultural, social, and economic importance. However, the importation of non-local subspecies reported by beekeepers has disrupted the natural geographical distribution area and the genetic diversity of the native honey bees. To assess the genetic diversity of A. m. intermissa and A. m. sahariensis, and their relationships with African and European subspecies, the COI-COII intergenic region was analyzed in 335 individuals, 68 sampled in Algeria, 71 in Europe, Madagascar, and the South West Indian Ocean archipelagos, and 196 sequences recovered from GenBank. The results show the presence of the A lineage exclusively in Algerian samples with the identification of 24 haplotypes of which 16 are described for the first time. These haplotypes were found to be shared by both subspecies, with A74 being the most common haplotype in the population studied. The sequence comparison indicates the existence of three polymorphisms of the COI-COII marker: P0Q, P0QQ, and P0QQQ. One new haplotype was identified in the M lineage in samples from France. No evidence of genetic introgression within the Algerian honey bee population was detected. These data enhance our knowledge of the genetic diversity and emphasize the importance of protecting these local subspecies.

Estimates of the vitality and performances of Apis mellifera mellifera and hybrid honey bee colonies in Siberia: a 13-year study

Honeybees display a great range of biological, behavioral, and economic traits, depending on their genetic origin and environmental factors. The high diversity of honeybees is the result of natural selection of specific phenotypes adapted to the local environment. Of particular interest is adaptation of local and non-local bee colonies to environmental conditions. To study the importance of genotype-environment interactions on the viability and productivity of local and non-local bee colonies, we analyzed the long-term dynamics of the main traits in dark forest bees (Apis mellifera mellifera) and hybrid colonies. From 2010 to 2022, a total of 64 colonies living in an apiary in Siberia, Russia, were monitored and tested to assess their biological, behavioral, and economic traits in a temperate continental climate. We detected significant correlations between the studied biological and behavioral indicators of the bee colony such as colony strength, overwintering ability, infection of colonies with diseases, hygienic behavior, and others. No relationships between the biological and economic (honey productivity) traits of bee colonies are shown. The overall result of our study is that local dark forest bee, A.m.mellifera, showed higher values for all analyzed traits than hybrid colonies. Compared to hybrids, dark forest bee colonies showed more gentleness, productivity, and survivorship. The results from our study indicate a specific local adaptation of the A.m.mellifera subspecies in a temperate continental climate. Siberia represents a unique region for the conservation of the dark forest bee. The creation of conservation areas is one way to protect local bee populations, well adapted to local environmental conditions, from uncontrolled importation of bee breeds from different regions.

Haplotype Diversity in mtDNA of Honeybee in the Czech Republic Confirms Complete Replacement of Autochthonous Population with the C Lineage

The study aimed to analyze the genetic diversity in the Czech population of Apis mellifera using mitochondrial DNA markers, tRNAleu-cox2 intergenic region and cox1 gene. A total of 308 samples of bees were collected from the entire Czech Republic (from colonies and flowers in 13 different regions). Following sequencing, several polymorphisms and haplotypes were identified. Analysis of tRNAleu-cox2 sequences revealed three DraI haplotypes (C, A1, and A4). The tRNAleu-cox2 region yielded 10 C lineage haplotypes, one of which is a newly described variant. Three A lineage haplotypes were identified, two of which were novel. A similar analysis of cox1 sequences yielded 16 distinct haplotypes (7 new) within the population. The most prevalent tRNAleu-cox2 haplotype identified was C1a, followed by C2a, C2c, C2l, and C2d. For the cox1 locus, the most frequent haplotypes were HpB02, HpB01, HpB03, and HpB04. The haplotype and nucleotide diversity indices were high in both loci, in tRNAleu-cox2 with values of 0.682 and 0.00172, respectively, and in cox1 0.789 and 0.00203, respectively. The Tajima's D values were negative and lower in tRNAleu-cox2 than in cox1. The most frequent haplotypes were uniformly distributed across all regions of the Czech Republic. No haplotype of the indigenous M lineage was identified. High diversity and the occurrence of rare haplotypes indicate population expansion and continuous import of tribal material of the C lineage.

A preliminary survey reveals that common viruses are found at low titers in a wild population of honey bees (Apis mellifera)

A major threat to honey bee (Apis mellifera Linnaeus, Hymenoptera: Apidae) health continues to be parasitism by the mite Varroa destructor, which has been linked to high colony losses worldwide. Besides feeding on developing and adult bees, Varroa is also a prolific vector of honey bee-associated viruses. Because they live in unmanaged conditions, wild honey bee colonies are not treated against Varroa, which has enabled the natural selection of more mite-tolerant bees. To date, few studies have explored the prevalence of viruses in unmanaged colonies. The Welder Wildlife Refuge (WWR) in Texas is a unique site to study the viral landscape of unmanaged honey bees in the United States. The goals of this study were to identify and quantify viruses in wild colonies at the WWR, to examine changes in the prevalence of viruses in these colonies over time, and to compare the presence and titers of viruses between wild colonies at the WWR and those from the nearest managed apiary. We collected bees from colonies at the WWR in 2013, 2016, and 2021, and analyzed selected viruses for their presence and titers via quantitative polymerase chain reaction. In 2021, we also sampled bees from the nearest managed apiary for comparison. We found low average virus titers in all wild colonies sampled, and no difference in virus titers between colonies at the WWR and those from the managed apiary. Our study indicates that virus titers in wild colonies at the WWR are similar to those found in nearby colonies, and that these titers fluctuate over time.

Haplotypes of the tRNAleu-COII mtDNA Region in Russian Apis mellifera Populations

Analysis of the mtDNA tRNAleu-COII locus is a widely used tool to establish belonging to a particular evolutionary lineage of Apis mellifera L. (lineages A, M, C, O, and Y). In Russia, most of the area was once inhabited by Apis mellifera mellifera from the M evolutionary lineage, but the introduction of bee subspecies from the southern regions of Russia (A. m. caucasica, A. m. carnica) and from abroad (A. m. carnica, A. m. ligustica) led to fragmentation of their native range. In this study, the results of assessing the haplotype number for the tRNAleu-COII locus of mtDNA in Russian Apis mellifera populations were presented. We analyzed 269 colonies from 19 regions of Russia. As a result, two evolutionary lineages were identified: the East European lineage C (26.4%) and the Northwestern European lineage M (73.6%). A total of 29 haplotypes were identified, 8 of them were already reported, and 21 were found to be novel. From the C lineage, haplotypes C1, C2, C2c, C2j, and C3 were predominant. All M lineage samples from Russia belong to the M17 and M4' haplogroups but have only minor variations in the form of nucleotide substitutions. An analysis of publications devoted to the tRNAleu-COII locus haplotypes, as well as an analysis of the available tRNAleu-COII sequences in GenBank, showed that there is still a problem with the haplotype nomenclature.

Honey Quality Control: Review of Methodologies for Determining Entomological Origin

Honey is a widely consumed natural product, and its entomological origin can significantly influence its market value. Therefore, traceability of the entomological origin of honey should also be considered in honey quality control protocols. Although several methods exist, such as physicochemical characterization and bioactivity profiling of honey of different entomological origins, the most promising three methods for entomological authentication of honey include protein-based identification, chemical profiling, and a DNA-based method. All of these methods can be applied for reliable identification of the entomological origin of honey. However, as the honey is a complex matrix, the inconsistency of the results obtained by these methods is a pragmatic challenge, and therefore, the use of each method in all the cases is questionable. Most of these methodologies can be used for authentication of newly harvested honey and it is worth understanding the possibility of using these methods for authentication of relatively old samples. Most probably, using DNA-based methods targeting small fragments of DNA can provide the best result in old samples, however, the species-specific primers targeting short fragments are limited and not available for all species. Therefore, using universal primers in combination with a DNA metabarcoding approach can be a good solution that requires further investigation. This present article describes the applications of different methods, their pros, and their cons to identify honey based on entomological origin.

Defining honeybee subspecies in an evolutionary context warrants strategized conservation

Despite the urgent need for conservation consideration, strategic action plans for the preservation of the Asian honeybee, Apis cerana Fabricius, 1793, remain lacking. Both the convergent and divergent adaptations of this widespread insect have led to confusing phenotypical traits and inconsistent infraspecific taxonomy. Unclear subspecies boundaries pose a significant challenge to honeybee conservation efforts, as it is difficult to effectively prioritize conservation targets without a clear understanding of subspecies identities. Here, we investigated genome variations in 362 worker bees representing almost all populations of mainland A. cerana to understand how evolution has shaped its population structure. Whole-genome single nucleotide polymorphisms (SNPs) based on nuclear sequences revealed eight putative subspecies, with all seven peripheral subspecies exhibiting mutually exclusive monophyly and distinct genetic divergence from the widespread central subspecies. Our results demonstrated that most classic morphological traits, including body size, were related to the climatic variables of the local habitats and did not reflect the true evolutionary history of the organism. Thus, such morphological traits were not suitable for subspecific delineation. Conversely, wing vein characters showed relative independence to the environment and supported the subspecies boundaries inferred from nuclear genomes. Mitochondrial phylogeny further indicated that the present subspecies structure was a result of multiple waves of population divergence from a common ancestor. Based on our findings, we propose that criteria for subspecies delineation should be based on evolutionary independence, trait distinction, and geographic isolation. We formally defined and described eight subspecies of mainland A. cerana. Elucidation of the evolutionary history and subspecies boundaries enables a customized conservation strategy for both widespread and endemic honeybee conservation units, guiding colony introduction and breeding.

Effects of local domestication warrant attention in honey bee population genetics

Honey bees, Apis mellifera, have for millennia been managed and exploited by humans and introduced into most suitable regions worldwide. However, given the lack of records for many introduction events, treating A. mellifera populations as native would predictably bias genetic studies regarding origin and evolution. Here, we used the Dongbei bee, a well-documented population, introduced beyond the natural distribution range approximately 100 years ago, to elucidate the effects of local domestication on animal population genetic analyses. Strong domestication pressure was detected in this population, and the genetic divergence between Dongbei bee and its ancestral subspecies was found to have occurred at the lineage level. Results of phylogenetic and time divergence analyses could consequently be misinterpreted. Proposing new subspecies or lineages and performing analyses of origin should thus strive to eliminate anthropogenic effects. We highlight the need for definitions of landrace and breed in honey bee sciences and make preliminary suggestions.

Meta-analysis of genetic diversity and intercolony relatedness among reproductives in commercial honey bee populations

Honey bee colonies are large kin groups, each with a single mother queen and thousands of female workers. Queen bees are highly polyandrous, each mating with an average of approximately 12 drones from other colonies. We used a meta-analysis approach to compare the pedigree relationships of honey bee reproductives (queens and their mates) across five different studies and to quantify the overall genetic diversity of breeding populations. We compared the inferred genotypes of queens and their mates from microsatellite analyses of worker offspring from a feral Africanized honey bee population (which served as a negative control for inbreeding), an experimentally derived population of sister queens (which served as a positive control for inbreeding), and three separate commercially managed populations. We then compared the relatedness of all drones mated to each queen (mate-mate), all queens within each population (queen-queen), each queen with each of her mates (queen-mate), and all drones within each population (drone-drone). We found, as expected, the lowest levels of genetic similarity in the outcrossed population and highest levels of genetic similarity in the inbred population. Levels of genetic similarity among the managed honey bee populations were intermediate but closer to that of the inbred population. Genetic structuring of the entire breeding population resulted in two major subpopulations, likely deriving from breeders on the east and west coast. The effects that these findings have on the overall population genetic diversity of managed honey bees is discussed.

Honey bee populations of the USA display restrictions in their mtDNA haplotype diversity

The genetic diversity of the USA honey bee (Apis mellifera L.) populations was examined through a molecular approach using two mitochondrial DNA (mtDNA) markers. A total of 1,063 samples were analyzed for the mtDNA intergenic region located between the cytochrome c oxidase I and II (COI-COII) and 401 samples were investigated for the NADH dehydrogenase 2 (ND2) coding gene. The samples represented 45 states, the District of Colombia and two territories of the USA. Nationwide, three maternal evolutionary lineages were identified: the North Mediterranean lineage C (93.79%), the West Mediterranean lineage M (3.2%) and the African lineage A (3.01%). A total of 27 haplotypes were identified, 13 of them (95.11%) were already reported and 14 others (4.87%) were found to be novel haplotypes exclusive to the USA. The number of haplotypes per state/territory ranged between two and eight and the haplotype diversity H ranged between 0.236-0.763, with a nationwide haplotype diversity of 0.597. Furthermore, the honey bee populations of the USA were shown to rely heavily (76.64%) on two single haplotypes (C1 = 38.76%, C2j = 37.62%) of the same lineage characterizing A. m. ligustica and A. m. carnica subspecies, respectively. Molecular-variance parsimony in COI-COII and ND2 confirmed this finding and underlined the central and ancestral position of C2d within the C lineage. Moreover, major haplotypes of A. m. mellifera (M3a, M7b, M7c) were recorded in six states (AL, AR, HI, MO, NM and WA). Four classic African haplotypes (A1e, A1v, A4, A4p) were also identified in nine states and Puerto Rico, with higher frequencies in southern states like LA, FL and TX. This data suggests the need to evaluate if a restricted mtDNA haplotype diversity in the US honey bee populations could have negative impacts on the beekeeping sustainability of this country.

Cutting corners: The impact of storage and DNA extraction on quality and quantity of DNA in honeybee (Apis mellifera) spermatheca

The purpose of our study was to investigate methods of short-term storage that allow preservation, transport and retrieval of genetic information contained in honeybee queen's spermatheca. Genotyping of the honeybee colony requires well ahead planned sample collection, depending on the type of data to be acquired. Sampling and genotyping of spermatheca's content instead of individual offspring is timesaving, allowing answers to the questions related to patriline composition immediately after mating. Such procedure is also cheaper and less error prone. For preservation either Allprotect Tissue Reagent (Qiagen) or absolute ethanol were used. Conditions during transportation were simulated by keeping samples 6-8 days at room temperature. Six different storing conditions of spermathecas were tested, complemented with two DNA extraction methods. We have analysed the concentration of DNA, RNA, and proteins in DNA extracts. We also analysed how strongly the DNA is subjected to fragmentation (through amplification of genetic markers ANT2 and tRNA^leu-COX2) and whether the quality of the extracted DNA is suitable for microsatellite (MS) analysis. Then, we tested the usage of spermatheca as a source of patriline composition in an experiment with three instrumentally inseminated virgin queens and performed MS analysis of the extracted DNA from each spermatheca, as well as queens' and drones' tissue. Our results show that median DNA concentration from spermathecas excised prior the storage, regardless of the storing condition and DNA extraction method, were generally lower than median DNA concentration obtained from spermathecas dissected from the whole queens after the storage. Despite the differences in DNA yield from the samples subjected to different storing conditions there was no significant effect of storage method or the DNA extraction method on the amplification success, although fewer samples stored in EtOH amplified successfully in comparison to ATR storing reagent. However, we recommend EtOH as a storing reagent due to its availability, low price, simplicity in usage in the field and in the laboratory, and capability of good preservation of the samples for DNA analysis during transport at room temperature.

A genotyping by sequencing approach can disclose Apis mellifera population genomic information contained in honey environmental DNA

Awareness has been raised over the last years on the genetic integrity of autochthonous honey bee subspecies. Genomic tools available in Apis mellifera can make it possible to measure this information by targeting individual honey bee DNA. Honey contains DNA traces from all organisms that contributed or were involved in its production steps, including the honey bees of the colony. In this study, we designed and tested a genotyping by sequencing (GBS) assay to analyse single nucleotide polymorphisms (SNPs) of A. mellifera nuclear genome using environmental DNA extracted from honey. A total of 121 SNPs (97 SNPs informative for honey bee subspecies identification and 24 SNPs associated with relevant traits of the colonies) were used in the assay to genotype honey DNA, which derives from thousands of honey bees. Results were integrated with information derived from previous studies and whole genome resequencing datasets. This GBS method is highly reliable in estimating honey bee SNP allele frequencies of the whole colony from which the honey derived. This assay can be used to identify the honey bee subspecies of the colony that produced the honey and, in turn, to authenticate the entomological origin of the honey.

Urban ecosystem drives genetic diversity in feral honey bee

Urbanization can change biodiversity in both directions, positive and negative, and despite the rising global trend of urban beekeeping, little is known about the impact of urbanization on the genetic diversity of honey bees. We investigate how urbanization affects the genetic variability of feral and managed honey bee colonies that are spread throughout the entire city, even in highly urban areas, through genetic analysis of 82 worker bees. We found convincing evidence of high genetic differentiation between these two groups. Additionally, by comparing city samples with 241 samples from 46 apiaries in rural parts of the country, variations in mitochondrial tRNAleu-cox2 intergenic region and microsatellite loci indicated that feral colonies have distinct patterns of genetic diversity. These results, with evidence that feral honey bees find niches within highly modified and human-dominated urban landscapes, lead us to conclude that urbanization is a driver of the genetic diversity of feral honey bees in the city.

Africanized honeybee population (Apis mellifera L.) in Nicaragua: Forewing length and mitotype lineages

Various subspecies of Apis mellifera L. were introduced to Central America since colonization 500 years ago. Hybridization increased with the entrance of the Africanized bee in Nicaragua in 1984. Rustic beekeeping activities and numerous feral swarms define the genetic pattern, reflected in phenotypic heterogeneity and remarkable differences in the behaviour of the bee colonies, especially the nest defence. Due to these facts, the question emerge about the degree of Africanization of honeybee colonies in Nicaragua. In this study, we identified Africanized honeybee colonies based on the single character "mean forewing length" and we corroborated our results by determining mitotypes using mtDNA analysis. Morphometric and genetic approaches were realized in three different geographical zones of Nicaragua and related to beehive characteristics and management. Worker bee samples were taken from the inside of 146 hives from 26 apiaries. Abdominal colour as phenotypic character was the first examination, followed by measurement of 1460 right forewings to determine corresponding probability of Africanization. More than 60% of the beehives showed phenotypic heterogeneity and mean forewing length of 8.74 mm (SD 0.16 mm) indicated a high degree of Africanization. Those results provided a selection of 96 worker bees to perform PCR of two worker bees per hive. For mitochondrial DNA analysis 14 samples from sentinel apiaries were added. Three from 61 beehives presented bees with different mtDNA. Throughout, three mitotypes of the African (A) lineage were detected; one mitotype is still unidentified. Mitotype A1 A. mellifera iberiensis was represented by 88 bees and mitotype A4 A. mellifera scutellata by 21 bees. Phylogenetic analysis confirmed the PCR findings. No associations were found between mitotypes, forewing length, beehive characteristics and management. A high degree of Africanization in A. mellifera colonies represented by two predominating mitotypes from the A lineage, prevail in Neotropical Nicaragua, with mitotype A4 predominating at higher altitudes.

The Carniolan Honeybee from Slovenia—A Complete and Annotated Mitochondrial Genome with Comparisons to Closely Related Apis mellifera Subspecies

Simple Summary

The western honeybee, Apis mellifera, is a globally distributed bee species with many recognised subspecies, one of which is Apis mellifera carnica, the Carniolan honeybee. Apis m. carnica is native to southern Central Europe and parts of the Balkans, with the locus classicus in Slovenia. It is also widely popular with beekeepers in parts of Central and Northern Europe and other parts of the world, including the USA, Canada, and even New Zealand. In Slovenia, A. m. carnica is protected, with measures to conserve the subspecies’ autochthonous domestic population in place. Such efforts depend heavily upon genomic and phylogenetic information. In this study, we sequenced and annotated the mitochondrial genome of a specimen from Slovenia and compared the obtained data with a previously published sample of the A. m. carnica from Austria and the closely related Italian honeybee A. m. ligustica. We found several features unique to the new mitochondrial genome. We also phylogenetically analyzed the relationship between our sequence and the selected available A. mellifera mitochondrial sequences. The acquired position of the sequenced A. m. carnica from Slovenia on the phylogenetic tree brings new evidence for close relationships among C and O lineages and reflects their recent historical matrilinear ancestry.

Abstract

The complete mitochondrial genome of the Carniolan honeybee (Apis mellifera carnica) from Slovenia, a homeland of this subspecies, was acquired in two contigs from WGS data and annotated. The newly obtained mitochondrial genome is a circular closed loop of 16,447 bp. It comprises 37 genes (13 protein coding genes, 22 tRNA genes, and 2 rRNA genes) and an AT-rich control region. The order of the tRNA genes resembles the order characteristic of A. mellifera. The mitogenomic sequence of A. m. carnica from Slovenia contains 44 uniquely coded sites in comparison to the closely related subspecies A. m. ligustica and to A. m. carnica from Austria. Furthermore, 24 differences were recognised in comparison between A. m. carnica and A. m. ligustica subspecies. Among them, there are three SNPs that affect translation in the nd2, nd4, and cox2 genes, respectively. The phylogenetic placement of A. m. carnica from Slovenia within C lineage deviates from the expected position and changes the perspective on relationship between C and O lineages. The results of this study represent a valuable addition to the information available in the phylogenomic studies of A. mellifera—a pollinator species of worldwide importance. Such genomic information is essential for this local subspecies’ conservation and preservation as well as its breeding and selection.

Assessment of genetic variation in Apis mellifera jemenitica (Hymenoptera: Apidae) based on mitochondrial Cytochrome Oxidase Subunit II and III

Morphometric and genetic characterization of many Apis mellifera subspecies are well-documented. A. m. jemenetica occurs naturally in Africa and Asia. In this study, genetic variation of mitochondrial Cytochrome Oxidase II (COII) and III (COIII) were analysed in 133 specimens of the endemic honeybee colonies within Saudi Arabia. The COII gene sequence length was 684 bp comprising nine synonymous (1.3%) and two non-synonymous single nucleotide polymorphisms (SNPs) (0.87%). Five variants of COII were not previously documented, one variant (MT755968) showed an extra restriction site when subjected to type II restriction endonuclease from Arthrobacter protophormiae (Apol) or to Haemophilus influenzae Rf (Hinf1). Changes in COII sequence separated samples into three haplogroups. Whereas, COIII gene sequence length was 780 bp, including 18 synonymous and five non-synonymous SNPs. Furthermore, variation in COII sequence was more informative based on restriction profiles and on amino acid changes compared with COIII gene sequence. Variants of COIII showed identical restriction sites when subjected to type II restriction endonuclease from Deinococcus radiophilus (DraI), and revealed high similarity to African subspecies. Results of this study are very useful in understanding genetic diversity and characterization of A. mellifera subspecies.

Further Evidence of Population Admixture in the Serbian Honey Bee Population

Socioeconomic interests and beekeeper preferences have often taken precedence over the conservation of locally native honey bee subspecies, leading to the predominance of admixture populations in human-dominated areas. To assess the genetic diversity of contemporary managed Serbian honey bee colonies, we used 14 microsatellite loci and analyzed 237 worker bees from 46 apiaries in eight localities of northern and southern Serbia. Furthermore, we compared data for nine microsatellite loci with 338 individuals from Italy, Hungary, Poland, and Spain. The standard parameters of genetic diversity in Serbian honey bee populations were in line with other analyses, although somewhat smaller. STRUCTURE analysis showed the existence of two equally distributed genetic clusters and Analysis of molecular variances could not confirm the presence of a geographically discrete population but showed local differences. Discriminant analysis of principal components showed overlapping of worker bees from different parts of Serbia. Clear genetic differentiation can be observed when comparing all populations between geographical regions and their corresponding subspecies. The absence of the A. m. macedonica subspecies from its historical distribution range in southern Serbia as well as the lack of distinctive geographical groups suggest that selective breeding, queen import, and migratory beekeeping practices strongly influence the genetic structure and diversity of honey bees, leading to the genetic uniformization and creation of the admixture population.

Population Structure and Diversity in European Honey Bees (Apismellifera L.)-An Empirical Comparison of Pool and Individual Whole-Genome Sequencing

BACKGROUND

Whole-genome sequencing has become routine for population genetic studies. Sequencing of individuals provides maximal data but is rather expensive and fewer samples can be studied. In contrast, sequencing a pool of samples (pool-seq) can provide sufficient data, while presenting less of an economic challenge. Few studies have compared the two approaches to infer population genetic structure and diversity in real datasets. Here, we apply individual sequencing (ind-seq) and pool-seq to the study of Western honey bees (Apis mellifera).

METHODS

We collected honey bee workers that belonged to 14 populations, including 13 subspecies, totaling 1347 colonies, who were individually (139 individuals) and pool-sequenced (14 pools). We compared allele frequencies, genetic diversity estimates, and population structure as inferred by the two approaches.

RESULTS

Pool-seq and ind-seq revealed near identical population structure and genetic diversities, albeit at different costs. While pool-seq provides genome-wide polymorphism data at considerably lower costs, ind-seq can provide additional information, including the identification of population substructures, hybridization, or individual outliers.

CONCLUSIONS

If costs are not the limiting factor, we recommend using ind-seq, as population genetic structure can be inferred similarly well, with the advantage gained from individual genetic information. Not least, it also significantly reduces the effort required for the collection of numerous samples and their further processing in the laboratory.

Thrice out of Asia and the adaptive radiation of the western honey bee

The origin of the western honey bee Apis mellifera has been intensely debated. Addressing this knowledge gap is essential for understanding the evolution and genetics of one of the world’s most important pollinators. By analyzing 251 genomes from 18 native subspecies, we found support for an Asian origin of honey bees with at least three expansions leading to African and European lineages. The adaptive radiation of honey bees involved selection on a few genomic “hotspots.” We found 145 genes with independent signatures of selection across all bee lineages, and these genes were highly associated with worker traits. Our results indicate that a core set of genes associated with worker and colony traits facilitated the adaptive radiation of honey bees across their vast distribution.

Hornets and Honey Bees: A Coevolutionary Arms Race between Ancient Adaptations and New Invasive Threats

Hornets and honey bees have a long history of coevolution resulting in a plethora of captivating adaptations and counteradaptations between predator and prey. From simple physiological mechanisms to complex behavioral strategies, some Vespa hornets have specialized in hunting honey bees, while the latter have put in place effective defenses to counteract their attack. Both hornets and honey bees have evolved the ability to detect the odors and the pheromones emitted by the other to locate the prey or to spot foraging predators. Hornets often rely on their bigger size, heavily armored body and destructive attacks, while honey bees differentiated collective defense responses finely coordinated to deter or kill the hornet menace. However, when new species of hornets and honey bees come into contact, the absence of coevolution can have a heavy impact on the defenseless bees. The evolutionary arms race between hornets and honey bees provides not only compelling examples of adaptations and counteradaptations between predator and prey, but could also represent a starting point for the development of effective and sustainable strategies to protect honey bees and beekeeping activities and to control invasive alien species of hornets.

Assessment of genetic variation in Apis mellifera jemenitica (Hymenoptera: Apidae) using Cytochrome Oxidase I gene sequences

The Arabian Honeybee Apis mellifera jemenitica is endemic to the Arabian Peninsula. It is highly adapted to temperature extremes and drought dominating the region. In this study, the mitochondrial Cytochrome Oxidase I (COI) was analyzed in 133 specimens of A. m. jemenitica from eight localities along the Red Sea cost of Saudi Arabia. Results revealed 33 synonymous, and 6 non-synonymous mutations within the COI sequences, resulting in change of 4 amino acids. Phylogenetic analysis based on either type of mutations revealed two main haplogroups accounting for 94% of the samples. In total Eighteen new haplotypes were identified and uploaded in the genebank, Fourteen of them are restricted to one/both haplogroups. All haplotypes identified in this study clustered with reference COI sequences of the sub-lineag Z (African Lineage). However one Haplotype (MW428270) represents high COI variability compared to other haplotypes and may resemble different evolutionary sub-lineage. Tajima's Neutrality Test (Ps = 0.025; D = -1.5) indicated population size expansion that took place after selective sweep and/or purifying selection.

MtDNA Analysis Indicates Human-Induced Temporal Changes of Serbian Honey Bees Diversity

Simple Summary

The western honey bee is one of the most economically and ecologically important species currently facing serious challenges in its whole area of distribution. The honey bee is a highly diverse species with about 30 subspecies that are adapted to regional climate factors, vegetation, pests and pathogens. The local populations of honey bees are rapidly changing and their diversity is constantly manipulated by beekeepers through the import of foreign queens, selection and migratory beekeeping. This manipulation may lead to such changes that honey bees lose their ability to thrive in the areas that were previously suitable for their wellbeing. To see how this human interference changed the genetic variability of native honey bee populations from Serbia, we sequenced part of the mitochondrial genome and compared them with published sequences. Our results suggest that human influence significantly changes the natural composition of honey bees in Serbia and that the presence of some previously reported subspecies could not be confirmed.

Abstract

Local populations of Apis mellifera are rapidly changing by modern beekeeping through the introduction of nonnative queens, selection and migratory beekeeping. To assess the genetic diversity of contemporary managed honey bees in Serbia, we sequenced mitochondrial tRNA^leu-cox2 intergenic region of 241 worker bees from 46 apiaries at eight localities. Nine haplotypes were observed in our samples, with C2d being the most common and widespread. To evaluate genetic diversity patterns, we compared our data with 1696 sequences from the NCBI GenBank from neighbouring countries and Serbia. All 32 detected haplotypes belonged to the Southeast Europe lineage C, with two newly described haplotypes from our sample. The most frequent haplotype was C2d, followed by C2c and C1a. To distinguish A. m. carnica from A. m. macedonica, both previously reported in Serbia, PCR-RFLP analysis on the COI gene segment of mtDNA was used, and the result showed only the presence of A.m. carnica subspecies. An MDS plot constructed on pairwise F_ST values showed significant geographical stratification. Our samples are grouped together, but distant from the Serbian dataset from the GenBank. This, with the absence of A. m. macedonica subspecies from its historic range of distribution in southern Serbia, indicates that honey bee populations are changing rapidly due to the anthropogenic influence.

Mitochondrial DNA Variation in Peruvian Honey Bee (Apis mellifera L.) Populations Using the tRNAleu-cox2 Intergenic Region

Simple Summary

Currently, the genetic diversity of Peruvian honey bee populations is unknown. Only two studies were carried out and suggest that many regions of Peru present Africanized honey bee colonies and a varied degree of Africanization. To molecularly characterize and know more about the genetic background of Peruvian honey bees, the highly polymorphic tRNAleu-cox2 was used. This study analyzed 512 colonies in three regions of Peru: Lima, Piura, and Junín. The results indicated that 65% colonies correspond to lineage A (African), 33.8% colonies to lineage C (Eastern European), and 1.2% colonies to lineage M (Western European). A total of 24 haplotypes were identified: 16 haplotypes belong to lineage A (sub-lineage A_I (13), sub-lineage A_III (03)), lineage C (06), and lineage M (02), and 15 of them are for the first time reported and represented by A1t, A1u, A1w, A4p, A4q, A4s, A4t, A4u, A4v, A4w, 30d, A30e, A65, M7b, and M7c. Piura and Lima presented higher proportions of African haplotypes and lower proportions of haplotypes from lineage C than Lima. Very few haplotypes of lineage M were identified, whose presence could be due to accidental purchases or traces of honey bee introductions from lineage M in the 19th century. Hence, studies about the diversity and genetic structure of Peruvian honey bee populations are necessary to promote adequate, sustainable management and establish conservation and breeding programs.

Abstract

Mitochondrial DNA variations of Peruvian honey bee populations were surveyed by using the tRNA^leu-cox2 intergenic region. Only two studies have characterized these populations, indicating the presence of Africanized honey bee colonies in different regions of Peru and varied levels of Africanization, but the current status of its genetic diversity is unknown. A total of 512 honey bee colonies were sampled from three regions to characterize them. Our results revealed the presence of European and African haplotypes: the African haplotypes identified belong to sub-lineage A_I (13) and sub-lineage A_III (03), and the European haplotypes to lineages C (06) and M (02). Of 24 haplotypes identified, 15 new sequences are reported here (11 sub-lineage A_I, 2 sub-lineage A_III, and 2 lineage M). Peruvian honey bee populations presented a higher proportion from African than European haplotypes. High proportions of African haplotype were reported for Piura and Junín, unlike Lima, which showed more European haplotypes from lineage C. Few colonies belonging to lineage M would represent accidental purchase or traces of the introduction into Peru in the 19th century.

Distribution of the Main Apis mellifera Mitochondrial DNA Lineages in Italy Assessed Using an Environmental DNA Approach

Growing interest has been emerging on the need to monitor the genetic integrity of the European Apis mellifera subspecies that could be threatened by the human-mediated dispersion of non-native populations and lines. Mitochondrial DNA (mtDNA) lineages can provide useful information for this purpose. In this study, we took advantage of the environmental DNA (eDNA) contained in the honey, which can be analyzed to detect the main groups of mitotypes of the honey bees that produced it. In this study, we applied this eDNA to produce a distribution map all over the Italian peninsula and the two major islands (Sicily and Sardinia) of the following three honey bee mtDNA lineages: A, C and M. A total of 607 georeferenced honey samples, produced in all Italian regions, was analyzed to detect these lineages. The A lineage was widespread in Sicily, as expected, considering that A. m. siciliana carries the African lineage. Surprisingly, this lineage was also reported in about 14% of all other samples produced in almost all continental regions, and in Sardinia. The applied method obtained an updated distribution map of honey bee mtDNA lineages that could be useful to design policies for the conservation of Italian honey bee genetic resources.

Genetic network analysis between Apis mellifera subspecies based on mtDNA argues the purity of specimens from North Africa, the Levant and Saudi Arabia

Objectives

This study aimed to analyze the genetic relationships between honey bee subspecies using reference specimens and recently collected specimens from different parts of the world. The purity of these specimens was discussed in light of the obtained results.

Methods

The genetic networks were constructed between 21 subspecies of honey bees, Apis mellifera L.: 9 in Africa, 7 in Europe and 5 in Asia. The analysis was performed using the mtDNA of these subspecies and the Population Analysis with Reticulate Trees software. Some subspecies were represented by more than two specimens based on the available online sequences.

Results and conclusions

The subspecies A. m. sahariensis from Africa showed unique characteristics and is genetically isolated than all other studied bee subspecies. Specimens collected from Saudi Arabia showed genetic relatedness to A. m. jemenitica, A. m. lamarckii, and some European subspecies, suggesting high degree of hybridization. The close genetic relationship between the Egyptian bees, A. m. lamarckii, and the Syrian bees, A. m. syriaca, were emphasized. The overall genetic network showed the presence of three distinct branches in relation to geographical locations. The high accurateness of the used analysis was confirmed by previous phylogenetic studies as well as the genetic relationships between hybrid bees of A. m. capensis and A. m. scutellata. The genetic networks showed the presence of bee subspecies from Africa in all branches including Europe and Asia. The study suggests the impurity of some specimens mostly due to the hybridization between subspecies. Specific recommendations for future conservation efforts of bees were presented in light of this study.

Mitochondrial DNA Suggests the Introduction of Honeybees of African Ancestry to East-Central Europe

Simple Summary

In Europe, a well-known threat to the conservation of honeybee diversity is the loss of genetic uniqueness of local populations due to beekeepers’ preference for a few genetic lineages. However, due to climate change and large-scale ongoing movement of breeding individuals, the expansion of bees of African origin could represent another threat. This issue has not yet been recognised in detail, although bees bearing African mitochondrial DNA occur in South-West and South Europe due to natural gene flow. Here, we determine the diversity of mitochondrial DNA in honey bees from East-Central Europe. We sequenced the COI-COII region in 427 bees sampled along two 900 km transects (17.5° N and 23° E). We found that 1.64% of bees (95% CI: 0.66–3.35%) had African mitochondrial DNA. It is unlikely that their presence in the area resulted from natural migration but instead human-driven introductions of hybrids of African ancestry. This expansion deserves more attention, as it may contribute to the dissemination of undesirable traits, parasites and diseases.

Abstract

In Europe, protecting the genetic diversity of Apis mellifera is usually perceived in the context of limiting the spread of the evolutionary C-lineage within the original range of the M-lineage. However, due to climate change and large-scale ongoing movement of breeding individuals, the expansion of bees from the African A-lineage could represent another threat. This issue has not yet been investigated in detail, although A-mitotypes occur in South-West and South Europe due to natural gene flow. Here, we determine the diversity of mtDNA in honey bees from East-Central Europe. We sequenced the COI-COII region in 427 bees sampled along two 900 km transects (17.5° N and 23° E). We found that 1.64% of bees (95% CI: 0.66–3.35 %) had A-mitotypes. It is unlikely that their presence in the area resulted from natural migration but instead human driven introductions of hybrids of African ancestry. This expansion deserves more attention, as it may contribute to the dissemination of undesirable traits, parasites and diseases.

Coexistence of honeybees with distinct mitochondrial haplotypes and hybridised nuclear genomes on the Comoros Islands

The honeybee, Apis mellifera, is a globally distributed species that has spread both naturally and by humans across the globe resulting in many natural and secondary contact zones. The geographic isolation of honeybees is likely to contribute to genetic differentiation. Secondary contact has resulted in hybridization at the nuclear genome, but replacement of mitochondrial. Here, we used a mitochondrial marker and 19 microsatellite markers to test for the variations in the mitochondrial and nuclear genomes of honeybee populations on the Comoros islands. We used samples of 160 workers for mtDNA analysis and 288 workers from 16 colonies spread across the three islands for microsatellite analyses. Our results showed that the wild honeybee populations of the Comoros Islands consist of coexisting mitochondrial haplotypes. One belongs to the typical African A-lineage, and the other, the newly described L-lineage, is closely related to Apis koschevnikovi, a honeybee species native to Southeast Asia. The nuclear genomes show complete hybridization, high genetic diversity, and strong differentiation according to the island of origin. Based on our results, we hypothesise that the Asian honeybee could have been transported from Southeast Asia to Madagascar and Comoros via the human migrations that occurred 6000 years ago, and has hybridised with African honeybees at the nuclear genome, but maternal ancestry still can be traced using the mtDNA markers. We conclude that mtDNA plays a pivotal role in adaptation to the local environment, with both haplotypes of the honeybees of Comoros contributing significantly to the mito-nuclear coadaptation resulting in maintenance at almost equal frequency.

Disentangling Ethiopian Honey Bee (Apis mellifera) Populations Based on Standard Morphometric and Genetic Analyses

Simple Summary

We conducted this population study of Ethiopian honey bees, using morphometric and genetic methods, to decipher their controversial classification. These honey bees are highly diverse and showed differentiation based on size and genetic information according to prevailing agro-ecological conditions, demonstrating morphological and molecular signatures of local adaptation. The results of both morphometric and genetic analyses suggest that Ethiopian honey bees differ from populations in the neighboring geographic regions and are characterized by extensive gene flow within the country, enhanced by honey bee colony trade. Consequently, future research that includes studying traits of vitality, behavior and colony performance of honey bees in remaining pocket areas of highland agro-ecological zones could contribute to the development of appropriate conservation management.

Abstract

The diversity and local differentiation of honey bees are subjects of broad general interest. In particular, the classification of Ethiopian honey bees has been a subject of debate for decades. Here, we conducted an integrated analysis based on classical morphometrics and a putative nuclear marker (denoted r7-frag) for elevational adaptation to classify and characterize these honey bees. Therefore, 660 worker bees were collected out of 66 colonies from highland, midland and lowland agro-ecological zones (AEZs) and were analyzed in reference to populations from neighboring countries. Multivariate morphometric analyses show that our Ethiopian samples are separate from Apis mellifera scutellata, A. m. jemenitica, A. m. litorea and A. m. monticola, but are closely related to A. m. simensis reference. Linear discriminant analysis showed differentiation according to AEZs in the form of highland, midland and lowland ecotypes. Moreover, size was positively correlated with elevation. Similarly, our Ethiopian samples were differentiated from A. m. monticola and A. m. scutellata based on r7-frag. There was a low tendency towards genetic differentiation between the Ethiopian samples, likely impacted by increased gene flow. However, the differentiation slightly increased with increasing elevational differences, demonstrated by the highland bees that showed higher differentiation from the lowland bees (F_ST = 0.024) compared to the midland bees (F_ST = 0.015). An allelic length polymorphism was detected (denoted as d) within r7-frag, showing a patterned distribution strongly associated with AEZ (X² = 11.84, p < 0.01) and found predominantly in highland and midland bees of some pocket areas. In conclusion, the Ethiopian honey bees represented in this study are characterized by high gene flow that suppresses differentiation.

Authoritative subspecies diagnosis tool for European honey bees based on ancestry informative SNPs

BACKGROUND

With numerous endemic subspecies representing four of its five evolutionary lineages, Europe holds a large fraction of Apis mellifera genetic diversity. This diversity and the natural distribution range have been altered by anthropogenic factors. The conservation of this natural heritage relies on the availability of accurate tools for subspecies diagnosis. Based on pool-sequence data from 2145 worker bees representing 22 populations sampled across Europe, we employed two highly discriminative approaches (PCA and F_ST) to select the most informative SNPs for ancestry inference.

RESULTS

Using a supervised machine learning (ML) approach and a set of 3896 genotyped individuals, we could show that the 4094 selected single nucleotide polymorphisms (SNPs) provide an accurate prediction of ancestry inference in European honey bees. The best ML model was Linear Support Vector Classifier (Linear SVC) which correctly assigned most individuals to one of the 14 subspecies or different genetic origins with a mean accuracy of 96.2% ± 0.8 SD. A total of 3.8% of test individuals were misclassified, most probably due to limited differentiation between the subspecies caused by close geographical proximity, or human interference of genetic integrity of reference subspecies, or a combination thereof.

CONCLUSIONS

The diagnostic tool presented here will contribute to a sustainable conservation and support breeding activities in order to preserve the genetic heritage of European honey bees.

Exploring the non-coding regions in the mtDNA of some honey bee species and subspecies

The sequence of the DNA contains coding and non-coding regions. The role of the non-coding regions is not known and is hypothesized to maintain the structure of the DNA. This study aimed to investigate the structure of the non-coding sequences in honey bees utilizing bioinformatics. The non-coding sequences of the mtDNA of three honey bee species Apis dorosata, Apis florea, Apis cerana, and ten subspecies of Apis mellifera were investigated. Different techniques were utilized to explore the non-coding regions of these bees including sequence analysis, phylogenetic relationships, enzymatic digestion, and statistical tests. Variations in size and sequences of nucleotides were detected in the studied species and subspecies, but with the same nucleotide abundance (i.e. nucleotides A were more than T and nucleotides G were less than C). The phylogenetic tree based on the non-coding regions was partially similar to the known phylogenetic relationships between these bees. The enzymatic digestion using four restriction enzymes confirmed the results of the phylogenetic relationships. The statistical analysis based on numerical codes for nucleotides showed the absence of significant variations between the studied bees in their sequences in a similar way to results of neutrality tests. This study suggests that the non-coding regions have the same functional role in all the studied bees regardless of the number of nucleotides, and not just to maintain the structure of the DNA. This is approximately the first study to shade lights on the non-coding regions of the mtDNA of honey bees.

Non-Invasive Genotyping of Honey Bee Queens Apis mellifera L.: Transition of the DraI mtDNA COI-COII Test to In Silico

Simple Summary

Each honey bee colony has a single queen which plays a crucial role in the survival and wellbeing of the entire hive. Honey bee genetic analysis and selection and breeding programs rely on destructive methods using worker bees; workers are numerous in a colony and can be quickly replaced. In this study, we tested and validated a fast and efficient non-destructive method to conduct genetic analysis directly on honey bee queens. We also describe a new method for the transition to in silico of a widely used honey bee genetic marker by reconciling both cleaved amplified polymorphic sequences and Sanger sequencing approaches. Both new approaches will provide significant service to honey bee breeding and selection programs, as well as facilitating and standardizing honey bee haplotype identification among research institutions.

Abstract

The honey bee Apis mellifera L. colony is headed by a single and indispensable queen, whose duty it is to ensure brood production and provide pheromonal stability within the colony. This study presents a non-invasive method that allows the identification of the queen maternal lineage and subspecies using the remaining tissue of her clipped wing. The DraI mtDNA COI-COII (DmCC) test was applied to various sizes of queen and worker wings and the results were compared with data obtained from other bee tissues. Furthermore, we propose a new method allowing in silico transition of the DmCC test and haplotype identification based on extended sequencing of the tRNAleu and COII genes. Our results show that DNA extracted by Chelex 10% from one-third of a queen’s wing is deemed adequate for a successful identification of her maternal evolutionary lineage, haplotype and subspecies. The in silico method proposed in this study fully adheres to the established guidelines of the DmCC, provides a universal standard for haplotype identification, and offers faster and more precise results by reconciling both cleaved amplified polymorphic sequences (CAPS) and Sanger sequencing approaches.

A revision of subspecies structure of western honey bee Apis mellifera

The taxonomy of honey bee A. mellifera contains a lot of issues due to the specificity of population structure, features of biology and resolutions of honey bee subspecies discrimination methods. There are a lot of transition zones between ranges of subspecies which led to the gradual changes of characteristics among neighbor subspecies. The modern taxonomic pattern of honey bee Apis mellifera is given in this paper. Thirty-three distinct honey bee subspecies are distributed across all Africa (11 subspecies), Western Asia and the Middle East (9 subspecies), and Europe (13 subspecies). All honey bee subspecies are subdivided into 5 evolutionary lineages: lineage A (10 subspecies) and its sublineage Z (3 subspecies), lineage M (3 subspecies), lineage C (10 subspecies), lineage O (3 subspecies), lineage Y (1 subspecies), lineage C or O (3 subspecies).

Mitochondrial genomes illuminate the evolutionary history of the Western honey bee (Apis mellifera)

Western honey bees (Apis mellifera) are one of the most important pollinators of agricultural crops and wild plants. Despite the growth in the availability of sequence data for honey bees, the phylogeny of the species remains a subject of controversy. Most notably, the geographic origin of honey bees is uncertain, as are the relationships among its constituent lineages and subspecies. We aim to infer the evolutionary and biogeographical history of the honey bee from mitochondrial genomes. Here we analyse the full mitochondrial genomes of 18 A. mellifera subspecies, belonging to all major lineages, using a range of gene sampling strategies and inference models to identify factors that may have contributed to the recovery of incongruent results in previous studies. Our analyses support a northern African or Middle Eastern origin of A. mellifera. We show that the previously suggested European and Afrotropical cradles of honey bees are the result of phylogenetic error. Monophyly of the M, C, and O lineages is strongly supported, but the A lineage appears paraphyletic. A. mellifera colonised Europe through at least two pathways, across the Strait of Gibraltar and via Asia Minor.

Genetic diversity and population structure of two subspecies of western honey bees (Apis mellifera L.) in the Republic of South Africa as revealed by microsatellite genotyping

Apis mellifera scutellata and Apis mellifera capensis, two native subspecies of western honey bees in the Republic of South Africa (RSA), are important to beekeepers in their native region because beekeepers use these bees for honey production and pollination purposes. Additionally, both bees are important invasive pests outside of their native ranges. Recently, whole mitogenome sequencing and single nucleotide polymorphisms were used to study their genetic diversity. To add to our knowledge of the molecular ecology of both bees, we tested the ability of microsatellites to be used as a tool to discriminate between A.m. capensis and A.m. scutellata. We analyzed the genetic variability and overall population structure of both bee subspecies and hybrids of the two by genotyping individuals collected from RSA (N = 813 bees from 75 apiaries) at 19 microsatellite DNA loci. Overall, populations averaged between 9.2 and 11.3 alleles per locus, with unbiased heterozygosity values ranging from 0.81 to 0.86 per population. Bayesian clustering analyses revealed two distinct evolutionary units, though the results did not match those of earlier morphometric and molecular analyses. This suggests that the microsatellites we tested were not sufficient for subspecies identification purposes, especially for Cape and hybrid bees. Nevertheless, the microsatellite data highlight the considerable genetic diversity within both populations and a larger-than-expected hybridization zone between the natural distributions of A.m. capensis and A.m. scutellata.

Genetic Differentiation of Eastern Honey Bee (Apis cerana) Populations Across Qinghai-Tibet Plateau-Valley Landforms

Many species of high-altitude plateaus tend to be narrowly distributed along river valleys at lower elevations due to a limitation of suitable habitats. The eastern honeybee (Apis cerana) is such a species and this study explored the effects of long and narrow geographic distributions on honeybee populations. Genetic differentiation and diversity were assessed across populations of the southeastern Qinghai-Tibet Plateau. A total of 492 honeybee samples from eight sampling sites in four valleys were analyzed for the genetic differentiation and diversity of 31 microsatellite loci and mitochondrial tRNA^leu-COII fragments. The following results were obtained: (1) Microsatellite genetic differentiation coefficients (F _ST) ranged from 0.06 to 0.16, and mitochondrial F _ST estimates ranged from 0.18 to 0.70 for different sampling sites in the same valley, indicating genetic differentiation. (2) Honeybees in adjacent valleys were also genetically differentiated. The F _ST of microsatellites and mitochondria were 0.04-0.29 and 0.06-0.76, respectively. (3) Likely a result of small population sizes, the observed genetic diversity was low. The observed impedance of honeybee gene flow among valleys increased both genetic differentiation and population numbers in the Qinghai-Tibet Plateau. This study contributes significantly to the current understanding of the mechanism underlying population genetic differentiation and highlights the potential effects of utilizing genetic resources that are subject to the ecological conditions of the long and narrow geographic distributions of plateau-valley landforms.