Next Article in Journal
Data on Herbivore Performance and Plant Herbivore Damage Identify the Same Plant Traits as the Key Drivers of Plant–Herbivore Interaction
Next Article in Special Issue
Geographical Distribution and Selection of European Honey Bees Resistant to Varroa destructor
Previous Article in Journal
Tick Preventive Behaviors and Practices Adopted by Medical Students from Poland, Germany, and Thailand in Relation to Socio-Demographic Conditions and Their Knowledge of Ticks and Tick-Borne Diseases
Previous Article in Special Issue
Substantial Genetic Progress in the International Apis mellifera carnica Population Since the Implementation of Genetic Evaluation
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Insects
Volume 11
Issue 12
10.3390/insects11120864
insects-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:
Communication

Selective Breeding for Low and High Varroa destructor Growth in Honey Bee (Apis mellifera) Colonies: Initial Results of Two Generations

by
Alvaro De la Mora
1,
Berna Emsen
1,2,
Nuria Morfin
1,
Daniel Borges
3,
Les Eccles
3,
Paul G. Kelly
1,
Paul H. Goodwin
1 and
Ernesto Guzman-Novoa
1,*
1
School of Environmental Sciences, University of Guelph, 50 Stone Road East, Guelph, ON N1G 2W1, Canada
2
Department of Animal Science, Agricultural Faculty, Ataturk University, Erzurum 25240, Turkey
3
Ontario Beekeepers’ Association Technology Transfer Program, 185, 5420 Hwy 6 N, Orchard Park Office, Guelph, ON N1H 6J2, Canada
*
Author to whom correspondence should be addressed.
Insects 2020, 11(12), 864; https://doi.org/10.3390/insects11120864
Submission received: 16 November 2020 / Revised: 27 November 2020 / Accepted: 3 December 2020 / Published: 4 December 2020
(This article belongs to the Special Issue Honeybee Breeding)
Browse Figure
Versions Notes

Abstract

:

Simple Summary

The mite Varroa destructor is considered the most damaging parasite of honey bees worldwide. Beekeepers use synthetic chemical products to control mite infestations in colonies, but the parasites soon develop resistance to them, which compromises their control. One alternative control strategy is the development of Varroa-resistant honey bees. Therefore, a breeding program was initiated to select for lower and higher rates of Varroa-population growth (LVG and HVG, respectively) and deformed wing virus (DWV) levels, which is transmitted by the mites. After two years of bidirectional selection, LVG colonies had Varroa population increases over the summer of 1.7 fold compared to 9.6 fold for HVG colonies. Additionally, HVG colonies had higher mite infestation rates in adult bees compared to LVG colonies. DWV presence and levels were higher in HVG colonies than in LVG colonies and winter mortality rates were 26% and 14% for the HVG and LVG bee types, respectively. The results of this study thus far indicate that selection for LVG may result in colonies with lower Varroa infestation rates, lower prevalence, and levels of DWV and higher colony winter survivorship. Future work will focus on determining mechanisms responsible for genetic differences and in identifying genes associated with Varroa-resistance in honey bees.

Abstract

After two years of bidirectional selection for low and high rates of Varroa destructor population growth (LVG and HVG, respectively) in honey bee (Apis mellifera) colonies in Ontario, Canada, significant differences between the two genotypes were observed. LVG colonies had V. destructor population increases over the summer of 1.7 fold compared to 9.6 fold for HVG colonies by Generation 2. Additionally, HVG colonies had significantly higher mite infestation rates in adult bees compared to LVG colonies for both selected generations. DWV prevalence and levels were significantly higher in HVG colonies than in LVG colonies in Generation 1 but not in Generation 2. Winter mortality rates of Generation 1 colonies were significantly different at 26% and 14% for the HVG and LVG genotypes, respectively. The results of this study thus far indicate that selection for LVG may result in colonies with lower V. destructor infestation rates, lower prevalence, and levels of DWV and higher colony winter survivorship. Future work will focus on determining what mechanisms are responsible for the genotypic differences, estimating genetic parameters, and molecular analyses of the genotypes to identify candidate genes associated with resistance to V. destructor and DWV that could potentially be used for marker-assisted selection.

1. Introduction

The mite Varroa destructor is considered the most damaging parasite of the western honey bee (Apis mellifera) worldwide. This mite has been frequently associated with the collapse of colonies in numerous countries [1,2,3,4,5]. V. destructor feeds upon the hemolymph and fat tissues of honey bees and transmits pathogenic viruses to them [6,7,8]. Additionally, V. destructor parasitism compromises immune responses, causes weight loss, shortens lifespan, and decreases honey yields of honey bees [9,10,11,12,13,14].
Another major factor associated with honey bee mortality is deformed wing virus (DWV), which is vectored by V. destructor [6,15]. DWV replicates in the mite, reaching high levels that are transmitted to the bees when parasitized by V. destructor [16]. DWV infections cause wing deformities and shorten the lifespan of bees [6,17]. Thus, both V. destructor and DWV are associated with honey bee mortality and colony losses. In locations with cold and long winters, such as Canada, over-wintering conditions favor the mites and DWV to cause more damage than in countries with less severe winters [18].
Most beekeepers use synthetic miticides to control V. destructor infestations in colonies, but the mites soon develop resistance to their active compounds, which has compromised their effectiveness [19]. Therefore, there is a need for alternative control strategies. One such strategy could be the development of Varroa-resistant strains of honey bees. Several attempts for breeding bees with increased resistance to V. destructor have been carried out in Europe and North America [20,21]. In North America, in particular, strains of bees showing some level of resistance to V. destructor have been produced in breeding programs to develop a Minnesota hygienic stock selected for increased hygienic behavior [22], a Russian honey bee stock selected for low mite infestation levels [23], a Varroa-sensitive hygiene (VSH) stock selected for removal of mite-infested brood [24,25], and an Indiana “mite-biter” honey bee stock selected for grooming behavior [26,27].
Developing a winter hardy, mite and virus tolerant stock would greatly contribute to a honey bee Integrated Pest Management (IPM) program, and is probably the most sustainable method of reducing colony losses due to the parasite. A previous study showed a strong correlation between lower V. destructor and lower DWV levels in honey bee colonies [28], suggesting that selecting for one trait may also select for the other. Therefore, a breeding program was initiated to select for lower rates of V. destructor population growth monitoring infection rates of DWV. This work was implemented in Ontario, Canada, with the collaboration of the Ontario Queen Breeders Association and the University of Guelph’s Honey Bee Research Center. Future work will examine the mechanisms that restrain V. destructor and DWV growth and the genes linked to increased resistance against these pathogens. This is a report of the initial results from two generations of bidirectional selection showing divergence in V. destructor population growth and DWV levels of a bidirectional selection study after two generations of selection.

2. Materials and Methods

2.1. Experimental Procedures

Experiments were conducted at apiaries of the Ontario Queen Breeders Association (OQBA) and the University of Guelph, ON, Canada. During the first year in 2018 (Generation 0), more than 300 colonies with queens of different genetic backgrounds equally represented (A. m. ligustica, A. m. carnica, and Buckfast strain) located in 15 apiaries were evaluated for fallen V. destructor mites twice [28]. The colonies were evaluated first in mid-spring (May) and then a second time 16 weeks later in late summer (September). Fallen mites were trapped on sticky papers composed of manila folders coated with Crisco® vegetable shortening placed underneath screened bottom boards of the hives. Averages of three-day fallen mite numbers were done twice in mid-spring and compared with those of late summer. The six colonies with the greatest proportional increase in mite numbers between the two evaluations were designated as high Varroa population growth (HVG) colonies, and the six colonies with the smallest proportional increase between the evaluations were designated as low Varroa population growth (LVG) colonies. These colonies were preselected to be used for larvae grafting the following year, expecting that at least three of each type would survive the winter. After assessments, all colonies were treated against V. destructor infestations using amitraz-impregnated plastic strips (Apivar®, Veto-Pharma, Palaiseau, France) as per the manufacturer.
In mid-spring of the following year (May 2019), three HVG and three LVG colonies selected in the previous year were used for larvae grafting to rear approx. 150 queens of each genotype. The six colonies used for larvae grafting were the three with the lowest and highest V. destructor population growth, respectively, in the previous season. Two original genetic backgrounds were represented in the queens used for grafting purposes. Two of the three HVG selected colonies were from A. m. carnica background and one from the Buckfast strain, whereas for the LVG selected colonies, two were from the Buckfast strain and one from A. m. carnica origin. The queens produced were allowed to open mate at a common mating yard isolated at least 5 km from other apiaries. Fifteen colonies each in 10 apiaries of the participating queen breeders were dequeened, and then each was divided into two hives. One HVG queen was introduced into one half of each divided hive and one LVG queen was introduced into the other half. This procedure ensured similar starting mite levels as well as colony and apiary conditions for both genotypes. In total, ca. 150 colonies were developed from queens of each genotype to produce bees of Generation 1. Fallen mite populations were determined with sticky papers as described above after establishing the colonies and then again in late summer. For Generation 1 in 2019, the six most HVG and six most LVG colonies were selected as per Generation 0. These procedures were repeated for a third year (Generation 2 in 2020). Winter mortality was visually assessed for colonies of Generation 1 the following spring after selection and evaluations (March 2020). Winter mortality for colonies of Generation 2 will be evaluated in the spring of 2021.

2.2. Varroa Destructor Infestation Rates in Bees

The selected colonies (the three HVG and three LVG colonies used for grafting) as well as 20 additional randomly chosen colonies of the two genotypes were sampled during late summer to determine mite infestation rates and DWV prevalence and levels. For mite infestation rates, approx. 300 adult bees were collected from each colony in jars containing 70% ethanol and the number of mites per 100 bees was calculated for each sample as per Dietemann et al. [29]. Adults and brood from the same colonies were collected and immediately frozen at −70 °C for future assessment of brood infestation as well as for detection and quantification of DWV.

2.3. DWV Identification and Quantification

Total RNA was extracted from samples of 10 bees per colony using One Step-RNA reagent (Bio Basic, Markham, ON, Canada) following the manufacturer’s instructions. For cDNA synthesis, 2 μg of total RNA was reverse-transcribed using the RevertAid H Minus First Strand cDNA Synthesis Kit (Fermentas Life Sciences, Burlington, ON, Canada), following the instructions of the manufacturer.
DWV was amplified and quantified in qRT-PCR reactions. The primer sequences and qRT-PCR protocols used to identify DWV type A were those by Di Prisco et al. [30]. PCR reactions were done with a BioRad CFX96™ thermocycler (Bio-Rad Laboratories, Mississauga, ON, Canada) and PowerUp™SYBRgreen™(2X) (Applied Biosystems, Foster City, CA, USA) on 96-well plates (Bio-Rad Laboratories, Mississauga, ON, Canada). Each qPCR reaction of 20 µL contained 2 µL of cDNA, 10 µL of SYBRgreen™, 0.4 µL of reverse and forward DWV primers (200 nM), and 7.2 µL of nuclease free H2O. A negative control was included in each run by adding 2 µL of nuclease free H2O instead of cDNA, and a positive control from previously identified DWV positive bee sample by qRT-PCR was also included. Each reaction consisted of one cycle at 48 °C for 15 min, one at 95 °C for 10 min, 40 cycles at 95 °C for 15 s and 60 °C for 60 s, followed by one cycle at 68 °C for 7 min. Calibration curves to convert Ct values to DWV gc were done using 300 bp gblocks® (Integrated DNA Technologies, Coralville, IA, USA) that included the sequence of the forward primer, amplicon, and reverse primer. The lyophilized gBlock® was diluted with 50 µL of ds H2O to obtain an initial concentration of 10 ng/µL that was used to make serial dilutions from 109 to 101 copies [14]. Using a plot of Ct values versus viral copy number (log10), a linear equation was used to calculate the DWV genome copy numbers. Ten-fold serial dilutions using gblocks® were also used to optimize the qRT-PCR reaction to obtain an efficiency of 95–105%. Three technical repetitions were done for each qRT-PCR run. Randomly selected amplicons of presumed DWV were sequenced at the University of Guelph Laboratory services to confirm identity.

2.4. Statistical Analyses

The data on fold changes in mite population growth were log-transformed as the data did not comply with normality based on the Shapiro Wilk test. The data on mite infestation rates were square-root transformed and the data on DWV copies were log-transformed for the same reasons. The transformed data were subjected to analyses of variance and Fisher LSD tests to separate means when significance was detected. Selection differentials and responses to selection were calculated. Differences in DWV prevalence were tested with χ2 tests. All statistical analyses were performed using the R statistical program [31].

3. Results

3.1. V. destructor Population Growth

Mite populations in LVG colonies that are presumably more resistant to V. destructor increased 1.7 fold by Generation 2 of selection, which was nearly six times lower than the 9.6-fold mite increase in HVG colonies that are presumed to be more susceptible to V. destructor (Figure 1). Significant differences were found for genotype (F1,646 = 59.3, p < 0.0001), generation (F2,646 = 3.2, p < 0.05), and genotype × generation interaction (F2,646 = 18.9, p < 0.0001). Selection differentials (SD) for the colonies used for larvae grafting were 45.5 ± 13.1 and 21.5 ± 11.5 for Generations 1 and 2, respectively, whereas the responses to selection were 6.6 ± 0.8 and 8.0 ± 1.3 for Generations 1 and 2, respectively. Additionally, HVG colonies had significantly higher mite infestation rates in adult bees compared to LVG colonies in the two selected generations (Table 1). There were significant effects of genotype (F1,76 = 65.2, p < 0.0001), generation (F1,76 = 6.7, p < 0.05), and genotype × generation interaction (F1,76 = 5.7, p < 0.05) for this trait.

3.2. DWV Prevalence and Levels

DWV prevalence was significantly higher in HVG colonies than in LVG colonies in Generation 1 (95% versus 65%, respectively, χ2 = 5.6, p < 0.05) but not in Generation 2 (55% versus 30%, respectively, χ2 = 2.6, p > 0.05). Similarly, the level of DWV was significantly higher in HVG colonies than in LVG colonies in Generation 1, but not in Generation 2 where there were no significant differences (Table 1). For DWV levels, significant differences were found for genotype (F1,45 = 5.7, p < 0.05) and generation (F1,45 = 34.4, p < 0.0001) but not for genotype × generation interaction (F1,45 = 2.2, p > 0.05).

3.3. Winter Colony Mortality

Winter mortality for Generation 1 colonies was 26% and 14% for the HVG and LVG genotypes, respectively, which differed significantly (χ2 = 4.5, p < 0.05).

4. Discussion

The initial results of this study clearly show that selecting for LVG resulted in significantly lower V. destructor infestation rates in colonies than selecting for HVG after two generations of selection. Furthermore, significantly higher winter colony survivorship rates were observed in colonies of the LVG than the HVG genotype. Additionally, LVG colonies had significantly lower prevalence and levels of DWV after one generation of selection than HVG colonies, but these did not significantly differ in the second generation of selection. It is possible that the DWV levels had been affected not only by the selection procedures on the maternal side, but also by the source of drones that mated with the selected queens. It is known that DWV can be transmitted through drone semen [32], which could affect the results of selection since the paternal contribution was not controlled in this study.
The above results support the hypothesis that a strategy of selecting for reduced V. destructor levels in colonies also may have reduced DWV levels and prevalence (in Generation 1), as well as decreased winter mortality of colonies. Thus, selecting for LVG colonies provides multiple benefits to honey bee colonies. These results agree with those of previous work showing decreased rates of V. destructor population growth in colonies selected for increased resistance to Varroa compared to a nonselected population of Italian bees [27], or bees selected for increased growth of mite populations [33].
Many variables influence V. destructor population growth, including climatic factors and genetic factors of the mites and host bees [34,35,36,37,38,39,40,41]. Additionally, infestation rates and V. destructor population growth are affected by reinfestation from drifting bees [42,43]. Therefore, future work should investigate the impact of such factors on the V. destructor populations in LVG and HVG colonies.
If the bee genotype is an important contributor to the genetic variability in mite population growth in the colonies, then it should be possible to achieve higher resistance to Varroa in selected populations of honey bees by breeding for LVG. Previous studies have shown that V. destructor population growth in colonies is a selectable and heritable trait. For example, Lodesani et al. [33] found an estimate of h2 = 0.84 for LVG in adult bees, whereas Harbo and Harris [38] found a h2 of 1.24 ± 0.49 for the proportion of mites in brood of honey bee colonies in which mite populations were allowed to grow for 10 weeks. More recently, Maucourt et al. [44] also estimated a h2 = 0.44 ± 0.56 for rates of V. destructor infestation in Canadian honey bee colonies. Future work will be done to determine the heritability estimates of the traits measured in this study. Furthermore, various statistical methods will be used to compare LVG and HVG genotypes when data from additional generations are available. It is possible that the standard errors calculated with ANOVA for the first generation of selection are continued into successive generations, which would then bias the estimation of breeding parameters.
The causes of different rates of infestation with V. destructor or infection with DWV between LVG and HVG bees remain to be discovered, and thus work is underway to determine if the divergence between the LVG and HVG genotypes for these traits continues and is sustained over additional generations of selection, and to determine the potential factors and mechanisms involved in causing the divergence in mite population growth between the two selected genotypes. Furthermore, other parameters associated with V. destructor resistance can be measured on the two genotypes, such as brood versus adult bee V. destructor infestation [38,41,44,45]. Finally, the molecular and cellular basis for lower V. destructor populations, DWV prevalence, and DWV levels between LVG and HVG colonies will be later examined, such as by genome sequencing to identify single nucleotide polymorphisms (SNPs) and high-throughput transcript sequencing to identify changes in gene expression related to resistance. This may also allow for the identification of targets for marker-assisted selection for LVG in the future.

5. Conclusions

The results of this study thus far are consistent with the idea that selection for LVG can result in colonies having lower V. destructor infestation rates, lower DWV prevalence and levels, and higher winter colony survivorship. In 2021, there will be further selection of a third generation of bees for LVG and HVG, as well as analysis of factors that could help explain differences in V. destructor population growth in the selected genotypes, as well as estimations of genetic parameters. There will also be a focus on molecular analysis of bees from the third generation of both bee genotypes to pinpoint candidate genes that could be associated with resistance to V. destructor and DWV that could be potentially used for marker-assisted selection.

Author Contributions

Conceptualization, E.G.-N. and P.H.G.; methodology, A.D.l.M., B.E., N.M., D.B., L.E. and P.G.K.; formal analyses, A.D.l.M. and E.G.-N.; investigation, A.D.l.M., B.E., N.M., D.B. and L.E.; resources, E.G.-N.; original draft preparation, A.D.l.M.; review and editing, A.D.l.M., N.M., E.G.-N. and P.H.G.; supervision, L.E., P.G.K., E.G.-N. and P.H.G.; project administration, E.G.-N.; funding acquisition, E.G.-N. and P.H.G. All authors have read and agreed to the published version of the manuscript.

Funding

This study was partially funded by the Ontario Ministry of Agriculture, Food and Rural Affairs (New Directions Grant No. 2017-3142).

Acknowledgments

We thank the participating queen breeders from the Ontario Bee Breeders Association for their assistance and use of their colonies. We also thank members of the Ontario Tech Transfer Program for their help with colony sampling and evaluations. Logistic and experimental work was also provided by Stephanie Otto, Catherine Vanderheyden, Wendy Shipsides, and Pia Marquardt-Salathe.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Van Engelsdorp, D.; Evans, J.D.; Sagerman, C.; Mullin, C.; Haubruge, E.; Nguyen, B.K.; Frazier, M.; Frazier, J.; Diana Cox-Foster, D.; Chen, Y.; et al. Colony collapse disorder: A descriptive study. PLoS ONE 2009, 4, e6481. [Google Scholar]
  2. Guzman-Novoa, E.; Eccles, L.; Calvete, Y.; McGowan, J.; Kelly, P.G.; Correa-Benítez, A. Varroa destructor is the main culprit for the death and reduced populations of overwintered honey bee (Apis mellifera) colonies in Ontario, Canada. Apidologie 2010, 41, 443–450. [Google Scholar] [CrossRef] [Green Version]
  3. Le Conte, Y.; Ellis, M.; Ritter, W. Varroa mites and honey bee health: Can Varroa explain part of the colony losses? Apidologie 2010, 41, 353–363. [Google Scholar] [CrossRef] [Green Version]
  4. Rosenkranz, P.; Aumeier, P.; Ziegelmann, B. Biology and control of Varroa destructor. J. Invertebr. Pathol. 2010, 103, 96–119. [Google Scholar] [CrossRef] [PubMed]
  5. Meixner, M.D.; Francis, R.M.; Gajda, A.; Kryger, P.; Andonov, S.; Uzunov, A.; Topolska, G.; Costa, C.; Amiri, E.; Berg, S.; et al. Occurrence of parasites and pathogens in honey bee colonies used in a European genotype-environment interactions experiment. J. Apic. Res. 2014, 53, 215–229. [Google Scholar] [CrossRef]
  6. Genersch, E.; Aubert, M. Emerging and re-emerging viruses of the honey bee (Apis mellifera L.). Vet. Res. 2010, 41, 54. [Google Scholar] [CrossRef] [Green Version]
  7. Ramsey, S.D.; Ochoa, R.; Bauchan, G.; Gulbronson, C.; Mowery, J.D.; Cohen, A.; Lim, D.; Joklik, J.; Cicero, J.M.; Ellis, J.D.; et al. Varroa destructor feeds primarily on honey bee fat body tissue and not hemolymph. Proc. Natl. Acad. Sci. USA 2019, 116, 1792–1801. [Google Scholar] [CrossRef] [Green Version]
  8. Anguiano-Baez, R.; Guzman-Novoa, E.; Hamiduzzaman, M.M.; Espinosa-Montaño, L.G.; Correa-Benítez, A. Varroa destructor (Mesostigmata: Varroidae) parasitism and climate differentially influence the prevalence, levels and overt infections of deformed wing virus in honey bees (Hymenoptera: Apidae). J. Insect Sci. 2016, 16, 44. [Google Scholar] [CrossRef] [Green Version]
  9. De Jong, D.; De Jong, P.H.; Goncalves, L.S. Weight loss and other damage to developing worker honey bees (Apis mellifera) due to infestation with Varroa jacobsoni. J. Apic. Res. 1982, 21, 165–167. [Google Scholar] [CrossRef]
  10. Gregory, P.G.; Evans, J.D.; Rinderer, T.; de Guzman, L. Conditional immune-gene suppression of honey bees parasitized by Varroa mites. J. Insect Sci. 2005, 5, 7–12. [Google Scholar] [CrossRef]
  11. 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]
  12. Murilhas, A. Varroa destructor infestation impact on Apis mellifera carnica capped brood production, bee population and honey storage in a Mediterranean climate. Apidologie 2002, 33, 271–281. [Google Scholar] [CrossRef]
  13. Koleoglu, G.; Goodwin, P.H.; Reyes-Quintana, M.; Hamiduzzaman, M.M.; Guzman-Novoa, E. Varroa destructor parasitism reduces hemocyte concentrations and prophenol oxidase gene expression in bees from two populations. Parasitol. Res. 2018, 117, 1175–1183. [Google Scholar] [CrossRef] [Green Version]
  14. Morfin, N.; Goodwin, P.H.; Guzman-Novoa, E. Interaction of Varroa destructor and sublethal clothianidin doses during the larval stage on subsequent adult honey bee (Apis mellifera L.) health, cellular immunity, deformed wing virus levels and differential gene expression. Microorganisms 2020, 8, 858. [Google Scholar] [CrossRef]
  15. Dainat, B.; Evans, J.D.; Chen, Y.P.; Gauthier, L.; Neumann, P. Predictive markers of honey bee colony collapse. PLoS ONE 2012, 7, e32151. [Google Scholar] [CrossRef]
  16. Gisder, S.; Aumeier, P.; Genersch, E. Deformed wing virus: Replication and viral load in mites (Varroa destructor). J. Gen. Virol. 2009, 90, 463–467. [Google Scholar] [CrossRef] [PubMed]
  17. Reyes-Quintana, M.; Espinoza-Montano, L.G.; Prieto, D.; Koleoglu, G.; Petukhova, T.; Correa-Benitez, A.; Guzman-Novoa, E. Impact of Varroa destructor and deformed wing virus on emergence, cellular immunity, wing integrity and survivorship of Africanized honey bees in Mexico. J. Invertebr. Pathol. 2019, 164, 43–48. [Google Scholar] [CrossRef] [PubMed]
  18. Currie, R.W.; Pernal, S.F.; Guzman-Novoa, E. Honey bee colony loses in Canada. J. Apic. Res. 2010, 49, 104–106. [Google Scholar] [CrossRef]
  19. Wallner, K. Varroacides and their residues in bee products. Apidologie 1999, 30, 235–248. [Google Scholar] [CrossRef] [Green Version]
  20. Büchler, R.; Berg, S.; Le Conte, Y. Breeding for resistance to Varroa destructor in Europe. Apidologie 2010, 41, 393–408. [Google Scholar] [CrossRef] [Green Version]
  21. Rinderer, T.E.; Harris, J.W.; Hunt, G.J.; de Guzman, L.I. Breeding for resistance to Varroa destructor in North America. Apidologie 2010, 41, 409–424. [Google Scholar] [CrossRef] [Green Version]
  22. Spivak, M.; Reuter, G.S.; Lee, K.; Ranum, B. The future of the MN Hygienic stock of bees in good hands. Am. Bee J. 2009, 149, 965–967. [Google Scholar]
  23. Rinderer, T.E.; de Guzman, L.I.; Delatte, G.T.; Stelzer, J.A.; Lancaster, V.A.; Kuznetsov, V.; Beaman, L.; Watts, R.; Harris, J.W. Resistance to the parasitic mite Varroa destructor in honey bees from far-eastern Russia. Apidologie 2001, 32, 381–394. [Google Scholar] [CrossRef] [Green Version]
  24. Harbo, J.R.; Harris, J.W. Suppressed mite reproduction explained by the behavior of adult bees. J. Apic. Res. 2005, 44, 21–23. [Google Scholar] [CrossRef]
  25. Danka, R.G.; Harris, J.W.; Dodds, G.E. Selection of VSH-derived “Pol-line” honey bees and evaluation of their Varroa-resistance characteristics. Apidologie 2016, 47, 483–490. [Google Scholar] [CrossRef]
  26. Hunt, G.; Given, J.K.; Tsuruda, J.M.; Andino, G.K. Breeding mite-biting bees to control Varroa. Bee Cult. 2016, 8, 41–47. [Google Scholar]
  27. Morfin, N.; Given, K.; Evans, M.; Guzman-Novoa, E.; Hunt, G.J. Grooming behavior and gene expression of the Indiana “mite-biter” honey bee stock. Apidologie 2019. [Google Scholar] [CrossRef] [Green Version]
  28. Emsen, B.; Hamiduzzaman, M.M.; Goodwin, P.H.; Guzman-Novoa, E. Lower virus infections in Varroa destructor-infested and uninfested brood and adult honey bees (Apis mellifera) of a low mite population growth colony compared to a high mite population growth colony. PLoS ONE 2015, 10, e0118885. [Google Scholar] [CrossRef] [Green Version]
  29. Dietemann, V.; Nazzi, F.; Martin, S.J.; Anderson, D.L.; Locke, B.; Delaplane, K.S.; Wauquiez, Q.; Tannahill, C.; Frey, E.; Ziegelmann, B.; et al. Standard methods for varroa research. J. Apic. Res. 2013, 52, 1–54. [Google Scholar] [CrossRef] [Green Version]
  30. Di Prisco, G.; Cavaliere, V.; Annoscia, D.; Varricchio, P.; Caprio, E.; Nazzi, F.; Gargiulo, G.; Pennacchio, F. Neonicotinoid clothianidin adversely affects insect immunity and promotes replication of a viral pathogen in honey bees. Proc. Natl. Acad. Sci. USA 2013, 110, 18466–18471. [Google Scholar] [CrossRef] [Green Version]
  31. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2020; Available online: http://www.R-project.org/ (accessed on 10 August 2020).
  32. 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] [PubMed]
  33. Lodesani, M.; Crailsheim, K.; Moritz, R.F.A. Effect of some characters on the population growth of mite Varroa jacobsoni in Apis mellifera L colonies and results of a bi-directional selection. J. Appl. Entomol. 2002, 126, 130–137. [Google Scholar] [CrossRef] [Green Version]
  34. Le Conte, Y.; Arnold, G.; Desenfant, P.H. Influence of brood temperature and hygrometry variation on the development of the honey bee ectoparasite Varroa jacobsoni (Mesostigmata: Varroidae). Environ. Entomol. 1990, 19, 1780–1785. [Google Scholar] [CrossRef]
  35. Guzman-Novoa, E.; Sánchez, A.; Page, R.E., Jr.; García, T. Susceptibility of European and Africanized honeybees (Apis mellifera L) and their hybrids to Varroa jacobsoni Oud. Apidologie 1996, 27, 93–103. [Google Scholar] [CrossRef] [Green Version]
  36. Guzman-Novoa, E.; Emsen, B.; Unger, P.; Espinosa-Montaño, L.G.; Petukhova, T. Genotypic variability and relationships between mite infestation levels, mite damage, grooming intensity, and removal of Varroa destructor mites in selected strains of worker honey bees (Apis mellifera L.). J. Invertebr. Pathol. 2012, 110, 314–320. [Google Scholar] [CrossRef]
  37. de Guzman, L.I.; Rinderer, T.E.; Stelzer, J.A. Ocurrence of two genotypes of Varroa jacobsoni Oud. in North America. Apidologie 1999, 30, 31–36. [Google Scholar] [CrossRef]
  38. Harbo, J.R.; Harris, J.W. Heritability in honey bees (Hymenoptera: Apidae) of characteristics associated with resistance to Varroa jacobsoni (Mesostigmata: Varroidae). J. Econ. Entomol. 1999, 92, 261–265. [Google Scholar] [CrossRef] [Green Version]
  39. Harris, J.W.; Harbo, J.R.; Villa, J.D.; Danka, R.G. Variable population growth of Varroa destructor (Msostigmata: Varroidae) in colonies of honey bees (Hymenoptera: Apidae) during a 10-year period. Environ. Entomol. 2003, 32, 1305–1312. [Google Scholar] [CrossRef]
  40. Medina-Flores, C.A.; Guzman-Novoa, E.; Hamiduzzaman, M.M.; Aréchiga-Flores, C.F.; López-Carlos, M.A. Africanized honey bees (Apis mellifera) have low infestation levels of the mite Varroa destructor in different ecological regions in Mexico. Gen. Mol. Res. 2014, 13, 7282–7293. [Google Scholar] [CrossRef]
  41. Russo, R.M.; Liendo, M.C.; Landi, L.; Pietronave, H.; Merke, J.; Fain, H.; Muntaabski, I.; Palacio, M.A.; Rodríguez, G.A.; Lanzavecchia, S.B.; et al. Grooming behavior in naturally Varroa-resistant Apis mellifera colonies from North-Central Argentina. Front. Ecol. Evol. 2020, 8, 590281. [Google Scholar] [CrossRef]
  42. DeGrandi-Hoffman, G.; Ahumada, F.; Zazueta, V.; Chambers, M.; Hidalgo, G.; Watkins deJong, E. Population growth of Varroa destructor (Acari: Varroidae) in honey bee colonies is affected by the number of foragers with mites. Exp. Appl. Acarol. 2016, 69, 21–34. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Nolan, M.P.; Delaplane, K.S. Distance between honey bee Apis mellifera colonies regulates populations of Varroa destructor at a landscape scale. Apidologie 2017, 48, 8–16. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Maucourt, S.; Fortin, F.; Robert, C.; Giovenazzo, P. Genetic parameters of honey bee colonies traits in a Canadian selection program. Insects 2020, 11, 587. [Google Scholar] [CrossRef] [PubMed]
  45. Büchler, R.; Kovacic, M.; Buchegger, M.; Puskadija, Z.; Hoppe, A.; Brascamp, E.W. Evaluation of traits for the selection of Apis mellifera against Varroa destructor. Insects 2020, 11, 618. [Google Scholar] [CrossRef] [PubMed]
Insects 11 00864 g001 550
Figure 1. Mean Varroa destructor population growth (fold change ± SE) per colony of two honey bee genotypes selected for high and low Varroa population growth (HVG and LVG, respectively) during two generations (n > 100 per genotype in each generation). Different letters indicate significant differences within and between genotypes.
Figure 1. Mean Varroa destructor population growth (fold change ± SE) per colony of two honey bee genotypes selected for high and low Varroa population growth (HVG and LVG, respectively) during two generations (n > 100 per genotype in each generation). Different letters indicate significant differences within and between genotypes.
Insects 11 00864 g001
Table
Table 1. Percent adult bee infestation with V. destructor (± SE) and mean log copies. of deformed wing virus (DWV) (±SE) found in HVG and LVG colonies after two generations of selection. Different letters after means and SE, indicate significant differences between genotypes.
Table 1. Percent adult bee infestation with V. destructor (± SE) and mean log copies. of deformed wing virus (DWV) (±SE) found in HVG and LVG colonies after two generations of selection. Different letters after means and SE, indicate significant differences between genotypes.
GenotypeGeneration% InfestationLog Copies of DWV
HVGOne9.61 ± 1.08 a8.14 ± 0.42 a
LVGOne5.15 ± 0.50 b6.23 ± 0.49 b
HVGTwo9.15 ± 0.91 a4.51 ± 0.37 c
LVGTwo2.62 ± 0.42 c4.07 ± 0.15 c
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

De la Mora, A.; Emsen, B.; Morfin, N.; Borges, D.; Eccles, L.; Kelly, P.G.; Goodwin, P.H.; Guzman-Novoa, E. Selective Breeding for Low and High Varroa destructor Growth in Honey Bee (Apis mellifera) Colonies: Initial Results of Two Generations. Insects 2020, 11, 864. https://doi.org/10.3390/insects11120864

AMA Style

De la Mora A, Emsen B, Morfin N, Borges D, Eccles L, Kelly PG, Goodwin PH, Guzman-Novoa E. Selective Breeding for Low and High Varroa destructor Growth in Honey Bee (Apis mellifera) Colonies: Initial Results of Two Generations. Insects. 2020; 11(12):864. https://doi.org/10.3390/insects11120864

Chicago/Turabian Style

De la Mora, Alvaro, Berna Emsen, Nuria Morfin, Daniel Borges, Les Eccles, Paul G. Kelly, Paul H. Goodwin, and Ernesto Guzman-Novoa. 2020. "Selective Breeding for Low and High Varroa destructor Growth in Honey Bee (Apis mellifera) Colonies: Initial Results of Two Generations" Insects 11, no. 12: 864. https://doi.org/10.3390/insects11120864

APA Style

De la Mora, A., Emsen, B., Morfin, N., Borges, D., Eccles, L., Kelly, P. G., Goodwin, P. H., & Guzman-Novoa, E. (2020). Selective Breeding for Low and High Varroa destructor Growth in Honey Bee (Apis mellifera) Colonies: Initial Results of Two Generations. Insects, 11(12), 864. https://doi.org/10.3390/insects11120864

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