Population genetic diversity and dynamics of the honey bee brood pathogen Melissococcus plutonius in a region with high prevalence.

European foulbrood (EFB) is a honey bee brood disease caused by the bacterium Melissococcus plutonius. Large-scale EFB outbreaks have been reported in several countries in recent decades, which entail costly sanitation measures of affected apiaries to restrict the spread of this contagious pathogen. To mitigate its impact, a better understanding of the population dynamics of the etiological agent is required. We here used multi-locus sequence typing (MLST) to infer the genetic diversity and geographical distribution of 160 M. plutonius isolates collected from EFB symptomatic honey bee colonies seven years apart. Isolates belonged to three clonal complexes (CCs) known worldwide and to 12 sequence types (STs), of which five were novel. Phylogenetic and clustering analyses showed that some of these novel sequence types have likely evolved locally during a period of outbreak, but most disappeared again. We further screened the isolates for melissotoxin A (mtxA), a putative virulence gene. The prevalence of STs in which mtxA was frequent increased over time, suggesting that this gene promotes spread. Despite the increased frequency of this gene in the population, the total number of cases decreased, which could be due to stricter control measures implemented before the second sampling period. Our results provide a better understanding of M. plutonius population dynamics and help identify knowledge gaps that limit efficient control of this emerging disease.

Quantitative and Qualitative Changes in the Deformed Wing Virus Population in Honey Bees Associated with the Introduction or Removal of Varroa destructor.

Varroa destructor is an ectoparasitic mite associated with significant losses of honeybee colonies globally. The mite vectors a range of pathogenic viruses, the most important of which is the Deformed wing virus (DWV). In the absence of Varroa, DWV exists as a low-level, highly diverse virus population. However, when transmitted by Varroa, certain variants become highly elevated, and may become near-clonal and cause symptomatic infections. Mite transmission between colonies can occur when parasitised workers drift from or rob adjacent hives. These activities can result in elevated mite levels, but the resulting change in the DWV population, the primary determinant of winter colony losses, has not been determined. In reciprocal studies, we investigated the influence of the removal of mites, or their acquisition, on the DWV population. When mites were removed from heavily infested colonies, there was a striking and rapid reduction in virus load. Conversely, siting Varroa-naïve colonies in a mite-infested apiary resulted in the acquisition of mites and concomitant changes in the virus population. We observed both near-clonal and highly divergent virus populations regardless of titre, suggesting changes were stochastic and colony-specific. Our findings have implications for the outcome of strategies in areas with total or patchy implementation of Varroa control plans.

Landscape-scale drivers of pollinator communities may depend on land-use configuration.

Research into pollinators in managed landscapes has recently combined approaches of pollination ecology and landscape ecology, because key stressors are likely to interact across wide areas. While laboratory and field experiments are valuable for furthering understanding, studies are required to investigate the interacting drivers of pollinator health and diversity across a broader range of landscapes and a wider array of taxa. Here, we use a network of 96 study landscapes in six topographically diverse regions of Britain, to test the combined importance of honeybee density, insecticide loadings, floral resource availability and habitat diversity to pollinator communities. We also explore the interactions between these drivers and the cover and proximity of semi-natural habitat. We found that among our four drivers, only honeybee density was positively related to wild pollinator abundance and diversity, and the positive association between abundance and floral resources depended on insecticide loadings and habitat diversity. By contrast, our exploratory models including habitat composition metrics revealed a complex suite of interactive effects. These results demonstrate that improving pollinator community composition and health is unlikely to be achieved with general resource enhancements only. Rather, local land-use context should be considered in fine-tuning pollinator management and conservation. This article is part of the theme issue 'Natural processes influencing pollinator health: from chemistry to landscapes'.

Identifying the climatic drivers of honey bee disease in England and Wales.

Honey bee colony health has received considerable attention in recent years, with many studies highlighting multifactorial issues contributing to colony losses. Disease and weather are consistently highlighted as primary drivers of colony loss, yet little is understood about how they interact. Here, we combined disease records from government honey bee health inspections with meteorological data from the CEDA to identify how weather impacts EFB, AFB, CBP, varroosis, chalkbrood and sacbrood. Using R-INLA, we determined how different meteorological variables influenced disease prevalence and disease risk. Temperature caused an increase in the risk of both varroosis and sacbrood, but overall, the weather had a varying effect on the six honey bee diseases. The risk of disease was also spatially varied and was impacted by the meteorological variables. These results are an important step in identifying the impacts of climate change on honey bees and honey bee diseases.

Transitions in symbiosis: evidence for environmental acquisition and social transmission within a clade of heritable symbionts.

A dynamic continuum exists from free-living environmental microbes to strict host-associated symbionts that are vertically inherited. However, knowledge of the forces that drive transitions in symbiotic lifestyle and transmission mode is lacking. Arsenophonus is a diverse clade of bacterial symbionts, comprising reproductive parasites to coevolving obligate mutualists, in which the predominant mode of transmission is vertical. We describe a symbiosis between a member of the genus Arsenophonus and the Western honey bee. The symbiont shares common genomic and predicted metabolic properties with the male-killing symbiont Arsenophonus nasoniae, however we present multiple lines of evidence that the bee Arsenophonus deviates from a heritable model of transmission. Field sampling uncovered spatial and seasonal dynamics in symbiont prevalence, and rapid infection loss events were observed in field colonies and laboratory individuals. Fluorescent in situ hybridisation showed Arsenophonus localised in the gut, and detection was rare in screens of early honey bee life stages. We directly show horizontal transmission of Arsenophonus between bees under varying social conditions. We conclude that honey bees acquire Arsenophonus through a combination of environmental exposure and social contacts. These findings uncover a key link in the Arsenophonus clades trajectory from free-living ancestral life to obligate mutualism, and provide a foundation for studying transitions in symbiotic lifestyle.

Chronic bee paralysis as a serious emerging threat to honey bees.

Chronic bee paralysis is a well-defined viral disease of honey bees with a global distribution that until recently caused rare but severe symptomatology including colony loss. Anecdotal evidence indicates a recent increase in virus incidence in several countries, but no mention of concomitant disease. We use government honey bee health inspection records from England and Wales to test whether chronic bee paralysis is an emerging infectious disease and investigate the spatiotemporal patterns of disease. The number of chronic bee paralysis cases increased exponentially between 2007 and 2017, demonstrating chronic bee paralysis as an emergent disease. Disease is highly clustered spatially within most years, suggesting local spread, but not between years, suggesting disease burnt out with periodic reintroduction. Apiary and county level risk factors are confirmed to include scale of beekeeping operation and the history of honey bee imports. Our findings offer epidemiological insight into this damaging emerging disease.

Characterisation of the British honey bee metagenome.

The European honey bee (Apis mellifera) plays a major role in pollination and food production. Honey bee health is a complex product of the environment, host genetics and associated microbes (commensal, opportunistic and pathogenic). Improved understanding of these factors will help manage modern challenges to bee health. Here we used DNA sequencing to characterise the genomes and metagenomes of 19 honey bee colonies from across Britain. Low heterozygosity was observed in many Scottish colonies which had high similarity to the native dark bee. Colonies exhibited high diversity in composition and relative abundance of individual microbiome taxa. Most non-bee sequences were derived from known honey bee commensal bacteria or pathogens. However, DNA was also detected from additional fungal, protozoan and metazoan species. To classify cobionts lacking genomic information, we developed a novel network analysis approach for clustering orphan DNA contigs. Our analyses shed light on microbial communities associated with honey bees and demonstrate the power of high-throughput, directed metagenomics for identifying novel biological threats in agroecosystems.

The invasion, provenance and diversity of Vespa velutina Lepeletier (Hymenoptera: Vespidae) in Great Britain.

The yellow-legged or Asian hornet (Vespa velutina colour form nigrithorax) was introduced into France from China over a decade ago. Vespa velutina has since spread rapidly across Europe, facilitated by suitable climatic conditions and the ability of a single nest to disperse many mated queens over a large area. Yellow-legged hornets are a major concern because of the potential impact they have on populations of many beneficial pollinators, most notably the western honey bee (Apis mellifera), which shows no effective defensive behaviours against this exotic predator. Here, we present the first report of this species in Great Britain. Actively foraging hornets were detected at two locations, the first around a single nest in Gloucestershire, and the second a single hornet trapped 54 km away in Somerset. The foraging activity observed in Gloucestershire was largely restricted to within 700 m of a single nest, suggesting highly localised movements. Genetic analyses of individuals from the Gloucestershire nest and the single hornet from Somerset suggest that these incursions represent an expansion of the European population, rather than a second incursion from Asia. The founding queen of the Gloucestershire nest mated with a single male, suggesting that sexual reproduction may have occurred in an area of low nest density. Whilst the nest contained diploid adult males, haploid 'true' males were only present at the egg stage, indicating that the nest was detected and removed before the production of queens. Members of the public reported additional dead hornets associated with camping equipment recently returned from France and imported timber products, highlighting possible pathways of incursion. The utility of microsatellites to inform surveillance during an incursion and the challenge of achieving eradication of this damaging pest are discussed.

Predicting the spread of the Asian hornet (Vespa velutina) following its incursion into Great Britain.

The yellow-legged or Asian hornet (Vespa velutina) is native to South-East Asia, and is a voracious predator of pollinating insects including honey bees. Since its accidental introduction into South-Western France in 2004, V. velutina has spread to much of western Europe. The presence of V. velutina in Great Britain was first confirmed in September 2016. The likely dynamics following an initial incursion are uncertain, especially the risk of continued spread, and the likely success of control measures. Here we extrapolate from the situation in France to quantify the potential invasion of V. velutina in Great Britain. We find that, without control, V. velutina could colonise the British mainland rapidly, depending upon how the Asian hornet responds to the colder climate in Britain compared to France. The announcement that a second hornet had been discovered in Somerset, increases the chance that the invasion first occurred before 2016. We therefore consider the likely site of first invasion and the probabilistic position of additional founding nests in late 2016 and early 2017. Given the potential dispersion of V. velutina, we conclude that vigilance is required over a large area to prevent the establishment of this threat to the pollinator population.

Invasion dynamics of Asian hornet, Vespa velutina (Hymenoptera: Vespidae): a case study of a commune in south-west France.

Asian hornet, Vespa velutina Lepeletier nests were discovered in 2007 in Andernos-les-Bains on the south-west coast of France, 3 years after the first reported sightings in France. The number of nests increased in the commune over the following 7 years, despite local authorities enacting a destruction policy. The nests existed in close proximity to one another leading to a high density of over 10 nests per square kilometre in urban areas. New information on the chosen habitat for nests is presented, and the differences between primary and secondary locations are evident, with primary nests mostly occupying buildings and man-made structures, while secondary nests were found on trees. Using Bayesian inference methods, we fit a basic model to the observational data, which allows us to estimate key demographic parameters. This model fit is highly informative for predicting V. velutina spread and colonisation of other at-risk regions, and suggests that local control has a limited impact on the spread of V. velutina once established within a region.

Efficient use of sentinel sites: detection of invasive honeybee pests and diseases in the UK.

Sentinel sites, where problems can be identified early or investigated in detail, form an important part of planning for exotic disease outbreaks in humans, livestock and plants. Key questions are: how many sentinels are required, where should they be positioned and how effective are they at rapidly identifying new invasions? The sentinel apiary system for invasive honeybee pests and diseases illustrates the costs and benefits of such approaches. Here, we address these issues with two mathematical modelling approaches. The first approach is generic and uses probabilistic arguments to calculate the average number of affected sites when an outbreak is first detected, providing rapid and general insights that we have applied to a range of infectious diseases. The second approach uses a computationally intensive, stochastic, spatial model to simulate multiple outbreaks and to determine appropriate sentinel locations for UK apiaries. Both models quantify the anticipated increase in success of sentinel sites as their number increases and as non-sentinel sites become worse at detection; however, unexpectedly sentinels perform relatively better for faster growing outbreaks. Additionally, the spatial model allows us to quantify the substantial role that carefully positioned sentinels can play in the rapid detection of exotic invasions.

A pan-European epidemiological study reveals honey bee colony survival depends on beekeeper education and disease control.

Reports of honey bee population decline has spurred many national efforts to understand the extent of the problem and to identify causative or associated factors. However, our collective understanding of the factors has been hampered by a lack of joined up trans-national effort. Moreover, the impacts of beekeeper knowledge and beekeeping management practices have often been overlooked, despite honey bees being a managed pollinator. Here, we established a standardised active monitoring network for 5 798 apiaries over two consecutive years to quantify honey bee colony mortality across 17 European countries. Our data demonstrate that overwinter losses ranged between 2% and 32%, and that high summer losses were likely to follow high winter losses. Multivariate Poisson regression models revealed that hobbyist beekeepers with small apiaries and little experience in beekeeping had double the winter mortality rate when compared to professional beekeepers. Furthermore, honey bees kept by professional beekeepers never showed signs of disease, unlike apiaries from hobbyist beekeepers that had symptoms of bacterial infection and heavy Varroa infestation. Our data highlight beekeeper background and apicultural practices as major drivers of honey bee colony losses. The benefits of conducting trans-national monitoring schemes and improving beekeeper training are discussed.

Microsporidia infection impacts the host cell's cycle and reduces host cell apoptosis.

Intracellular parasites can alter the cellular machinery of host cells to create a safe haven for their survival. In this regard, microsporidia are obligate intracellular fungal parasites with extremely reduced genomes and hence, they are strongly dependent on their host for energy and resources. To date, there are few studies into host cell manipulation by microsporidia, most of which have focused on morphological aspects. The microsporidia Nosema apis and Nosema ceranae are worldwide parasites of honey bees, infecting their ventricular epithelial cells. In this work, quantitative gene expression and histology were studied to investigate how these two parasites manipulate their host's cells at the molecular level. Both these microsporidia provoke infection-induced regulation of genes involved in apoptosis and the cell cycle. The up-regulation of buffy (which encodes a pro-survival protein) and BIRC5 (belonging to the Inhibitor Apoptosis protein family) was observed after infection, shedding light on the pathways that these pathogens use to inhibit host cell apoptosis. Curiously, different routes related to cell cycle were modified after infection by each microsporidia. In the case of N. apis, cyclin B1, dacapo and E2F2 were up-regulated, whereas only cyclin E was up-regulated by N. ceranae, in both cases promoting the G1/S phase transition. This is the first report describing molecular pathways related to parasite-host interactions that are probably intended to ensure the parasite's survival within the cell.

Variability of residue concentrations of ciprofloxacin in honey from treated hives.

Honey bees (Apis mellifera L.) were treated with a model veterinary drug compound (ciprofloxacin) in a 3-year study (2012-14) to investigate the variability of residue concentration in honey. Sucrose solution containing ciprofloxacin was administered to 45 hives (1 g of ciprofloxacin per hive) at the beginning of the honey flow in late May/mid-June 2012, 2013 and 2014. Buckfast honey bees (A. mellifera - hybrid) were used in years 2012 and 2013. Carniolan honey bees (A. mellifera carnica) were used instead of the Buckfast honey bees as a replacement due to unforeseen circumstances in the final year of the study (2014). Honey was collected over nine scheduled time points from May/June till late October each year. Up to five hives were removed and their honey analysed per time point. Honey samples were analysed by liquid chromatography-tandem mass spectrometry (LC-MS/MS) to determine ciprofloxacin concentration. Statistical assessment of the data shows that the inter-hive variation of ciprofloxacin concentrations in 2012/13 is very different compared with that of 2014 with relative standard deviations (RSDs) of 138% and 61%, respectively. The average ciprofloxacin concentration for 2014 at the last time point was more than 10 times the concentration compared with samples from 2012/13 at the same time point. The difference between the 2012/13 data compared with the 2014 data is likely due to the different type of honey bees used in this study (2012/13 Buckfast versus 2014 Carniolan). Uncertainty estimates for honey with high ciprofloxacin concentration (upper 95th percentile) across all hives for 55-day withdrawal samples gave residual standard errors (RSEs) of 22%, 20% and 11% for 2012, 2013 and 2014, respectively. If the number of hives were to be reduced for future studies, RSEs were estimated to be 52% (2012), 54% (2013) and 26% (2014) for one hive per time point (nine total hives).

Identifying bacterial predictors of honey bee health.

Non-targeted approaches are useful tools to identify new or emerging issues in bee health. Here, we utilise next generation sequencing to highlight bacteria associated with healthy and unhealthy honey bee colonies, and then use targeted methods to screen a wider pool of colonies with known health status. Our results provide the first evidence that bacteria from the genus Arsenophonus are associated with poor health in honey bee colonies. We also discovered Lactobacillus and Leuconostoc spp. were associated with healthier honey bee colonies. Our results highlight the importance of understanding how the wider microbial population relates to honey bee colony health.

Transcriptome analysis of the synganglion from the honey bee mite, Varroa destructor and RNAi knockdown of neural peptide targets.

Varroa mites (Varroa destructor) and the viruses that they transmit are one of the major contributing factors to the global honey bee crisis. Gene products within the nervous system are the targets of all the insecticides currently used to control Varroa but there is a paucity of transcriptomic data available for Varroa neural tissues. A cDNA library from the synganglia ("brains") of adult female Varroa was constructed and 600 ESTs sequenced and analysed revealing several current and potential druggable targets. Contigs coding for the deformed wing virus (DWV) variants V. destructor virus-1 (VDV-1) and the recombinant (VDV-1DVD) were present in the synganglion library. Negative-sense RNA-specific PCR indicated that VDV-1 replicates in the Varroa synganglion and all other tissues tested, but we could not detect DWV replicating in any Varroa tissue. Two neuropeptides were identified in the synganlion EST library: a B-type allatostatin and a member of the crustacean hyperglycaemic hormone (CHH) superfamily. Knockdown of the allatostatin or the CHH-like gene by double-stranded RNA-interference (dsRNAi) resulted in 85% and 55% mortality, respectively, of Varroa. Here, we present the first transcriptomic survey in Varroa and demonstrate that neural genes can be targeted by dsRNAi either for genetic validation of putative targets during drug discovery programmes or as a potential control measure in itself.

Honey Bee Colonies Headed by Hyperpolyandrous Queens Have Improved Brood Rearing Efficiency and Lower Infestation Rates of Parasitic Varroa Mites.

A honey bee queen mates on wing with an average of 12 males and stores their sperm to produce progeny of mixed paternity. The degree of a queen's polyandry is positively associated with measures of her colony's fitness, and observed distributions of mating number are evolutionary optima balancing risks of mating flights against benefits to the colony. Effective mating numbers as high as 40 have been documented, begging the question of the upper bounds of this behavior that can be expected to confer colony benefit. In this study we used instrumental insemination to create three classes of queens with exaggerated range of polyandry--15, 30, or 60 drones. Colonies headed by queens inseminated with 30 or 60 drones produced more brood per bee and had a lower proportion of samples positive for Varroa destructor mites than colonies whose queens were inseminated with 15 drones, suggesting benefits of polyandry at rates higher than those normally obtaining in nature. Our results are consistent with two hypotheses that posit conditions that reward such high expressions of polyandry: (1) a queen may mate with many males in order to promote beneficial non-additive genetic interactions among subfamilies, and (2) a queen may mate with many males in order to capture a large number of rare alleles that regulate resistance to pathogens and parasites in a breeding population. Our results are unique for identifying the highest levels of polyandry yet detected that confer colony-level benefit and for showing a benefit of polyandry in particular toward the parasitic mite V. destructor.

Pathogens as Predictors of Honey Bee Colony Strength in England and Wales.

Inspectors with the UK National Bee Unit were asked for 2007-2008 to target problem apiaries in England and Wales for pathogen screening and colony strength measures. Healthy colonies were included in the sampling to provide a continuum of health conditions. A total of 406 adult bee samples was screened and yielded 7 viral, 1 bacterial, and 2 microsporidial pathogens and 1 ectoparasite (Acarapis woodi). In addition, 108 samples of brood were screened and yielded 4 honey bee viruses. Virus prevalence varied from common (deformed wing virus, black queen cell virus) to complete absence (Israeli acute paralysis virus). When colonies were forced into one of two classes, strong or weak, the weak colonies contained more pathogens in adult bees. Among observed pathogens, only deformed wing virus was able to predict colony strength. The effect was negative such that colonies testing positive for deformed wing virus were likely to have fewer combs of bees or brood. This study constitutes the first record for Nosema ceranae in Great Britain. These results contribute to the growing body of evidence linking pathogens to poor honey bee health.

Biogeography of Paenibacillus larvae, the causative agent of American foulbrood, using a new multilocus sequence typing scheme.

American foulbrood is the most destructive brood disease of honeybees (Apis mellifera) globally. The absence of a repeatable, universal typing scheme for the causative bacterium Paenibacillus larvae has restricted our understanding of disease epidemiology. We have created the first multilocus sequence typing scheme (MLST) for P. larvae, which largely confirms the previous enterobacterial repetitive intergenic consensus (ERIC)-polymerase chain reaction (PCR)-based typing scheme's divisions while providing added resolution and improved repeatability. We have used the new scheme to determine the distribution and biogeography of 294 samples of P. larvae from across six continents. We found that of the two most epidemiologically important ERIC types, ERIC I was more diverse than ERIC II. Analysis of the fixation index (FST ) by distance suggested a significant relationship between genetic and geographic distance, suggesting that population structure exists in populations of P. larvae. Interestingly, this effect was only observed within the native range of the host and was absent in areas where international trade has moved honeybees and their disease. Correspondence analysis demonstrated similar sequence type (ST) distributions between native and non-native countries and that ERIC I and II STs mainly have differing distributions. The new typing scheme facilitates epidemiological study of this costly disease of a key pollinator.

Parasite pressures on feral honey bees (Apis mellifera sp.).

Feral honey bee populations have been reported to be in decline due to the spread of Varroa destructor, an ectoparasitic mite that when left uncontrolled leads to virus build-up and colony death. While pests and diseases are known causes of large-scale managed honey bee colony losses, no studies to date have considered the wider pathogen burden in feral colonies, primarily due to the difficulty in locating and sampling colonies, which often nest in inaccessible locations such as church spires and tree tops. In addition, little is known about the provenance of feral colonies and whether they represent a reservoir of Varroa tolerant material that could be used in apiculture. Samples of forager bees were collected from paired feral and managed honey bee colonies and screened for the presence of ten honey bee pathogens and pests using qPCR. Prevalence and quantity was similar between the two groups for the majority of pathogens, however feral honey bees contained a significantly higher level of deformed wing virus than managed honey bee colonies. An assessment of the honey bee race was completed for each colony using three measures of wing venation. There were no apparent differences in wing morphometry between feral and managed colonies, suggesting feral colonies could simply be escapees from the managed population. Interestingly, managed honey bee colonies not treated for Varroa showed similar, potentially lethal levels of deformed wing virus to that of feral colonies. The potential for such findings to explain the large fall in the feral population and the wider context of the importance of feral colonies as potential pathogen reservoirs is discussed.