Ontogeny of immunity and natural viral infection in Apis mellifera drones and workers.

The most common viral diseases affecting honey bees (Apis mellifera) in Israel include deformed wing viruses (DWV-A and DWV-B) and acute paralysis viruses (ABPV and IAPV). These viruses are transmitted within and between colonies, both horizontally and vertically. All members of the colony contribute to this transmission, on the other hand individual and social immunity, particularly hygienic behaviour, may affect the outcome of the process. In this study, we evaluated the ontogeny of natural infections of DWV-A, DWV-B, ABPV and IAPV, their prevalence and loads, in workers and drones from high (H) and low (L) hygienic colonies. In parallel, we evaluated the expression of two immune genes: peptidoglycan recognition protein S2(PGRP-S2) and hymenoptaecin. The prevalence of DWV-B and IAPV increased with age and was higher in workers than in drones. ABPV was not detected in drones. The expression of both immune genes was significantly affected by age and sex. Drones from H colonies had higher expression of these genes. The increased expression of immune genes with drones' age, particularly in hygienic colonies, suggest additional value of honey bee breeding for hygienic behaviour for sustainable beekeeping.

Epidemiology of the Microsporidium Nosema ceranae in Four Mediterranean Countries.

Nosema ceranae is a highly prevalent intracellular parasite of honey bees' midgut worldwide. This Microsporidium was monitored during a long-term study to evaluate the infection at apiary and intra-colony levels in six apiaries in four Mediterranean countries (France, Israel, Portugal, and Spain). Parameters on colony strength, honey production, beekeeping management, and climate were also recorded. Except for São Miguel (Azores, Portugal), all apiaries were positive for N. ceranae, with the lowest prevalence in mainland France and the highest intra-colony infection in Israel. A negative correlation between intra-colony infection and colony strength was observed in Spain and mainland Portugal. In these two apiaries, the queen replacement also influenced the infection levels. The highest colony losses occurred in mainland France and Spain, although they did not correlate with the Nosema infection levels, as parasitism was low in France and high in Spain. These results suggest that both the effects and the level of N. ceranae infection depends on location and beekeeping conditions. Further studies on host-parasite coevolution, and perhaps the interactions with other pathogens and the role of honey bee genetics, could assist in understanding the difference between nosemosis disease and infection, to develop appropriate strategies for its control.

Multiple benefits of breeding honey bees for hygienic behavior.

Honey bee colonies are prone to invasion by pests and pathogens. The combination of the parasitic mite Varroa destructor (Varroa) and the multiple viruses it vectors, is a major driver of colony losses. Breeding for hygienic behavior to reduce Varroa populations is considered a sustainable way to reduce the impact of Varroa on honey bee health. However, hygienic behavior may have a cost to the health of individual bees, both in terms of viral infection risk and immune function. To determine whether selection for hygienic behavior at the colony level is associated with trade-offs in honey bee viral infection and immune function, we compared Varroa populations, viral loads, and individual immune function between honey bee colonies that were bred for high and low hygienic behavior. Specifically, we measured Varroa infestation, Deformed wing virus DWV-A, DWV-B, Acute bee paralysis virus (ABPV), and Israeli acute paralysis virus IAPV viral genome levels in bee samples from artificially inseminated queens in our bi-directional selection program for hygienic behavior in Israel. In addition, we evaluated the expression of 12 genes from the Jak-STAT, Toll, IMD and RNAi immune pathways. We found significantly lower Varroa infestation and DWV loads in highly hygienic colonies than in colonies exhibiting low hygienic behavior. In addition, workers of the hygienic colonies had significantly higher expression of the immune genes PGRP-S2 and hymenoptaecin compared to workers from low hygienic colonies. These results indicate no trade-offs in breeding for hygienic behavior. Hygienic honey bees were associated with reduced Varroa populations and reduced DWV prevalence or load at the colony level. Individual immunity of hygienic bees was increased, which could also contribute to lower virus levels, although lower Varroa levels due to social immunity presumably contributed as well. In sum, we demonstrate multiple health benefits of breeding for honey bee hygiene.

Metatranscriptome Analysis of Sympatric Bee Species Identifies Bee Virus Variants and a New Virus, Andrena-Associated Bee Virus-1.

Bees are important plant pollinators in agricultural and natural ecosystems. High average annual losses of honey bee (Apis mellifera) colonies in some parts of the world, and regional population declines of some mining bee species (Andrena spp.), are attributed to multiple factors including habitat loss, lack of quality forage, insecticide exposure, and pathogens, including viruses. While research has primarily focused on viruses in honey bees, many of these viruses have a broad host range. It is therefore important to apply a community level approach in studying the epidemiology of bee viruses. We utilized high-throughput sequencing to evaluate viral diversity and viral sharing in sympatric, co-foraging bees in the context of habitat type. Variants of four common viruses (i.e., black queen cell virus, deformed wing virus, Lake Sinai virus 2, and Lake Sinai virus NE) were identified in honey bee and mining bee samples, and the high degree of nucleotide identity in the virus consensus sequences obtained from both taxa indicates virus sharing. We discovered a unique bipartite + ssRNA Tombo-like virus, Andrena-associated bee virus-1 (AnBV-1). AnBV-1 infects mining bees, honey bees, and primary honey bee pupal cells maintained in culture. AnBV-1 prevalence and abundance was greater in mining bees than in honey bees. Statistical modeling that examined the roles of ecological factors, including floral diversity and abundance, indicated that AnBV-1 infection prevalence in honey bees was greater in habitats with low floral diversity and abundance, and that interspecific virus transmission is strongly modulated by the floral community in the habitat. These results suggest that land management strategies that aim to enhance floral diversity and abundance may reduce AnBV-1 spread between co-foraging bees.

Infection of a Lepidopteran Cell Line with Deformed Wing Virus.

Many attempts to develop a reliable cell cultured-based system to study honey bee virus infections have encountered substantial difficulties. We investigated the ability of a cell line from a heterologous insect to sustain infection by a honey bee virus. For this purpose, we infected the Lepidopteran hemocytic cell line (P1) with Deformed wing virus (DWV). The genomic copies of DWV increased upon infection, as monitored by quantitative RT-PCR. Moreover, a tagged-primer-based RT-PCR analysis showed the presence of DWV negative-sense RNA in the cells, indicating virus replication. However, the DWV from infected cells was mildly infectious to P1 cells. Similar results were obtained when the virus was injected into Apis mellifera pupae. Thus, though the virus yields from the infected cells appeared to be very low, we show for the first time that DWV can replicate in a heterologous cell line. Given the availability of many other insect cell lines, our study paves the way for future exploration in this direction. In the absence of adequate A. mellifera cell lines, exploring the ability of alternative cell lines to enable honey bee virus infections could provide the means to study and understand the viral infectious cycle at the cellular level and facilitate obtaining purified isolates of these viruses.

Natural Variation in Flower Color and Scent in Populations of Eruca sativa (Brassicaceae) Affects Pollination Behavior of Honey Bees.

Plants of Eruca sativa Mill. (Brassicaceae) from desert and Mediterranean populations in Israel differ in flower color and size. In the desert habitat, the population has higher abundance of flowers with cream color and longer petals, whereas in the Mediterranean habitat, the population has higher abundance of flowers with yellow and shorter petals. Choice experiments with honey bee foragers (Apis mellifera Linn., Apidae, Hymenoptera), the main pollinator in the natural habitat in Israel, confirmed that they are more attracted to the yellow flower morph than to the cream one. A proboscis extension response test indicated that honey bees are able to discriminate between flower scents of the desert and Mediterranean populations. Considering the advantage of plants of the yellow morph in attracting pollinators, we further tested in a common garden experiment whether these possess higher fitness than plants of the desert population. Indeed, a significant association was found between flower color and fruit set, and seed mass. In general, our results provide evidence for ecotypic differentiation between populations imposed by pollinators. The advantage of the yellow color morph in attracting pollinators may explain its dominance among plants of the Mediterranean population. We discuss why the cream color morph may be dominant in the desert habitat, considering the possibility of different pollinators, tradeoffs between traits, or pleiotropy.

New Viruses from the Ectoparasite Mite Varroa destructor Infesting Apis mellifera and Apis cerana.

Varroa destructor is an ectoparasitic mite of Asian or Eastern honeybees Apis cerana(A. cerana) which has become a serious threat to European subspecies of Western honeybees Apis mellifera (A. mellifera) within the last century. V.destructor and its vectored honeybee viruses became serious threats for colony survival. This is a short period for pathogen- and host-populations to adapt. To look for possible variation in the composition of viral populations we performed RNA metagenomic analysis of the Western honeybee subspecies A. m. ligustica, A. m.syriaca, A. m. intermissa, and A. cerana and their respective V. destructor mites. The analysis revealed two novel viruses: Varroa orthomyxovirus-1 (VOV-1) in A. mellifera and V. destructor and a Hubei like-virga virus-14 homolog in V. destructor. VOV-1 was more prevalent in V. destructor than in A. mellifera and we found evidence for viral replication in both hosts. Interestingly, we found differences in viral loads of A. cerana and their V. destructor, A. m. intermissa, and its V. destructor showed partial similarity, while A. m.ligustica and A. m.syriaca and their varroa where very similar. Deformed wing virus exhibited 82.20%, 99.20%, 97.90%, and 0.76% of total viral reads in A. m. ligustica, A. m. syriaca, A. m. intermissa, and A. cerana, respectively. This is the first report of a complete segmented-single-stranded negative-sense RNA virus genome in honeybees and V. destructor mites.

Presence of Apis Rhabdovirus-1 in Populations of Pollinators and Their Parasites from Two Continents.

The viral ecology of bee communities is complex, where viruses are readily shared among co-foraging bee species. Additionally, in honey bees (Apis mellifera), many viruses are transmitted - and their impacts exacerbated - by the parasitic Varroa destructor mite. Thus far, the viruses found to be shared across bee species and transmitted by V. destructor mites are positive-sense single-stranded RNA viruses. Recently, a negative-sense RNA enveloped virus, Apis rhabdovirus-1 (ARV-1), was found in A. mellifera honey bees in Africa, Europe, and islands in the Pacific. Here, we describe the identification - using a metagenomics approach - of ARV-1 in two bee species (A. mellifera and Bombus impatiens) and in V. destructor mites from populations collected in the United States and Israel. We confirmed the presence of ARV-1 in pools of A. mellifera, B. impatiens, and V. destructor from Israeli and U.S. populations by RT-PCR and found that it can reach high titers in individual honey bees and mites (107-108 viral genomic copies per individual). To estimate the prevalence of ARV-1 in honey bee populations, we screened 104 honey bee colonies across Israel, with 21 testing ARV-1-positive. Tagged-primer-mediated RT-PCR analysis detected the presence of the positive-sense ARV-1 RNA in A. mellifera and V. destructor, indicating that ARV-1 replicates in both hosts. This is the first report of the presence of ARV-1 in B. impatiens and of the replication of a rhabdovirus in A. mellifera and V. destructor. Our data suggest that Varroa mites could act as an ARV-1 vector; however, the presence of ARV-1 in B. impatiens (which are not parasitized by Varroa) suggests that it may not require the mite for transmission and ARV-1 may be shared among co-foraging bee species. Given that ARV-1 is found in non-Apis bee species, and because "ARV" is used for the Adelaide River virus, we propose that this virus should be called bee rhabdovirus 1 and abbreviated BRV-1. These results greatly expand our understanding of the diversity of viruses that can infect bee communities, though further analysis is required to determine how infection with this virus impacts these different hosts.