A blood test to monitor bee health across a European network of agricultural sites of different land-use by MALDI BeeTyping mass spectrometry

There are substantial concerns about impaired honey bee health and colony losses due to several poorly understood factors. We used MALDI profiling (MALDI BeeTyping®) analysis to investigate how some environmental and management factors under field conditions across Europe affected the honey bee haemolymph peptidome (all peptides in the circulatory fluid), as a profile of molecular markers representing the immune status of Apis mellifera. Honey bees were exposed to a range of environmental stressors in 128 agricultural sites across eight European countries in four biogeographic zones, with each country contributing eight sites each for two different cropping systems: oilseed rape (OSR) and apple (APP). The full haemolymph peptide profiles, including the presence and levels of three key immunity markers, namely the antimicrobial peptides (AMPs) Apidaecin, Abaecin and Defensin-1, allowed the honey bee responses to environmental variables to be discriminated by country, crop type and site. When considering just the AMPs, it was not possible to distinguish between countries by the prevalence of each AMP in the samples. However, it was possible to discriminate between countries on the amounts of the AMPs, with the Swedish samples in particular expressing high amounts of all AMPs. A machine learning model was developed to discriminate the haemolymphs of bees from APP and OSR sites. The model was 90.6 % accurate in identifying the crop type from the samples used to build the model. Overall, MALDI BeeTyping® of bee haemolymph represents a promising and cost-effective "blood test" for simultaneously monitoring dozens of peptide markers affected by environmental stressors at the landscape scale, thus providing policymakers with new diagnostic and regulatory tools for monitoring bee health.

Pesticide use negatively affects bumble bees across European landscapes

Sustainable agriculture requires balancing crop yields with the effects of pesticides on non-target organisms, such as bees and other crop pollinators. Field studies demonstrated that agricultural use of neonicotinoid insecticides can negatively affect wild bee species1,2, leading to restrictions on these compounds3. However, besides neonicotinoids, field-based evidence of the effects of landscape pesticide exposure on wild bees is lacking. Bees encounter many pesticides in agricultural landscapes4-9 and the effects of this landscape exposure on colony growth and development of any bee species remains unknown. Here we show that the many pesticides found in bumble bee-collected pollen are associated with reduced colony performance during crop bloom, especially in simplified landscapes with intensive agricultural practices. Our results from 316 Bombus terrestris colonies at 106 agricultural sites across eight European countries confirm that the regulatory system fails to sufficiently prevent pesticide-related impacts on non-target organisms, even for a eusocial pollinator species in which colony size may buffer against such impacts10,11. These findings support the need for postapproval monitoring of both pesticide exposure and effects to confirm that the regulatory process is sufficiently protective in limiting the collateral environmental damage of agricultural pesticide use.

Distribution of infectious and parasitic agents among three sentinel bee species across European agricultural landscapes

Infectious and parasitic agents (IPAs) and their associated diseases are major environmental stressors that jeopardize bee health, both alone and in interaction with other stressors. Their impact on pollinator communities can be assessed by studying multiple sentinel bee species. Here, we analysed the field exposure of three sentinel managed bee species (Apis mellifera, Bombus terrestris and Osmia bicornis) to 11 IPAs (six RNA viruses, two bacteria, three microsporidia). The sentinel bees were deployed at 128 sites in eight European countries adjacent to either oilseed rape fields or apple orchards during crop bloom. Adult bees of each species were sampled before their placement and after crop bloom. The IPAs were detected and quantified using a harmonised, high-throughput and semi-automatized qPCR workflow. We describe differences among bee species in IPA profiles (richness, diversity, detection frequencies, loads and their change upon field exposure, and exposure risk), with no clear patterns related to the country or focal crop. Our results suggest that the most frequent IPAs in adult bees are more appropriate for assessing the bees' IPA exposure risk. We also report positive correlations of IPA loads supporting the potential IPA transmission among sentinels, suggesting careful consideration should be taken when introducing managed pollinators in ecologically sensitive environments.

Sulfoxaflor effects depend on the interaction with other pesticides and Nosema ceranae infection in the honey bee (Apis mellifera)

Honey bees health is compromised by many factors such as the use of agrochemicals in agriculture and the various diseases that can affect them. Multiple studies have shown that these factors can interact, producing a synergistic effect that can compromise the viability of honey bees. This study analyses the interactions between different pesticides and the microsporidium Nosema ceranae and their effect on immune and detoxification gene expression, sugar consumption and mortality in the Iberian western honey bee (Apis mellifera iberiensis). For this purpose, workers were infected with N. ceranae and subjected to a sugar-water diet with field concentrations of the pesticides sulfoxaflor, azoxystrobin and glyphosate. Increased sugar intake and altered immune and cytochrome P450 gene expression were observed in workers exposed to sulfoxaflor and infected with N. ceranae. None of the pesticides affected Nosema spore production in honey bee gut. Of the three pesticides tested (alone or in combination) only sulfoxaflor increased mortality in honey bees. Taken together, our results suggest that the effects of sulfoxaflor were attenuated in contact with other pesticides, and that Nosema infection leads to increase sugar intake in sulfoxaflor-exposed bees. Overall, this underlines the importance of studying the interaction between different stressors to understand their overall impact not only on honey bee but also on wild bees health.

Evolutionarily diverse origins of deformed wing viruses in western honey bees

Novel transmission routes can allow infectious diseases to spread, often with devastating consequences. Ectoparasitic varroa mites vector a diversity of RNA viruses, having switched hosts from the eastern to western honey bees (Apis cerana to Apis mellifera). They provide an opportunity to explore how novel transmission routes shape disease epidemiology. As the principal driver of the spread of deformed wing viruses (mainly DWV-A and DWV-B), varroa infestation has also driven global honey bee health declines. The more virulent DWV-B strain has been replacing the original DWV-A strain in many regions over the past two decades. Yet, how these viruses originated and spread remains poorly understood. Here, we use a phylogeographic analysis based on whole-genome data to reconstruct the origins and demography of DWV spread. We found that, rather than reemerging in western honey bees after varroa switched hosts, as suggested by previous work, DWV-A most likely originated in East Asia and spread in the mid-20th century. It also showed a massive population size expansion following the varroa host switch. By contrast, DWV-B was most likely acquired more recently from a source outside East Asia and appears absent from the original varroa host. These results highlight the dynamic nature of viral adaptation, whereby a vector's host switch can give rise to competing and increasingly virulent disease pandemics. The evolutionary novelty and rapid global spread of these host-virus interactions, together with observed spillover into other species, illustrate how increasing globalization poses urgent threats to biodiversity and food security.

Bee Trypanosomatids: First Steps in the Analysis of the Genetic Variation and Population Structure of Lotmaria passim, Crithidia bombi and Crithidia mellificae

Trypanosomatids are among the most prevalent parasites in bees but, despite the fact that their impact on the colonies can be quite important and that their infectivity may potentially depend on their genotypes, little is known about the population diversity of these pathogens. Here we cloned and sequenced three non-repetitive single copy loci (DNA topoisomerase II, glyceraldehyde-3-phosphate dehydrogenase and RNA polymerase II large subunit, RPB1) to produce new genetic data from Crithidia bombi, C. mellificae and Lotmaria passim isolated from honeybees and bumblebees. These were analysed by applying population genetic tools in order to quantify and compare their variability within and between species, and to obtain information on their demography and population structure. The general pattern for the three species was that (1) they were subject to the action of purifying selection on nonsynonymous variants, (2) the levels of within species diversity were similar irrespective of the host, (3) there was evidence of recombination among haplotypes and (4) they showed no haplotype structuring according to the host. C. bombi exhibited the lowest levels of synonymous variation (π_S= 0.06 ± 0.04 %) - and a mutation frequency distribution compatible with a population expansion after a bottleneck - that contrasted with the extensive polymorphism displayed by C. mellificae (π_S= 2.24 ± 1.00 %), which likely has a more ancient origin. L. passim showed intermediate values (π_S= 0.40 ± 0.28 %) and an excess of variants a low frequencies probably linked to the spread of this species to new geographical areas.

“Migratory beekeeping and its influence on the prevalence and dispersal of pathogens to managed and wild bees”

Demand for food is growing along with the human population, leading to an increase in plant production. Many crops are pollinated by insects, so the global demand for managed pollinators is also increasing. The honey bee has traditionally been considered the main provider of crop pollination services. For providing it beekeepers seasonally transport hives to different locations after the flowering of different crops. These movements could be detrimental to pollinators by: i) stressing honey bees, making them more susceptible to pathogens and parasites; ii) spreading bee parasites and pathogens across locations; iii) increasing the transmission of parasites and pathogens between managed and wild pollinators and vice versa (spillover and spillback, respectively). To understand the impact of migratory beekeeping on bee health, we conducted a systematic review to identify the main trends and provide a complete picture of existing knowledge on the subject. We found 52 studies analysing pathogen-related impacts of migratory beekeeping on honey bees. However, only 16 investigations tested the effect of migratory practices on the prevalence and spread of pathogens and parasites. We found no studies that assessed the impact of migratory beekeeping on the occurrence and spread of pests and diseases in wild bees. In general, migratory beekeeping tends to increase the prevalence of pathogens and parasites in honey bee colonies. However, the results were very heterogeneous, probably due to several uncontrolled underlying factors such as management, biological and geographical factors, and the interactions between them. In conclusion, there is an urgent need for studies to assess the impact of migratory beekeeping on bee health, given the current global bee decline and the expected increase in migratory beekeeping due to climate change and crop pollination demand.

•

Migratory beekeeping may affect the health of both managed and wild bees.

•

Only 16 studies compare the prevalence of pathogens between migratory and stationary hives.

•

Results are heterogeneous due to the complex interaction of the underlying factors.

•

Overall, the prevalence of pathogens and parasites increased in migratory honey bees.

•

The impact of migratory beekeeping on the health of wild bees is unknown.

Detection of Microsporidia in Pollinator Communities of a Mediterranean Biodiversity Hotspot for Wild Bees

Insect pollination is crucial for the maintenance of natural and managed ecosystems but the functioning of this ecosystem service is threatened by a worldwide decline of pollinators. Key factors in this situation include the spread and interspecific transmission of pathogens worldwide through the movement of managed pollinators. Research on this field has been mainly conducted in some particular species, while studies assessing the interspecific transmission of pathogens at a community level are scarce. However, this information is pivotal to design strategies to protect pollinators. Herein, we analysed the prevalence of two common microsporidia pathogens of managed honey bees (Nosema ceranae and N. apis) in bee communities of semiarid Mediterranean areas from the Southeast of the Iberian Peninsula. Our results confirm the ability of N. ceranae to disperse across wild bee communities in semiarid Mediterranean ecosystems since it was detected in 36 Apoidea species (39% of the sampling; for the first time in nine genera). The prevalence of the pathogen did not show any phylogenetic signal which suggests a superfamily host range of the pathogen or that wild bees may be acting only as vectors of N. ceranae. In addition, N. apis was detected in an Eucera species, which is the second time it has been detected by molecular techniques in a host other than the honey bee. Our study represents the primary assessment of the prevalence of microsporidia at community level in Mediterranean areas and provides outstanding results on the ability of Nosema pathogens to spread across the landscape.

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.

The Effect of Migratory Beekeeping on the Infestation Rate of Parasites in Honey Bee (Apis mellifera) Colonies and on Their Genetic Variability

Migratory beekeeping is a widely extended practice aimed at increasing the yield of products and pollination services of honey bee colonies. However, it represents a stress factor, as it facilitates the dissemination of diseases and may compromise the genetic identity of the colonies involved. To analyze the extent of these effects, pathogens infestation rate and genetic composition were monitored in a field experiment comparing stationary and migratory colonies sharing the same environmental conditions but differing in management (stationary vs. migratory) and genetic background. We studied the pathogens infestation rate (Varroa destructor, Nosema spp., and Deformed Wing Virus (DWV)) at four different times: before migratory operation, two weeks later, at the end of the migratory period, and two weeks after the return of the migratory hives. An increased incidence of V. destructor and Nosema ceranae and a lower DWV viral load were found in migratory colonies. Temporary changes in genetic diversity were detected regardless of colony type, suggesting that stressors other than management affect the genetic diversity of the colonies. Our study demonstrates that migratory practices have variable effects on the health and genetic diversity of honey bee colonies, which should be taken into account for the development of sustainable beekeeping.

Wide diversity of parasites in Bombus terrestris (Linnaeus, 1758) revealed by a high-throughput sequencing approach

Assessing the extent of parasite diversity requires the application of appropriate molecular tools, especially given the growing evidence of multiple parasite co-occurrence. Here, we compared the performance of a next-generation sequencing technology (Ion PGM ™ System) in 12 Bombus terrestris specimens that were PCR-identified as positive for trypanosomatids (Leishmaniinae) in a previous study. These bumblebees were also screened for the occurrence of Nosematidae and Neogregarinorida parasites using both classical protocols (either specific PCR amplification or amplification with broad-range primers plus Sanger sequencing) and Ion PGM sequencing. The latter revealed higher parasite diversity within individuals, especially among Leishmaniinae (which were present as a combination of Lotmaria passim, Crithidia mellificae and Crithidia bombi), and the occurrence of taxa never reported in these hosts: Crithidia acanthocephali and a novel neogregarinorida species. Furthermore, the complementary results produced by the different sets of primers highlighted the convenience of using multiple markers to minimize the chance of some target organisms going unnoticed. Altogether, the deep sequencing methodology offered a more comprehensive way to investigate parasite diversity than the usual identification methods and provided new insights whose importance for bumblebee health should be further analysed.

Intercontinental long-distance seed dispersal across the Mediterranean Basin explains population genetic structure of a bird-dispersed shrub

Long-distance dispersal (LDD) is a pivotal process for plants determining their range of distribution and promoting gene flow among distant populations. Most fleshy-fruited species rely on frugivorous vertebrates to disperse their seeds across the landscape. While LDD events are difficult to record, a few ecological studies have shown that birds move a sizeable number of ingested seeds across geographic barriers, such as sea straits. The foraging movements of migrant frugivores across distant populations, including those separated by geographic barriers, creates a constant flow of propagules that in turn shapes the spatial distributions of the genetic variation in populations. Here, we have analysed the genetic diversity and structure of 74 populations of Pistacia lentiscus, a fleshy-fruited shrub widely distributed in the Mediterranean Basin, to elucidate whether the Mediterranean Sea acts as a geographic barrier or alternatively whether migratory frugivorous birds promote gene flow among populations located on both sides of the sea. Our results show reduced genetic distances among populations, including intercontinental populations, and they show a significant genetic structure across an eastern-western axis. These findings are consistent with known bird migratory routes that connect the European and African continents following a north-southwards direction during the fruiting season of many fleshy-fruited plants. Further, approximate Bayesian analysis failed to explain the observed patterns as a result of historical population migrations at the end of Last Glacial Maximum. Therefore, anthropic and/or climatic changes that would disrupt the migratory routes of frugivorous birds might have genetic consequences for the plant species they feed upon.

The toxic unit approach as a risk indicator in honey bees surveillance programmes: A case of study in Apis mellifera iberiensis

The influence of genetic diversity and exposure to xenobiotics on the prevalence of pathogens was studied within the context of a voluntary epidemiological study in Spanish apiaries of Apis mellifera iberiensis, carried out during the spring season of years 2014 and 2015. As such, the evolutionary lineages of the honey bee colonies were identified, a multiresidue analysis of xenobiotics was carried out in beebread and worker bee samples, and the Toxic Unit (TUm) was estimated for each sampled apiary. The relationship between lineages and the most prevalent pathogens (Nosema ceranae, Varroa destructor, trypanosomatids, Black Queen Cell Virus; and Deformed Wing Virus) was analysed with contingency tables, and the possible relationships between TUm and the prevalence of these pathogens were studied by using a factor analysis. The statistical analysis supported the associations between V. destructor and Deformed Wing Virus (DWV), and between N. ceranae and Black Queen Cell Virus (BQCV), but the association between these pathogens and trypanosomatids was not observed. TUm values varied between 5.5 × 10^-6 and 3.65 × 10^-1. When TUm < 3.35 × 10^-4, it was mainly determined by coumaphos, tau-fluvalinate and/or chlorfenvinphos. At higher values, other insecticides also contributed to TUm, although a clear predominance was not seen up to TUm ≥ 1.83 × 10^-2, when it was mainly defined by acrinathrin, spinosad and/or imidacloprid. The possible cumulative effect from the joint action of xenobiotics was >10% in the 63% of the cases. The prevalence of pathogens did not appear to be influenced by the distribution of evolutionary lineages and, while the prevalence of V. destructor was not found to be determined by TUm, there was a trend towards an increasing prevalence of N. ceranae when TUm ≥ 23 10^-4. This study is an example of using TUm approach beyond the field of the ecotoxicology.

An integrative approach to discern the seed dispersal role of frugivorous guilds in a Mediterranean semiarid priority habitat

Seed dispersal is an essential process to maintain the viability of plant populations, and understanding this ecological process allows management strategies to be developed to conserve ecosystems. European Union priority habitat 5220* is defined as “Mediterranean arborescent shrubland with Ziziphus lotus” and it represents a favorable microclimate within the severe climatic conditions typical of the semiarid south-eastern region of the Iberian Peninsula. Therefore, the study of seed dispersal in this priority habitat by different frugivorous guilds, is a challenge for its conservation. In this study, we have characterized a mutualistic network of seed dispersal that is mediated by vertebrates (mammals and birds) in the protected habitat 5220*. The aims of this study were to: (i) identify the seed disperser community; (ii) analyze the relative role of key species in the dispersal process; and (iii) compare the functional ecology of the seed dispersal process between mammals and birds. As such, we collected animal faeces to determine seed dispersers taxonomy, identifying the mammals through the visual aspect of the faeces and the birds by DNA barcoding. In the case of birds, we also collected regurgitated seeds in which the disperser species was also identified through molecular techniques. This allowed us to build-up a mutualistic network and to identify the relative role of these animals in seed dispersal. Our results showed that mammals and birds fulfilled complementary roles in seed dispersal, with birds representing the main dispersers of key plants within the 5220* habitat, and mammals the main dispersers of human-cultivated plants. Herein, we provide a useful approach with relevant information that can be used to propose management policies that focus on restoring the threatened 5220* habitat, promoting the role of birds to disperse key species that structure plant communities of this priority habitat.

Coexistence of genetically different Varroa destructor in Apis mellifera colonies

The aim of this study was to investigate the genetic diversity of Varroa destructor parasitizing Apis mellifera colonies and to test for possible host-parasite association at the mitochondrial DNA (mtDNA) level. Six A. mellifera haplotypes (including a novel C2aa) and five haplotypes of V. destructor were detected in 29 analyzed colonies from eight sampling sites in Serbia. We revealed the presence of the K and S1 haplotypes as well as KS1 and KP1 heteroplasmic mite individuals in all localities, while the P1 haplotype was only found in four sampling sites. Significant differences in V. destructor genetic diversity were found at both apiary and colony levels, with mite haplotypes coexisting in almost all tested colonies. In addition, a significant correlation between the number of analyzed mites per colony and the number of identified V. destructor haplotypes was observed. However, no significant host-parasite relationship was found, suggesting that mites bearing different haplotypes as well as those heteroplasmic individuals are well adapted to the host, A. mellifera, independently of the identified haplotype present in each colony. Our results will contribute to future population and biogeographic studies concerning V. destructor infesting A. mellifera, as well as to better understanding their host-parasite relationship.

El código de barras de ADN confirma la distribución de <em>Bombus magnus</em> (Vogt, 1911) (Hymenoptera: Apidae) en la península Ibérica

Bombus magnus (Vogt, 1911) (Hymenoptera: Apidae) es una de las tres especies crípticas pertenecientes al complejo lucorum junto con B. lucorum (Linnaeus, 1761) y B. cryptarum (Fabricius, 1775). En la península Ibérica solo se encuentran B. lurocum y B. magnus pero la presencia de esta última no ha sido confirmada al sur de los Pirineos. Dada su similitud morfológica, usamos la región del código de barras de ADN para identificar 113 individuos de este complejo de especies en un muestreo ibérico. Los resultados confirman la presencia de B. magnus en los Pirineos y amplían su distribución actual hacia la meseta Norte ibérica. Dados estos resultados, sugerimos que ha de revisarse su distribución y el estado de conservación de esta especie en la península Ibérica.

Beekeeping practices and geographic distance, not land use, drive gene flow across tropical bees

Across the globe, wild bees are threatened by ongoing natural habitat loss, risking the maintenance of plant biodiversity and agricultural production. Despite the ecological and economic importance of wild bees and the fact that several species are now managed for pollination services worldwide, little is known about how land use and beekeeping practices jointly influence gene flow. Using stingless bees as a model system, containing wild and managed species that are presumed to be particularly susceptible to habitat degradation, here we examine the main drivers of tropical bee gene flow. We employ a novel landscape genetic approach to analyse data from 135 populations of 17 stingless bee species distributed across diverse tropical biomes within the Americas. Our work has important methodological implications, as we illustrate how a maximum-likelihood approach can be applied in a meta-analysis framework to account for multiple factors, and weight estimates by sample size. In contrast to previously held beliefs, gene flow was not related to body size or deforestation, and isolation by geographic distance (IBD) was significantly affected by management, with managed species exhibiting a weaker IBD than wild ones. Our study thus reveals the critical importance of beekeeping practices in shaping the patterns of genetic differentiation across bee species. Additionally, our results show that many stingless bee species maintain high gene flow across heterogeneous landscapes. We suggest that future efforts to preserve wild tropical bees should focus on regulating beekeeping practices to maintain natural gene flow and enhancing pollinator-friendly habitats, prioritizing species showing a limited dispersal ability.

Population Genetics of Nosema apis and Nosema ceranae: One Host (Apis mellifera) and Two Different Histories

Two microsporidians are known to infect honey bees: Nosema apis and Nosema ceranae. Whereas population genetics data for the latter have been released in the last few years, such information is still missing for N. apis. Here we analyze the patterns of nucleotide polymorphism at three single-copy loci (PTP2, PTP3 and RPB1) in a collection of Apis mellifera isolates from all over the world, naturally infected either with N. apis (N = 22) or N. ceranae (N = 23), to provide new insights into the genetic diversity, demography and evolution of N. apis, as well as to compare them with evidence from N. ceranae. Neutral variation in N. apis and N. ceranae is of the order of 1%. This amount of diversity suggests that there is no substantial differentiation between the genetic content of the two nuclei present in these parasites, and evidence for genetic recombination provides a putative mechanism for the flow of genetic information between chromosomes. The analysis of the frequency spectrum of neutral variants reveals a significant surplus of low frequency variants, particularly in N. ceranae, and suggests that the populations of the two pathogens are not in mutation-drift equilibrium and that they have experienced a population expansion. Most of the variation in both species occurs within honey bee colonies (between 62%-90% of the total genetic variance), although in N. apis there is evidence for differentiation between parasites isolated from distinct A. mellifera lineages (20%-34% of the total variance), specifically between those collected from lineages A and C (or M). This scenario is consistent with a long-term host-parasite relationship and contrasts with the lack of differentiation observed among host-lineages in N. ceranae (< 4% of the variance), which suggests that the spread of this emergent pathogen throughout the A. mellifera worldwide population is a recent event.

Stable genetic diversity despite parasite and pathogen spread in honey bee colonies

In the last decades, the rapid spread of diseases, such as varroosis and nosemosis, associated with massive honey bee colonies mortality around the world has significantly decreased the number and size of honey bee populations and possibly their genetic diversity. Here, we compare the genetic diversity of Iberian honey bee colonies in two samplings performed in 2006 and 2010 in relation to the presence of the pathogenic agents Nosema apis, Nosema ceranae, and Varroa destructor in order to determine whether parasite and pathogen spread in honey bee colonies reflects changes in genetic diversity. We found that the genetic diversity remained similar, while the incidence of N. ceranae increased and the incidence of N. apis and V. destructor decreased slightly. These results indicate that the genetic diversity was not affected by the presence of these pathogenic agents in the analyzed period. However, the two groups of colonies with and without Nosema/Varroa detected showed significant genetic differentiation (G test). A detailed analysis of the allelic segregation of microsatellite loci in Nosema/Varroa-negative colonies and parasitized ones revealed two outlier loci related to genes involved in immune response.

A worldwide survey of genome sequence variation provides insight into the evolutionary history of the honeybee Apis mellifera

The honeybee Apis mellifera has major ecological and economic importance. We analyze patterns of genetic variation at 8.3 million SNPs, identified by sequencing 140 honeybee genomes from a worldwide sample of 14 populations at a combined total depth of 634×. These data provide insight into the evolutionary history and genetic basis of local adaptation in this species. We find evidence that population sizes have fluctuated greatly, mirroring historical fluctuations in climate, although contemporary populations have high genetic diversity, indicating the absence of domestication bottlenecks. Levels of genetic variation are strongly shaped by natural selection and are highly correlated with patterns of gene expression and DNA methylation. We identify genomic signatures of local adaptation, which are enriched in genes expressed in workers and in immune system- and sperm motility-related genes that might underlie geographic variation in reproduction, dispersal and disease resistance. This study provides a framework for future investigations into responses to pathogens and climate change in honeybees.

Signatures of selection in the Iberian honey bee (Apis mellifera iberiensis) revealed by a genome scan analysis of single nucleotide polymorphisms

Understanding the genetic mechanisms of adaptive population divergence is one of the most fundamental endeavours in evolutionary biology and is becoming increasingly important as it will allow predictions about how organisms will respond to global environmental crisis. This is particularly important for the honey bee, a species of unquestionable ecological and economical importance that has been exposed to increasing human-mediated selection pressures. Here, we conducted a single nucleotide polymorphism (SNP)-based genome scan in honey bees collected across an environmental gradient in Iberia and used four FST -based outlier tests to identify genomic regions exhibiting signatures of selection. Additionally, we analysed associations between genetic and environmental data for the identification of factors that might be correlated or act as selective pressures. With these approaches, 4.4% (17 of 383) of outlier loci were cross-validated by four FST -based methods, and 8.9% (34 of 383) were cross-validated by at least three methods. Of the 34 outliers, 15 were found to be strongly associated with one or more environmental variables. Further support for selection, provided by functional genomic information, was particularly compelling for SNP outliers mapped to different genes putatively involved in the same function such as vision, xenobiotic detoxification and innate immune response. This study enabled a more rigorous consideration of selection as the underlying cause of diversity patterns in Iberian honey bees, representing an important first step towards the identification of polymorphisms implicated in local adaptation and possibly in response to recent human-mediated environmental changes.

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.

Genetic structure of Balearic honeybee populations based on microsatellite polymorphism

The genetic variation of honeybee colonies collected in 22 localities on the Balearic Islands (Spain) was analysed using eight polymorphic microsatellite loci. Previous studies have demonstrated that these colonies belong either to the African or west European evolutionary lineages. These populations display low variability estimated from both the number of alleles and heterozygosity values, as expected for the honeybee island populations. Although genetic differentiation within the islands is low, significant heterozygote deficiency is present, indicating a subpopulation genetic structure. According to the genetic differentiation test, the honeybee populations of the Balearic Islands cluster into two groups: Gimnesias (Mallorca and Menorca) and Pitiusas (Ibiza and Formentera), which agrees with the biogeography postulated for this archipelago. The phylogenetic analysis suggests an Iberian origin of the Balearic honeybees, thus confirming the postulated evolutionary scenario for Apis mellifera in the Mediterranean basin. The microsatellite data from Formentera, Ibiza and Menorca show that ancestral populations are threatened by queen importations, indicating that adequate conservation measures should be developed for protecting Balearic bees.

Genetic structure of Balearic honeybee populations based on microsatellite polymorphism

The genetic variation of honeybee colonies collected in 22 localities on the Balearic Islands (Spain) was analysed using eight polymorphic microsatellite loci. Previous studies have demonstrated that these colonies belong either to the African or west European evolutionary lineages. These populations display low variability estimated from both the number of alleles and heterozygosity values, as expected for the honeybee island populations. Although genetic differentiation within the islands is low, significant heterozygote deficiency is present, indicating a subpopulation genetic structure. According to the genetic differentiation test, the honeybee populations of the Balearic Islands cluster into two groups: Gimnesias (Mallorca and Menorca) and Pitiusas (Ibiza and Formentera), which agrees with the biogeography postulated for this archipelago. The phylogenetic analysis suggests an Iberian origin of the Balearic honeybees, thus confirming the postulated evolutionary scenario for Apis mellifera in the Mediterranean basin. The microsatellite data from Formentera, Ibiza and Menorca show that ancestral populations are threatened by queen importations, indicating that adequate conservation measures should be developed for protecting Balearic bees.