Genetic diversity within honey bee colonies affects pathogen load and relative virus levels in honey bees, Apis mellifera L

Neonicotinoid exposure increases Varroa destructor (Mesostigmata: Varroidae) mite parasitism severity in honey bee colonies and is not mitigated by increased colony genetic diversity

Agrochemical exposure is a major contributor to ecological declines worldwide, including the loss of crucial pollinator species. In addition to direct toxicity, field-relevant doses of pesticides can increase species' vulnerabilities to other stressors, including parasites. Experimental field demonstrations of potential interactive effects of pesticides and additional stressors are rare, as are tests of mechanisms via which pollinators tolerate pesticides. Here, we controlled honey bee colony exposure to field-relevant concentrations of 2 neonicotinoid insecticides (clothianidin and thiamethoxam) in pollen and simultaneously manipulated intracolony genetic heterogeneity. We showed that exposure increased rates of Varroa destructor (Anderson and Trueman) parasitism and that while increased genetic heterogeneity overall improved survivability, it did not reduce the negative effect size of neonicotinoid exposure. This study is, to our knowledge, the first experimental field demonstration of how neonicotinoid exposure can increase V. destructor populations in honey bees and also demonstrates that colony genetic diversity cannot mitigate the effects of neonicotinoid pesticides.

Seasonal trends of the ABPV, KBV, and IAPV complex in Italian managed honey bee (Apis mellifera L.) colonies

Acute bee paralysis virus (ABPV), Kashmir bee virus (KBV), and Israeli acute paralysis virus (IAPV) usually persist as covert infections in honey bee colonies. They can cause rapid bee mortality in cases of severe infection, often associated with high Varroa destructor infestation, by which they are transmitted. In various countries, these viruses have been associated with colony collapse. Despite their potential danger, these viruses are often disregarded, and little information is available on their occurrence in many countries, including Italy. In 2021, 370 apiaries representing all of the Italian regions were investigated in four different months (June, September, November, and March) for the presence of ABPV, KBV, and IAPV. IAPV was not found in any of the apiaries investigated, whereas 16.45% and 0.67% of the samples tested positive for ABPV and KBV, respectively. Most ABPV cases occurred in late summer-autumn in both northern and southern regions. We observed a scattered pattern of KBV-positive colonies that did not allow any seasonal or regional trends to be discerned. Differences observed among regions and months were potentially related to the dynamics of varroa infestation, viral genetic variations, and different climatic conditions resulting in variations in bee behaviour. This study improves our understanding of the circulation of bee viruses and will contribute to better disease prevention and preservation of bee health.

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

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

Viral Co-Infections and Antiviral Immunity in Honey Bees

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

Observation of Genetic Gain with Instrumental Insemination of Honeybee Queens

Controlling mating in the honeybee (Apis mellifera) is part of one of the greatest challenges for the beekeeping industry's genetic selection programs due to specific characteristics of their reproduction. Several techniques for supervising honeybee mating with relative effective control have been developed over the years to allow honeybee selection. As part of this project, we compared the genetic gains for several colony performance traits, obtained using the BLUP-animal method, according to the selection pressure applied in controlled reproduction (directed fertilization versus instrumental insemination). Our results show similar genetic gains for hygienic behavior and honey production between colonies whether queens were fertilized naturally or via instrumental insemination, as well as similar or lower genetic gains for colonies with queens inseminated for spring development. In addition, we noticed greater fragility in queens following insemination. These findings show that instrumental insemination is an effective tool for reproductive control in genetic selection and for estimating breeding values more precisely. However, this technique does not result in queens of superior genetic quality for commercial purposes.

Meta-analysis of genetic diversity and intercolony relatedness among reproductives in commercial honey bee populations

Honey bee colonies are large kin groups, each with a single mother queen and thousands of female workers. Queen bees are highly polyandrous, each mating with an average of approximately 12 drones from other colonies. We used a meta-analysis approach to compare the pedigree relationships of honey bee reproductives (queens and their mates) across five different studies and to quantify the overall genetic diversity of breeding populations. We compared the inferred genotypes of queens and their mates from microsatellite analyses of worker offspring from a feral Africanized honey bee population (which served as a negative control for inbreeding), an experimentally derived population of sister queens (which served as a positive control for inbreeding), and three separate commercially managed populations. We then compared the relatedness of all drones mated to each queen (mate-mate), all queens within each population (queen-queen), each queen with each of her mates (queen-mate), and all drones within each population (drone-drone). We found, as expected, the lowest levels of genetic similarity in the outcrossed population and highest levels of genetic similarity in the inbred population. Levels of genetic similarity among the managed honey bee populations were intermediate but closer to that of the inbred population. Genetic structuring of the entire breeding population resulted in two major subpopulations, likely deriving from breeders on the east and west coast. The effects that these findings have on the overall population genetic diversity of managed honey bees is discussed.

Morphometric and genetic characterization as tools for selection of Apis mellifera (Hymenoptera: Apidae) stocks in an area of natural hybridization in Argentina

Beekeepers around the world select bees' characteristics that facilitate and favor production. In regions where hybridization among lineages is taking place, this selection is a challenge, given that these regions are "natural laboratories", where the action of evolutionary processes of a population or species occurs in real time. A natural honeybee (Apis mellifera) hybrid zone exists in Argentina between 28° and 35° South, where Africanized (AHB) and European (EHB) populations converge. In this zone, beekeepers use selected genetic resources of European origin mostly, since the local Africanized bees show a higher defensive behavior, which is not desirable for management. Although EHB colonies have many advantages for honey production, they are not fully adapted to the subtropical climate and are susceptible to certain parasitosis such as varroosis. In addition, both AHB and EHB mate in drone congregation areas (DCAs), where males and virgin queens fly to meet, resulting in variability in the desired characteristics. In this study, we explored the degree of hybridization within a DCA and its reference apiary, located in the province of Entre Ríos, by applying two complementary techniques. First, morphotypes with different degrees of hybridization between European and African subspecies were observed in the reference apiary, indicating a high sensitivity of this morphometric approach to detect hybridization in these populations. Second, a genetic analysis revealed haplotypes of both origins for drones in DCAs, with a higher prevalence of European haplotypes, while all the colonies from the reference apiary exhibited European haplotypes. Overall, our results are in line with the strong impact that commercial beekeeping has on the genetics of DCAs. We show how wing morphometry may be used to monitor hybridization between European and African subspecies, a tool that may be evaluated in other regions of the world where hybridization occurs.

Genetic slippage after sex maintains diversity for parasite resistance in a natural host population

Although parasite-mediated selection is a major driver of host evolution, its influence on genetic variation for parasite resistance is not yet well understood. We monitored resistance in a large population of the planktonic crustacean Daphnia magna over 8 years, as it underwent yearly epidemics of the bacterial pathogen Pasteuria ramosa. We observed cyclic dynamics of resistance: Resistance increased throughout the epidemics, but susceptibility was restored each spring when hosts hatched from sexual resting stages. Host resting stages collected across the year showed that largely resistant host populations can produce susceptible sexual offspring. A genetic model of resistance developed for this host-parasite system, based on multiple loci and strong epistasis, is in partial agreement with our findings. Our results reveal that, despite strong selection for resistance in a natural host population, genetic slippage after sexual reproduction can be a strong factor for the maintenance of genetic diversity of host resistance.

The ontogeny of selection on genetic diversity in harvester ants

Selection may favour traits throughout an individual's lifetime or at a particular life stage. In many species of social insects, established colonies that are more genetically diverse outperform less diverse colonies with respect to a variety of traits that contribute to fitness, but whether selection favours high diversity in small colonies is unknown. We tested the hypothesis that selection favours genetically diverse colonies during the juvenile period using a multi-year field experiment with the harvester ant, Pogonomyrmex occidentalis. We used controlled matings to generate colonies that varied in genetic diversity and transplanted them into the field. We monitored their survival for seven (the 2015 cohort, n = 149) and six (the 2016 cohort, n = 157) years. Genetically more diverse colonies had greater survival, resulting in significant viability selection. However, in both cohorts survival was not influenced by genetic diversity until colonies were three years old. We suggest that changes in their internal organization enabled colonies to use the benefits of multiple genotypes, and discuss possible mechanisms that can generate this pattern.

Phenomic analysis of the honey bee pathogen-web and its dynamics on colony productivity, health and social immunity behaviors

Many pathogens and parasites have evolved to overwhelm and suppress their host's immune system. Nevertheless, the interactive effects of these agents on colony productivity and wintering success have been relatively unexplored, particularly in large-scale phenomic studies. As a defense mechanism, honey bees have evolved remarkable social behaviors to defend against pathogen and parasite challenges, which reduce the impact of disease and improve colony health. To investigate the complex role of pathogens, parasites and social immunity behaviors in relation to colony productivity and outcomes, we extensively studied colonies at several locations across Canada for two years. In 2016 and 2017, colonies founded with 1-year-old queens of diverse genetic origin were evaluated, which represented a generalized subset of the Canadian bee population. During each experimental year (May through April), we collected phenotypic data and sampled colonies for pathogen analysis in a standardized manner. Measures included: colony size and productivity (colony weight, cluster size, honey production, and sealed brood population), social immunity traits (hygienic behavior, instantaneous mite population growth rate, and grooming behavior), as well as quantification of gut parasites (Nosema spp., and Lotmaria passim), viruses (DWV-A, DWV-B, BQCV and SBV) and external parasites (Varroa destructor). Our goal was to examine: 1) correlations between pathogens and colony phenotypes; 2) the dynamics of pathogens and parasites on colony phenotypes and productivity traits; and 3) the effects of social immunity behaviors on colony pathogen load. Our results show that colonies expressing high levels of some social immunity behaviors were associated with low levels of pathogens/parasites, including viruses, Nosema spp., and V. destructor. In addition, we determined that elevated viral and Nosema spp. levels were associated with low levels of colony productivity, and that five out of six pathogenic factors measured were negatively associated with colony size and weight in both fall and spring periods. Finally, this study also provides information about the incidence and abundance of pathogens, colony phenotypes, and further disentangles their inter-correlation, so as to better understand drivers of honey bee colony health and productivity.

Direct evidence for increased disease resistance in polyandrous broods exists only in eusocial Hymenoptera

Background

The ‘genetic diversity’ hypothesis posits that polyandry evolved as a mechanism to increase genetic diversity within broods. One extension of this hypothesis is the ‘genetic diversity for disease resistance’ hypothesis (GDDRH). Originally designed for eusocial Hymenoptera, GDDRH states that polyandry will evolve as an effect of lower parasite prevalence in genetically variable broods. However, this hypothesis has been broadly applied to several other taxa. It is unclear how much empirical evidence supports GDDRH specifically, especially outside eusocial Hymenoptera.

Results

This question was addressed by conducting a literature review and posteriorly conducting meta-analyses on the data available using Hedges’s g. The literature review found 10 direct and 32 indirect studies with both having a strong publication bias towards Hymenoptera. Two meta-analyses were conducted and both found increased polyandry (direct tests; n = 8, g = 0.2283, p =  < 0.0001) and genetic diversity generated by other mechanisms (indirect tests; n = 10, g = 0.21, p =  < 0.0001) reduced parasite load. A subsequent moderator analysis revealed that there were no differences among Orders, indicating there may be applicability outside of Hymenoptera. However, due to publication bias and low sample size we must exercise caution with these results.

Conclusion

Despite the fact that the GDDRH was developed for Hymenoptera, it is frequently applied to other taxa. This study highlights the low amount of direct evidence supporting GDDRH, particularly outside of eusocial Hymenoptera. It calls for future research to address species that have high dispersal rates and contain mixes of solitary and communal nesting.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01925-3.

Living with relatives offsets the harm caused by pathogens in natural populations

Living with relatives can be highly beneficial, enhancing reproduction and survival. High relatedness can, however, increase susceptibility to pathogens. Here, we examine whether the benefits of living with relatives offset the harm caused by pathogens, and if this depends on whether species typically live with kin. Using comparative meta-analysis of plants, animals, and a bacterium (nspecies = 56), we show that high within-group relatedness increases mortality when pathogens are present. In contrast, mortality decreased with relatedness when pathogens were rare, particularly in species that live with kin. Furthermore, across groups variation in mortality was lower when relatedness was high, but abundances of pathogens were more variable. The effects of within-group relatedness were only evident when pathogens were experimentally manipulated, suggesting that the harm caused by pathogens is masked by the benefits of living with relatives in nature. These results highlight the importance of kin selection for understanding disease spread in natural populations.

Does genetic diversity protect host populations from parasites? A meta-analysis across natural and agricultural systems

If parasites transmit more readily between closely related hosts, then parasite burdens should decrease with increased genetic diversity of host populations. This important hypothesis is often accepted at face value—notorious epidemics of crop monocultures testify to the vulnerability of host populations that have been purged of diversity. Yet the relationship between genetic diversity and parasitism likely varies across contexts, differing between crop and noncrop hosts and between experimental and natural host populations. Here, we used a meta‐analytic approach to ask if host diversity confers protection against parasites over the range of contexts in which it has been tested. We synthesized the results of 102 studies, comprising 2004 effect sizes representing a diversity of approaches and host‐parasite systems. Our results validate a protective effect of genetic diversity, while revealing significant variation in its strength across biological and empirical contexts. In experimental host populations, genetic diversity reduces parasitism by ∼20% for noncrop hosts and by ∼50% for crop hosts. In contrast, observational studies of natural host populations show no consistent relationship between genetic diversity and parasitism, with both strong negative and positive correlations reported. This result supports the idea that, if parasites preferentially attack close relatives, the correlation of genetic diversity with parasitism could be positive or negative depending upon the potential for host populations to evolve in response to parasite selection. Taken together, these results reinforce genetic diversity as a priority for both conservation and agriculture and emphasize the challenges inherent to drawing comparisons between controlled experimental populations and dynamic natural populations.

Digging into the Genomic Past of Swiss Honey Bees by Whole-Genome Sequencing Museum Specimens

Historical specimens in museum collections provide opportunities to gain insights into the genomic past. For the Western honey bee, Apis mellifera L., this is particularly important because its populations are currently under threat worldwide and have experienced many changes in management and environment over the last century. Using Swiss Apis mellifera mellifera as a case study, our research provides important insights into the genetic diversity of native honey bees prior to the industrial-scale introductions and trade of non-native stocks during the 20th century—the onset of intensive commercial breeding and the decline of wild honey bees following the arrival of Varroa destructor. We sequenced whole-genomes of 22 honey bees from the Natural History Museum in Bern collected in Switzerland, including the oldest A. mellifera sample ever sequenced. We identify both, a historic and a recent migrant, natural or human-mediated, which corroborates with the population history of honey bees in Switzerland. Contrary to what we expected, we find no evidence for a significant genetic bottleneck in Swiss honey bees, and find that genetic diversity is not only maintained, but even slightly increased, most probably due to modern apicultural practices. Finally, we identify signals of selection between historic and modern honey bee populations associated with genes enriched in functions linked to xenobiotics, suggesting a possible selective pressure from the increasing use and diversity of chemicals used in agriculture and apiculture over the last century.

Host genetic diversity limits parasite success beyond agricultural systems: a meta-analysis

There is evidence that human activities are reducing the population genetic diversity of species worldwide. Given the prediction that parasites better exploit genetically homogeneous host populations, many species could be vulnerable to disease outbreaks. While agricultural studies have shown the devastating effects of infectious disease in crop monocultures, the widespread nature of this diversity-disease relationship remains unclear in natural systems. Here, we provide broad support that high population genetic diversity can protect against infectious disease by conducting a meta-analysis of 23 studies, with a total of 67 effect sizes. We found that parasite functional group (micro- or macroparasite) affects the presence of the effect and study setting (field or laboratory-based environment) influences the magnitude. Our study also suggests that host genetic diversity is overall a robust defence against infection regardless of host reproduction, parasite host range, parasite diversity, virulence and the method by which parasite success was recorded. Combined, these results highlight the importance of monitoring declines of host population genetic diversity as shifts in parasite distributions could have devastating effects on at-risk populations in nature.

Detection of replicative Kashmir Bee Virus and Black Queen Cell Virus in Asian hornet Vespa velutina (Lepelieter 1836) in Italy

Information concerning the pathogenic role of honey bee viruses in invasive species are still scarce. The aim of this investigation was to assess the presence of several honey bee viruses, such as Black Queen Cell Virus (BQCV), Kashmir Bee Virus (KBV), Slow Paralysis Virus (SPV), Sac Brood Virus (SBV), Israeli Acute Paralysis Virus (IAPV), Acute Bee Paralysis Virus (ABPV), Chronic Bee Paralysis Virus (CBPV), in Vespa velutina specimens collected in Italy during 2017. Results of this investigation indicate that among pathogens, replicative form of KBV and BQCV were detected, assessing the spillover effect of both these viruses from managed honey bees to hornets.

Tracking the Genetic Stability of a Honey Bee (Hymenoptera: Apidae) Breeding Program With Genetic Markers

A genetic stock identification (GSI) assay was developed in 2008 to distinguish Russian honey bees from other honey bee stocks that are commercially produced in the United States. Probability of assignment (POA) values have been collected and maintained since the stock release in 2008 to the Russian Honey Bee Breeders Association. These data were used to assess stability of the breeding program and the diversity levels of the contemporary breeding stock through comparison of POA values and genetic diversity parameters from the initial release to current values. POA values fluctuated throughout 2010-2016, but have recovered to statistically similar levels in 2016 (POA(2010) = 0.82, POA(2016) = 0.74; P = 0.33). Genetic diversity parameters (i.e., allelic richness and gene diversity) in 2016 also remained at similar levels when compared to those in 2010. Estimates of genetic structure revealed stability (FST(2009/2016) = 0.0058) with a small increase in the estimate of the inbreeding coefficient (FIS(2010) = 0.078, FIS(2016) = 0.149). The relationship among breeding lines, based on genetic distance measurement, was similar in 2008 and 2016 populations, but with increased homogeneity among lines (i.e., decreased genetic distance). This was expected based on the closed breeding system used for Russian honey bees. The successful application of the GSI assay in a commercial breeding program demonstrates the utility and stability of such technology to contribute to and monitor the genetic integrity of a breeding stock of an insect species.

Parasites and Their Social Hosts

The study of parasitism in socially living organisms shows that social group size correlates with the risk of infection, but group structure - and thus differences in contact networks - is generally more important. Also, genetic makeup or environmental conditions have effects. 'Social immunity' focuses on defence against parasites that are particular to social living. Recently, the role of socially transmitted microbiota for defence has become a focus, too. But whether and how parasites adapt to social organisms - beyond adaptation to solitary hosts - is poorly understood. Genomic and proteomic methods, as well as network analysis, will be tools that hold promise for many unsolved questions, but to expand our concepts in the first place is a much needed agenda.

Early gut colonizers shape parasite susceptibility and microbiota composition in honey bee workers

Microbial symbionts living within animal guts are largely composed of resident bacterial species, forming communities that often provide benefits to the host. Gut microbiomes of adult honey bees (Apis mellifera) include core residents such as the betaproteobacterium Snodgrassella alvi, alongside transient parasites such as the protozoan Lotmaria passim To test how these species affect microbiome composition and host physiology, we administered S alvi and/or L passim inocula to newly emerged worker bees from four genetic backgrounds (GH) and reared them in normal (within hives) or stressed (protein-deficient, asocial) conditions. Microbiota acquired by normal bees were abundant but quantitatively differed across treatments, indicating treatment-associated dysbiosis. Pretreatment with S. alvi made normal bees more susceptible to L. passim and altered developmental and detoxification gene expression. Stressed bees were more susceptible to L. passim and were depauperate in core microbiota, yet supplementation with S. alvi did not alter this susceptibility. Microbiomes were generally more variable by GH in stressed bees, which also showed opposing and comparatively reduced modulation of gene expression responses to treatments compared with normal bees. These data provide experimental support for a link between altered gut microbiota and increased parasite and pathogen prevalence, as observed from honey bee colony collapse disorder.