Heritability estimates of the novel trait 'suppressed in ovo virus infection' in honey bees (Apis mellifera)

New insights into honey bee viral and bacterial seasonal infection patterns using third-generation nanopore sequencing on honey bee haemolymph

Honey bees are rapidly declining, which poses a significant threat to our environment and agriculture industry. These vital insects face a disease complex believed to be caused by a combination of parasites, viruses, pesticides, and nutritional deficiencies. However, the real aetiology is still enigmatic. Due to the conventional analysis methods, we still lack complete insights into the honey bee virome and the presence of pathogenic bacteria. To fill this knowledge gap, we employed third-generation nanopore metagenomic sequencing on honey bee haemolymph to monitor the presence of pathogens over almost a year. This study provides valuable insights into the changes in bacterial and viral loads within honey bee colonies. We identified different pathogens in the honey bee haemolymph, which are not included in honey bee screenings. These pathogens comprise the Apis mellifera filamentous virus, Apis rhabdoviruses, and various bacteria such as Frischella sp. and Arsenophonus sp. Furthermore, a sharp contrast was observed between young and old bees. Our research proposes that transgenerational immune priming may play a role in shaping infection patterns in honey bees. We observed a significant increase in pathogen loads in the spring, followed by a notable decrease in pathogen presence during the summer and autumn months. However, certain pathogens seem to be able to evade this priming effect, making them particularly intriguing as potential factors contributing to mortality. In the future, we aim to expand our research on honey bee transgenerational immune priming and investigate its potential in natural settings. This knowledge will ultimately enhance honey bee health and decrease colony mortality.

Eggs sampling as an effective tool for identifying the incidence of viruses in honey bees involved in artificial queen rearing

The Carniolan honey bee (Apis mellifera carnica) plays an essential role in crop pollination, environment diversity, and the production of honey bee products. However, the health of individual honey bees and their colonies is under pressure due to multiple stressors, including viruses as a significant threat to bees. Monitoring various virus infections could be a crucial selection tool during queen rearing. In the present study, samples from all developmental stages (eggs, larvae, pupae, and queens) were screened for the incidence of seven viruses during queen rearing in Slovenia. The screening of a total of 108 samples from five queen breeders was performed by the RT-qPCR assays. The results showed that the highest incidence was observed for black queen cell virus (BQCV), Lake Sinai virus 3 (LSV3), deformed wing virus B (DWV-B), and sacbrood virus (SBV). The highest viral load was detected in queens (6.07 log10 copies/queen) and larvae (5.50 log10 copies/larva) for BQCV, followed by SBV in larvae (5.47 log10 copies/larva). When comparing all the honey bee developmental stages, the eggs exhibited general screening for virus incidence and load in queen mother colonies. The results suggest that analyzing eggs is a good indicator of resilience to virus infection during queen development.

Quantitative trait loci mapping for survival of virus infection and virus levels in honey bees

Israeli acute paralysis virus (IAPV) is a highly virulent, Varroa-vectored virus that is of global concern for honey bee health. Little is known about the genetic basis of honey bees to withstand infection with IAPV or other viruses. We set up and analyzed a backcross between preselected honey bee colonies of low and high IAPV susceptibility to identify quantitative trait loci (QTL) associated with IAPV susceptibility. Experimentally inoculated adult worker bees were surveyed for survival and selectively sampled for QTL analysis based on SNPs identified by whole-genome resequencing and composite interval mapping. Additionally, natural titers of other viruses were quantified in the abdomen of these workers via qPCR and also used for QTL mapping. In addition to the full dataset, we analyzed distinct subpopulations of susceptible and non-susceptible workers separately. These subpopulations are distinguished by a single, suggestive QTL on chromosome 6, but we identified numerous other QTL for different abdominal virus titers, particularly in the subpopulation that was not susceptible to IAPV. The pronounced QTL differences between the susceptible and non-susceptible subpopulations indicate either an interaction between IAPV infection and the bees' interaction with other viruses or heterogeneity among workers of a single cohort that manifests itself as IAPV susceptibility and results in distinct subgroups that differ in their interaction with other viruses. Furthermore, our results indicate that low susceptibility of honey bees to viruses can be caused by both, virus tolerance and virus resistance. QTL were partially overlapping among different viruses, indicating a mixture of shared and specific processes that control viruses. Some functional candidate genes are located in the QTL intervals, but their genomic co-localization with numerous genes of unknown function delegates any definite characterization of the underlying molecular mechanisms to future studies.

Influence of virus abundances in donor colonies and nurse hives on queens of Apis mellifera during the rearing process

Background

Honeybees are one of the three most important animals for mankind. In order to be safe and increase number of bee colonies for pollination, the breeding of queens is necessary. For several decades, bees were selected on economic and behavioral aspects. With the appearance of the neozootic mite Varroa destructor beekeepers were forced to adapt their methods. Varroa destructor can act as a vector for many different bee pathogenic viruses and by this potentiates its devastating impact.

Aim

Methods of rearing queens were not evaluated since the mites' appearance. Besides scientific approaches, viruses received too little attention in regard to the rearing process of honeybee queens. Herein, we present a detailed analysis of virus abundances [Aparavirus, acute bee paralysis virus (ABPV); Triatovirus, black queen cell virus (BQCV); Cripavirus, chronic bee paralysis virus (CBPV); and Iflaviruses, deformed wings virus (DWV), Sacbrood virus (SBV), VDV-1] in breeding hives, donating first instar larvae, hives that are nursing these larvae until the pupa stage, and on queens of Apis mellifera in a breeding apiary.

Methods

Nurse and donor colonies of the queen-rearing process were sampled in the year 2020 and analyzed by RT qPCR. Virus quantifications were correlated with queen mortalities and seasonal effects.

Results

Virus detections increased in reared queens, however, the elevated virus titers did not increase the mortality of the queens until their exclosure. Moreover, we observed a lower interrelation between virus abundance in queens and their original donor colonies, than between nurse hives and their nursed queens.

Conclusion

The bee pathogenic viruses ABPV, BQCV, CBPV, DWV, SBV, and VDV-1 do not influence the mortality of bee queens during the rearing process. Whether respective virus loads result in sublethal or long-term effects remains to be elucidated.

Viral Co-Infections and Antiviral Immunity in Honey Bees

Over the past few decades, honey bees have been facing an increasing number of stressors. Beyond individual stress factors, the synergies between them have been identified as a key factor in the observed increase in colony mortality. However, these interactions are numerous and complex and call for further research. Here, in line with our need for a systemic understanding of the threats that they pose to bee health, we review the interactions between honey bee viruses. As viruses are obligate parasites, the interactions between them not only depend on the viruses themselves but also on the immune responses of honey bees. Thus, we first summarise our current knowledge of the antiviral immunity of honey bees. We then review the interactions between specific pathogenic viruses and their interactions with their host. Finally, we draw hypotheses from the current literature and suggest directions for future research.

Expression of Molecular Markers of Resilience against Varroa destructor and Bee Viruses in Ethiopian Honey Bees (Apis mellifera simensis) Focussing on Olfactory Sensing and the RNA Interference Machinery

Varroa destructor mites and the viruses it vectors are two major factors leading to high losses of honey bees (Apis mellifera) colonies worldwide. However, honey bees in some African countries show resilience to varroa infestation and/or virus infections, although little is known about the mechanisms underlying this resilience. In this study, we investigated the expression profiles of some key molecular markers involved in olfactory sensing and RNA interference, as these processes may contribute to the bees' resilience to varroa infestation and virus infection, respectively. We found significantly higher gene expression of the odorant binding protein, OBP14, in the antennae of Ethiopian bees compared to Belgian bees. This result suggests the potential of OBP14 as a molecular marker of resilience to mite infestation. Scanning electron microscopy showed no significant differences in the antennal sensilla occurrence and distribution, suggesting that resilience arises from molecular processes rather than morphological adaptations. In addition, seven RNAi genes were upregulated in the Ethiopian honey bees and three of them-Dicer-Drosha, Argonaute 2, and TRBP2-were positively correlated with the viral load. We can conclude that the antiviral immune response was triggered when bees were experiencing severe viral infection and that this might contribute to the bees' resilience to viruses.

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.

Consequences of incorrect genetic parameter estimates for single-trait and multi-trait genetic evaluations in honeybees

Genetic and residual variances of traits are important input parameters for best linear unbiased prediction (BLUP) breeding value estimation. In honeybees, estimates of these variances are often associated with large standard errors, entailing a risk to perform genetic evaluations under wrong premises. The consequences hereof have not been sufficiently studied. In particular, there are no adequate investigations on this topic accounting for multi-trait selection or genetic peculiarities of the honeybee. We performed simulation studies and explored the consequences of selection for honeybee populations with a broad range of true and assumed genetic parameters. We found that in single-trait evaluations, the response to selection was barely compromised by assuming erroneous parameters, so that reductions in genetic progress after 20 years never exceeded 21%. Phenotypic selection appeared inferior to BLUP selection, particularly under low heritabilities. Parameter choices for genetic evaluation had great effects on inbreeding development. By wrongly assuming high heritabilities, inbreeding rates were reduced by up to 74%. When parallel selection was performed for two traits, the right choice of genetic parameters appeared considerably more crucial as several incorrect premises yielded inadvertent negative selection for one of the traits. This phenomenon occurred in multiple constellations in which the selection traits expressed a negative genetic correlation. It was not reflected in the estimated breeding values. Our results indicate that breeding efforts heavily rely on detailed knowledge on genetic parameters, particularly when multi-trait selection is performed. Thus, considerable effort should be invested into precise parameter estimations.

Influence of model selection and data structure on the estimation of genetic parameters in honeybee populations

Estimating genetic parameters of quantitative traits is a prerequisite for animal breeding. In honeybees, the genetic variance separates into queen and worker effects. However, under data paucity, parameter estimations that account for this peculiarity often yield implausible results. Consequently, simplified models that attribute all genetic contributions to either the queen (queen model) or the workers (worker model) are often used to estimate variance components in honeybees. However, the causes for estimations with the complete model (colony model) to fail and the consequences of simplified models for variance estimates are little understood. We newly developed the necessary theory to compare parameter estimates that were achieved by the colony model with those of the queen and worker models. Furthermore, we performed computer simulations to quantify the influence of model choice, estimation algorithm, true genetic parameters, rates of controlled mating, apiary sizes, and phenotype data completeness on the success of genetic parameter estimations. We found that successful estimations with the colony model were only possible if at least some of the queens mated controlled on mating stations. In that case, estimates were largely unbiased if more than 20% of the colonies had phenotype records. The simplified queen and worker models proved more stable and yielded plausible parameter estimates for almost all settings. Results obtained from these models were unbiased when mating was uncontrolled, but with controlled mating, the simplified models consistently overestimated heritabilities. This study elucidates the requirements for variance component estimation in honeybees and provides the theoretical groundwork for simplified honeybee models.

Understanding social resilience in honeybee colonies

Honeybee colonies experience high losses, induced by several stressors that can result in the collapse of colonies. Experiments show what effects stressors, such as parasites, pathogens and pesticides, can have on individual honeybees as well as colonies. Although individuals may die, colonies do not always collapse from such disturbances. As a superorganism, the colony can maintain or return back to homeostasis through colony mechanisms. This capacity is defined as social resilience. When the colony faces a high stress load, this may lead to breakdown in mechanisms, loss in resilience and eventually colony collapse. Before social resilience can be measured in honeybees, we need to examine the mechanisms in colonies that allow recovery and maintenance after stressor exposure. Here, we discuss some of these mechanisms and how they affect the social resilience of honeybee colonies. Understanding social resilience in honeybees is essential to managing colony health and loss prevention.

Beeporter: Tools for high-throughput analyses of pollinator-virus infections

Pollinators are in decline thanks to the combined stresses of disease, pesticides, habitat loss, and climate. Honey bees face numerous pests and pathogens but arguably none are as devastating as Deformed wing virus (DWV). Understanding host-pathogen interactions and virulence of DWV in honey bees is slowed by the lack of cost-effective high-throughput screening methods for viral infection. Currently, analysis of virus infection in bees and their colonies is tedious, requiring a well-equipped molecular biology laboratory and the use of hazardous chemicals. Here we describe virus clones tagged with green fluorescent protein (GFP) or nanoluciferase (nLuc) that provide high-throughput detection and quantification of virus infections. GFP fluorescence is measured noninvasively in living bees via commonly available long-wave UV light sources and a smartphone camera, or a standard ultraviolet transilluminator gel imaging system. Nonlethal monitoring with GFP allows continuous screening of virus growth and serves as a direct breeding tool for identifying honey bee parents with increased antiviral resistance. Expression using the nLuc reporter strongly correlates with virus infection levels and is especially sensitive. Using multiple reporters, it is also possible to visualize competition, differential virulence, and host tissue targeting by co-occuring pathogens. Finally, it is possible to directly assess the risk of cross-species "spillover" from honey bees to other pollinators and vice versa.

Developmental environment shapes honeybee worker response to virus infection

The consequences of early-life experiences are far reaching. In particular, the social and nutritional environments that developing animals experience can shape their adult phenotypes. In honeybees, larval nutrition determines the eventual social roles of adults as reproductive queens or sterile workers. However, little is known about the effects of developmental nutrition on important adult worker phenotypes such as disease resilience. In this study, we manipulated worker developmental nutrition in two distinct ways under semi-natural field conditions. In the first experiment, we restricted access to nutrition via social isolation by temporarily preventing alloparental care. In the second experiment, we altered the diet quality experienced by the entire colony, leading to adult bees that had developed entirely in a nutritionally restricted environment. When bees from these two experiments reached the adult stage, we challenged them with a common bee virus, Israeli acute paralysis virus (IAPV) and compared mortality, body condition, and the expression of immune genes across diet and viral inoculation treatments. Our findings show that both forms of early life nutritional stress, whether induced by lack of alloparental care or diet quality restriction, significantly reduced bees' resilience to virus infection and affected the expression of several key genes related to immune function. These results extend our understanding of how early life nutritional environment can affect phenotypes relevant to health and highlight the importance of considering how nutritional stress can be profound even when filtered through a social group. These results also provide important insights into how nutritional stress can affect honeybee health on a longer time scale and its potential to interact with other forms of stress (i.e. disease).

Breeding for Virus Resistance and Its Effects on Deformed Wing Virus Infection Patterns in Honey Bee Queens

Viruses, and in particular the deformed wing virus (DWV), are considered as one of the main antagonists of honey bee health. The 'suppressed in ovo virus infection' trait (SOV) described for the first time that control of a virus infection can be achieved from genetically inherited traits and that the virus state of the eggs is indicative for this. This research aims to explore the effect of the SOV trait on DWV infections in queens descending from both SOV-positive (QDS+) and SOV-negative (QDS-) queens. Twenty QDS+ and QDS- were reared from each time four queens in the same starter-finisher colony. From each queen the head, thorax, ovaries, spermatheca, guts and eviscerated abdomen were dissected and screened for the presence of the DWV-A and DWV-B genotype using qRT-PCR. Queens descending from SOV-positive queens showed significant lower infection loads for DWV-A and DWV-B as well as a lower number of infected tissues for DWV-A. Surprisingly, differences were less expressed in the reproductive tissues, the ovaries and spermatheca. These results confirm that selection on the SOV trait is associated with increased virus resistance across viral genotypes and that this selection drives DWV towards an increased tissue specificity for the reproductive tissues. Further research is needed to explore the mechanisms underlying the interaction between the antiviral response and DWV.

Intra-Colonial Viral Infections in Western Honey Bees (Apis Mellifera)

RNA viruses play a significant role in the current high losses of pollinators. Although many studies have focused on the epidemiology of western honey bee (Apis mellifera) viruses at the colony level, the dynamics of virus infection within colonies remains poorly explored. In this study, the two main variants of the ubiquitous honey bee virus DWV as well as three major honey bee viruses (SBV, ABPV and BQCV) were analyzed from Varroa-destructor-parasitized pupae. More precisely, RT-qPCR was used to quantify and compare virus genome copies across honey bee pupae at the individual and subfamily levels (i.e., patrilines, sharing the same mother queen but with different drones as fathers). Additionally, virus genome copies were compared in cells parasitized by reproducing and non-reproducing mite foundresses to assess the role of this vector. Only DWV was detected in the samples, and the two variants of this virus significantly differed when comparing the sampling period, colonies and patrilines. Moreover, DWV-A and DWV-B exhibited different infection patterns, reflecting contrasting dynamics. Altogether, these results provide new insight into honey bee diseases and stress the need for more studies about the mechanisms of intra-colonial disease variation in social insects.

Comparing Survival of Israeli Acute Paralysis Virus Infection among Stocks of U.S. Honey Bees

Among numerous viruses that infect honey bees (Apis mellifera), Israeli acute paralysis virus (IAPV) can be linked to severe honey bee health problems. Breeding for virus resistance may improve honey bee health. To evaluate the potential for this approach, we compared the survival of IAPV infection among stocks from the U.S. We complemented the survival analysis with a survey of existing viruses in these stocks and assessing constitutive and induced expression of immune genes. Worker offspring from selected queens in a common apiary were inoculated with IAPV by topical applications after emergence to assess subsequent survival. Differences among stocks were small compared to variation within stocks, indicating the potential for improving honey bee survival of virus infections in all stocks. A positive relation between worker survival and virus load among stocks further suggested that honey bees may be able to adapt to better cope with viruses, while our molecular studies indicate that toll-6 may be related to survival differences among virus-infected worker bees. Together, these findings highlight the importance of viruses in queen breeding operations and provide a promising starting point for the quest to improve honey bee health by selectively breeding stock to be better able to survive virus infections.