Shift in virus composition in honeybees (Apis mellifera) following worldwide invasion by the parasitic mite and virus vector Varroa destructor.

Invasive vectors can induce dramatic changes in disease epidemiology. While viral emergence following geographical range expansion of a vector is well known, the influence a vector can have at the level of the host's pathobiome is less well understood. Taking advantage of the formerly heterogeneous spatial distribution of the ectoparasitic mite Varroa destructor that acts as potent virus vector among honeybees Apis mellifera, we investigated the impact of its recent global spread on the viral community of honeybees in a retrospective study of historical samples. We hypothesized that the vector has had an effect on the epidemiology of several bee viruses, potentially altering their transmissibility and/or virulence, and consequently their prevalence, abundance, or both. To test this, we quantified the prevalence and loads of 14 viruses from honeybee samples collected in mite-free and mite-infested populations in four independent geographical regions. The presence of the mite dramatically increased the prevalence and load of deformed wing virus, a cause of unsustainably high colony losses. In addition, several other viruses became more prevalent or were found at higher load in mite-infested areas, including viruses not known to be actively varroa-transmitted, but which may increase opportunistically in varroa-parasitized bees.

Virus infections in honeybee colonies naturally surviving ectoparasitic mite vectors.

Western honeybee populations, Apis mellifera, in Europe have been known to survive infestations of the ectoparasitic mite, Varroa destructor, by means of natural selection. Proposed mechanisms in literature have been focused on the management of this parasite, however literature remains scare on the differences in viral ecology between colonies that have adapted to V. destructor and those that are consistently treated for it. Samples were collected from both a mite-surviving and a sympatric mite-susceptible honeybee population in Norway. The prevalence and abundances of 10 viruses, vectored by the parasite or not, were investigated in adult host workers and pupae as well as in V. destructor mites. Here we show that the mite-vectored Deformed wing virus (DWV-A) is often lower in both abundance and prevalence in the mite-surviving population in tandem with lower phoretic mite infestations compared to the mite susceptible population. However, the non-mite-vectored Black queen cell virus (BQCV), had both a higher abundance and prevalence in the mite-surviving population compared to the susceptible population. The data therefore suggest that general adaptations to virus infections may be unlikely to explain colony survival. Instead, mechanisms suppressing mite reproduction and therefore the impact seem to be more important.

Host brood traits, independent of adult behaviours, reduce Varroa destructor mite reproduction in resistant honeybee populations.

The ectoparasitic mite Varroa destructor is an invasive species of Western honey bees (Apis mellifera) and the largest pathogenic threat to their health world-wide. Its successful invasion and expansion is related to its ability to exploit the worker brood for reproduction, which results in an exponential population growth rate in the new host. With invasion of the mite, wild honeybee populations have been nearly eradicated from Europe and North America, and the survival of managed honeybee populations relies on mite population control treatments. However, there are a few documented honeybee populations surviving extended periods without control treatments due to adapted host traits that directly impact Varroa mite fitness. The aim of this study was to investigate if Varroa mite reproductive success was affected by traits of adult bee behaviours or by traits of the worker brood, in three mite-resistant honey bee populations from Sweden, France and Norway. The mite's reproductive success was measured and compared in broods that were either exposed to, or excluded from, adult bee access. Mite-resistant bee populations were also compared with a local mite-susceptible population, as a control group. Our results show that mite reproductive success rates and mite fecundity in the three mite-resistant populations were significantly different from the control population, with the French and Swedish populations having significantly lower reproductive rates than the Norwegian population. When comparing mite reproduction in exposed or excluded brood treatments, no differences were observed, regardless of population. This result clearly demonstrates that Varroa mite reproductive success can be suppressed by traits of the brood, independent of adult worker bees.

Cutting corners: The impact of storage and DNA extraction on quality and quantity of DNA in honeybee (Apis mellifera) spermatheca.

The purpose of our study was to investigate methods of short-term storage that allow preservation, transport and retrieval of genetic information contained in honeybee queen's spermatheca. Genotyping of the honeybee colony requires well ahead planned sample collection, depending on the type of data to be acquired. Sampling and genotyping of spermatheca's content instead of individual offspring is timesaving, allowing answers to the questions related to patriline composition immediately after mating. Such procedure is also cheaper and less error prone. For preservation either Allprotect Tissue Reagent (Qiagen) or absolute ethanol were used. Conditions during transportation were simulated by keeping samples 6-8 days at room temperature. Six different storing conditions of spermathecas were tested, complemented with two DNA extraction methods. We have analysed the concentration of DNA, RNA, and proteins in DNA extracts. We also analysed how strongly the DNA is subjected to fragmentation (through amplification of genetic markers ANT2 and tRNAleu-COX2) and whether the quality of the extracted DNA is suitable for microsatellite (MS) analysis. Then, we tested the usage of spermatheca as a source of patriline composition in an experiment with three instrumentally inseminated virgin queens and performed MS analysis of the extracted DNA from each spermatheca, as well as queens' and drones' tissue. Our results show that median DNA concentration from spermathecas excised prior the storage, regardless of the storing condition and DNA extraction method, were generally lower than median DNA concentration obtained from spermathecas dissected from the whole queens after the storage. Despite the differences in DNA yield from the samples subjected to different storing conditions there was no significant effect of storage method or the DNA extraction method on the amplification success, although fewer samples stored in EtOH amplified successfully in comparison to ATR storing reagent. However, we recommend EtOH as a storing reagent due to its availability, low price, simplicity in usage in the field and in the laboratory, and capability of good preservation of the samples for DNA analysis during transport at room temperature.

Virus Prevalence in Egg Samples Collected from Naturally Selected and Traditionally Managed Honey Bee Colonies across Europe.

Monitoring virus infections can be an important selection tool in honey bee breeding. A recent study pointed towards an association between the virus-free status of eggs and an increased virus resistance to deformed wing virus (DWV) at the colony level. In this study, eggs from both naturally surviving and traditionally managed colonies from across Europe were screened for the prevalence of different viruses. Screenings were performed using the phenotyping protocol of the 'suppressed in ovo virus infection' trait but with qPCR instead of end-point PCR and a primer set that covers all DWV genotypes. Of the 213 screened samples, 109 were infected with DWV, 54 were infected with black queen cell virus (BQCV), 3 were infected with the sacbrood virus, and 2 were infected with the acute bee paralyses virus. It was demonstrated that incidences of the vertical transmission of DWV were more frequent in naturally surviving than in traditionally managed colonies, although the virus loads in the eggs remained the same. When comparing virus infections with queen age, older queens showed significantly lower infection loads of DWV in both traditionally managed and naturally surviving colonies, as well as reduced DWV infection frequencies in traditionally managed colonies. We determined that the detection frequencies of DWV and BQCV in honey bee eggs were lower in samples obtained in the spring than in those collected in the summer, indicating that vertical transmission may be lower in spring. Together, these patterns in vertical transmission show that honey bee queens have the potential to reduce the degree of vertical transmission over time.

Identification of 121 variants of honey bee Vitellogenin protein sequences with structural differences at functional sites.

Proteins are under selection to maintain central functions and to accommodate needs that arise in ever-changing environments. The positive selection and neutral drift that preserve functions result in a diversity of protein variants. The amount of diversity differs between proteins: multifunctional or disease-related proteins tend to have fewer variants than proteins involved in some aspects of immunity. Our work focuses on the extensively studied protein Vitellogenin (Vg), which in honey bees (Apis mellifera) is multifunctional and highly expressed and plays roles in immunity. Yet, almost nothing is known about the natural variation in the coding sequences of this protein or how amino acid-altering variants might impact structure-function relationships. Here, we map out allelic variation in honey bee Vg using biological samples from 15 countries. The successful barcoded amplicon Nanopore sequencing of 543 bees revealed 121 protein variants, indicating a high level of diversity in Vg. We find that the distribution of non-synonymous single nucleotide polymorphisms (nsSNPs) differs between protein regions with different functions; domains involved in DNA and protein-protein interactions contain fewer nsSNPs than the protein's lipid binding cavities. We outline how the central functions of the protein can be maintained in different variants and how the variation pattern may inform about selection from pathogens and nutrition.

CSI Pollen: Diversity of Honey Bee Collected Pollen Studied by Citizen Scientists.

A diverse supply of pollen is an important factor for honey bee health, but information about the pollen diversity available to colonies at the landscape scale is largely missing. In this COLOSS study, beekeeper citizen scientists sampled and analyzed the diversity of pollen collected by honey bee colonies. As a simple measure of diversity, beekeepers determined the number of colors found in pollen samples that were collected in a coordinated and standardized way. Altogether, 750 beekeepers from 28 different regions from 24 countries participated in the two-year study and collected and analyzed almost 18,000 pollen samples. Pollen samples contained approximately six different colors in total throughout the sampling period, of which four colors were abundant. We ran generalized linear mixed models to test for possible effects of diverse factors such as collection, i.e., whether a minimum amount of pollen was collected or not, and habitat type on the number of colors found in pollen samples. To identify habitat effects on pollen diversity, beekeepers' descriptions of the surrounding landscape and CORINE land cover classes were investigated in two different models, which both showed that both the total number and the rare number of colors in pollen samples were positively affected by 'urban' habitats or 'artificial surfaces', respectively. This citizen science study underlines the importance of the habitat for pollen diversity for bees and suggests higher diversity in urban areas.

Using Citizen Science to Scout Honey Bee Colonies That Naturally Survive Varroa destructor Infestations.

Citizen Science contributes significantly to the conservation of biodiversity, but its application to honey bee research has remained minimal. Even though certain European honey bee (Apis mellifera) populations are known to naturally survive Varroa destructor infestations, it is unclear how widespread or common such populations are. Such colonies are highly valuable for investigating the mechanisms enabling colony survival, as well as for tracking the conservation status of free-living honey bees. Here, we use targeted Citizen Science to identify potentially new cases of managed or free-living A. mellifera populations that survive V. destructor without mite control strategies. In 2018, a survey containing 20 questions was developed, translated into 13 languages, and promoted at beekeeping conferences and online. After three years, 305 reports were collected from 28 countries: 241 from managed colonies and 64 from free-living colonies. The collected data suggest that there could be twice as many naturally surviving colonies worldwide than are currently known. Further, online and personal promotion seem to be key for successful recruitment of participants. Although the survivor status of these colonies still needs to be confirmed, the volume of reports and responses already illustrate how effectively Citizen Science can contribute to bee research by massively increasing generated data, broadening opportunities for comparative research, and fostering collaboration between scientists, beekeepers, and citizens. The success of this survey spurred the development of a more advanced Citizen Science platform, Honey Bee Watch, that will enable a more accurate reporting, confirmation, and monitoring of surviving colonies, and strengthen the ties between science, stakeholders, and citizens to foster the protection of both free-living and managed honey bees.

Adaptive population structure shifts in invasive parasitic mites, Varroa destructor.

Comparative studies of genetic diversity and population structure can shed light on the ecological and evolutionary factors governing host-parasite interactions. Even though invasive parasites are considered of major biological importance, little is known about their adaptative potential when infesting the new hosts. Here, the genetic diversification of Varroa destructor, a novel parasite of Apis mellifera originating from Asia, was investigated using population genetics to determine how the genetic structure of the parasite changed in distinct European populations of its new host. To do so, mites infesting two categories of hosts in four European regions were compared: (a) adapted hosts surviving through means of natural selection, thereby expected to impose strong selective pressure on the mites, and (b) treated host populations, surviving mite infestations because acaricides are applied, therefore characterized by a relaxed selection imposed by the host on the mites. Significant genetic divergence was found across regions, partially reflecting the invasion pattern of V. destructor throughout Europe and indicating local adaptation of the mite to the host populations. Additionally, varying degrees of genotypic changes were found between mites from adapted and treated colonies. Altogether, these results indicate that V. destructor managed to overcome the genetic bottlenecks following its introduction in Europe and that host-mediated selection fostered changes in the genetic structure of this mite at diverse geographic scales. These findings highlight the potential of parasites to adapt to their local host populations and confirm that adaptations developed within coevolutionary dynamics are a major determinant of population genetic changes.

Adapted tolerance to virus infections in four geographically distinct Varroa destructor-resistant honeybee populations.

The ectoparasitic mite, Varroa destructor, is unarguably the leading cause of honeybee (Apis mellifera) mortality worldwide through its role as a vector for lethal viruses, in particular, strains of the Deformed wing virus (DWV) and Acute bee paralysis virus (ABPV) complexes. Several honeybee populations across Europe have well-documented adaptations of mite-resistant traits but little is known about host adaptations towards the virus infections vectored by the mite. The aim of this study was to assess and compare the possible contribution of adapted virus tolerance and/or resistance to the enhanced survival of four well-documented mite-resistant honeybee populations from Norway, Sweden, The Netherlands and France, in relation to unselected mite-susceptible honeybees. Caged adult bees and laboratory reared larvae, from colonies of these four populations, were inoculated with DWV and ABPV in a series of feeding infection experiments, while control groups received virus-free food. Virus infections were monitored using RT-qPCR assays in individuals sampled over a time course. In both adults and larvae the DWV and ABPV infection dynamics were nearly identical in all groups, but all mite-resistant honeybee populations had significantly higher survival rates compared to the mite-susceptible honeybees. These results suggest that adapted virus tolerance is an important component of survival mechanisms.

Reproductive success of the parasitic mite (Varroa destructor) is lower in honeybee colonies that target infested cells with recapping.

Cell recapping is a behavioural trait of honeybees (Apis mellifera) where cells with developing pupae are uncapped, inspected, and then recapped, without removing the pupae. The ectoparasitic mite Varroa destructor, unarguably the most destructive pest in apiculture world-wide, invades the cells of developing pupae to feed and reproduce. Honeybees that target mite infested cells with this behaviour may disrupt the reproductive cycle of the mite. Hence, cell recapping has been associated with colony-level declines in mite reproduction. In this study we compared the colony-level efficacy of cell recapping (how often infested cells are recapped) to the average mite fecundity in A. mellifera. Our study populations, known to be adapted to V. destructor, were from Avignon, France, Gotland, Sweden, and Oslo, Norway, and were compared to geographically similar, treated control colonies. The results show that colonies with a higher recapping efficacy also have a lower average mite reproductive success. This pattern was likely driven by the adapted populations as they had the largest proportion of highly-targeted cell recapping. The consistent presence of this trait in mite-resistant and mite-susceptible colonies with varying degrees of expression may make it a good proxy trait for selective breeding on a large scale.

Local Variation in Recombination Rates of the Honey Bee (Apis mellifera) Genome among Samples from Six Disparate Populations.

Meiotic recombination is an essential component of eukaryotic sexual reproduction but its frequency varies within and between genomes. Although it is well-established that honey bees have a high recombination rate with about 20 cM/Mbp, the proximate and ultimate causes of this exceptional rate are poorly understood. Here, we describe six linkage maps of the Western Honey Bee Apis mellifera that were produced with consistent methodology from samples from distinct parts of the species' near global distribution. We compared the genome-wide rates and distribution of meiotic crossovers among the six maps and found considerable differences. Overall similarity of local recombination rates among our samples was unrelated to geographic or phylogenetic distance of the populations that our samples were derived from. However, the limited sampling constrains the interpretation of our results because it is unclear how representative these samples are. In contrast to previous studies, we found only in two datasets a significant relation between local recombination rate and GC content. Focusing on regions of particularly increased or decreased recombination in specific maps, we identified several enriched gene ontologies in these regions and speculate about their local adaptive relevance. These data are contributing to an increasing comparative effort to gain an understanding of the intra-specific variability of recombination rates and their evolutionary role in honey bees and other social insects.

Evaluation of Suppressed Mite Reproduction (SMR) Reveals Potential for Varroa Resistance in European Honey Bees (Apis mellifera L.).

In the fight against the Varroa destructor mite, selective breeding of honey bee (Apis mellifera L.) populations that are resistant to the parasitic mite stands as a sustainable solution. Selection initiatives indicate that using the suppressed mite reproduction (SMR) trait as a selection criterion is a suitable tool to breed such resistant bee populations. We conducted a large European experiment to evaluate the SMR trait in different populations of honey bees spread over 13 different countries, and representing different honey bee genotypes with their local mite parasites. The first goal was to standardize and validate the SMR evaluation method, and then to compare the SMR trait between the different populations. Simulation results indicate that it is necessary to examine at least 35 single-infested cells to reliably estimate the SMR score of any given colony. Several colonies from our dataset display high SMR scores indicating that this trait is present within the European honey bee populations. The trait is highly variable between colonies and some countries, but no major differences could be identified between countries for a given genotype, or between genotypes in different countries. This study shows the potential to increase selective breeding efforts of V. destructor resistant populations.

Honey bee predisposition of resistance to ubiquitous mite infestations.

Host-parasite co-evolution history is lacking when parasites switch to novel hosts. This was the case for Western honey bees (Apis mellifera) when the ectoparasitic mite, Varroa destructor, switched hosts from Eastern honey bees (Apis cerana). This mite has since become the most severe biological threat to A. mellifera worldwide. However, some A. mellifera populations are known to survive infestations, largely by suppressing mite population growth. One known mechanism is suppressed mite reproduction (SMR), but the underlying genetics are poorly understood. Here, we take advantage of haploid drones, originating from one queen from the Netherlands that developed Varroa-resistance, whole exome sequencing and elastic-net regression to identify genetic variants associated with SMR in resistant honeybees. An eight variants model predicted 88% of the phenotypes correctly and identified six risk and two protective variants. Reproducing and non-reproducing mites could not be distinguished using DNA microsatellites, which is in agreement with the hypothesis that it is not the parasite but the host that adapted itself. Our results suggest that the brood pheromone-dependent mite oogenesis is disrupted in resistant hosts. The identified genetic markers have a considerable potential to contribute to a sustainable global apiculture.

Reduced Postcapping Period in Honey Bees Surviving Varroa destructor by Means of Natural Selection.

The ectoparasitic mite Varroa destructor is a key factor for colony losses in European honey bee subspecies (Apis mellifera), but it is also known that some host populations have adapted to the mite by means of natural selection. The role of a shorter host brood postcapping period in reducing mite reproductive success has been investigated in other surviving subspecies, however its role in the adaptation of European honey bee populations has not been addressed. Here, we use a common garden approach to compare the length of the worker brood postcapping period in a Norwegian surviving honey bee population with the postcapping period of a local susceptible population. The data show a significantly shorter postcapping period in the surviving population for ~10% of the brood. Since even small differences in postcapping period can significantly reduce mite reproductive success, this mechanism may well contribute to natural colony survival. It appears most likely that several mechanisms acting together produce the full mite-surviving colony phenotype.

Rapid parallel evolution overcomes global honey bee parasite.

In eusocial insect colonies nestmates cooperate to combat parasites, a trait called social immunity. However, social immunity failed for Western honey bees (Apis mellifera) when the ectoparasitic mite Varroa destructor switched hosts from Eastern honey bees (Apis cerana). This mite has since become the most severe threat to A. mellifera world-wide. Despite this, some isolated A. mellifera populations are known to survive infestations by means of natural selection, largely by supressing mite reproduction, but the underlying mechanisms of this are poorly understood. Here, we show that a cost-effective social immunity mechanism has evolved rapidly and independently in four naturally V. destructor-surviving A. mellifera populations. Worker bees of all four 'surviving' populations uncapped/recapped worker brood cells more frequently and targeted mite-infested cells more effectively than workers in local susceptible colonies. Direct experiments confirmed the ability of uncapping/recapping to reduce mite reproductive success without sacrificing nestmates. Our results provide striking evidence that honey bees can overcome exotic parasites with simple qualitative and quantitative adaptive shifts in behaviour. Due to rapid, parallel evolution in four host populations this appears to be a key mechanism explaining survival of mite infested colonies.

Molecular methods indicate lack of spread of Acarapis woodi introduced to honey bees in western Norway.

The tracheal mite, Acarapis woodi, may be one of many factors contributing to the decline in honey bee (Apis mellifera) populations. Databases on the widespread distribution of A. woodi exist, but the data seem patchy. Norway is not listed as being infested, although there have been at least two separate introductions of the parasite. Investigations in 2003, 2006, and 2009 using standard microscopy methods indicated persistence of A. woodi in honey bees in this region. In 2013, we conducted another survey. Samples were sent in from 335 beehives belonging to 39 apiaries, and all were asked to complete a questionnaire. Analysis for A. woodi in the submitted samples was by PCR, with sequencing of positive results. The results described in this article indicate that this parasite still persists in some apiaries in this region, but at a low, and possibly decreasing, level, with positive results obtained from just two (5.1%) of the apiaries. Of the 17 beekeepers that answered the questionnaire, none reported symptoms of infestation with A. woodi. Sequencing of PCR products indicated a difference between the two A. woodi isolates. Our results were generally encouraging regarding the apparent lack of spread of A. woodi, within the County. Furthermore, the sequencing results may indicate two separate introductions rather than spread. Nevertheless, the scarcity of data, the vulnerability of honey bee populations globally, and the potential contribution of this parasite to reduced survival, indicate that the situation should be continued to be monitored. In addition, Norwegian beekeepers should be made aware of, and follow, restrictions regarding import and transport of bees, both nationally and internationally.

Norwegian honey bees surviving Varroa destructor mite infestations by means of natural selection.

BACKGROUND: Managed, feral and wild populations of European honey bee subspecies, Apis mellifera, are currently facing severe colony losses globally. There is consensus that the ectoparasitic mite Varroa destructor, that switched hosts from the Eastern honey bee Apis cerana to the Western honey bee A. mellifera, is a key factor driving these losses. For >20 years, breeding efforts have not produced European honey bee colonies that can survive infestations without the need for mite control. However, at least three populations of European honey bees have developed this ability by means of natural selection and have been surviving for >10 years without mite treatments. Reduced mite reproductive success has been suggested as a key factor explaining this natural survival. Here, we report a managed A. mellifera population in Norway, that has been naturally surviving consistent V. destructor infestations for >17 years. METHODS: Surviving colonies and local susceptible controls were evaluated for mite infestation levels, mite reproductive success and two potential mechanisms explaining colony survival: grooming of adult worker bees and Varroa Sensitive Hygiene (VSH): adult workers specifically detecting and removing mite-infested brood. RESULTS: Mite infestation levels were significantly lower in surviving colonies and mite reproductive success was reduced by 30% when compared to the controls. No significant differences were found between surviving and control colonies for either grooming or VSH. DISCUSSION: Our data confirm that reduced mite reproductive success seems to be a key factor for natural survival of infested A. mellifera colonies. However, neither grooming nor VSH seem to explain colony survival. Instead, other behaviors of the adult bees seem to be sufficient to hinder mite reproductive success, because brood for this experiment was taken from susceptible donor colonies only. To mitigate the global impact of V. destructor, we suggest learning more from nature, i.e., identifying the obviously efficient mechanisms favored by natural selection.

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.

The effects of primiparity on reproductive performance in the brown bear.

We studied the effects of primiparity on litter size, offspring size, and cub loss in brown bears (Ursus arctos) in two study areas (north, south) in Sweden from 1987 to 2006. Sexually selected infanticide (SSI) has been suggested previously as a mortality factor in our study populations. Females in the south became primiparous earlier than females in the north. Primiparous females had significantly smaller litters of cubs than multiparous females. We found no evidence that primiparity was costly in terms of the interlitter interval. Primiparous mothers had a higher probability of cub loss than multiparous mothers. The probability of cub loss was analyzed separately for the pre-mating and the mating season. The probability of cub loss by primiparous females in the pre-mating season increased with both increasing population density and deteriorating food conditions, whereas the probability of cub loss during the mating season decreased with increasing age of primiparity and increased with male turnover (a variable predicting SSI). The temporal patterns of cub loss by primiparous females suggested that the critical times for reproductive success by primiparous females were the pre-mating season (from birth to shortly after leaving the den) and the mating season. Cub loss in these periods was independent and caused by different factors. Cub loss before the mating season seemed to be most influenced by food conditions, whereas that during the mating season appeared to be caused by SSI.