Next Article in Journal
Metagenomics Analysis Reveals an Extraordinary Inner Bacterial Diversity in Anisakids (Nematoda: Anisakidae) L3 Larvae
Next Article in Special Issue
Whole Genome Sequence Analysis of a Novel Apilactobacillus Species from Giant Honeybee (Apis dorsata) Gut Reveals Occurrence of Genetic Elements Coding Prebiotic and Probiotic Traits
Previous Article in Journal
Nicotine Increases Macrophage Survival through α7nAChR/NF-κB Pathway in Mycobacterium avium paratuberculosis Infection
Previous Article in Special Issue
Seed Meals from Brassica nigra and Eruca sativa Control Artificial Nosema ceranae Infections in Apis mellifera
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Microorganisms
Volume 9
Issue 5
10.3390/microorganisms9051087
microorganisms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Loreley Castelli
1,†,
María Laura Genchi García
2,3,4,†,
Anne Dalmon
5,
Daniela Arredondo
1,
Karina Antúnez
1,
Ciro Invernizzi
6,
Francisco José Reynaldi
2,4,
Yves Le Conte
5 and
Alexis Beaurepaire
5,7,*
1
Laboratorio de Microbiología y Salud de las Abejas, Departamento de Microbiología, Instituto de Investigaciones Biológicas Clemente Estable, Montevideo 11600, Uruguay
2
Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), La Plata 1900, Buenos Aires, Argentina
3
Instituto Multidisciplinario de Biología Celular (IMBICE), La Plata 1900, Buenos Aires, Argentina
4
Laboratorio de Virología, Facultad de Ciencias Veterinarias, Universidad Nacional de La Plata (LAVIR-FCV-UNLP), La Plata 1900, Buenos Aires, Argentina
5
Abeilles et Environnement, INRAE, 84000 Avignon, France
6
Sección Etología, Facultad de Ciencias, Montevideo 11400, Uruguay
7
Institute of Bee Health, University of Bern, 3003 Bern, Switzerland
*
Author to whom correspondence should be addressed.
†
These authors have contributed equally to this work.
Microorganisms 2021, 9(5), 1087; https://doi.org/10.3390/microorganisms9051087
Submission received: 31 March 2021 / Revised: 1 May 2021 / Accepted: 11 May 2021 / Published: 18 May 2021
(This article belongs to the Special Issue Microorganisms in Pollinators: Interactions with Other Factors)
Browse Figures
Versions Notes

Abstract

:
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.

1. Introduction

RNA viruses stand amongst the most diverse and widespread disease-causing agents [1,2]. In these pathogens, genetic variation is enhanced by the accumulation of mutations during replication and their rearrangement by genetic recombination [3,4]. This high mutation rate allows RNA viruses to generate collections of closely related viral genomes around one or more master genotypes [5,6,7]. Virus transmission can occur in two different ways: horizontally, i.e., among individuals of the same generation (directly through contact between individuals, contaminated colony food or food sources, or indirectly through biological vectors), and vertically, i.e., from parents to offspring [8]. Notably, the horizontal transmission of viruses is widespread in social species, where barriers between individuals are significantly reduced due to diverse factors, such as high densities and high rates of interaction [9,10].
RNA viruses play a major role in the pollinator losses currently observed across the globe [11,12]. Many studies have focused on their impact on the health of the western honey bee (Apis mellifera) [13,14]. In this species, the impact of RNA viruses is greatly enhanced by the invasive ectoparasite mite Varroa destructor. In fact, while feeding on honey bees, Varroa can act as a vector for highly virulent viruses [15,16]. V. destructor also influences the immune status of its host, triggering uncontrolled replication through immunosuppression [17,18]. As a consequence, after the introduction of this parasite into novel ecosystems, honey bee viral landscapes significantly change [19,20,21,22]. Among the 70 different viruses that have been characterized in honey bees, a great proportion were reported to be transmitted by V. destructor, including the ubiquitous Deformed wing virus (DWV), Acute bee paralysis virus (ABPV), Kashmir bee virus (KBV) and Israeli acute bee paralysis virus (IAPV) (reviewed in Beaurepaire et al. [23]). Due to the multiple routes of infection, resulting in frequent co-infections [24], Black queen cell virus (BQCV) and Sacbrood virus (SBV) are also often detected in A. mellifera colonies [25,26].
In honey bee epidemiology, the colony is often considered the basic unit, suggesting that viral loads are uniform among individuals from the same hive. However, due to their polyandrous mating system, A. mellifera colonies display genetic heterogeneity across individuals emerging from the same queen, which may be of high relevance when studying pathogens [27]. Thus, within a colony, workers can share the same parents, or their fathers can differ, resulting in half-sister groups, also called ‘subfamilies’ or ‘patrilines’. According to Sherman et al. [28], queens benefit from mating with more than one male to reduce the probability of having an offspring mainly composed of genotypes susceptible to a pathogen. Polyandry may indeed provide beneficial non-additive genetic interactions among subfamilies by capturing rare alleles that regulate resistance to pathogens and parasites in a population. Consequently, genetic variation within colonies may increase the chance of surviving diseases [29,30]. This theory predicts that different subfamilies found in a colony express distinct resistance levels to pathogens with which they co-evolved [31]. However, empirical studies testing this theoretical prediction remain scarce.
The aim of this study was to analyze the patterns of infection of some of the major honey bee RNA viruses at the intra-colonial level. This was achieved by comparing virus infection levels between individual A. mellifera pupae of known patrilines across distinct periods and under different V. destructor infestation statuses. The results documented here reveal new insights into the ecology of honey bee viruses and stress the need to further consider intra-colonial heterogeneity in disease infection for social insect research.

2. Materials and Methods

Colonies of A. mellifera were sampled in the INRAE institute’s apiary (Avignon, France) from the end of August to the end of October 2018 to study the relationships between mite infestation and worker subfamilies (see [32]). In brief, worker brood cells were screened in colonies not treated against V. destructor for over a year before running the experiments. Reproductive characteristics of V. destructor mites were analyzed by opening brood cells and examining their content (pupae age, recapping behavior, foundress number and all other Varroa stages).
The phenotyped individuals were dissected directly upon sampling. DNA from 486 collected pupae (purple eye stage) was isolated using a hind leg per individual according to the Chelex extraction protocol described by Walsh et al. [33]. The patrilines of the pupae were reconstructed using microsatellite markers, and the relationship between mite infestation traits and brood genotypes was tested. The rest of the collected pupae were stored on dry ice during sampling and transferred at −80 °C for further analysis directly after the phenotyping.
In the current study, individuals belonging to the most prevalent subfamilies in the three colonies genotyped from Beaurepaire et al. [32] (mean = 9.82 pupae per subfamilies) were analyzed. These samples, all infected by one reproducing or non-reproducing foundress V. destructor, originated from three colonies with a total of 36 samples per colony. Accordingly, the 36 selected pupae per colony belonged to 3–4 distinct subfamilies (Table 1). These 108 individuals were used to study patterns of infection with four honey bee viruses, including two DWV variants (DWV-A, DWV-B, ABPV, SBV and BQCV) (Table 2).
Each pupa was crushed in 800 µL of QIAzol (Qiagen, Courtaboeuf, France) containing a 9 mm silver bead using a TissueLyser (Qiagen) at 30 Hz for 30 s with 30 s intervals and repeated four times, according to the protocol described in Dalmon et al. [34]. The tube was centrifuged for 1 min at 12,000× g and 4 °C. RNA extraction was performed using 0.5 mL of crude extract of each pupa. An RNeasy Plus Universal mini kit (Qiagen) was used according to the manufacturer’s protocol.
The RNA was quantified using a Nanodrop® spectrophotometer (Thermo Scientific, Illkirch, France) and normalized at 500 ng/µL. Reverse transcription was performed with the high-capacity RNA to cDNA kit (Applied Biosystems, Saint Aubin, France) using 1 µg RNA in 20 µL total reaction volume containing random primers according to the manufacturer’s protocol.
The number of virus genome copies of each virus was determined by quantitative PCR using a StepOne-Plus Real-Time PCR System (Applied Biosystems®, Thermo Scientific, Illkirch, France). Three microliters of ten-fold diluted cDNA were added to 5 µL of SYBR Green master mix (Applied Biosystems). Five reactions were conducted independently, each containing 10 µM of forward and reverse virus and variant-specific primers [35,36] (Table 2). Amplification was performed using the following program: 10 min at 95 °C, 40 cycles of 15 s at 95 °C, 1 min at 60 °C and a melting curve was generated from 60 to 95 °C. Absolute reference to the corresponding synthetic DNA fragments (Eurofins Genomics, Ebersberg, Germany, Table S1) was used to calculate the mean number of virus genome copies per pupa. Detection thresholds were set up at Cq = 31. Quantification was replicated twice per sample and included negative controls for reverse transcription and qPCR. All samples except negative controls tested positive at least for one virus.
Table
Table 2. Set of primers for quantitative PCR. The name, primer sequence, fragment size of the PCR product, annealing temperature (Tm) and primer references are provided. Corresponding synthetic DNA fragments used for quantification are listed in Table S1.
Table 2. Set of primers for quantitative PCR. The name, primer sequence, fragment size of the PCR product, annealing temperature (Tm) and primer references are provided. Corresponding synthetic DNA fragments used for quantification are listed in Table S1.
NameSequence (5′–3′)Size (bp)TmRef.
DWVnew-F1TACTAGTGCTGGTTTTCCTTT15679.5 °C[37]
DWVA-R1CTCATTAACTGTGTCGTTGAT
DWVB-R1CTCATTAACTGAGTTGTTGTC
ABPV-F6548TCATACCTGCCGATCAAG19782.1 °C[36]
KIABPV-B6707CTGAATAATACTGTGCGTATC
SBV-qF3164TTGGAACTACGCATTCTCTG33578.5 °C[36]
SBV-qB3461CTCTAACCTCGCATCAAC
BQCV-qF7893AGTGGCGGAGATGTATGC29481.6 °C[36]
BQCV-qB8150GGAGGTGAAGTGGCTATATC
To assess the intra-colonial patterns of viral infections across the sampled pupae, variations of honey bee viruses across V. destructor infestation levels and reproduction type, and among brood age categories, sampling periods, colonies and subfamilies, were compared. All statistical tests described below were conducted with R v. 4.0.4 [37], and the levels of significance used were corrected for multiple comparisons using the Bonferroni method to adjust p-values when required.
Given that pupae of different ages (i.e., 5 to 9 days old) were collected during the sampling, virus genome copies were compared between the different age categories using a linear regression to test whether virus infection slopes differed significantly from zero. These tests were non-significant, thereby enabling the grouping of samples of different ages for subsequent analyses (Figure S1). As the samples were collected during two distinct sampling periods (I: end of August until mid-September and II: October), levels of virus infection at these two periods were compared using Wilcoxon–Mann–Whitney tests for each colony independently. In cases of significant time differences, samples were grouped by period for subsequent analyses. One colony (colony B) was removed from this temporal analysis as only three samples were collected in this colony during the second sampling period (Table 1).
Differences of virus genome copies among the three colonies and among the patrilines within each colony were tested using non-parametric tests: Wilcoxon–Mann–Whitney tests were used when comparing two groups, Kruskal–Wallis tests were used to compare >2 groups together and Dunn tests were used in the case of significant differences for subsequent pairwise comparisons. Notably, only groups with at least six samples were considered for these statistical analyses.
Finally, the number of virus genome copies of the colonies was compared to the mite infestation level in the brood in the same matching period. Moreover, as all pupae chosen for the experiment were infested with a single foundress mite, successfully reproducing or otherwise, additional Wilcoxon–Mann–Whitney tests were applied to assess the effect of mite reproduction on the virus genome copies.

3. Results

3.1. Overall Virus Quantification

Of the four viruses investigated (DWV, ABPV, SBV and BQCV), only one (DWV) was found to be over the detection thresholds (Appendix A). All samples were infected by the two variants of this virus (DWV-A and DWV-B). Preliminary observations of the distribution of DWV variants in the pupae showed that most samples were either moderately (<1010 genome copies/pupa) or strongly (>1012 genome copies/pupa) infected (Figure 1). Interestingly, when observing the number of virus genome copies with the two variants, four groups could be distinguished: (i) a high number of genome copies of both variants, (ii) a low number of genome copies of both variants, (iii) a high number of genome copies of DWV-A but a low number of genome copies of DWV-B and (iv) a low number of genome copies of DWV-A but a high number of genome copies of DWV-B (Figure 2).

3.2. Virus Genome Copies across Time and Colonies

The comparison of the DWV variants genome copies across the two collection periods revealed contrasting temporal patterns in the two colonies tested. Significant increases in DWV-A were detected in both colonies, while the levels of DWV-B did not change significantly over time (Table 3, Figure 3). Given these temporal patterns, the different periods were considered independently for the subsequent analyses for both variants. Comparing the virus genome copies across the different pupae age stages (day 5 to 8 post-capping) did not reveal any significant differences, although slight decreases in virus genome copies during pupal development could be observed (Figure S1).
When analyzing the DWV-A and DWV-B infection levels across the colonies, significant differences were found during the first sampling period between colony B and the other two colonies (Figure 3).

3.3. Virus Genome Copies across Subfamilies

Comparing the levels of DWV variants between the different subfamilies within colonies and sampling periods revealed some significant differences across the groups tested. More precisely, significant differences were found when comparing two and three patrilines for their infection levels with DWV-A and DWV-B, respectively (Table S2, Figure 4). Interestingly, these results varied for the two DWV variants investigated, i.e., some patrilines were differentially infected by one variant, but not with the other.

3.4. Virus Genome Copies and Mite Infestation

When comparing the level of mite infestations of the colonies with the number of virus genome copies of DWV-A and DWV-B, similar patterns were found across the colonies and time periods. More precisely, virus genome copies of both variants appeared to be positively correlated to levels of mite infestations in the brood for two colonies (B and E), but not for the third colony (D) (Figure S2). Finally, comparing pupae infected by a cell with a reproducing vs. a non-reproducing mite did not reveal any statistically significant differences between the two groups (Wilcoxon–Mann–Whitney tests, all p > 0.05, Table S3).

4. Discussion

In this study, patterns of infections by honey bee viruses were investigated within A. mellifera colonies. Although four common viruses were screened for, only one (DWV) was detected. Investigating two variants of this virus (DWV-A and DWV-B) revealed novel insights into the biology of this major bee pathogen by unravelling new features of its epidemiology across time, colonies, subfamilies and individuals.
To date, the natural infection of honey bee viruses has been principally studied in adults, and very rarely in the brood [24]. In the present study, 108 five to nine-day-old pupae were sampled and analyzed for four major honey bee viruses. However, three out of the four viruses analyzed were below the detection thresholds. Notably, low infection with BQCV, ABPV and SBV in A. mellifera pupae sampled in fall in France was also documented in another study [38]. Given the variety of factors (e.g., time, colonies and patrilines) playing a significant role in DWV infection patterns documented here, it is important to conduct follow-up studies including higher sample sizes per colony to perform more powerful statistical analyses (i.e., models) in order to determine the contribution of each factor individually, as well as interactions between variables. Another important aspect to consider is that, although the DWV-A and DWV-B are known to recombine readily [34,39], recombinants between these two variants could not be analyzed with the methods used here. Nevertheless, the results presented in this study bring several new insights into the dynamics of the two major DWV variants.
First, the analysis of DWV variant infection levels overall, as well as between the analyzed colonies, revealed that the pupae were either moderately (<1010 genome copies/pupae) or strongly (>1012 genome copies/pupae) infected. Furthermore, investigations of the relationships between the two variants of DWV at the individual level revealed four distinct groups of pupae: high infections with both variants, moderate infections with both variants and both cases of high infection with one strain and moderate infection with the other. Such puzzling patterns have been described in previous studies, with DWV-A and DWV-B showing a bimodal distribution between symptomatic and asymptomatic adult honey bees [40] and delimitations between groups of about 106–8 DWV-A copies/bee and 109–10 DWV-B copies/bee [41,42]. Here, the virus genome copies data indicate that the majority of the pupae were infected to levels corresponding to future symptomatic adults [43]. Such high virus infections may result from the fact that all cells analyzed were parasitized by one V. destructor foundress and from the absence of treatments against this DWV-vector for more than one year in the colonies prior to sampling.
Two of the colonies analyzed here (colonies D and E) originated from a Varroa-surviving population [44], while the third one (colony B) was from a local susceptible stock that requires frequent mite treatment to survive. This might explain the higher number of genome copies of DWV variants observed in colony B compared to the other two colonies. Notably, some populations of naturally surviving A. mellifera colonies (i.e., that do not require Varroa treatment) display lower virus infection levels in the brood [45]. However, more colonies from both the surviving and susceptible stocks used here should be tested to confirm that this surviving A. mellifera population displays tolerance or resistance to virus infections.
As DWV variant levels differed across colonies, the temporal dynamics of virus genome copies were analyzed at the colony level. This analysis revealed a variation of temporal infection between the two DWV variants, with significant differences in virus genome copies detected between the two sampling periods when comparing the levels of DWV-A, though not of DWV-B. As similar patterns were detected when comparing the genome copies of the two variants to the colony-level brood infestation with V. destructor, the vector was probably not responsible for the differences observed. This suggests that other factors governing the dynamics of the two DWV variants might explain the temporal differences observed when comparing DWV-A and DWV-B over the sampling periods. For instance, differential mortality associated with the two variants [46,47] might have played a role, since more lethal pathogens are generally less likely to transmit [8].
The results of this study also showed that patterns of DWV-A and DWV-B infections can differ significantly across honey bee pupae, even when these hosts belong to the same colony and share the same environment (i.e., the pupae were all infested by a V. destructor foundress and were collected on the same frame). Moreover, when individuals were grouped by subfamilies, significant differences in DWV-A and/or DWV-B genome copies were found. Notably, the significant variation of virus genome copies across subfamilies is in line with previous work showing differential infection in honey bee subfamilies with two other pathogens: Paenibacillus larvae [48] and Ascosphaera apis [49]. Three non-exclusive hypotheses come to mind when trying to explain the results obtained here. First, DWV, as well as the other viruses that were screened for in this study (ABPV, SBV and BQCV), can be transmitted venereally [50,51,52]. Thus, it is possible that the different levels of the two DWV variants observed were caused by venereal transmission, i.e., from differentially infected drones to their offspring. Additionally, differences in DWV genome copies may have also resulted from interactions with other microorganisms, such as pathogens [53,54] or the host’s gut microbiota [55]. Finally, differential resilience towards virus infection exhibited by the pupae subfamilies could explain these results [56,57]. In this case, resistance mechanisms inherited by the pupae from their sires may have caused the patterns observed. A recent study documenting the ‘suppressed in ovo virus infection’ showed that virus resistance has a strong genetic component in A. mellifera [58], thereby supporting the latter hypothesis. From this statement, we can infer that virus infections may vary in the different pupae subfamilies within a colony due to the differential expression of resistance traits across pupae originating from different sires. However, although significant differences were found when comparing subfamilies, important variations were observed within patrilines as well. This suggests that other factors might play a role in determining DWV infection levels in A. mellifera pupae and stresses the need for more studies to identify the source of this variation.
In this study, significant variations of infection between patrilines were observed for DWV-A and DWV-B. DWV-B is known to have proliferated in Europe more recently than DWV-A [47,59]. Despite the differences of exposure of the host to the variants, the variability across honey bee subfamilies documented here suggests that the interactions between the two DWV variants and their hosts have both led to the evolution of adaptations in A. mellifera and/or the viruses. However, the fact that the differences of DWV-A and DWV-B infection were not consistent across patrilines suggests that distinct biological mechanisms might be at play in the interactions among A. mellifera, V. destructor and these two variants.

5. Conclusions

The results documented here show novel patterns of intra-colonial variation in virus infection that reveal an underestimated complexity in the spatiotemporal dynamics of infection by honey bee viruses. These findings stress that intra-colonial heterogeneity should be more carefully considered when conducting projects about honey bee diseases, as well as other social insects with similar life-history traits (i.e., polyandry or polygyny). Colonies are often considered an epidemiological unit, but this might fail to account for the variability described here, leading to erroneous conclusions. In this study, potential co-evolutionary interactions between hosts and virus genotypes may have caused the differential interactions between the two variants of DWV and the individual and subfamily levels of infections. If parental genetic components are indeed responsible for variations in virus infection levels, new avenues for the breeding and selection of more disease-resilient honey bee colonies can be explored. This could contribute to ensuring more sustainable beekeeping globally, as well as providing new knowledge on host–pathogen interactions and the importance of genetic diversity in disease epidemiology.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/microorganisms9051087/s1, Figure S1: Comparison of virus levels across brood stages; Figure S2: Viral genome copies and mite infestation levels, Table S1: Standard Gene Information; Table S2: Result of statistical tests comparing the levels of DWV across patrilines; Table S3: Results of the Wilcoxon–Mann–Whitney tests to compare cells with reproducing and non-reproducing mites.

Author Contributions

Conceptualization, A.D., Y.L.C. and A.B.; methodology, A.D., K.A., F.J.R., C.I., Y.L.C. and A.B.; formal analysis, L.C., M.L.G.G., A.D. and A.B.; resources, A.D. and Y.L.C.; data curation, L.C., M.L.G.G., D.A. and A.B.; writing—original draft preparation, L.C., M.L.G.G. and A.B.; writing—review and editing, A.D., D.A., K.A., C.I., F.J.R. and Y.L.C.; visualization, D.A. and A.B.; supervision, A.D., K.A., C.I., F.J.R. and Y.L.C.; funding acquisition, A.D., F.J.R., C.I., K.A. and Y.L.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by The Persephone Charitable and Environmental Trust and the ECOS-SUD program.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All additional data can be obtained from the corresponding author upon reasonable request.

Acknowledgments

The authors want to thank Didier Crauser for his help with the honey bee colony management, as well as the ECOS-Sud program and The Persephone Charitable and Environmental Trust for funding this research. The authors would also like to thank the Vinetum foundation and the University of Bern for their kind support.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Full dataset including cq values for each sample analyzed.

References

  1. Mihara, T.; Koyano, H.; Hingamp, P.; Grimsley, N.; Goto, S.; Ogata, H. Taxon Richness of “Megaviridae” Exceeds those of Bacteria and Archaea in the Ocean. Microbes Environ. 2018, 33, 162–171. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Shi, M.; Lin, X.-D.; Tian, J.-H.; Chen, L.-J.; Chen, X.; Li, C.-X.; Qin, X.-C.; Li, J.; Cao, J.-P.; Eden, J.-S.; et al. Redefining the invertebrate RNA virosphere. Nat. Cell Biol. 2016, 540, 539–543. [Google Scholar] [CrossRef]
  3. Froissart, R.; Roze, D.; Uzest, M.; Galibert, L.; Blanc, S.; Michalakis, Y. Recombination Every Day: Abundant Recombination in a Virus during a Single Multi-Cellular Host Infection. PLoS Biol. 2005, 3, e89. [Google Scholar] [CrossRef] [Green Version]
  4. Moya, A.; Holmes, E.C.; González-Candelas, F. The population genetics and evolutionary epidemiology of RNA viruses. Nat. Rev. Genet. 2004, 2, 279–288. [Google Scholar] [CrossRef]
  5. Domingo, E.; Sheldon, J.; Perales, C. Viral Quasispecies Evolution. Microbiol. Mol. Biol. Rev. 2012, 76, 159–216. [Google Scholar] [CrossRef] [Green Version]
  6. Domingo, E.; Perales, C. Viral quasispecies. PLoS Genet. 2019, 15, e1008271. [Google Scholar] [CrossRef] [Green Version]
  7. Dolan, P.T.; Whitfield, Z.J.; Andino, R. Mechanisms and Concepts in RNA Virus Population Dynamics and Evolution. Annu. Rev. Virol. 2018, 5, 69–92. [Google Scholar] [CrossRef]
  8. Fries, I.; Camazine, S. Implications of horizontal and vertical pathogen transmission for honey bee epidemiology. Apidologie 2001, 32, 199–214. [Google Scholar] [CrossRef] [Green Version]
  9. Kappeler, P.M.; Cremer, S.; Nunn, C.L. Sociality and health: Impacts of sociality on disease susceptibility and transmission in animal and human societies. Philos. Trans. R. Soc. B Biol. Sci. 2015, 370, 20140116. [Google Scholar] [CrossRef] [Green Version]
  10. Caillaud, D.; Levréro, F.; Cristescu, R.; Gatti, S.; Dewas, M.; Douadi, M.; Gautier-Hion, A.; Raymond, M.; Ménard, N. Gorilla susceptibility to Ebola virus: The cost of sociality. Curr. Biol. 2006, 16, R489–R491. [Google Scholar] [CrossRef] [Green Version]
  11. McMenamin, A.J.; Flenniken, M.L. Recently identified bee viruses and their impact on bee pollinators. Curr. Opin. Insect Sci. 2018, 26, 120–129. [Google Scholar] [CrossRef]
  12. Gisder, S.; Genersch, E. Viruses of commercialized insect pollinators. J. Invertebr. Pathol. 2017, 147, 51–59. [Google Scholar] [CrossRef]
  13. Chen, Y.P.; Pettis, J.S.; Collins, A.; Feldlaufer, M.F. Prevalence and Transmission of Honeybee Viruses. Appl. Environ. Microbiol. 2006, 72, 606–611. [Google Scholar] [CrossRef] [Green Version]
  14. McMenamin, A.J.; Genersch, E. Honey bee colony losses and associated viruses. Curr. Opin. Insect Sci. 2015, 8, 121–129. [Google Scholar] [CrossRef]
  15. Sumpter, D.J.T.; Martin, S.J. The dynamics of virus epidemics inVarroa-infested honey bee colonies. J. Anim. Ecol. 2004, 73, 51–63. [Google Scholar] [CrossRef]
  16. Nazzi, F.; Brown, S.P.; Annoscia, D.; Del Piccolo, F.; Di Prisco, G.; Varricchio, P.; Della Vedova, G.; Cattonaro, F.; Caprio, E.; Pennacchio, F. Synergistic Parasite-Pathogen Interactions Mediated by Host Immunity Can Drive the Collapse of Honeybee Colonies. PLOS Pathog. 2012, 8, e1002735. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Yang, X.; Cox-Foster, D.L. Impact of an ectoparasite on the immunity and pathology of an invertebrate: Evidence for host immunosuppression and viral amplification. Proc. Natl. Acad. Sci. USA 2005, 102, 7470–7475. [Google Scholar] [CrossRef] [Green Version]
  18. Di Prisco, G.; Annoscia, D.; Margiotta, M.; Ferrara, R.; Varricchio, P.; Zanni, V.; Caprio, E.; Nazzi, F.; Pennacchio, F. A mutualistic symbiosis between a parasitic mite and a pathogenic virus undermines honey bee immunity and health. Proc. Natl. Acad. Sci. USA 2016, 113, 3203–3208. [Google Scholar] [CrossRef] [Green Version]
  19. Wilfert, L.; Long, G.; Leggett, H.C.; Schmid-Hempel, P.; Butlin, R.; Martin, S.J.M.; Boots, M. Deformed wing virus is a recent global epidemic in honeybees driven by Varroa mites. Science 2016, 351, 594–597. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Mondet, F.; De Miranda, J.R.; Kretzschmar, A.; Le Conte, Y.; Mercer, A.R. On the Front Line: Quantitative Virus Dynamics in Honeybee (Apis mellifera L.) Colonies along a New Expansion Front of the Parasite Varroa destructor. PLoS Pathog. 2014, 10, e1004323. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Martin, S.J.; Highfield, A.C.; Brettell, L.; Villalobos, E.M.; Budge, G.E.; Powell, M.; Nikaido, S.; Schroeder, D.C. Global Honey Bee Viral Landscape Altered by a Parasitic Mite. Science 2012, 336, 1304–1306. [Google Scholar] [CrossRef]
  22. Remnant, E.J.; Shi, M.; Buchmann, G.; Blacquière, T.; Holmes, E.C.; Beekman, M.; Ashe, A. A Diverse Range of Novel RNA Viruses in Geographically Distinct Honey Bee Populations. J. Virol. 2017, 91, e0015-17. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Beaurepaire, A.; Piot, N.; Doublet, V.; Antunez, K.; Campbell, E.; Chantawannakul, P.; Chejanovsky, N.; Gajda, A.; Heerman, M.; Panziera, D.; et al. Diversity and Global Distribution of Viruses of the Western Honey Bee, Apis mellifera. Insects 2020, 11, 239. [Google Scholar] [CrossRef] [PubMed]
  24. Yañez, O.; Piot, N.; Dalmon, A.; De Miranda, J.R.; Chantawannakul, P.; Panziera, D.; Amiri, E.; Smagghe, G.; Schroeder, D.; Chejanovsky, N. Bee Viruses: Routes of Infection in Hymenoptera. Front. Microbiol. 2020, 11, 943. [Google Scholar] [CrossRef] [PubMed]
  25. Al Naggar, Y.; Paxton, R.J. Mode of Transmission Determines the Virulence of Black Queen Cell Virus in Adult Honey Bees, Posing a Future Threat to Bees and Apiculture. Viruses 2020, 12, 535. [Google Scholar] [CrossRef]
  26. Dubois, E.; Dardouri, M.; Schurr, F.; Cougoule, N.; Sircoulomb, F.; Thiéry, R. Outcomes of honeybee pupae inoculated with deformed wing virus genotypes A and B. Apidologie 2019, 51, 18–34. [Google Scholar] [CrossRef] [Green Version]
  27. Van Baalen, M.; Beekman, M. The Costs and Benefits of Genetic Heterogeneity in Resistance against Parasites in Social Insects. Am. Nat. 2006, 167, 568–577. [Google Scholar] [CrossRef]
  28. Sherman, P.W.; Seeley, T.D.; Reeve, H.K. Parasites, Pathogens, and Polyandry in Social Hymenoptera. Am. Nat. 1988, 131, 602–610. [Google Scholar] [CrossRef]
  29. Delaplane, K.S.; Pietravalle, S.; Brown, M.A.; Budge, G.E. Honey Bee Colonies Headed by Hyperpolyandrous Queens Have Improved Brood Rearing Efficiency and Lower Infestation Rates of Parasitic Varroa Mites. PLoS ONE 2015, 10, 0142985. [Google Scholar] [CrossRef] [Green Version]
  30. Simone-Finstrom, M.; Walz, M.; Tarpy, D.R. Genetic diversity confers colony-level benefits due to individual immunity. Biol. Lett. 2016, 12, 20151007. [Google Scholar] [CrossRef]
  31. Kraus, F.B.; Page, R.E. Parasites, Pathogens and Polyandry in Social Insects. Am. Nat. 1998, 151, 383–391. [Google Scholar] [CrossRef] [PubMed]
  32. Beaurepaire, A.; Sann, C.; Arredondo, D.; Mondet, F.; Le Conte, Y. Behavioral Genetics of the Interactions between Apis mellifera and Varroa destructor. Insects 2019, 10, 299. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Dorman, N. Citations. Biotechnique 2012, 53, 69. [Google Scholar] [CrossRef] [Green Version]
  34. Dalmon, A.; Desbiez, C.; Coulon, M.; Thomasson, M.; Le Conte, Y.; Alaux, C.; Vallon, J.; Moury, B. Evidence for positive selection and recombination hotspots in Deformed wing virus (DWV). Sci. Rep. 2017, 7, srep41045. [Google Scholar] [CrossRef]
  35. Locke, B.; Forsgren, E.; Fries, I.; De Miranda, J.R. Acaricide Treatment Affects Viral Dynamics in Varroa destructor-Infested Honey Bee Colonies via both Host Physiology and Mite Control. Appl. Environ. Microbiol. 2011, 78, 227–235. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Kevill, J.L.; Highfield, A.; Mordecai, G.J.; Martin, S.J.; Schroeder, D.C. ABC Assay: Method Development and Application to Quantify the Role of Three DWV Master Variants in Overwinter Colony Losses of European Honey Bees. Viruses 2017, 9, 314. [Google Scholar] [CrossRef] [Green Version]
  37. R Core Team R: A Language and Environment for Statistical Computing. 2016. Available online: https://www.r-project.org/ (accessed on 17 May 2021).
  38. Tentcheva, D.; Gauthier, L.; Zappulla, N.; Dainat, B.; Cousserans, F.; Colin, M.E.; Bergoin, M. Prevalence and Seasonal Variations of Six Bee Viruses in Apis mellifera L. and Varroa destructor Mite Populations in France. Appl. Environ. Microbiol. 2004, 70, 7185–7191. [Google Scholar] [CrossRef] [Green Version]
  39. Moore, J.; Jironkin, A.; Chandler, D.; Burroughs, N.; Evans, D.J.; Ryabov, E.V. Recombinants between Deformed wing virus and Varroa destructor virus-1 may prevail in Varroa destructor-infested honeybee colonies. J. Gen. Virol. 2010, 92, 156–161. [Google Scholar] [CrossRef]
  40. Schurr, F.; Tison, A.; Militano, L.; Cheviron, N.; Sircoulomb, F.; Rivière, M.-P.; Ribière-Chabert, M.; Thiéry, R.; Dubois, E. Validation of quantitative real-time RT-PCR assays for the detection of six honeybee viruses. J. Virol. Methods 2019, 270, 70–78. [Google Scholar] [CrossRef]
  41. Ryabov, E.V.; Wood, G.R.; Fannon, J.M.; Moore, J.D.; Bull, J.C.; Chandler, D.; Mead, A.; Burroughs, N.; Evans, D.J. A Virulent Strain of Deformed Wing Virus (DWV) of Honeybees (Apis mellifera) Prevails after Varroa destructor-Mediated, or In Vitro, Transmission. PLoS Pathog. 2014, 10, 1004230. [Google Scholar] [CrossRef] [Green Version]
  42. Highfield, A.C.; El Nagar, A.; Mackinder, L.C.M.; Noël, L.M.-L.J.; Hall, M.J.; Martin, S.J.; Schroeder, D.C. Deformed Wing Virus Implicated in Overwintering Honeybee Colony Losses. Appl. Environ. Microbiol. 2009, 75, 7212–7220. [Google Scholar] [CrossRef] [Green Version]
  43. Dainat, B.; Evans, J.D.; Chen, Y.P.; Gauthier, L.; Neumann, P. Dead or Alive: Deformed Wing Virus and Varroa destructor Reduce the Life Span of Winter Honeybees. Appl. Environ. Microbiol. 2011, 78, 981–987. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Le Conte, Y.; De Vaublanc, G.; Crauser, D.; Jeanne, F.; Rousselle, J.-C.; Bécard, J.-M. Honey bee colonies that have survivedVarroa destructor. Apidologie 2007, 38, 566–572. [Google Scholar] [CrossRef] [Green Version]
  45. Oddie, M.A.Y.; Lanz, S.; Dahle, B.; Yañez, O.; Neumann, P. Virus infections in Varroa destructor-resistant honeybees. bioRxiv 2020. [Google Scholar] [CrossRef]
  46. McMahon, D.P.; Natsopoulou, M.E.; Doublet, V.; Fürst, M.; Weging, S.; Brown, M.J.F.; Gogol-Döring, A.; Paxton, R.J. Elevated virulence of an emerging viral genotype as a driver of honeybee loss. Proc. R. Soc. B Boil. Sci. 2016, 283, 20160811. [Google Scholar] [CrossRef] [PubMed]
  47. Kevill, J.L.; De Souza, F.S.; Sharples, C.; Oliver, R.; Schroeder, D.C.; Martin, S.J. DWV-A Lethal to Honey Bees (Apis mellifera): A Colony Level Survey of DWV Variants (A, B, and C) in England, Wales, and 32 States across the US. Viruses 2019, 11, 426. [Google Scholar] [CrossRef] [Green Version]
  48. Oldroyd, B.P.; Palmer, K.A. Evidence for intra-colonial genetic variance in resistance to American foulbrood of honey bees (Apis mellifera): Further support for the parasite/pathogen hypothesis for the evolution of polyandry. Naturwissenschaften 2003, 90, 265–268. [Google Scholar] [CrossRef] [PubMed]
  49. Invernizzi, C.; Peñagaricano, F.; Tomasco, I.H. Intracolonial genetic variability in honeybee larval resistance to the chalkbrood and American foulbrood parasites. Insectes Sociaux 2009, 56, 233–240. [Google Scholar] [CrossRef]
  50. Yañez, O.; Jaffé, R.; Jarosch, A.; Fries, I.; Moritz, R.F.A.; Paxton, R.J.; De Miranda, J.R. Deformed wing virus and drone mating flights in the honey bee (Apis mellifera): Implications for sexual transmission of a major honey bee virus. Apidologie 2012, 43, 17–30. [Google Scholar] [CrossRef] [Green Version]
  51. Amiri, E.; Meixner, M.D.; Kryger, P. Deformed wing virus can be transmitted during natural mating in honey bees and infect the queens. Sci. Rep. 2016, 6, 33065. [Google Scholar] [CrossRef]
  52. Prodělalová, J.; Moutelíková, R.; Titěra, D. Multiple Virus Infections in Western Honeybee (Apis mellifera L.) Ejaculate Used for Instrumental Insemination. Viruses 2019, 11, 306. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Eiri, D.M.; Suwannapong, G.; Endler, M.; Nieh, J.C. Nosema ceranae Can Infect Honey Bee Larvae and Reduces Subsequent Adult Longevity. PLoS ONE 2015, 10, e0126330. [Google Scholar] [CrossRef] [Green Version]
  54. Urbieta-Magro, A.; Higes, M.; Meana, A.; Barrios, L.; Martín-Hernández, R. Age and Method of Inoculation Influence the Infection of Worker Honey Bees (Apis mellifera) by Nosema ceranae. Insects 2019, 10, 417. [Google Scholar] [CrossRef] [Green Version]
  55. Dosch, C.; Manigk, A.; Streicher, T.; Tehel, A.; Paxton, R.; Tragust, S. The Gut Microbiota Can Provide Viral Tolerance in the Honey Bee. Microorganisms 2021, 9, 871. [Google Scholar] [CrossRef] [PubMed]
  56. Locke, B.; Forsgren, E.; De Miranda, J.R. Increased Tolerance and Resistance to Virus Infections: A Possible Factor in the Survival of Varroa destructor-Resistant Honey Bees (Apis mellifera). PLoS ONE 2014, 9, 99998. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Thaduri, S.; Locke, B.; Granberg, F.; De Miranda, J.R. Temporal changes in the viromes of Swedish Varroa-resistant and Varroa-susceptible honeybee populations. PLoS ONE 2018, 13, e0206938. [Google Scholar] [CrossRef] [PubMed]
  58. De Graaf, D.C.; Laget, D.; De Smet, L.; Boúúaert, D.C.; Brunain, M.; Veerkamp, R.F.; Brascamp, E.W. Heritability estimates of the novel trait ‘suppressed in ovo virus infection’ in honey bees (Apis mellifera). Sci. Rep. 2020, 10, 1–10. [Google Scholar] [CrossRef] [PubMed]
  59. Norton, A.M.; Remnant, E.J.; Buchmann, G.; Beekman, M. Accumulation and Competition Amongst Deformed Wing Virus Genotypes in Naïve Australian Honeybees Provides Insight Into the Increasing Global Prevalence of Genotype, B. Front. Microbiol. 2020, 11, 620. [Google Scholar] [CrossRef] [PubMed]
Microorganisms 09 01087 g001 550
Figure 1. Overall frequency of DWV-A and DWV-B. The distribution of infections with DWV variants (Log10) is shown as the proportion of pupae with between 107 and 1015 DWV genome copies/nymph. (a) Histogram for DWV-A; (b) histogram for DWV-B.
Figure 1. Overall frequency of DWV-A and DWV-B. The distribution of infections with DWV variants (Log10) is shown as the proportion of pupae with between 107 and 1015 DWV genome copies/nymph. (a) Histogram for DWV-A; (b) histogram for DWV-B.
Microorganisms 09 01087 g001
Microorganisms 09 01087 g002 550
Figure 2. Relations between DWV variants. The relationships between the DWV variant (DWV-A and DWV-B) levels in the samples are represented by plotting the infection with both variants for each pupa.
Figure 2. Relations between DWV variants. The relationships between the DWV variant (DWV-A and DWV-B) levels in the samples are represented by plotting the infection with both variants for each pupa.
Microorganisms 09 01087 g002
Microorganisms 09 01087 g003 550
Figure 3. DWV variant genome copies across colonies and sampling periods. The graphs represent the levels of DWV-A and DWV-B across colonies and sampling periods. (a): Sampling period I (end of August to mid-September); (b): Sampling period II (October). The asterisks represent statistical significance after Bonferroni corrections (Period I: Kruskal–Wallis and Dunn tests; Period II: Wilcoxon–Mann–Whitney tests). Only groups with a sample size ≥ 14 were included (sample sizes: Period I: colony B: 33; Colony D: 14; Colony E: 22. Period II: Colony D: 18; Colony E: 18). Asterisks reflect significance levels.
Figure 3. DWV variant genome copies across colonies and sampling periods. The graphs represent the levels of DWV-A and DWV-B across colonies and sampling periods. (a): Sampling period I (end of August to mid-September); (b): Sampling period II (October). The asterisks represent statistical significance after Bonferroni corrections (Period I: Kruskal–Wallis and Dunn tests; Period II: Wilcoxon–Mann–Whitney tests). Only groups with a sample size ≥ 14 were included (sample sizes: Period I: colony B: 33; Colony D: 14; Colony E: 22. Period II: Colony D: 18; Colony E: 18). Asterisks reflect significance levels.
Microorganisms 09 01087 g003
Microorganisms 09 01087 g004 550
Figure 4. DWV variant genome copies across patrilines. Boxplots representing the level of infection across pupae grouped by patrilines within each colony and sampling periods (Col B: period I; Col D: period II; Col E period I and II). The asterisks indicate levels of significance of the non-parametric tests conducted (Wilcoxon–Mann–Whitney or Dunn tests corrected with Bonferroni). (a): DWV-A genome copies; (b): DWV-B genome copies. Error bars represent 95% confidence intervals.
Figure 4. DWV variant genome copies across patrilines. Boxplots representing the level of infection across pupae grouped by patrilines within each colony and sampling periods (Col B: period I; Col D: period II; Col E period I and II). The asterisks indicate levels of significance of the non-parametric tests conducted (Wilcoxon–Mann–Whitney or Dunn tests corrected with Bonferroni). (a): DWV-A genome copies; (b): DWV-B genome copies. Error bars represent 95% confidence intervals.
Microorganisms 09 01087 g004
Table
Table 1. Number of samples and subfamilies per A. mellifera colony screened for RNA viruses. The colony of origin (Colony); subfamily as described in Beaurepaire et al. [32] (Subfamily = patriline); sample size per subfamily (N); per period (NPeriod: I for end of August to mid-September and II for October) and per foundress reproduction groups (NReproduction: R for reproducing foundress and NR for non-reproducing foundress).
Table 1. Number of samples and subfamilies per A. mellifera colony screened for RNA viruses. The colony of origin (Colony); subfamily as described in Beaurepaire et al. [32] (Subfamily = patriline); sample size per subfamily (N); per period (NPeriod: I for end of August to mid-September and II for October) and per foundress reproduction groups (NReproduction: R for reproducing foundress and NR for non-reproducing foundress).
ColonySubfamilyNNPeriodNReproduction
BA12I: 12; II: 0R: 5; NR: 7
B10I: 8; II: 2R: 2; NR: 8
C6I: 5; II: 1R: 4; NR: 2
E8I: 8; II: 0R: 8; NR: 0
DG7I: 1; II: 6R: 5; NR: 2
J12I: 6; II: 6R: 9; NR: 3
K10I: 4; II: 6R: 8; NR: 2
L7I: 3; II: 4R: 5; NR: 2
EQ12I: 6; II: 6R: 10; NR: 2
R12I: 8; II: 4R: 10; NR: 2
S12I: 4; II: 8R: 10; NR: 2
Table
Table 3. Statistical analysis of temporal changes in DWV genome copies. The table indicates the result of Wilcoxon–Mann–Whitney tests conducted to study the differences of DWV, DWV-A and DWV-B between the two sampling periods (I: end of August to mid-September and II: mid- to late October). Asterisks reflect significance levels. NA: not available (sample sizes too low).
Table 3. Statistical analysis of temporal changes in DWV genome copies. The table indicates the result of Wilcoxon–Mann–Whitney tests conducted to study the differences of DWV, DWV-A and DWV-B between the two sampling periods (I: end of August to mid-September and II: mid- to late October). Asterisks reflect significance levels. NA: not available (sample sizes too low).
ColonyDWV-ADWV-B
Col BNANA
Col DW = 231.5, p = 0.012 *W = 147, p = 0.8328
Col EW = 68, p = 0.003 **W = 116, p = 0.1496
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Castelli, L.; Genchi García, M.L.; Dalmon, A.; Arredondo, D.; Antúnez, K.; Invernizzi, C.; Reynaldi, F.J.; Le Conte, Y.; Beaurepaire, A. Intra-Colonial Viral Infections in Western Honey Bees (Apis Mellifera). Microorganisms 2021, 9, 1087. https://doi.org/10.3390/microorganisms9051087

AMA Style

Castelli L, Genchi García ML, Dalmon A, Arredondo D, Antúnez K, Invernizzi C, Reynaldi FJ, Le Conte Y, Beaurepaire A. Intra-Colonial Viral Infections in Western Honey Bees (Apis Mellifera). Microorganisms. 2021; 9(5):1087. https://doi.org/10.3390/microorganisms9051087

Chicago/Turabian Style

Castelli, Loreley, María Laura Genchi García, Anne Dalmon, Daniela Arredondo, Karina Antúnez, Ciro Invernizzi, Francisco José Reynaldi, Yves Le Conte, and Alexis Beaurepaire. 2021. "Intra-Colonial Viral Infections in Western Honey Bees (Apis Mellifera)" Microorganisms 9, no. 5: 1087. https://doi.org/10.3390/microorganisms9051087

APA Style

Castelli, L., Genchi García, M. L., Dalmon, A., Arredondo, D., Antúnez, K., Invernizzi, C., Reynaldi, F. J., Le Conte, Y., & Beaurepaire, A. (2021). Intra-Colonial Viral Infections in Western Honey Bees (Apis Mellifera). Microorganisms, 9(5), 1087. https://doi.org/10.3390/microorganisms9051087

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop