The genome of the blind bee louse fly reveals deep convergences with its social host and illuminates Drosophila origins.

Social insects' nests harbor intruders known as inquilines,1 which are usually related to their hosts.2,3 However, distant non-social inquilines may also show convergences with their hosts,4,5 although the underlying genomic changes remain unclear. We analyzed the genome of the wingless and blind bee louse fly Braula coeca, an inquiline kleptoparasite of the western honey bee, Apis mellifera.6,7 Using large phylogenomic data, we confirmed recent accounts that the bee louse fly is a drosophilid8,9 and showed that it had likely evolved from a sap-breeder ancestor associated with honeydew and scale insects' wax. Unlike many parasites, the bee louse fly genome did not show significant erosion or strict reliance on an endosymbiont, likely due to a relatively recent age of inquilinism. However, we observed a horizontal transfer of a transposon and a striking parallel evolution in a set of gene families between the honey bee and the bee louse fly. Convergences included genes potentially involved in metabolism and immunity and the loss of nearly all bitter-tasting gustatory receptors, in agreement with life in a protective nest and a diet of honey, pollen, and beeswax. Vision and odorant receptor genes also exhibited rapid losses. Only genes whose orthologs in the closely related Drosophila melanogaster respond to honey bee pheromone components or floral aroma were retained, whereas the losses included orthologous receptors responsive to the anti-ovarian honey bee queen pheromones. Hence, deep genomic convergences can underlie major phenotypic transitions during the evolution of inquilinism between non-social parasites and their social hosts.

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.

Genotypic trade-off between appetitive and aversive capacities in honeybees.

Honey bees can learn both appetitive and aversive associations, using two olfactory conditioning protocols. Appetitive conditioning of the proboscis extension response (PER) involves associating an odor, the conditioned stimulus (CS) with a sucrose solution, the unconditioned stimulus (US). Conversely, aversive conditioning of the sting extension response (SER) involves associating the odor CS with an electric or thermal shock US. Here, we investigated the relationship between bees' appetitive and aversive learning capacities at the individual level and the influence of bees' genotype. As learning performance was shown to depend on an individuals' sensitivity to the US, we systematically measured four different traits in each individual bee: sensitivity to sucrose, PER learning performance with a sucrose US, sensitivity to temperature, SER learning with a temperature US. First, we confirmed for both conditioning types that learning performance correlates with US responsiveness. Second, we found a trade-off between appetitive and aversive learning performances: bees that were better appetitive learners (and had a lower sucrose US threshold) learned less efficiently in the aversive conditioning (and had a higher temperature US threshold). Because the honey bee queen typically mates with 15-20 males, the workers from a honey bee hive belong to as many different patrilines, allowing for the search of the genetic determinism of cognitive abilities. Using microsatellite analysis, we show that a genetic determinism underlies the trade-off between appetitive and aversive capacities, with appetitively vs aversively biased patrilines. The honey bee hive thus appears as a genetically structured cognitive community.

The Conservation of Native Honey Bees Is Crucial.

Recent studies have emphasized the role of the western honey bee, Apis mellifera, as a managed agricultural species worldwide, but also as a potential threat to endangered wild pollinators. This has resulted in the suggestion that honey bees should be regulated in natural areas to conserve wild pollinators. We argue that this perspective fails to appreciate the multifaceted nature of honey bees as native or introduced species with either managed or wild colonies. Wild populations of A. mellifera are currently imperiled, and natural areas are critical for the conservation of local subspecies and genotypes. We propose that a differentiation between managed and wild populations is required and encourage integrated conservation planning for all endangered wild bees, including A. mellifera.

Hygroregulation, a key ability for eusocial insects: Native Western European honeybees as a case study.

Sociality has brought many advantages to various hymenoptera species, including their ability of regulating physical factors in their nest (e.g., temperature). Although less studied, humidity is known to be important for egg, larval and pupal development, and also for nectar concentration. Two subspecies of Apis mellifera of the M evolutionary lineage were used as models to test the ability of a superorganism (i.e. honeybee colony) to regulate the humidity in its nest (i.e. "hygroregulation hypothesis") in four conservation centers: two in France (A. m. mellifera) and two in Portugal (A. m. iberiensis). We investigated the ability of both subspecies to regulate the humidity in hives daily, but also during the seasons for one complete year. Our data and statistical analysis demonstrated the capacity of the bees to regulate humidity in their hive, regardless of the day, season or subspecies. Furthermore, the study showed that humidity in beehives is stable even during winter, when brood is absent, and when temperature is known to be less stable in the beehives. These results suggest that humidity is important for honeybees at every life stage, maybe because of the 'imprint' of the evolutionary history of this hymenopteran lineage.

High sample throughput genotyping for estimating C-lineage introgression in the dark honeybee: an accurate and cost-effective SNP-based tool.

The natural distribution of the honeybee (Apis mellifera L.) has been changed by humans in recent decades to such an extent that the formerly widest-spread European subspecies, Apis mellifera mellifera, is threatened by extinction through introgression from highly divergent commercial strains in large tracts of its range. Conservation efforts for A. m. mellifera are underway in multiple European countries requiring reliable and cost-efficient molecular tools to identify purebred colonies. Here, we developed four ancestry-informative SNP assays for high sample throughput genotyping using the iPLEX Mass Array system. Our customized assays were tested on DNA from individual and pooled, haploid and diploid honeybee samples extracted from different tissues using a diverse range of protocols. The assays had a high genotyping success rate and yielded accurate genotypes. Performance assessed against whole-genome data showed that individual assays behaved well, although the most accurate introgression estimates were obtained for the four assays combined (117 SNPs). The best compromise between accuracy and genotyping costs was achieved when combining two assays (62 SNPs). We provide a ready-to-use cost-effective tool for accurate molecular identification and estimation of introgression levels to more effectively monitor and manage A. m. mellifera conservatories.

Genetic diversity and differentiation among insular honey bee populations in the southwest Indian Ocean likely reflect old geographical isolation and modern introductions.

With globalization the Western honey bee has become a nearly cosmopolitan species, but it was originally restricted to the Old World. This renowned model of biodiversity has diverged into five evolutionary lineages and several geographic "subspecies." If Apis mellifera unicolor is indubitably an African subspecies endemic to Madagascar, its relationship with honey bees from three archipelagos in the southwest Indian Ocean (SWIO) hotspot of biodiversity is misunderstood. We compared recent mtDNA diversity data to an original characterization of the nuclear diversity from honey bees in the Mascarenes and Comoros archipelagos, using 14 microsatellites, but also additional mtDNA tRNALeu-cox2 analysis. Our sampling offers the most comprehensive dataset for the SWIO populations with a total of 3,270 colonies from 10 islands compared with 855 samples from Madagascar, 113 from Africa, and 138 from Europe. Comprehensive mitochondrial screening confirmed that honey bees from La Réunion, Mauritius, and Comoros archipelagos are mainly of African origin (88.1% out of 2,746 colonies) and that coexistence with European lineages occurs only in the Mascarenes. PCA, Bayesian, and genetic differentiation analysis showed that African colonies are not significantly distinct on each island, but have diversified among islands and archipelagos. FST levels progressively decreased in significance from European and African continental populations, to SWIO insular and continental populations, and finally among islands from the same archipelago. Among African populations, Madagascar shared a nuclear background with and was most closely related to SWIO island populations (except Rodrigues). Only Mauritius Island presented clear cytoplasmic disequilibrium and genetic structure characteristic of an admixed population undergoing hybridization, in this case, between A. m. unicolor and A. m. ligustica, A. m. carnica and A. m. mellifera-like individuals. Finally, global genetic clustering analysis helped to better depict the colonization and introduction pattern of honey bee populations in these archipelagos.

Genetic characterization of the honeybee ectoparasitic mite Varroa destructor from Benin (West Africa) using mitochondrial and microsatellite markers.

Varroa destructor is one of the scourges of global beekeeping. It was detected for the first time in Benin in 2011 on the honeybee Apis mellifera adansonii. The aim of this study was to identify the strain of Varroa sp. found and study its genetic diversity. In total 183 Varroa mites were sampled in 21 municipalities in Benin. The COI intergenic region of each mite mtDNA was amplified by PCR. The SacI restriction enzyme was used to determine the strains of Varroa sp. Only the Korean (K) haplotype, identical to the most prevalent strain in Africa, was detected. Analysis of the genetic diversity of Varroa mites with eight microsatellite loci (Simple Sequence Repeats) indicated a very low diversity of genotypes. Thus, V. destructor populations from Benin appear to make up a single group. Their clonal wealth ranges from 0.00 to 0.47. This study is an important step forward in the monitoring of the infestation of V. destructor.

Molecular characterization and functional expression of the Apis mellifera voltage-dependent Ca2+ channels.

Voltage-gated Ca(2+) channels allow the influx of Ca(2+) ions from the extracellular space upon membrane depolarization and thus serve as a transducer between membrane potential and cellular events initiated by Ca(2+) transients. Most insects are predicted to possess three genes encoding Cavα, the main subunit of Ca(2+) channels, and several genes encoding the two auxiliary subunits, Cavß and Cavα2δ; however very few of these genes have been cloned so far. Here, we cloned three full-length cDNAs encoding the three Cavα subunits (AmelCav1a, AmelCav2a and AmelCav3a), a cDNA encoding a novel variant of the Cavß subunit (AmelCavßc), and three full-length cDNAs encoding three Cavα2δ subunits (AmelCavα2δ1 to 3) of the honeybee Apis mellifera. We identified several alternative or mutually exclusive exons in the sequence of the AmelCav2 and AmelCav3 genes. Moreover, we detected a stretch of glutamine residues in the C-terminus of the AmelCav1 subunit that is reminiscent of the motif found in the human Cav2.1 subunit of patients with Spinocerebellar Ataxia type 6. All these subunits contain structural domains that have been identified as functionally important in their mammalian homologues. For the first time, we could express three insect Cavα subunits in Xenopus oocytes and we show that AmelCav1a, 2a and 3a form Ca(2+) channels with distinctive properties. Notably, the co-expression of AmelCav1a or AmelCav2a with AmelCavßc and AmCavα2δ1 produces High Voltage-Activated Ca(2+) channels. On the other hand, expression of AmelCav3a alone leads to Low Voltage-Activated Ca(2+) channels.

MtDNA COI-COII marker and drone congregation area: an efficient method to establish and monitor honeybee (Apis mellifera L.) conservation centres.

Honeybee subspecies have been affected by human activities in Europe over the past few decades. One such example is the importation of nonlocal subspecies of bees which has had an adverse impact on the geographical repartition and subsequently on the genetic diversity of the black honeybee Apis mellifera mellifera. To restore the original diversity of this local honeybee subspecies, different conservation centres were set up in Europe. In this study, we established a black honeybee conservation centre Conservatoire de l'Abeille Noire d'Ile de France (CANIF) in the region of Ile-de-France, France. CANIF's honeybee colonies were intensively studied over a 3-year period. This study included a drone congregation area (DCA) located in the conservation centre. MtDNA COI-COII marker was used to evaluate the genetic diversity of CANIF's honeybee populations and the drones found and collected from the DCA. The same marker (mtDNA) was used to estimate the interactions and the haplotype frequency between CANIF's honeybee populations and 10 surrounding honeybee apiaries located outside of the CANIF. Our results indicate that the colonies of the conservation centre and the drones of the DCA show similar stable profiles compared to the surrounding populations with lower level of introgression. The mtDNA marker used on both DCA and colonies of the conservation centre seems to be an efficient approach to monitor and maintain the genetic diversity of the protected honeybee populations.

Genotypic influence on aversive conditioning in honeybees, using a novel thermal reinforcement procedure.

In Pavlovian conditioning, animals learn to associate initially neutral stimuli with positive or negative outcomes, leading to appetitive and aversive learning respectively. The honeybee (Apis mellifera) is a prominent invertebrate model for studying both versions of olfactory learning and for unraveling the influence of genotype. As a queen bee mates with about 15 males, her worker offspring belong to as many, genetically-different patrilines. While the genetic dependency of appetitive learning is well established in bees, it is not the case for aversive learning, as a robust protocol was only developed recently. In the original conditioning of the sting extension response (SER), bees learn to associate an odor (conditioned stimulus - CS) with an electric shock (unconditioned stimulus - US). This US is however not a natural stimulus for bees, which may represent a potential caveat for dissecting the genetics underlying aversive learning. We thus first tested heat as a potential new US for SER conditioning. We show that thermal stimulation of several sensory structures on the bee's body triggers the SER, in a temperature-dependent manner. Moreover, heat applied to the antennae, mouthparts or legs is an efficient US for SER conditioning. Then, using microsatellite analysis, we analyzed heat sensitivity and aversive learning performances in ten worker patrilines issued from a naturally inseminated queen. We demonstrate a strong influence of genotype on aversive learning, possibly indicating the existence of a genetic determinism of this capacity. Such determinism could be instrumental for efficient task partitioning within the hive.

A fifth major genetic group among honeybees revealed in Syria.

BACKGROUND: Apiculture has been practiced in North Africa and the Middle-East from antiquity. Several thousand years of selective breeding have left a mosaic of Apis mellifera subspecies in the Middle-East, many uniquely adapted and survived to local environmental conditions. In this study we explore the genetic diversity of A. mellifera from Syria (n = 1258), Lebanon (n = 169) and Iraq (n = 35) based on 14 short tandem repeat (STR) loci in the context of reference populations from throughout the Old World (n = 732). RESULTS: Our data suggest that the Syrian honeybee Apis mellifera syriaca occurs in both Syrian and Lebanese territories, with no significant genetic variability between respective populations from Syria and Lebanon. All studied populations clustered within a new fifth independent nuclear cluster, congruent with an mtDNA Z haplotype identified in a previous study. Syrian honeybee populations are not associated with Oriental lineage O, except for sporadic introgression into some populations close to the Turkish and Iraqi borders. Southern Syrian and Lebanese populations demonstrated high levels of genetic diversity compared to the northern populations. CONCLUSION: This study revealed the effects of foreign queen importations on Syrian bee populations, especially for the region of Tartus, where extensive introgression of A. m. anatolica and/or A. m. caucasica alleles were identified. The policy of creating genetic conservation centers for the Syrian subspecies should take into consideration the influence of the oriental lineage O from the northern Syrian border and the large population of genetically divergent indigenous honeybees located in southern Syria.

Comparative susceptibility of three Western honeybee taxa to the microsporidian parasite Nosema ceranae.

Genetic diversity of a host species is a key factor to counter infection by parasites. Since two separation events and the beginning of beekeeping, the Western honeybee, Apis mellifera, has diverged in many phylogenetically-related taxa that share common traits but also show specific physiological, behavioural and morphological traits. In this study, we tested the hypothesis that A. mellifera taxa living in a same habitat should respond differently to parasites like Nosema ceranae, a microsporidia living in host's midgut. We used the Poulin and Combes' concept of virulence to compare the susceptibility of three A. mellifera taxa to N. ceranae infection. Three criteria were measured 10 days post-infection (dpi): the host mortality, the host sugar consumption and the development success of the parasite (i.e. number of spores produced). Interestingly, we showed that the observed variation in susceptibility to infection by N. ceranae is not linked to honeybee taxa but results from the variability between colonies, and that those differences are probably linked to genetic variations. The use of these three criteria allows us to conclude that the differences in susceptibility are mediated by a genetic variability in honeybee workers from resistance to tolerance. Finally, we discuss the consequences of our findings for beekeeping management.

THE ORIGIN OF WEST EUROPEAN SUBSPECIES OF HONEYBEES (APIS MELLIFERA): NEW INSIGHTS FROM MICROSATELLITE AND MITOCHONDRIAL DATA.

Apis mellifera is composed of three evolutionary branches including mainly African (branch A), western and northern European (branch M), and southeastern European (branch C) populations. The existence of morphological clines extending from the equator to the Polar Circle through Morocco and Spain raised the hypothesis that the branch M originated in Africa. Mitochondrial DNA analysis revealed that branches A and M were characterized by highly diverged lineages implying very remote links between both branches. It also revealed that mtDNA haplotypes from lineages A coexisted with haplotypes M in the Iberian Peninsula and formed a south-north frequency cline, suggesting that this area could be a secondary contact zone between the two branches. By analyzing 11 populations sampled along a France-Spain/Portugal-Morocco-Guinea transect at 8 microsatellite loci and the DraI RFLP of the COI-COII mtDNA marker, we show that Iberian populations do not present any trace of "africanization" and are very similar to French populations when considering microsatellite markers. Therefore, the Iberian Peninsula is not a transition area. The higher haplotype A variability observed in Spanish and Portuguese samples compared to that found in Africa is explained by a higher mutation rate and multiple and recent introductions. Selection appears to be the best explanation to the morphological and allozymic clines and to the diffusion and maintenance of African haplotypes in Spain and Portugal.