Polymorphism of locus COI-COII of mitochondrial DNA in the honeybee Apis mellifera L. from southern Ural region

In spite of high biodiversity within the honeybee species Apis mellifera L., only one geographical race, the dark-colored forest honeybee A. m. mellifera, is uniquely adapted to severe environmental conditions of Eurasian forest and forest-steppe zones. Within the vast range of this race, only single isolates remain, where the dark forest honeybee is purebred. The Bashkir population is supposed to be one of these isolates. Molecular-genetic assessment of the state of the gene pool of this population revealed that southern honeybee races were introduced into the Bashkortostan Republic with great intensity, which was above the assimilation capacity of the population. The main part of the former range of the Bashkir population represents a hybrid zone with approximately equal ratio between gene pools of local and introduced honeybees. Our studies provide the possibility to single out one extant reserve of A. m. mellifera, Burzyanskii raion, in which the proportion of local bees in the gene pool is 0.98.     

Molecular genetic analysis of five extant reserves of black honeybee <Emphasis Type="Italic">Apis melifera melifera</Emphasis> in the Urals and the Volga region

Local populations of the black honeybee Apis mellifera mellifera from the Urals and the Volga region were examined in comparison with local populations of southern honeybee subspecies A. m. caucasica and A. m. carpatica from the Caucasus and the Carpathians. Genetic analysis was performed on the basis of the polymorphism of nine microsatellite loci of nuclear DNA and the mtDNA COI–COII locus. On the territory of the Urals and the Volga region, five extant populations (reserves) of the black honeybee A. m. mellifera were identified, including the Burzyanskaya, Tatyshlinskaya, Yuzhno-Prikamskaya, Visherskaya, and Kambarskaya populations. These five populations are the basis of the modern gene pool of the black honeybee A. m. mellifera from the Urals and the Volga region. The greatest proportion of the remaining indigenous gene pool of A. m. mellifera (the core of the gene pool of the population of A. m. mellifera) is distributed over the entire territory of Perm krai and the north of the Republic of Bashkortostan. For the population of A. m. mellifera from the Urals and the Volga region, the genetic standards were calculated, which will be useful for future population studies of honeybees.     

New approach to the mitotype classification in black honeybee <Emphasis Type="Italic">Apis mellifera mellifera</Emphasis> and Iberian honeybee <Emphasis Type="Italic">Apis mellifera iberiensis</Emphasis>

The black honeybee Apis mellifera mellifera L. is today the only subspecies of honeybee which is suitable for commercial breeding in the climatic conditions of Northern Europe with long cold winters. The main problem of the black honeybee in Russia and European countries is the preservation of the indigenous gene pool purity, which is lost as a result of hybridization with subspecies, A. m. caucasica, A. m. carnica, A. m. carpatica, and A. m. armeniaca, introduced from southern regions. Genetic identification of the subspecies will reduce the extent of hybridization and provide the gene pool conservation of the black honeybee. Modern classification of the honeybee mitotypes is mainly based on the combined use of the DraI restriction endonuclease recognition site polymorphism and sequence polymorphism of the mtDNA COI–COII region. We performed a comparative analysis of the mtDNA COI–COII region sequence polymorphism in the honeybees of the evolutionary lineage M from Ural and West European populations of black honeybee A. m. mellifera and Spanish bee A. m. iberiensis. A new approach to the classification of the honeybee M mitotypes was suggested. Using this approach and on the basis of the seven most informative SNPs of the mtDNA COI–COII region, eight honeybee mitotype groups were identified. In addition, it is suggested that this approach will simplify the previously proposed complicated mitotype classification and will make it possible to assess the level of the mitotype diversity and to identify the mitotypes that are the most valuable for the honeybee breeding and rearing.     

Association of Varroa destructor females in multiply infested cells of the honeybee Apis mellifera

The genetic diversity of Varroa destructor (Anderson &Trueman)is limited outside its natural range due to population bottlenecks and its propensity to inbreed.In light of the arms race between V.destructor and its honeybee (Apis mellifera L.)host, any mechanism enhancing population admixture of the mite may be favored.One way that admixture can occur is when two genetically dissimilar mites coinvade a brood cell, with the progeny of the foundresses admixing.We determined the relatedness of 393 pairs of V.destructor foundresses,each pair collected from a single bee brood cell (n =five colonies).We used six microsatellites to identify the genotypes of mites coinvading a cell and calculated the frequency of pairs with different or the same genotypes.We found no deviation from random coinvasion,but the frequency of cells infested by mites with different genotypes was high.This rate of recombination,coupled with a high transmission rate of mites,homogenized the allelic pool of mites within the apiary.     

Local honeybee (<Emphasis Type="Italic">Apis mellifera mellifera</Emphasis> L.) populations in the Urals

The COI-COII intergenic region of mitochondrial DNA (mtDNA) was studied in local honeybee (Apis mellifera mellifera) L. populations from the Middle and Southern Urals. Analysis of bee colonies in these regions revealed apiaries enriched in families descending from A. m. mellifera in the maternal lineage. These results confirm the suggestion of preservation of A. m. mellifera refuges in the Urals and provide grounds for work on the preservation of the gene pool of this bee variety, valuable for all Russia.     

Gene flow in admixed populations and implications for the conservation of the Western honeybee, Apis mellifera

Anthropogenic activity, especially modern apiculture, has considerable impact on the natural distribution of the Western honeybee, Apis mellifera, leading to the spread, replacement and fragmentation of many subspecies. This creates demand for the conservation of some subspecies, in particular, Apis mellifera mellifera, which once was widely distributed in Western Europe and nowadays is endangered through habitat loss and fragmentation. Moreover, A. m. mellifera may be further endangered by hybridisation in populations that now occur in artificial sympatry with other subspecies. Here, we quantify and compare individual hybridisation between sympatric and allopatric honeybee populations of A. m. mellifera and A. m. carnica using microsatellite markers and a Bayesian model-based approach. We had a special focus on pure breeding populations, which are a major tool in honeybee conservation. Our results demonstrate that subspecies are still highly differentiated, but gene flow is not prevented by the current management strategies, creating urgent demand for an improved conservation management of A. m. mellifera. However, the occurrence of a high number of pure individuals might suggest that some sort of hybrid barrier acts against the complete admixture of the two subspecies.     

Population structure of common juniper in the Cis-Urals and Southern Urals

The population structure of common juniper (Juniperus communis L.) growing in the Cis-Urals and Southern Urals has been studied using 17 morphological traits of generative and vegetative organs. A multivariate analysis of ten coenopopulations has recognized three phenotypically different local populations: Cis-Ural, forest Cis-Ural forest-steppe, and Southern Ural mountain populations. The Cis-Ural forest population is strictly associated with lowland pine forests of the northwestern part of the Bashkir Cis-Urals. The Cis-Ural forest-steppe population is located in the northwestern part of the Bashkir Cis-Urals and the southeastern part of the Udmurt Cis-Urals. The Southern Ural mountain population is located in the central part of the Southern Urals and is associated mainly with mountain pine and dark coniferous forests. The last population is divided into forest and forest-edge subpopulations; the first one is represented by typical undergrowth locations, whereas the second is associated with open steppelike slopes and forest edges. In general, based on morphological traits of generative organs, the revealed local subpopulations hold an intermediate position between the Eastern European and Siberian populations of common juniper. Based on the morphological traits of vegetative organs, Cis-Ural populations are considered related to the populations of the European part of Russia, whereas the mountain Southern Ural population resemble Siberian populations. Concerning morphological traits of generative organs, the intrapopulation phenotypic diversity of common juniper is higher for mountain habitats; in the case of vegetative organs, the maximum diversity is observed for lowland habitats. The character of phenotypic differentiation determines the need to conserve the gene pool of common juniper of the Cis-Urals and southern Urals on a population basis.     

SNPs selected by information content outperform randomly selected microsatellite loci for delineating genetic identification and introgression in the endangered dark European honeybee (Apis mellifera mellifera)

The honeybee (Apis mellifera) has been threatened by multiple factors including pests and pathogens, pesticides and loss of locally adapted gene complexes due to replacement and introgression. In western Europe, the genetic integrity of the native A. m. mellifera (M‐lineage) is endangered due to trading and intensive queen breeding with commercial subspecies of eastern European ancestry (C‐lineage). Effective conservation actions require reliable molecular tools to identify pure‐bred A. m. mellifera colonies. Microsatellites have been preferred for identification of A. m. mellifera stocks across conservation centres. However, owing to high throughput, easy transferability between laboratories and low genotyping error, SNPs promise to become popular. Here, we compared the resolving power of a widely utilized microsatellite set to detect structure and introgression with that of different sets that combine a variable number of SNPs selected for their information content and genomic proximity to the microsatellite loci. Contrary to every SNP data set, microsatellites did not discriminate between the two lineages in the PCA space. Mean introgression proportions were identical across the two marker types, although at the individual level, microsatellites' performance was relatively poor at the upper range of Q‐values, a result reflected by their lower precision. Our results suggest that SNPs are more accurate and powerful than microsatellites for identification of A. m. mellifera colonies, especially when they are selected by information content.     

Multiple host shifts by the emerging honeybee parasite,Varroa jacobsoni

Host shifts are a key mechanism of parasite evolution and responsible for the emergence of many economically important pathogens. Varroa destructor has been a major factor in global honeybee (Apis mellifera) declines since shifting hosts from the Asian honeybee (Apis cerana) > 50 years ago. Until recently, only two haplotypes of V. destructor (Korea and Japan) had successfully host shifted to A. mellifera. In 2008, the sister species V. jacobsoni was found for the first time parasitizing A. mellifera in Papua New Guinea (PNG). This recent host shift presents a serious threat to world apiculture but also provides the opportunity to examine host shifting in this system. We used 12 microsatellites to compare genetic variation of V. jacobsoni on A. mellifera in PNG with mites on A. cerana in both PNG and surrounding regions. We identified two distinct lineages of V. jacobsoni reproducing on A. mellifera in PNG. Our analysis indicated independent host shift events have occurred through small numbers of mites shifting from local A. cerana populations. Additional lineages were found in the neighbouring Papua and Solomon Islands that had partially host shifted to A. mellifera, that is producing immature offspring on drone brood only. These mites were likely in transition to full colonization of A. mellifera. Significant population structure between mites on the different hosts suggested host shifted V. jacobsoni populations may not still reproduce on A. cerana, although limited gene flow may exist. Our studies provide further insight into parasite host shift evolution and help characterize this new Varroa mite threat to A. mellifera worldwide.     

The role of epistatic interactions underpinning resistance to parasitic Varroa mites in haploid honey bee (Apis mellifera) drones

The Red Queen hypothesis predicts that host–parasite coevolutionary dynamics can select for host resistance through increased genetic diversity, recombination and evolutionary rates. However, in haplodiploid organisms such as the honeybee (Apis mellifera), models suggest the selective pressure is weaker than in diploids. Haplodiploid sex determination, found in A. mellifera, can allow deleterious recessive alleles to persist in the population through the diploid sex with negative effects predominantly expressed in the haploid sex. To overcome these negative effects in haploid genomes, epistatic interactions have been hypothesized to play an important role. Here, we use the interaction between A. mellifera and the parasitic mite Varroa destructor to test epistasis in the expression of resistance, through the inhibition of parasite reproduction, in haploid drones. We find novel loci on three chromosomes which explain over 45% of the resistance phenotype. Two of these loci interact only additively, suggesting their expression is independent of each other, but both loci interact epistatically with the third locus. With drone offspring inheriting only one copy of the queen's chromosomes, the drones will only possess one of two queen alleles throughout the years‐long lifetime of the honeybee colony. Varroa, in comparison, completes its highly inbred reproductive cycle in a matter of weeks, allowing it to rapidly evolve resistance. Faced with the rapidly evolving Varroa, a diversity of pathways and epistatic interactions for the inhibition of Varroa reproduction could therefore provide a selective advantage to the high levels of recombination seen in A. mellifera. This allows for the remixing of phenotypes despite a fixed queen genotype.     

Discrimination of Iranian honeybee populations (Apis mellifera meda) from commercial subspecies of Apis mellifera L. using morphometric and genetic methods

The morphological characters of honeybees have an important role for discriminating honeybee subspecies. In the present research, Iranian populations of honeybee (Apis mellifera) were collected from 19 areas in Iran. The samples were collected from stationary beekeeping sites. Moreover, pictures of honeybee forewings held in the Bee Data Bank in Oberursel were compared with Iranian honeybee populations. 19 morphological characters were measured for each forewing of worker honeybee to evaluate differentiation of Iranian honeybee populations from the commercial honeybee subspecies A. m. mellifera, A. m. carnica, A. m. caucasica and A. m. ligustica. Additionally, part of the tRNA^leu gene, an intergenic region and part of COII was used to confirm differentiation of the commercial subspecies from Iranian honeybee populations. Results of the cluster analyses showed that 19 morphological characters of forewings differentiated Iranian populations from the commercial subspecies. Moreover, the phylogenetic tree of part of the tRNA^leu gene, an intergenic region and part of COII differentiated the commercial subspecies from Iranian honeybee populations. Results of the discriminant function analyses (DFA) indicated that the references samples of A. m. meda overlapped with Iranian populations.     

Fine spatial structure of allozyme genotypes in isolated population of pedunculate oak <Emphasis Type="Italic">Quercus robur</Emphasis> L. (Fagaceae)

The extent and spatial pattern of genetic variation at polymorphic allozyme loci in a population of pedunculate oak Quercus robur from the Bashkir Transural region was investigated using autocorrelation analysis. In the plantation examined, statistically significant local concentration of most of the alleles in two-dimensional space was identified. The measures for protection of this small population located outside of the western border of the species range, in the mountain-steppe habitat, and characterized by specific gene pool, are suggested.     

The pheromones of laying workers in two honeybee sister species: <Emphasis Type="Italic">Apis cerana</Emphasis> and <Emphasis Type="Italic">Apis mellifera</Emphasis>

When a honeybee colony loses its queen, workers activate their ovaries and begin to lay eggs. This is accompanied by a shift in their pheromonal bouquet, which becomes more queen like. Workers of the Asian hive bee Apis cerana show unusually high levels of ovary activation and this can be interpreted as evidence for a recent evolutionary arms race between queens and workers over worker reproduction in this species. To further explore this, we compared the rate of pheromonal bouquet change between two honeybee sister species of Apis cerana and Apis mellifera under queenright and queenless conditions. We show that in both species, the pheromonal components HOB, 9-ODA, HVA, 9-HDA, 10-HDAA and 10-HDA have significantly higher amounts in laying workers than in non-laying workers. In the queenright colonies of A. mellifera and A. cerana, the ratios (9-ODA)/(9-ODA + 9-HDA + 10-HDAA + 10-HDA) are not significantly different between the two species, but in queenless A. cerana colonies the ratio is significant higher than in A. mellifera, suggesting that in A. cerana, the workers’ pheromonal bouquet is dominated by the queen compound, 9-ODA. The amount of 9-ODA in laying A. cerana workers increased by over 585% compared with the non-laying workers, that is 6.75 times higher than in A. mellifera where laying workers only had 86% more 9-ODA compared with non-laying workers.     

The impact of apiculture on the genetic structure of wild honeybee populations (<Emphasis Type="Italic">Apis mellifera</Emphasis>) in Sudan

Apiculture often relies on the importation of non-native honeybees (Apis mellifera) and large distance migratory beekeeping. These activities can cause biodiversity conflicts with the conservation of wild endemic honeybee subspecies. We studied the impact of large scale honeybee imports on managed and wild honeybee populations in Sudan, a centre of biodiversity of A. mellifera, using as set of linked microsatellite DNA loci and mitochondrial DNA markers. The mitochondrial DNA analyses showed that all wild honey bees exclusively belonged to African haplotypes. However, we revealed non-native haplotypes in managed colonies on apiaries reflecting unambiguous evidence of imports from European stock. Moreover, we found significantly higher linkage disequilibria for microsatellite markers in wild populations in regions with apiculture compared to wild populations which had no contact to beekeeping. Introgression of imported honeybees was only measurable at the population level in close vicinity to apicultural activities but not in wild populations which represent the vast majority of honeybees in Sudan.     

Conserving genetic diversity in the honeybee: Comments on Harpur et al. (2012)

The article by Harpur et al. (2012) ‘Management increases genetic diversity of honey bees via admixture’ concludes that ‘…honey bees do not suffer from reduced genetic diversity caused by management and, consequently, that reduced genetic diversity is probably not contributing to declines of managed Apis mellifera populations’. In the light of current honeybee and beekeeping declines and their consequences for honeybee conservation and the pollination services they provide, we would like to express our concern about the conclusions drawn from the results of Harpur et al. (2012). While many honeybee management practices do not imply admixture, we are convinced that the large-scale genetic homogenization of admixed populations could drive the loss of valuable local adaptations. We also point out that the authors did not account for the extensive gene flow that occurs between managed and wild/feral honeybee populations and raise concerns about the data set used. Finally, we caution against underestimating the importance of genetic diversity for honeybee colonies and highlight the importance of promoting the use of endemic honeybee subspecies in apiculture.     

Deciding on the wing: in-flight decision making and search space sampling in the red dwarf honeybee Apis florea

During reproductive swarming and seasonal migration, a honeybee swarm needs to locate and move to a new, suitable nest site. While the nest-site selection process in cavity-nesting species such as the European honeybee Apis mellifera is very precise with the swarm carefully selecting a single site, open-nesting species, such as Apis florea, lack such precision. These differences in precision in the nest-site selection process are thought to arise from the differing nest-site requirements of open- and cavity-nesting species. While A. florea can nest on almost any tree, A. mellifera is constrained by the scarcity of suitable nest sites. Here we show that imprecision in the nest-site selection process allows swarms to quickly reach a decision when many nest sites are available. In contrast, a very precise nest-site selection process slows down the decision-making process when nest sites are abundant.     

10,000 years in isolation? Honeybees (<Emphasis Type="Italic">Apis mellifera</Emphasis>) in Saharan oases

After the transition from a savannah to a desert about 10,000 years ago the isolated Saharan oases offer a unique case for studying the effect of population fragmentation and isolation over a period of many thousand years. We use the honeybee, Apis mellifera, as a test system because they are an abundant wild species in the African dry savannahs but are particularly sensitive to drift and bottlenecks in small isolated populations due to the small effective size resulting from male haploidy, the sex determination system and sociality. We compared the non-fragmented coastal population with the oases of Brak and Kufra using 15 polymorphic microsatellite loci assessing the mating frequency, colony density, gene diversity, and population differentiation. We found that the honeybee population of the remote oasis of Kufra is well isolated whereas those of the oasis of Brak and the coastal regions show genetic foot prints of introgression by commercial beekeeping. The isolated Kufra population showed no indications of inbreeding suggesting that the endemic population size is sufficient to ensure sustainable local survival.     

Genetic and morphological variation of bee-parasitic <Emphasis Type="Italic">Tropilaelaps</Emphasis> mites (Acari: Laelapidae): new and re-defined species

Mites in the genus Tropilaelaps are parasites of social honeybees. Two species, Tropilaelaps clareae and T. koenigerum, have been recorded and their primary hosts are presumed to be the giant honeybees of Asia, Apis dorsata and A. laboriosa. The most common species, T. clareae, is also an economically important pest of the introduced Western honeybee (A. mellifera) throughout Asia and is considered an emerging threat to world apiculture. In the studies reported here, genetic (mtDNA CO-I and nuclear ITS1-5.8S-ITS2 gene sequence) and morphological variation and host associations were examined among Tropilaelaps isolates collected from A. dorsata, A. laboriosa and A. mellifera throughout Asia and neighbouring regions. The results clearly indicate that the genus contains at least four species. Tropilaelaps clareae, previously assumed to be ubiquitous in Asia, was found to be two species, and it is here redefined as encompassing haplotypes (mites with distinct mtDNA gene sequences) that parasitise native A. dorsata breviligula and introduced A. mellifera in the Philippines and also native A. d. binghami on Sulawesi Island in Indonesia. Tropilaelaps mercedesae n. sp., which until now has been mistaken for T. clareae, encompasses haplotypes that, together with haplotypes of T. koenigerum, parasitise native A. d. dorsata in mainland Asia and Indonesia (except Sulawesi Island). It also parasitises introduced A. mellifera in these and surrounding regions and, with another new species, T. thaii n. sp., also parasitises A. laboriosa in mountainous Himalayan regions. Methods are described for identifying each species. These studies help to clarify the emerging threat of Tropilaelaps to world apiculture and will necessitate a revision of quarantine protocols for countries that import and export honeybees.     

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.     

Gene duplication and the evolution of phenotypic diversity in insect societies

Gene duplication is an important evolutionary process thought to facilitate the evolution of phenotypic diversity. We investigated if gene duplication was associated with the evolution of phenotypic differences in a highly social insect, the honeybee Apis mellifera. We hypothesized that the genetic redundancy provided by gene duplication could promote the evolution of social and sexual phenotypes associated with advanced societies. We found a positive correlation between sociality and rate of gene duplications across the Apoidea, indicating that gene duplication may be associated with sociality. We also discovered that genes showing biased expression between A. mellifera alternative phenotypes tended to be found more frequently than expected among duplicated genes than singletons. Moreover, duplicated genes had higher levels of caste‐, sex‐, behavior‐, and tissue‐biased expression compared to singletons, as expected if gene duplication facilitated phenotypic differentiation. We also found that duplicated genes were maintained in the A. mellifera genome through the processes of conservation, neofunctionalization, and specialization, but not subfunctionalization. Overall, we conclude that gene duplication may have facilitated the evolution of social and sexual phenotypes, as well as tissue differentiation. Thus this study further supports the idea that gene duplication allows species to evolve an increased range of phenotypic diversity.