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

Loreley Castelli, María Laura Genchi García, Anne Dalmon, Daniela Arredondo, Karina Antúnez, Ciro Invernizzi, Francisco José Reynaldi, Yves Le Conte, Alexis Beaurepaire

<figure><img src="image_1.png" /></figure>

► To cite this version:

Loreley Castelli, María Laura Genchi García, Anne Dalmon, Daniela Arredondo, Karina Antúnez, et al. Intra-Colonial Viral Infections in Western Honey Bees (Apis Mellifera). Microorganisms, 2021, 9 (5), pp.1087. 10.3390/microorganisms9051087. [https://hal-03252369](https://hal.inrae.fr/hal-03252369)

HAL Id: hal-03252369

https://hal.inrae.fr/hal-03252369

Submitted on 7 Jun 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

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

Loreley Castelli <fnref n="1" />,†, María Laura Genchi García <fnref n="2" />,<fnref n="3" />,<fnref n="4" />,†, Anne Dalmon <fnref n="5" />, Daniela Arredondo <fnref n="1" />, Karina Antúnez <fnref n="1" />, Ciro Invernizzi <fnref n="6" />, Francisco José Reynaldi <fnref n="2" />,<fnref n="4" />, Yves Le Conte <fnref n="5" /> and Alexis Beaurepaire <fnref n="5" />,<fnref n="7" />,*

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

Academic Editor: Dániel Cadar

Received: 31 March 2021

Accepted: 11 May 2021

Published: 18 May 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2021 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

**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., patriline, 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 patriline. 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.

**Keywords:** evolutionary biology; host-pathogen interactions; population genetics; viruses; pollinators

## 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-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

<footnote n="1">¹ Laboratorio de Microbiología y Salud de las Abejas, Departamento de Microbiología, Instituto de Investigaciones Biológicas Clemente Estable, Montevideo 11600, Uruguay; castelli.loreley@gmail.com (L.C.); danielarpapiol@gmail.com (D.A.); kantunez03@gmail.com (K.A.)</footnote>

<footnote n="2">² Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), La Plata 1900, Buenos Aires, Argentina; ml.genchigarcia@gmail.com (M.L.G.G.); freynaldi@yahoo.com (F.J.R.)</footnote>

<footnote n="3">³ Instituto Multidisciplinario de Biología Celular (IMBICE), La Plata 1900, Buenos Aires, Argentina</footnote>

<footnote n="4">⁴ Laboratorio de Virología, Facultad de Ciencias Veterinarias, Universidad Nacional de La Plata (LAVIR-FCV-UNLP), La Plata 1900, Buenos Aires, Argentina</footnote>

<footnote n="5">5 Abeilles et Environnement, INRAE, 84000 Avignon, France; anne.dalmon@inrae.fr (A.D.); yves.le-conte@inrae.fr (Y.L.C.)</footnote>

<footnote n="6">⁶ Sección Etología, Facultad de Ciencias, Montevideo 11400, Uruguay; cirobee@gmail.com</footnote>

<footnote n="7">⁷ Institute of Bee Health, University of Bern, 3003 Bern, Switzerland</footnote>

<footnote>* Correspondence: Alexis.Beaurepaire@vetsuisse.unibe.ch; Tel.: +41-31-631-57-66</footnote>

<footnote>† These authors have contributed equally to this work.</footnote>
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–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).

<table><caption>**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 (N<sub>Period</sub>: I for end of August to mid-September and II for October) and per foundress reproduction groups (N<sub>Reproduction</sub>: R for reproducing foundress and NR for non-reproducing foundress).</caption><thead><tr><th>Colony</th><th>Subfamily</th><th>N</th><th>N<sub>Period</sub></th><th>N<sub>Reproduction</sub></th></tr></thead><tbody><tr><td rowspan="5">B</td><td>A</td><td>12</td><td>I: 12; II: 0</td><td>R: 5; NR: 7</td></tr><tr><td>B</td><td>10</td><td>I: 8; II: 2</td><td>R: 2; NR: 8</td></tr><tr><td>C</td><td>6</td><td>I: 5; II: 1</td><td>R: 4; NR: 2</td></tr><tr><td>E</td><td>8</td><td>I: 8; II: 0</td><td>R: 8; NR: 0</td></tr><tr><td>G</td><td>7</td><td>I: 1; II: 6</td><td>R: 5; NR: 2</td></tr><tr><td rowspan="4">D</td><td>J</td><td>12</td><td>I: 6; II: 6</td><td>R: 9; NR: 3</td></tr><tr><td>K</td><td>10</td><td>I: 4; II: 6</td><td>R: 8; NR: 2</td></tr><tr><td>L</td><td>7</td><td>I: 3; II: 4</td><td>R: 5; NR: 2</td></tr><tr><td>Q</td><td>12</td><td>I: 6; II: 6</td><td>R: 10; NR: 2</td></tr><tr><td rowspan="2">E</td><td>R</td><td>12</td><td>I: 8; II: 4</td><td>R: 10; NR: 2</td></tr><tr><td>S</td><td>12</td><td>I: 4; II: 8</td><td>R: 10; NR: 2</td></tr></tbody></table>

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><caption>**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.</caption><thead><tr><th>Name</th><th>Sequence (5'-3')</th><th>Size (bp)</th><th>Tm</th><th>Ref.</th></tr></thead><tbody><tr><td>DWVnew-F1</td><td>TACTAGTGCTGGTTTTCCTTT</td><td></td><td></td><td></td></tr><tr><td>DWVA-R1</td><td>CTCATTAACTGTGTCGTTGAT</td><td>156</td><td>79.5 °C</td><td>[37]</td></tr><tr><td>DWVB-R1</td><td>CTCATTAACTGAGTTGTTGTC</td><td></td><td></td><td></td></tr><tr><td>ABPV-F6548</td><td>TCATACCTGCCGATCAAG</td><td></td><td></td><td></td></tr><tr><td>KIABPV-B6707</td><td>CTGAATAATACTGTGCGTATC</td><td>197</td><td>82.1 °C</td><td>[36]</td></tr><tr><td>SBV-qF3164</td><td>TTGGAACTACGCATTCTCTG</td><td></td><td></td><td></td></tr><tr><td>SBV-qB3461</td><td>CTCTAACCTCGCATCAAC</td><td>335</td><td>78.5 °C</td><td>[36]</td></tr><tr><td>BQCV-qF7893</td><td>AGTGGCGGAGATGTATGC</td><td></td><td></td><td></td></tr><tr><td>BQCV-qB8150</td><td>GGAGGTGAAGTGGCTATATC</td><td>294</td><td>81.6 °C</td><td>[36]</td></tr></tbody></table>

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 patrilineages 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 (<10¹⁰ genome copies/pupa) or strongly (>10¹² 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).

<figure><img src="image_4.png" /><figcaption>**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 10⁷ and 10¹⁵ DWV genome copies/nymph. (a) Histogram for DWV-A; (b) histogram for DWV-B.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**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.</figcaption></figure>

## 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).
<figure><img src="image_6.png" /><figcaption>**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.</figcaption></figure>

<table><caption>**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).</caption><thead><tr><th>Colony</th><th>DWV-A</th><th>DWV-B</th></tr></thead><tbody><tr><td>Col B</td><td>NA</td><td>NA</td></tr><tr><td>Col D</td><td>W = 231.5, p = 0.012 *</td><td>W = 147, p = 0.8328</td></tr><tr><td>Col E</td><td>W = 68, p = 0.003 **</td><td>W = 116, p = 0.1496</td></tr></tbody></table>

### 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.
<figure><img src="image_7.png" /><figcaption>**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.</figcaption></figure>

### 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 (<10¹⁰ genome copies/pupae) or strongly (>10¹² 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 10⁶⁻⁸ DWV-A copies/bee and 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–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 www.mdpi.com/2076-2607/9/5/1087/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.

**Conflicts of Interest:** The authors declare no conflicts 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, doi:10.1264/jsme2.me17203.

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, doi:10.1038/nature20167.

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, doi:10.1371/journal.pbio.0030089.

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, doi:10.1038/nrmicro863.

5. Domingo, E.; Sheldon, J.; Perales, C. Viral Quasispecies Evolution. *Microbiol. Mol. Biol. Rev.* **2012**, *76*, 159–216, doi:10.1128/mmbr.05023-11.

6. Domingo, E.; Perales, C. Viral quasispecies. *PLoS Genet.* **2019**, *15*, e1008271, doi:10.1371/journal.pgen.1008271.

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, doi:10.1146/annurev-virology-101416-041718.

8. Fries, I.; Camazine, S. Implications of horizontal and vertical pathogen transmission for honey bee epidemiology. *Apidologie* **2001**, *32*, 199–214, doi:10.1051/apido:2001122.

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, doi:10.1098/rstb.2014.0116.

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, doi:10.1016/j.cub.2006.06.017.

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, doi:10.1016/j.cois.2018.02.009.

12. Gisder, S.; Genersch, E. Viruses of commercialized insect pollinators. *J. Invertebr. Pathol.* **2017**, *147*, 51–59, doi:10.1016/j.jip.2016.07.010.

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, doi:10.1128/aem.72.1.606-611.2006.

14. McMenamin, A.J.; Genersch, E. Honey bee colony losses and associated viruses. *Curr. Opin. Insect Sci.* **2015**, *8*, 121–129, doi:10.1016/j.cois.2015.01.015.

15. Sumpter, D.J.T.; Martin, S.J. The dynamics of virus epidemics in Varroa-infested honey bee colonies. *J. Anim. Ecol.* **2004**, *73*, 51–63, doi:10.1111/j.1365-2656.2004.00776.x.

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, doi:10.1371/journal.ppat.1002735.

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, doi:10.1073/pnas.0501860102.

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, doi:10.1073/pnas.1523515113.

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, doi:10.1126/science.aac9976.

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, doi:10.1371/journal.ppat.1004323.

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, doi:10.1126/science.1220941.

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, doi:10.1128/jvi.00158-17.

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, doi:10.3390/insects11040239.
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, doi:10.3389/fmicb.2020.00943.

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, doi:10.3390/v12050535.

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, doi:10.1007/s13592-019-00701-z.

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, doi:10.1086/501169.

28. Sherman, P.W.; Seeley, T.D.; Reeve, H.K. Parasites, Pathogens, and Polyandry in Social Hymenoptera. *Am. Nat.* **1988**, *131*, 602–610, doi:10.1086/284809.

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*, e0142985, doi:10.1371/journal.pone.0142985.

30. Simone-Finstrom, M.; Walz, M.; Tarpy, D.R. Genetic diversity confers colony-level benefits due to individual immunity. *Biol. Lett.* **2016**, *12*, 20151007, doi:10.1098/rsbl.2015.1007.

31. Kraus, F.B.; Page, R.E. Parasites, Pathogens and Polyandry in Social Insects. *Am. Nat.* **1998**, *151*, 383–391.

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, doi:10.3390/insects10090299.

33. Dorman, N. Citations. *Biotechnique* **2012**, *53*, 69, doi:10.2144/000113897.

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, doi:10.1038/srep41045.

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, doi:10.1128/aem.06094-11.

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, doi:10.3390/v9110314.

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, doi:10.1128/aem.70.12.7185-7191.2004.

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, doi:10.1099/vir.0.025965-0.

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, doi:10.1016/j.jviromet.2019.04.020.

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*, e1004230, doi:10.1371/journal.ppat.1004230.

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, doi:10.1128/aem.02227-09.

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, doi:10.1128/aem.06537-11.

44. Le Conte, Y.; De Vaublanc, G.; Crauser, D.; Jeanne, F.; Rousselle, J.-C.; Bécard, J.-M. Honey bee colonies that have survived Varroa destructor. *Apidologie* **2007**, *38*, 566–572, doi:10.1051/apido:2007040.

45. Oddie, M.A.Y.; Lanz, S.; Dahle, B.; Yañez, O.; Neumann, P. Virus infections in Varroa destructor-resistant honeybees. *bioRxiv* **2020**, doi:10.1101/2020.10.16.342246.

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, doi:10.1098/rspb.2016.0811.

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, doi:10.3390/v11050426.

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, doi:10.1007/s00114-003-0418-3.

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, doi:10.1007/s00040-009-0016-2.
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, doi:10.1007/s13592-011-0088-7.

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, doi:10.1038/srep33065.

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, doi:10.3390/v11040306.

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, doi:10.1371/journal.pone.0126330.

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, doi:10.3390/insects10120417.

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, doi:10.3390/microorganisms9040871.

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*, e99998, doi:10.1371/journal.pone.0099998.

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, doi:10.1371/journal.pone.0206938.

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, doi:10.1038/s41598-020-71388-x.

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, doi:10.3389/fmicb.2020.00620.