Mitochondrial DNA Diversity of Honey Bees (Apis mellifera) from Unmanaged Colonies and Swarms in the United States

To study the genetic diversity of honey bees (Apis mellifera L.) from unmanaged colonies in the United States, we sequenced a portion of the mitochondrial DNA COI–COII region. From the 530 to 1,230 bp amplicon, we observed 23 haplotypes from 247 samples collected from 12 states, representing three of the four A. mellifera lineages known to have been imported into the United States (C, M, and O). Six of the 13 C lineage haplotypes were not found in previous queen breeder studies in the United States. The O lineage accounted for 9% of unmanaged colonies which have not yet been reported in queen breeder studies. The M lineage accounted for a larger portion of unmanaged samples (7%) than queen breeder samples (3%). Based on our mitochondrial DNA data, the genetic diversity of unmanaged honey bees in the United States differs significantly from that of queen breeder populations (p < 0.00001). The detection of genetically distinct maternal lineages of unmanaged honey bees suggests that these haplotypes may have existed outside the managed honey bee population for a long period.     

Admixture in Africanized honey bees (Apis mellifera) from Panamá to San Diego,California (U.S.A.)

The Africanized honey bee (AHB) is a New World amalgamation of several subspecies of the western honey bee (Apis mellifera), a diverse taxon historically grouped into four major biogeographic lineages: A (African), M (Western European), C (Eastern European), and O (Middle Eastern). In 1956, accidental release of experimentally bred “Africanized” hybrids from a research apiary in Sao Paulo, Brazil initiated a hybrid species expansion that now extends from northern Argentina to northern California (U.S.A.). Here, we assess nuclear admixture and mitochondrial ancestry in 60 bees from four countries (Panamá; Costa Rica, Mexico; U.S.A) across this expansive range to assess ancestry of AHB several decades following initial introduction and test the prediction that African ancestry decreases with increasing latitude. We find that AHB nuclear genomes from Central America and Mexico have predominately African genomes (76%–89%) with smaller contributions from Western and Eastern European lineages. Similarly, nearly all honey bees from Central America and Mexico possess mitochondrial ancestry from the African lineage with few individuals having European mitochondria. In contrast, AHB from San Diego (CA) shows markedly lower African ancestry (38%) with substantial genomic contributions from all four major honey bee lineages and mitochondrial ancestry from all four clades as well. Genetic diversity measures from all New World populations equal or exceed those of ancestral populations. Interestingly, the feral honey bee population of San Diego emerges as a reservoir of diverse admixture and high genetic diversity, making it a potentially rich source of genetic material for honey bee breeding.     

Introgression in native populations of <Emphasis Type="Italic">Apis mellifera mellifera</Emphasis> L: implications for conservation

Hybridisation and introgression can have negative impacts on regional biodiversity through the potential erosion of locally adapted lineages. The honey bee (Apis mellifera L.) occurs in twenty-seven subspecies across Europe, is an extremely economically important insect, yet threatened by multifarious impacts. Transhumance of the most commercially appealing varieties threatens native honey bee diversity by introgression and subsequent loss of locally adapted traits, or even by complete removal of some subspecies from parts of the range. Here levels of admixture and introgression are examined in UK honey bees suspected to be from hives of the dark European honey bee (Apis mellifera mellifera). Microsatellite DNA and STRUCTURE analyses reveal that the studied populations are generally admixed, and discriminant analysis of principal components shows them to be intermediate between A. m. mellifera and Apis mellifera carnica populations. However, examining mitochondrial haplotype data (COI-COII intergenic spacer region) and nuclear DNA reveal that some hives are relatively pure (from 4 to 15 hives, depending on the Q-value threshold). Genetic diversity is relatively high in comparison with other European populations. Implications for conservation and management are discussed.     

The Honey Bee Pathosphere of Mongolia: European Viruses in Central Asia

Parasites and pathogens are apparent key factors for the detrimental health of managed European honey bee subspecies, Apis mellifera. Apicultural trade is arguably the main factor for the almost global distribution of most honey bee diseases, thereby increasing chances for multiple infestations/infections of regions, apiaries, colonies and even individual bees. This imposes difficulties to evaluate the effects of pathogens in isolation, thereby creating demand to survey remote areas. Here, we conducted the first comprehensive survey for 14 honey bee pathogens in Mongolia (N = 3 regions, N = 9 locations, N = 151 colonies), where honey bee colonies depend on humans to overwinter. In Mongolia, honey bees, Apis spp., are not native and colonies of European A. mellifera subspecies have been introduced ~60 years ago. Despite the high detection power and large sample size across Mongolian regions with beekeeping, the mite Acarapis woodi, the bacteria Melissococcus plutonius and Paenibacillus larvae, the microsporidian Nosema apis, Acute bee paralysis virus, Kashmir bee virus, Israeli acute paralysis virus and Lake Sinai virus strain 2 were not detected, suggesting that they are either very rare or absent. The mite Varroa destructor, Nosema ceranae and four viruses (Sacbrood virus, Black queen cell virus, Deformed wing virus (DWV) and Chronic bee paralysis virus) were found with different prevalence. Despite the positive correlation between the prevalence of V. destructor mites and DWV, some areas had only mites, but not DWV, which is most likely due to the exceptional isolation of apiaries (up to 600 km). Phylogenetic analyses of the detected viruses reveal their clustering and European origin, thereby supporting the role of trade for pathogen spread and the isolation of Mongolia from South-Asian countries. In conclusion, this survey reveals the distinctive honey bee pathosphere of Mongolia, which offers opportunities for exciting future research.     

Colonization history and population differentiation of the Honey Bees (Apis mellifera L.) in Puerto Rico

Honey bees (Apis mellifera L.) are the primary commercial pollinators across the world. The subspecies A. m. scutellata originated in Africa and was introduced to the Americas in 1956. For the last 60 years, it hybridized successfully with European subspecies, previous residents in the area. The result of this hybridization was called Africanized honey bee (AHB). AHB has spread since then, arriving to Puerto Rico (PR) in 1994. The honey bee population on the island acquired a mosaic of features from AHB or the European honey bee (EHB). AHB in Puerto Rico shows a major distinctive characteristic, docile behavior, and is called gentle Africanized honey bees (gAHB). We used 917 SNPs to examine the population structure, genetic differentiation, origin, and history of range expansion and colonization of gAHB in PR. We compared gAHB to populations that span the current distribution of A. mellifera worldwide. The gAHB population is shown to be a single population that differs genetically from the examined populations of AHB. Texas and PR groups are the closest genetically. Our results support the hypothesis that the Texas AHB population is the source of gAHB in Puerto Rico.     

Population structure of honey bees in the Carpathian Basin (Hungary) confirms introgression from surrounding subspecies

Carniolan honey bees (Apis mellifera carnica) are considered as an indigenous subspecies in Hungary adapted to most of the ecological and climatic conditions in this area. However, during the last decades Hungarian beekeepers have recognized morphological signs of the Italian honey bee (Apis mellifera ligustica). As the natural distribution of the honey bee subspecies can be affected by the importation of honey bee queens or by natural gene flow, we aimed at determining the genetic structure and characteristics of the local honey bee population using molecular markers. All together, 48 Hungarian and 84 foreign (Italian, Polish, Spanish, Liberian) pupae and/or workers were used for mitochondrial DNA analysis. Additionally, 53 sequences corresponding to 10 subspecies and the Buckfast hybrid were downloaded from GenBank. For the nuclear analysis, 236 Hungarian and 106 foreign honey bees were genotyped using nine microsatellites. Heterozygosity values, population‐specific alleles, FST values, principal coordinate analysis, assignment tests, structure analysis, and dendrograms were calculated. Haplotype and nucleotide diversity values showed moderate values. We found that one haplotype (H9) was dominant in Hungary. The presence of the black honey bee (Apis mellifera mellifera) was negligible, but a few individuals resembling other subspecies were identified. We proved that the Hungarian honey bee population is nearly homogeneous but also demonstrated introgression from the foreign subspecies. Both mitochondrial DNA and microsatellite analyses corroborated the observations of the beekeepers. Molecular analyses suggested that Carniolan honey bee in Hungary is slightly affected by Italian and black honey bee introgression. Genetic differences were detected between Polish and Hungarian Carniolan honey bee populations, suggesting the existence of at least two different gene pools within A. m. carnica.     

A Cell Line Resource Derived from Honey Bee (Apis mellifera) Embryonic Tissues

A major hindrance to the study of honey bee pathogens or the effects of pesticides and nutritional deficiencies is the lack of controlled in vitro culture systems comprised of honey bee cells. Such systems are important to determine the impact of these stress factors on the developmental and cell biology of honey bees. We have developed a method incorporating established insect cell culture techniques that supports sustained growth of honey bee cells in vitro. We used honey bee eggs mid to late in their embryogenesis to establish primary cultures, as these eggs contain cells that are progressively dividing. Primary cultures were initiated in modified Leibovitz’s L15 medium and incubated at 32^°C. Serial transfer of material from several primary cultures was maintained and has led to the isolation of young cell lines. A cell line (AmE-711) has been established that is composed mainly of fibroblast-type cells that form an adherent monolayer. Most cells in the line are diploid (2n = 32) and have the Apis mellifera karyotype as revealed by Giemsa stain. The partial sequence for the mitochondrial-encoded cytochrome c oxidase subunit I (Cox 1) gene in the cell line is identical to those from honey bee tissues and a consensus sequence for A. mellifera. The population doubling time is approximately 4 days. Importantly, the cell line is continuously subcultured every 10–14 days when split at a 1:3 ratio and is cryopreserved in liquid nitrogen. The cell culture system we have developed has potential application for studies aimed at honey bee development, genetics, pathogenesis, transgenesis, and toxicology.     

Insight into probiotic properties of lactic acid bacterial endosymbionts of Apis mellifera L. derived from the Polish apiary

Taking into account that fructophilic lactic acid bacteria (FLAB) can play an important role in the health of honey bees and can be used as probiotics, phenotypic properties of probiotic interest of Lactobacillus kunkeei (12 strains) and Fructobacillus fructossus bacteria (2 strains), isolated from Apis mellifera gastrointestinal tract, have been studied. We have evaluated survival of tested FLAB in honey bee gut, their susceptibility to antibiotics (ampicillin, erythromycin, tylosin), cell surface hydrophobicity, auto-aggregation ability, co-aggregation with model pathogenic bacteria, biofilm formation capacity, and effect of studied FLAB, added to sucrose syrup bee diet, on longevity of honey bees. The tested FLAB exhibited good gastrointestinal tract tolerance and high antibiotic susceptibility, which are important criteria in the screening of probiotic candidates. It was also found that all FLAB studied have high cell surface hydrophobicity and fulfil next selection criterion for their use as probiotics. Symbionts of A. mellifera showed also auto- and co-aggregation capacities regarded as valuable features for biofilm formation and inhibition of pathogens adhesion to the bee gut cells. Biofilm-development ability is a desired characteristic of probiotic lactic acid bacteria. As indicated by quantitative crystal violet-stained microplate assay and confocal laser scanning microscopy imaging, all studied A. mellifera gut isolates exhibit a biofilm positive phenotype. Moreover, it was also documented, on honey bees kept in cages, that supplementation of A. mellifera sucrose diet with FLAB decreases mortality and improves significantly longevity of honey bees. Presented research showed that A. mellifera FLAB symbionts are good candidates for application as probiotics.     

Reduced SNP Panels for Genetic Identification and Introgression Analysis in the Dark Honey Bee (Apis mellifera mellifera)

Beekeeping activities, especially queen trading, have shaped the distribution of honey bee (Apis mellifera) subspecies in Europe, and have resulted in extensive introductions of two eastern European C-lineage subspecies (A. m. ligustica and A. m. carnica) into the native range of the M-lineage A. m. mellifera subspecies in Western Europe. As a consequence, replacement and gene flow between native and commercial populations have occurred at varying levels across western European populations. Genetic identification and introgression analysis using molecular markers is an important tool for management and conservation of honey bee subspecies. Previous studies have monitored introgression by using microsatellite, PCR-RFLP markers and most recently, high density assays using single nucleotide polymorphism (SNP) markers. While the latter are almost prohibitively expensive, the information gained to date can be exploited to create a reduced panel containing the most ancestry-informative markers (AIMs) for those purposes with very little loss of information. The objective of this study was to design reduced panels of AIMs to verify the origin of A. m. mellifera individuals and to provide accurate estimates of the level of C-lineage introgression into their genome. The discriminant power of the SNPs using a variety of metrics and approaches including the Weir & Cockerham’s F_ST, an F_ST-based outlier test, Delta, informativeness (I_n), and PCA was evaluated. This study shows that reduced AIMs panels assign individuals to the correct origin and calculates the admixture level with a high degree of accuracy. These panels provide an essential tool in Europe for genetic stock identification and estimation of admixture levels which can assist management strategies and monitor honey bee conservation programs.     

Admixture increases diversity in managed honey bees: Reply to De la Rúa et al. (2013)

De la Rúa et al. (2013) express some concerns about the conclusions of our recent study showing that management increases genetic diversity of honey bees (Apis mellifera) by promoting admixture (Harpur et al. 2012). We provide a brief review of the literature on the population genetics of A. mellifera and show that we utilized appropriate sampling methods to estimate genetic diversity in the focal populations. Our finding of higher genetic diversity in two managed A. mellifera populations on two different continents is expected to be the norm given the large number of studies documenting admixture in honey bees. Our study focused on elucidating how management affects genetic diversity in honey bees, not on how to best manage bee colonies. We do not endorse the intentional admixture of honey bee populations, and we agree with De la Rúa et al. (2013) that native honey bee subspecies should be conserved.     

Thinking beyond Western commercial honeybee hives: towards improved conservation of honey bee diversity

A decline of wild pollinators, along with a decline of bee diversity, has been a cause of concern among academics and governmental organizations. According to IPBES, a lack of wild pollinator data contributes to difficulties in comprehensively analyzing the regional status of wild pollinators in Africa, Latin America, Asia and Oceania. It may have also contributed to the prevailing lack of awareness of the diversity of honey bees, of which the managed Apis mellifera is often considered as “the (only) honey bee,” despite the fact that there are eight other honey bee species extant in Asia. A survey of 100 journal articles published in 2016 shows that 57% of the studies still identified A. mellifera as “the honey bee.” In total, 80% of studies were conducted solely on A. mellifera. This focus on A. mellifera has also caused the honey standard of Codex Alimentarius and the European Union to be based solely on A. mellifera, causing improper evaluation of honeys from other species. We recommend adapting current standards to reflect the diversity of honey bees and in the process correct failures in the honey market and pave the way towards improved protection of honey bee species and their habitats.     

Genome-wide characterization of long intergenic non-coding RNAs (lincRNAs) provides new insight into viral diseases in honey bees Apis cerana and Apis mellifera

Background

Long non-coding RNAs (lncRNAs) are a class of RNAs that do not encode proteins. Recently, lncRNAs have gained special attention for their roles in various biological process and diseases.

Results

In an attempt to identify long intergenic non-coding RNAs (lincRNAs) and their possible involvement in honey bee development and diseases, we analyzed RNA-seq datasets generated from Asian honey bee (Apis cerana) and western honey bee (Apis mellifera). We identified 2470 lincRNAs with an average length of 1011 bp from A. cerana and 1514 lincRNAs with an average length of 790 bp in A. mellifera. Comparative analysis revealed that 5 % of the total lincRNAs derived from both species are unique in each species. Our comparative digital gene expression analysis revealed a high degree of tissue-specific expression among the seven major tissues of honey bee, different from mRNA expression patterns. A total of 863 (57 %) and 464 (18 %) lincRNAs showed tissue-dependent expression in A. mellifera and A. cerana, respectively, most preferentially in ovary and fat body tissues. Importantly, we identified 11 lincRNAs that are specifically regulated upon viral infection in honey bees, and 10 of them appear to play roles during infection with various viruses.

Conclusions

This study provides the first comprehensive set of lincRNAs for honey bees and opens the door to discover lincRNAs associated with biological and hormone signaling pathways as well as various diseases of honey bee.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1868-7) contains supplementary material, which is available to authorized users.     

Pyrosequencing analysis of the bacterial communities in the guts of honey bees Apis cerana and Apis mellifera in Korea

The bacterial communities in the guts of the adults and larvae of the Asian honey bee Apis cerana and the European honey bee Apis mellifera were surveyed by pyrosequencing the 16S rRNA genes. Most of the gut bacterial 16S rRNA gene sequences were highly similar to the known honey bee-specific ones and affiliated with Pasteurellaceae or lactic acid bacteria (LAB). The numbers of operational taxonomic units (OTUs, defined at 97% similarity) were lower in the larval guts (6 or 9) than in the adult guts (18 or 20), and the frequencies of Pasteurellaceae-related OTUs were higher in the larval guts while those of LAB-related OTUs in the adult guts. The frequencies of Lactococcus, Bartonella, Spiroplasma, Enterobacteriaceae, and Flavobacteriaceae-related OTUs were much higher in A. cerana guts while Bifidobacterium and Lachnospiraceae-related OTUs were more abundant in A. mellfera guts. The bacterial community structures in the midguts and hindguts of the adult honey bees were not different for A. cerana, but significantly different for A. mellifera. The above results substantiated the previous observation that honey bee guts are dominated by several specific bacterial groups, and also showed that the relative abundances of OTUs could be markedly changed depending on the developmental stage, the location within the gut, and the honey bee species. The possibility of using the gut bacterial community as an indicator of honey bee health was discussed.     

Structural and functional differences in the antennal olfactory system of worker honey bees of Apis mellifera and Apis cerana

Olfactory cues are important sensory modalities on individual discrimination, perception, and efficient orientation to food sources in most insects. In honey bees, which are well known as eusocial insects, olfactory cues are mainly used to maintain a colony. Although much research has been reported on olfactory systems in honey bee olfaction, little is known about the differences between two major honey bee species, the European honey bee Apis mellifera and the Asian honey bee Apis cerana. In order to understand the differences of olfactory characteristics in the two species, we compared the distribution of sensory hairs on the antennae and antennal olfactory responses, using electron microscopy, electrophysiological recording and molecular expression level of odorant receptors. Our present study demonstrated that the antennae of A. cerana have more olfactory sensilla than A. mellifera, responding more strongly to various floral volatile compounds. At the molecular level, olfactory co-receptor (Orco), which makes heterodimers with other conventional olfactory receptors, is more abundantly expressed in the antenna of A. cerana than of A. mellifera. These findings extend our understanding of the olfactory systems and behavioral responses to various ecological and biological signals in two closely related honey bee species.     

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

ObjectivesThis 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.MethodsThe 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 conclusionsThe 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.     

Infestation of Japanese Native Honey Bees by Tracheal Mite and Virus from Non-native European Honey Bees in Japan

Invasion of alien species has been shown to cause detrimental effects on habitats of native species. Insect pollinators represent such examples; the introduction of commercial bumble bee species for crop pollination has resulted in competition for an ecological niche with native species, genetic disturbance caused by mating with native species, and pathogen spillover to native species. The European honey bee, Apis mellifera, was first introduced into Japan for apiculture in 1877, and queen bees have been imported from several countries for many years. However, its effects on Japanese native honey bee, Apis cerana japonica, have never been addressed. We thus conducted the survey of honey bee viruses and Acarapis mites using both A. mellifera and A. c. japonica colonies to examine their infestation in native and non-native honey bee species in Japan. Honey bee viruses, Deformed wing virus (DWV), Black queen cell virus (BQCV), Israeli acute paralysis virus (IAPV), and Sacbrood virus (SBV), were found in both A. mellifera and A. c. japonica colonies; however, the infection frequency of viruses in A. c. japonica was lower than that in A. mellifera colonies. Based on the phylogenies of DWV, BQCV, and SBV isolates from A. mellifera and A. c. japonica, DWV and BQCV may infect both honey bee species; meanwhile, SBV has a clear species barrier. For the first time in Japan, tracheal mite (Acarapis woodi) was specifically found in the dead honey bees from collapsing A. c. japonica colonies. This paper thus provides further evidence that tracheal-mite-infested honey bee colonies can die during cool winters with no other disease present. These results demonstrate the infestation of native honey bees by parasite and pathogens of non-native honey bees that are traded globally.     

csd alleles in the red dwarf honey bee (Apis florea,Hymenoptera: Apidae) show exceptionally high nucleotide diversity

Abstract The single locus complementary sex determination (sl‐csd) gene is the primary gene determining the gender of honey bees (Apis spp.). While the csd gene has been well studied in the Western honey bee (Apis mellifera), and comparable data exist in both the Eastern honey bee (Apis cerana) and the giant honey bee (Apis dorsata), no studies have been conducted in the red dwarf honey bee, Apis florea. In this study we cloned the genomic region 3 of the A. florea csd gene from 60 workers, and identified 12 csd alleles. Analysis showed that similar to A. mellifera, region 3 of the csd gene contains a RS domain at the N terminal, a proline‐rich domain at the C terminal, and a hypervariable region in the middle. However, the A. florea csd gene possessed a much higher level of nucleotide diversity, compared to A. mellifera, A. cerana and Apis dorsata. We also show that similar to the other three Apis species, in A. florea, nonsynonymous mutations in the csd gene are selectively favored in young alleles.     

Olfactory responses of Aethina tumida murray (Coleoptera: Nitidulidae) to some major volatile compounds from hive materials and workers of Apis mellifera

The small hive beetle (Aethina tumida Murray) is an invasive pest affecting honey bee colonies. The beetles are known to be attracted to volatiles from hive products and honey bees like Apis mellifera L. Previously we reported the presence of five major compounds from the volatile extracts of hive materials; ethyl linolenate and ethyl palmitate from pollen dough, oleamide and tetracosane in fermenting honey, and oleamide and 5-methyl-2-phenyl-1H-indole from A. mellifera worker bees. This study tested the attractiveness of the aforementioned five volatile organic compounds to small hive beetles (SHB) by Y-tube olfactometric bioassay. Ethyl linolenate was highly attractive to both male and female adults of SHB. Ethyl palmitate was attractive to SHB only at higher concentration (0.01–01 mg/ml). Interestingly, tetracosane, 5-methyl-2-phenyl-1H-indole and oleamide were repellent for SHB of both sexes, but ethyl linolenate and ethyl palmitate as components of honey bee brood pheromone attracted SHB. The results highlight that SHB differentially utilizes volatile chemicals from hive materials and honey bees as cues to locate honey bee hives.     

PCR-RFLP of mitochondrial DNA reveals two origins of Apis mellifera in Taiwan

Beekeeping has been a highly valued industry in Taiwan. As a result, many subspecies of Apis mellifera have been introduced to Taiwan since 1911, leading to the hybridization of different subspecies. In order to know the matrilineal origins of Taiwan A. mellifera, a total of 280 samples collected from 33 apiaries throughout the island were examined. Using PCR-RFLP of four mitochondrial gene fragments, i.e., the non-coding region between tRNA^leu and cytochrome c oxidase subunit II (intergenic tRNA^leu-COII), cytochrome b (Cyt b), large subunit rRNA (Ls rRNA) and cytochrome c oxidase subunit I (COI), we only found two haplotypes exist in 280 samples. Haplotypes ababa and bbbaa account for 87% of these Western bees belonged to the Eastern European (C) lineage and 13% belonged to the Middle East (Z) lineage, respectively, with the latter being totally absent in northern Taiwan. African (A) and Mellifera (M) lineages, officially imported once in 1990s and 1930s respectively, were not detected. The identification of subspecies of A. mellifera and survey of their distribution on the island are expected to facilitate efficient breeding programs and establish a more booming beekeeping industry.