American foulbrood in a honeybee colony: spore-symptom relationship and feedbacks between disease and colony development
BMC Ecology volume 20, Article number: 15 (2020)
The most severe bacterial disease of honeybees is American foulbrood (AFB). The epidemiology of AFB is driven by the extreme spore resilience, the difficulty of bees to remove these spores, and the considerable incidence of undetected spore-producing colonies. The honeybee collective defence mechanisms and their feedback on colony development, which involves a division of labour at multiple levels of colony organization, are difficult to model. To better predict disease outbreaks we need to understand the feedback between colony development and disease progression within the colony. We therefore developed Bayesian models with data from forty AFB-diseased colonies monitored over an entire foraging season to (i) investigate the relationship between spore production and symptoms, (ii) disentangle the feedback loops between AFB epidemiology and natural colony development, and (iii) discuss whether larger insect societies promote or limit within-colony disease transmission.
Rather than identifying a fixed spore count threshold for clinical symptoms, we estimated the probabilities around the relationship between spore counts and symptoms, taking into account modulators such as brood amount/number of bees and time post infection. We identified a decrease over time in the bees-to-brood ratio related to disease development, which should ultimately induce colony collapse. Lastly, two contrasting theories predict that larger colonies could promote either higher (classical epidemiological SIR-model) or lower (increasing spatial nest segregation and more effective pathogen removal) disease prevalence.
AFB followed the predictions of the SIR-model, partly because disease prevalence and brood removal are decoupled, with worker bees acting more as disease vectors, infecting new brood, than as agents of social immunity, by removing infected brood. We therefore established a direct link between disease prevalence and social group size for a eusocial insect. We furthermore provide a probabilistic description of the relationship between AFB spore counts and symptoms, and how disease development and colony strength over a season modulate this relationship. These results help to better understand disease development within honeybee colonies, provide important estimates for further epidemiological modelling, and gained important insights into the optimal sampling strategy for practical beekeeping and honeybee research.
Honeybees are important pollinators in agricultural  and natural habitats . The demand for managed pollinators in agriculture has steadily increased during the past decades due to changing diets and an alarming decrease in natural pollinators in cultivated landscapes . At the same time, beekeepers worldwide are experiencing increased winter and seasonal colony losses [4,5,6,7]. Such losses stem from a combination of parasites and diseases, poor nutrition, inadequate beekeeping management practices and pesticide exposure; both individually and synergistically [8,9,10,11,12,13,14,15,16].
One of the major threats to colony health and beekeeping viability is American foulbrood (AFB); a contagious, lethal bacterial disease of honeybee brood that is widely distributed across the world. The disease causes great economic losses during outbreaks due to reduced productivity and material turnover [17,18,19,20]. American foulbrood is caused by Paenibacillus larvae, a spore-forming bacterium that produces extremely resilient spores which can remain viable for decades . Within a colony, P. larvae spores are spread by nurse bees performing in-hive tasks, such as cleaning, but especially through the feeding of young larvae with spore-contaminated food . Billions of spores are produced in the dying larvae . The dried larval remains (scales) are difficult to remove by workers and are a continuous source of infection for new brood. The lethality and epidemiology of AFB are driven by the resilience of the spores and the fact that the removal of diseased brood, a communal bee hygienic behaviour [23, 24], is not sufficient to remove this source of infection [20, 25]. The spores are distributed between colonies by swarming, robbing and in particular by beekeepers moving contaminated material between colonies [26, 27].
One major problem for the control of American foulbrood is that even though clinical symptoms are highly characteristic for the disease, they tend to appear late during the epidemic, when the colony’s hygienic behaviour to remove infected larvae before they produce spores, can no longer keep up with the epidemic. Estimations have been made that as much as 25% of spore-producing colonies remain undetected . Infections are therefore enzootic, since they remain in the population without external inputs , and occult, since they are present but largely visually undetected . Colonies can produce large amounts of infectious P. larvae spores with relatively few cases of symptomatic brood, thus escaping detection during routine beekeeper inspections while continuing to be a source of infection both within a beekeeping operation, between beekeepers (through sale of bees and equipment) and to feral and managed colonies within flight range through drifting and robbing .
One way to address the risk of pre-clinical infectious colonies for epidemic spread at multiple scales (local, regional, national) is to determine this risk directly from P. larvae spore levels in material sampled from the colony, thereby unambiguously identifying all infectious colonies rather than just those presenting detectable symptoms. It has been shown that adult bees provide most reliable samples for relating P. larvae spore levels to AFB symptoms, superior to either colony debris or honey samples [28, 31,32,33]. The spore load of individual bees is positively correlated to the likelihood of clinical symptoms [31, 34]. Since the spores are heavily concentrated in the brood frames and hive material, beekeepers can remove much of the colony-level spore burden by shaking the adult bees and queen into new, clean hives and frames with fresh wax foundation . Attempts have been made previously to predict clinical symptoms from the number of spores in a colony [36, 37]. However these attempts lacked two important elements for improving the reliability of such a calibration curve, namely an estimation of uncertainty and the usage of the recommended standardised grading of the severity of AFB symptoms.
We therefore included standardised AFB symptom grading scales  and multilevel Bayesian linear models to rectify these deficiencies, in order to identify more reliably the probability of AFB symptoms given a particular spore count. The results obtained are therefore directly applicable to practical beekeeping, as well as to research and epidemiological modelling. For example, the results can be used to parameterize other Bayesian models, by using the posterior probability predicted by these models as prior probability estimates for other Bayesian analyses, e.g. for predicting or analysing AFB transmission, infectivity, epidemiology, or symptoms in various actual or theoretical scenarios. Predicting AFB symptoms from spore counts follows a causational logic. However, from a practical perspective it would also be useful to explore the reverse relationship, i.e. to predict the colony spore levels from observed symptoms since the primary data obtained from colony inspections is the presence and severity of symptoms, which are then followed up with laboratory spore analyses of adult bee samples. This furthermore also serves as a quality control of the standardised AFB symptom grading system.
Our second objective was to disentangle the natural colony development and the colony-level disease development over a season. While both the development of P. larvae infection in larvae , and AFB disease epidemiology between colonies using colony infection data [18, 27] have been described, the factors shaping AFB virulence at the colony level are still largely unknown. We focused on the onset and the development of the disease, which are the most relevant disease stages from a practical as well as an epidemiological perspective. AFB can kill a colony within a single season, which in temperate regions of the northern hemisphere ranges from the beginning of April to end of September. Although one previous attempt at modelling AFB development suggests that the onset of symptoms is rather sudden , there has been no controlled study to evaluate the time course of AFB progression during a full season. In addition, honeybee colonies are complex entities (super-organisms) where most of the colony dynamics are driven by social interactions and decision making between its members, based on sensory input from within the colony and the environment, which makes it very difficult to predict the course of any disease . Particularly collective defence components, such as hygienic behaviour, and its feedback on colony development are difficult to model reliably, since it involves both brood and adult bees and affects role allocation and decision making at multiple levels in colony organization [24, 40].
Our last objective concerned whether larger honeybee colonies limit or accelerate the epidemiology of AFB. For example colonies that have naturally adapted to survive uncontrolled Varroa destructor infestation display unique and characteristic colony development traits, including smaller overall size, reduced drone brood production and lower brood-to-adult bee ratios, all of which are predicted to limit the reproductive potential for this mite . The classical deterministic epidemiological model, based on Susceptible, Infected and Resistant Hosts (the SIR-model ) would predict that larger colonies (more brood) should increase the spread of disease. However, living in social groups is clearly beneficial with regard to predation risk  and social hygiene, which may outweigh the higher infection risks . For eusocial insect societies in particular, larger colonies may be more able than smaller colonies to deploy counter measures to epidemic disease spread, such as removing diseased brood (hygienic behaviour) . Additionally, the increased spatial separation in larger colonies may delay the spread of infection . Considerable research has been conducted to understand how larger social groups cope with higher infection risk, and factors such as genetic diversity , group size [47, 48] and the structure of the social network  have been shown to be important for affecting the colony-level effects of infection and disease. Most of these studies aim at the mechanisms social animals use to minimize infection risk (e.g. the ability to detect and remove fungal diseases  or the effectiveness of cuticular antimicrobial defences in relation to colony size ). However, such an approach is contingent on the assumptions related to these biological mechanisms. Here we approach the relationships between disease severity and colony size directly, and therefore independent of mechanistic assumptions. For example the variation in the number of susceptible hosts (i.e. the amount of larval brood, for AFB) over a season depends on many intrinsic and external factors [52, 53]. Since American foulbrood kills the brood, a shift in the brood-to-worker ratio can be expected, but it is unclear whether larger honeybee colonies are better than smaller colonies at handling such a shift. In this study we will therefore attempt to describe how AFB affects colony strength parameters (amount of brood and number of adult bees) and how these parameters feed back to the epidemiology in the colony, in order to understand the relevance of colony size itself, independent of colony-level mechanisms, on disease spread.
The results describe the importance of a particular predictor for predicting a response variable, the direction of the predicted effect on the response variable, and the extent to which the primary relationship is modulated by secondary and tertiary predictors. Within each figure, the central subfigure shows the relationship between the primary predictor and the response variable if the influence of the other two secondary predictors (top and right) is neutral. Because these secondary predictors mostly co-vary (i.e. Fig. 1: Brood 80 at Time 20 is more likely than Brood 200 at Time 20) the three subfigures on the leading diagonal are the most likely outcomes, while the remaining six subfigures are primarily for explaining and understanding the interactions. A summary of the original data can be found in Additional file 1: (Figure S1 and Table S1). Throughout the results variables have capital letters.
Predicting clinical symptoms from spore counts
All response variables were relevant for predicting Symptoms (Table 2: P[effect > 0] = 100% for all predictors). Spores were a strong predictor of the Symptoms in all four models, both directly and through its interaction with the other variables (Table 1: M1–M4), and because subtracting the effect of secondary predictors from the Spores predictor left a high probability of an effect size larger than zero (Table 2: Spores-Brood P[effect > 0] = 91.3%; Spores-Time P[effect > 0] = 94.3%). Brood and Time were equally predictive for Symptoms (Table 2: Brood-Time P[effect > 0] = 55.6%), though much less so than Spores. The Symptoms increased first slightly and then strongly with increasing Spores, regardless of the time and brood (Fig. 1: Time 50/Brood 130). Symptoms increased over Time, both in absolute terms and in relation to a given spore count level (see also Additional file 2: Figure S2), and Symptoms also increased with increasing Brood (see also Additional file 2: Figure S3).
Brood was a stronger modulator of the Spores-Symptoms relationship than Time since the Spores × Brood interaction received more weight than the Spores × Time interaction (Table 1: M2 includes Spores × Brood and not Spores × Time). This is also illustrated by a stronger change of the Spores-Symptoms relationship along Brood than along Time (Fig. 1).
Lastly we calculated the probability of encountering Symptoms if no Spores are detected in an adult bee sample (Spores = 0), which resulted in a probability around 0.22.
Predicting spore counts from clinical Symptoms
All response variables were relevant for predicting Spores (Table 2). Symptoms were the strongest predictors of Spores, both directly and through its many significant interactions with other predictors; three out of the four selected models (Table 1), and since subtracting the effects of secondary predictors from symptoms left a high probability of an effect size larger than zero (Table 2). Bees seemed more important than Time for predicting Spores (Table 2).
Regardless of Time and Bees, Spores increased with increasing clinical symptoms (Fig. 2: Time 50/Bees 9; see also Additional file 2: Figure S4 for the full range of spore counts). In general though, Spores decreased over Time, as is illustrated by the decrease over Time for any given level of Symptoms (see also Additional file 2: Figure S5). Spores also decreased with increasing Bees irrespective of Time or Symptoms (Additional file 2: Figure S6). However, early in the season Spores increased with increasing number of bees.
Time was a stronger modulator of the Symptoms-Spores relationship than number of bees, since the model including only the Symptoms × Time interaction received 50% of the Akaike weight (Table 1). This is illustrated by a greater change in the slope of the Symptoms-Spores relationship in relation to Time than in relation to Bees (Fig. 2, Additional file 2: Figure S4).
We furthermore estimated the number of expected spores in an adult bee sample if symptoms are at level zero or level one, which resulted in around 158 and 228 spores, respectively. Lastly, we calculated the posterior of the difference for extreme and likely values of Spores in order to investigate a dilution effect of sampling 100 bees while colony size differed (see Additional file 2).
Predicting the effect of AFB disease on colony strength
All response variables were relevant for predicting the first variable of colony size, the number of adult bees (Table 2). Brood was the strongest predictor of Bees, as revealed by the importance of its interactions with the other variables in all four selected models (Table 1). This importance is further illustrated by the high probability of an effect larger than zero after subtracting the effect of Spores or Time as a co-predictor. The time post-infection seemed to be a much more important co-predictor of the number of bees than the number of spores (Table 2).
Regardless of Time or Brood the number of bees increased with increasing spore numbers (Fig. 3: Time 50/Brood 130). Similarly, regardless of Spores and Brood, Bees increased with Time (Fig. 3: see at zero spore count and in Additional file 2: Figure S7 at Spores 850) and Bees increase with increasing Brood if the other predictors are held at their mean/median values (Fig. 3: zero Time point; Additional file 2: Figure S8: Spores 850).
Time was a stronger modulator of the Bees-Spores relationship in the colonies than Brood as the model with the Spores × Time interaction was ranked higher than the model including the Spores × Brood interaction (Table 1). This is illustrated by the strong change over Time in the slope, from a positive to a negative relationship, while the changes with respect to Brood are less prominent (Fig. 3).
Regarding the second measure of colony strength, i.e. the amount of brood, we saw that all response variables were relevant for predicting brood amount (Table 2). The number of bees was by far the strongest predictor of brood amount in the first model, which included only the Symptoms × Bees interaction, and received 99% of all the predictors’ weight (Table 1). Therefore, here we only used M13 for the predictions instead of a weighted combination of all 4 models.
The importance of Bees for predicting the amount of brood is further illustrated by the high portability of an effect larger than zero after subtracting the co-predictor symptoms or time. Symptoms seemed more important than Time for predicting Brood (Table 2).
Regardless of Time and Bees, the Brood increased with increasing Symptoms (Fig. 3: Time 50/Bees 9). Similarly, Brood decreased with Time irrespective of the Symptoms and Bees (Fig. 4 and Additional file 2: Figure S9, any symptom score). Brood also increases with Bees (Fig. 4), although at very high symptoms, more bees did not lead to more brood anymore.
Bees was an overwhelmingly stronger modulator of the Brood-Symptoms relationship relative to Time since Symptoms × Bees was the only relevant interaction (Table 1). This is illustrated by the strong change in the slope of the brood-symptoms relationship in relation to the number of Bees, relative to the slope in relation to Time (Fig. 4).
The primary objective of our study was to determine the probability distribution for encountering colony-level AFB symptoms from P. larvae spore counts, as a superior approach to identifying a threshold. By taking a probabilistic approach, other factors affecting the relationship between spore counts and symptoms, such as colony size, brood availability and stage of the infection process were explicitly accounted for, something which is not possible with simple threshold values.
While we could describe the relationship between spore levels and symptom independent of the effects of time of the season and the amount of brood (by keeping these at their mean value), questions remain about the generality of the results. Other possible factors that may change the spore-symptom relationship are variability in bacterial virulence (an innate genetic property of different P. larvae strains ), colony genetics , resistance against P. larvae , and hygienic behaviour . All colonies in the study were placed in the same isolated apiary and were experimentally infected with sufficient Paenibacillus larvae spores to precipitate AFB disease. All colonies were therefore under the same infection pressure. This means that any additional contagion contributed by bees drifting between the colonies will be miniscule compared to the contagion developed within each colony, and will not have affected the results. Although proximity and connectedness (beekeeping and geolocation) are important determinants of the disease pressure in a colony  this is not applicable to the current study which was conducted in complete isolation from other beekeeping operations, as required by sanitary regulations.
From an applied perspective, it may be more interesting to get an idea of the spore counts given a symptom score. We found even for colony symptom scores of 0 (i.e. no disease), there is a 100% probability that the spore count is larger than zero, averaging around 158 Spores per bee (Table 2: Symptoms = 0). Symptomatic colonies with an AFB score of 1 corresponded to around 228 Spores per bee (Table 2: Symptoms = 1). This is considerably lower than the previously estimated threshold of 3000 spores per bee for AFB symptoms , but higher than the estimate of Lindström .
Scoring clinical symptoms may be biased given that symptoms may not be visible at an early stage of infection  and that symptoms in larger colonies with many brood frames to inspect may be underestimated by human eye [57, 58]. Regardless of the season or the spore levels, clinical symptoms increased with brood size in our study. It remains unclear if, and to what extent, the precision of the prediction of AFB symptoms suffers from increasing colony size. However, the AFB disease scoring is very sensitive at the lower range and any number of diseased cells above 100 corresponds the highest AFB symptom score, which would probably guard against such overestimation. Another result of our probabilistic approach is evidence that the AFB scoring method seems not to produce false positives, since for a spore count of 0, the probabilistic estimate for AFB symptoms does not reach 1, which is the minimum AFB-positive score (Table 2: Spores = 0: 0.22 ± 0.05). This confirms the accuracy of the scoring scale and previous findings that high spore levels will be detected in symptomatic colonies .
Our second objective was to clarify the interactive relationship between colony development and disease development. Brood amount and the time post infection were similarly important secondary predictive factors affecting the relationship between spore counts and symptoms. Both are important for the epidemiology of the disease: the amount of brood representing new, uninfected hosts and time being an obvious important factor in any epidemiological disease progression. Symptoms increased slowly with increasing brood size and over time. Including data from symptomatic colonies only would have shown a faster increase. However, we were also interested in predicting from the spore counts the probability of symptoms developing in colonies that passed visual AFB inspection. The presented models therefore investigated the disease development in infected colonies, rather than just in symptomatic colonies.
AFB kills progressively more brood as the epidemic intensifies, and the consequent shortage of new adult bees leads to progressive dwindling and eventual demise of the colony. This study describes this process in greater detail, including the behavioural responses of the colony. During normal colony development, and as long as there is forage available, more brood leads to more adult bees, which then leads to more brood and so on. In these experiments however, the amount of brood decreased throughout the entire season, rather than being restricted to the autumn (which is when the colonies normally transition to broodless winter colonies ). Admittedly we did not had uninfected control colonies the development of the brood could be compared to, which would have enabled us to estimate the brood loss over the season entirely attributed to AFB. Nevertheless, since no swarming occurred in these colonies during the season, we believe this decrease could therefore be (mostly) attributed to AFB. However, early in the season the colonies responded to AFB symptoms by increasing their brood production, shown by a greater increase in smaller colonies than in larger colonies (Fig. 4: steeper slope at Bees 6 than Bees 12). Nevertheless, later in the season this compensation attempt failed and the number of bees decreased with increasing bacterial spore levels (Fig. 3). This decrease may have been partly overestimated, due to the dilution effect in larger colonies, where the spore count per adult bee sample is shared between larger numbers of adult bees (see Additional file 2 for more explanations and implications for sampling strategy). The compensation failure can be seen more clearly by slower increase in brood amount with increasing symptoms severity (Fig. 4). Larger amounts of brood fail to hatch into adult workers, which is most clearly illustrated by the brood-to-bee ratio. At later time points, the same brood amount predicts a larger number of bees (Additional file 2: Figure S8) and the same number of bees predicts less brood (Additional file 2: Figure S10).
The observed increase in the amount of brood in response to a stressor has previously been described for Varroa destructor infested colonies , where larger colonies in the autumn were more likely to die the following year. A bee colony is an adaptable unit that uses brood rearing as one of the mechanisms to respond to external and internal stimuli, through both positive and negative feedback loops . A (perceived) deficit in either adults or healthy brood can be one of such internal stimuli, resulting in an elevated brood rearing effort in smaller colonies, and thus an elevated brood-to-adult ratio. We observed such a compensation attempt, at the most relevant time point of the life cycle for the colony development (spring), and we could show that this effect subsequently carried over to the adult (worker) bee stage. Although external factors such as foraging availability and quality, determined by the surrounding landscape, are also highly influential for brood rearing and overall colony strength [61, 62], all colonies in this study were located in the same apiary and therefore exposed to the same landscape and foraging conditions. Since the number of bees declined due to the disease we can also expect a feedback of lower food intake to further increase disease prevalence that creates a cycle of stress .
Our final objective was to investigate the epidemiological aspect of AFB. In these experiments we verify the expectation of the SIR model of epidemiology, but also identified peculiarities for the spread of AFB within a colony. The epidemic potential (given by the reproductive number R0) of a disease increases with increasing transmission rates and number of susceptible hosts . The particularities of disease epidemiology in social animals was reviewed recently . One theoretical prediction is that, contrary to the SIR model, disease prevalence may decrease with increasing group size, if the behavioural responses limiting disease prevalence or transmission become more effective with increasing group size, such as grooming behaviour in termites . The hygienic behaviour of honeybees that involves detecting and removing infected and asymptomatic brood would also affect the SIR model , since it systematically reduces the amount of infectious material in the colony, while the disappearance of brood would of itself act as a stimulus for rearing new (uninfected) brood, both of which are important parameters for the SIR model. In these experiments, clinical symptoms always increased with brood size (Additional file 2: Figure S3), thus favouring the traditional SIR model of epidemiology as explanation. Although small colony size due to poor nutrition may amplify disease susceptibility , we do expect that the disease will be more severe and increase more over the season in larger colonies. This would mean larger colonies are not more resilient against AFB and will decrease in size stronger than smaller ones. Smaller colonies have lower expected overwintering survival  which could lead to an additional colony loss in the next spring.
Contrary to the positive effects of hygienic grooming behaviour on disease in termites, our study found a positive relationship between colony size and disease symptoms, implying the hygienic behaviour of brood removal is perhaps ineffective at breaking this relationship. In fact, the adults carry bacterial spores and serve as vectors infecting new brood. The broader implication here is whether group size facilitates or hinders disease transmission in social animals will depend on what life stage is affected by the disease and how this effect translates to the other life stages or affect the task allocation  within eusocial insects.
We provide a novel, and potentially more reliable method for quantifying the relationship between P. larvae spore counts and AFB symptoms. Furthermore, we showed how AFB-caused brood mortality led to progressively fewer adult worker bees, eventually tipping the colony into a deadly negative spiral from which it could not escape. We identified that AFB disease epidemiology in honeybees follows the more traditional SIR model of epidemiology. We found little evidence of any beneficial effects of the hygienic behaviour of brood removal on containing the epidemic, especially since adult bees simultaneously also act as vectors of the disease. We extend the discussion further to larger eusocial societies exhibiting stronger social immunity by showing this seems not to apply for American foulbrood in honeybees as workers are removing diseased brood but also vectoring the disease. The study therefore emphasizes to consider how certain defence strategies will manifest themselves in other life stages of the eusocial society and shows the direct feedbacks between the epidemic over a season and the colony size.
On March 24th 2014, forty honeybee colonies located in an isolated apiary with a history of AFB in Beltsville, MD, USA were selected for the experiment (colonies owned by USDA-ARS Bee Research Lab). The experiment was originally designed to test the efficacy of a commercial honeybee specific lactic acid bacteria preparation against AFB relative to two negative controls (a placebo preparation and a no treatment control) and a positive control (the antibiotic Tylosin) (Additional file 3, ). The colonies were arranged in four rows of ten colonies each, with 1.5 m distance between individual colonies in each row and 1.5 m distance between rows and all entrances facing the same direction . All colonies were experimentally inoculated with the same dose of Paenibacillus larvae spores at to precipitate AFB epidemics with the four treatment groups distributed randomly among the 40 colonies. Although no active measures were taken to prevent drifting between the 40 colonies, the uniform inoculation of all colonies in the apiary and the spatial randomization of the treatment groups means that the effect of drifting bees on the AFB development in colonies is both minimal, and randomly distributed between the treatment groups. The randomization also means that any potential variability in the data caused by the treatment groups can be accounted for statistically in our modelling (Additional file 2). The colonies were assessed and adult bees sampled on April 23rd, and then 21, 37, 51, 79, and 105 days after the first assessment. On each sampling occasion approximately 200 adult bees were collected from the brood chamber per colony and the samples were stored at − 20 °C until spore estimation in the lab.
Colony assessments, AFB scoring, and spore counting
On each of the six sampling occasions, total colony size, the amount of brood and the severity of AFB symptoms were evaluated using standard protocols. Colony size (hereafter: Bees; with capital first letter) was estimated by a cumulative score of the proportion of each frame side that was occupied by adult bees . The amount of brood in the colony (hereafter: Brood) was estimated by a cumulative score of the number of 5 × 5 cm squares on each frame that were occupied by brood . The colony-level severity of AFB (hereafter: Symptoms) was estimated by a cumulative score of the visual inspection of each brood frame for signs of the disease [38, 71]. Each frame was rated using the recommended scale of 0 (no visible signs), 1 (fewer than 10 diseased cells), 2 (11–100 diseased cells), and 3 (more than 100 diseased cells). Samples of diseased cells were tested in the laboratory to confirm the diagnoses. The spore levels (hereafter: Spores) were determined from samples of 100 adult worker bees, as described previously [32, 38]. The raw data consisted of P. larvae colony forming units (CFU) and the data are presented as CFU per bee (see also Additional file 3). The effect of colony or disease development during the season is represented in the models by Time.
The data was obtained from 40 colonies, sampled once a month for 6 months during a single bee season (see Additional file 1 for original data overview). A Bayesian approach was used for the statistical modelling and analyses [72,73,74]. All variables were continuous counts and the analysis is similar to a multiple regression. The models were constructed in two steps (see Additional file 3: Model building and validation). First, two similar models with different random structures were compared. In step two, we compared eight models for each of the four response variables. Each model included the three main effects and all combinations of their interactions. The four most important models from step two were then used in the analysis by weighting the predictions in order to include modulations of one predictor by the other two predictors. To understand the effect of each predictor we calculated the posterior of the response variable along the full observed range of one explanatory variable while keeping the remaining two explanatory variables constant, conventionally at their mean/median value. For a better understanding of the complex interaction between the three continuous predictors two additional values were selected for each of the 4 models to investigate interactive effects. For example approximately the mean (132.6), 1st (78) and 3rd (191) quantile were used for brood (hereafter: Brood 130, Brood 80, Brood 200, respectively; see also Additional file 3: Model building and validation). Furthermore we selected specific values of a predictor and summarized the posterior of the response variable in order to answer specific questions (e.g.: Table 2: Symptoms at Spores level zero: Spores = 0). We also calculated three scenarios (Extreme dilution, Likely dilution 1, and Likely dilution 2) by subtracting the posterior of one set of values from the posterior of another set. The resulting posteriors can be seen as pairwise comparisons among these sets. In order to further investigate the predictiveness of each main effect compared to the other we subtracted the posteriors from each other . The posterior of each main effect was weighted the same way as for the predictions and the smaller was subtracted from the larger in order to calculate the posterior of the difference. In all models we used minimal informative priors and the posterior was generated as a Monte Carlo sample (2000 iterations; Hamilton Monte Carlo; 1000 warm up, 1000 sampling the chains) using STAN  handled from R  using function from McElreath .
Gallai N, Salles JM, Settele J, Vaissière BE. Economic valuation of the vulnerability of world agriculture confronted with pollinator decline. Ecol Econ. 2009;68:810–21.
Hung K-LJ, Kingston JM, Albrecht M, Holway DA, Kohn JR. The worldwide importance of honey bees as pollinators in natural habitats. Proc R Soc B Biol Sci. 2018. https://doi.org/10.1098/rspb.2017.2140.
Aizen MA, Harder LD. The global stock of domesticated honey bees is growing slower than agricultural demand for pollination. Curr Biol. 2009;19:915–8.
Steinhauer N, Kulhanek K, Antúnez K, Human H, Chantawannakul P, Chauzat MP, et al. Drivers of colony losses. Curr Opin Insect Sci. 2018;26:142–8.
McMenamin AJ, Genersch E. Honey bee colony losses and associated viruses. Curr Opin Insect Sci. 2015;8:121–9.
vanEngelsdorp D, Caron D, Hayes J, Underwood R, Henson M, Rennich K, et al. A national survey of managed honey bee 2010–2011 winter colony losses in the USA: results from the Bee Informed Partnership. J Apic Res. 2012;51:115–24. https://doi.org/10.3896/IBRA.184.108.40.206.
Gray A, Brodschneider R, Adjlane N, Ballis A, Brusbardis V, Charrière J-D, et al. Loss rates of honey bee colonies during winter 2017/2018 in 36 countries participating in the COLOSS survey, including effects of forage sources. J Apic Res. 2019;58:479–85. https://doi.org/10.1080/00218839.2019.1615661.
Tosi S, Nieh JC, Sgolastra F, Cabbri R, Medrzycki P. Neonicotinoid pesticides and nutritional stress synergistically reduce survival in honey bees. Proc R Soc B Biol Sci. 2017;284:20171711. https://doi.org/10.1098/rspb.2017.1711.
Doublet V, Labarussias M, de Miranda JR, Moritz RFA, Paxton RJ. Bees under stress: sublethal doses of a neonicotinoid pesticide and pathogens interact to elevate honey bee mortality across the life cycle. Environ Microbiol. 2015;17:969–83.
López JH, Krainer S, Engert A, Schuehly W, Riessberger-Gallé U, Crailsheim K. Sublethal pesticide doses negatively affect survival and the cellular responses in American foulbrood-infected honeybee larvae. Sci Rep. 2017. https://doi.org/10.1038/srep40853.
Klein S, Cabirol A, Devaud J-M, Barron AB, Lihoreau M. Why bees are so vulnerable to environmental stressors. Trends Ecol Evol. 2016;32:268–78.
Becher MA, Osborne JL, Thorbek P, Kennedy PJ, Grimm V. Towards a systems approach for understanding honeybee decline: a stocktaking and synthesis of existing models. J Appl Ecol. 2013;50:868–80.
Vanengelsdorp D, Meixner MD. A historical review of managed honey bee populations in Europe and the United States and the factors that may affect them. J Invertebr Pathol. 2010;103:80–95.
van Dooremalen C, Cornelissen B, Poleij-Hok-Ahin C, Blacquière T. Single and interactive effects of Varroa destructor, Nosema spp., and imidacloprid on honey bee colonies (Apis mellifera). Ecosphere. 2018. https://doi.org/10.1002/ecs2.2378.
Straub L, Williams GR, Vidondo B, Khongphinitbunjong K, Retschnig G, Schneeberger A, et al. Neonicotinoids and ectoparasitic mites synergistically impact honeybees. Sci Rep. 2019;9:8159. https://doi.org/10.1038/s41598-019-44207-1.
Genersch E, von der Ohe W, Kaatz H, Schroeder A, Otten C, Büchler R, et al. The German bee monitoring project: a long term study to understand periodically high winter losses of honey bee colonies. Apidologie. 2010;41:332–52. https://doi.org/10.1051/apido/2010014.
Chauzat MP, Jacques A, EPILOBEE consortium, Laurent M, Bougeard S, Hendrikx P, et al. Risk indicators affecting honeybee colony survival in Europe: one year of surveillance. Apidologie. 2016;47:348–78.
Mill AC, Rushton SP, Shirley MDF, Smith GC, Mason P, Brown MA, et al. Clustering, persistence and control of a pollinator brood disease: epidemiology of American foulbrood. Environ Microbiol. 2014;16:3753–63.
Ebeling J, Knispel H, Hertlein G, Fünfhaus A, Genersch E. Biology of Paenibacillus larvae, a deadly pathogen of honey bee larvae. Appl Microbiol Biotechnol. 2016;100:7387–95. https://doi.org/10.1007/s00253-016-7716-0.
Genersch E. American Foulbrood in honeybees and its causative agent, Paenibacillus larvae. J Invertebr Pathol. 2010;103(SUPPL. 1):S10–9. https://doi.org/10.1016/j.jip.2009.06.015.
Hasemann L. How long can spores of American foulbrood live? Am Bee J. 1961;101:298–9.
Lindström A, Korpela S, Fries I. The distribution of Paenibacillus larvae spores in adult bees and honey and larval mortality, following the addition of American foulbrood diseased brood or spore-contaminated honey in honey bee (Apis mellifera) colonies. J Invertebr Pathol. 2008;99:82–6.
Spivak M, Reuter GS. Resistance to American foulbrood disease by honey bee colonies Apis mellifera bred for hygienic behavior. Apidologie. 2001;32:555–65. https://doi.org/10.1051/apido:2001103.
Wilson-Rich N, Spivak M, Fefferman NH, Starks PT. Genetic, individual, and group facilitation of disease resistance in insect societies. Annu Rev Entomol. 2009;54:405–23. https://doi.org/10.1146/annurev.ento.53.103106.093301.
Bailey L, Ball B. Honey bee pathology. 2nd ed. Amsterdam: Elsevier; 1991.
Fries I, Camazine S. Implications of horizontal and vertical pathogen transmission for honey bee epidemiology. Apidologie. 2001;32:199–214. https://doi.org/10.1051/apido:2001122.
Datta S, Bull JC, Budge GE, Keeling MJ. Modelling the spread of American foulbrood in honeybees. J R Soc Interface. 2013;10:20130650–20130650. https://doi.org/10.1098/rsif.2013.0650.
Gillard M, Charriere JD, Belloy L. Distribution of Paenibacillus larvae spores inside honey bee colonies and its relevance for diagnosis. J Invertebr Pathol. 2008;99:92–5.
Erban T, Ledvinka O, Kamler M, Nesvorna M, Hortova B, Tyl J, et al. Honeybee (Apis mellifera)-associated bacterial community affected by American foulbrood: detection of Paenibacillus larvae via microbiome analysis/631/158/855/631/326/2565/855/38/23/38/22/38/47 article. Sci Rep. 2017;7:1–10. https://doi.org/10.1038/s41598-017-05076-8.
Fries I, Lindström A, Korpela S. Vertical transmission of American foulbrood (Paenibacillus larvae) in honey bees (Apis mellifera). Vet Microbiol. 2006;114:269–74.
Lindström A, Fries I. Sampling of adult bees for detection of American foulbrood (Paenibacillus larvae subsp larvae) spores in honey bee (Apis mellifera) colonies. J Apic Res. 2005;44:82–6. https://doi.org/10.1080/00218839.2005.11101154.
Forsgren E, Laugen AT. Prognostic value of using bee and hive debris samples for the detection of American foulbrood disease in honey bee colonies. Apidologie. 2014;45:10–20.
Nordström S, Forsgren E, Fries I. Comparative diagnosis of American foulbrood using samples of adult honey bees and honey. Apic Sci. 2002;46:5–13.
Goodwin RM, Perry JH, Haine HM. A study on the presence of Bacillus larvae spores carried by adult honey bees to identify colonies with clinical symptoms of American foulbrood disease. J Apic Res. 1996;35:118–20.
Del Hoyo ML, Basualdo M, Lorenzo A, Palacio MA, Rodriguez EM, Bedascarrasbure E. Effect of shaking honey bee colonies affected by American foulbrood on Paenibacillus larvae larvae spore loads. J Apic Res. 2001;40:65–9.
Lindström A. Distribution of Paenibacillus larvae spores among adult honey bees (Apis mellifera) and the relationship with clinical symptoms of American foulbrood. Microb Ecol. 2008;56:253–9. https://doi.org/10.1007/s00248-007-9342-y.
Gende L, Satta A, Ligios V, Ruiu L, Buffa F, Fernandez N, et al. Searching for an American foulbrood early detection threshold by the determination of paenibacillus larvae spore load in worker honey bees. Bull Insectol. 2011;64:229–33.
de Graaf DC, Alippi AM, Antúnez K, Aronstein KA, Budge G, De Koker D, et al. Standard methods for American foulbrood research. J Apic Res. 2013;52:1–28. https://doi.org/10.3896/IBRA.220.127.116.11.
Jatulan EO, Rabajante JF, Banaay CGB, Fajardo AC, Jose EC. A mathematical model of intra-colony spread of American foulbrood in European honeybees (Apis mellifera L.). PLoS ONE. 2015;10:1–13.
Pie MR, Rosengaus RB, Traniello JFA. Nest architecture, activity pattern, worker density and the dynamics of disease transmission in social insects. J Theor Biol. 2004;226:45–51.
Locke B, Fries I. Characteristics of honey bee colonies (Apis mellifera) in Sweden surviving Varroa destructor infestation. Apidologie. 2011;42:533–42. https://doi.org/10.1007/s13592-011-0029-5.
Anderson RM, May RM. Coevolution of hosts and parasites. Parasitology. 1982;85:411. https://doi.org/10.1017/S0031182000055360.
Stephan JG, Low M, Stenberg JA, Björkman C. Predator hunting mode and host plant quality shape attack-abatement patterns of predation risk in an insect herbivore. Ecosphere. 2016;7:e01541. https://doi.org/10.1002/ecs2.1541.
Ezenwa VO, Worsley-Tonks KEL. Social living simultaneously increases infection risk and decreases the cost of infection. Proc R Soc B Biol Sci. 2018;285:20182142. https://doi.org/10.1098/rspb.2018.2142.
Schmid-Hempel P. Parasites and their social hosts. Trends Parasitol. 2017;33:453–62.
Tarpy DR. Genetic diversity within honeybee colonies prevents severe infections and promotes colony growth. Proc R Soc B Biol Sci. 2003;270:99–103. https://doi.org/10.1098/rspb.2002.2199.
Simone-Finstrom M, Walz M, Tarpy DR. Genetic diversity confers colony-level benefits due to individual immunity. Biol Lett. 2016;12:20151007. https://doi.org/10.1098/rsbl.2015.1007.
Cremer S, Armitage SAO, Schmid-Hempel P. Social Immunity. Curr Biol. 2007;17:R693–702. https://doi.org/10.1016/j.cub.2007.06.008.
Nunn CL, Jordan F, McCabe CM, Verdolin JL, Fewell JH. Infectious disease and group size: more than just a numbers game. Philos Trans R Soc B Biol Sci. 2015;370:20140111–20140111. https://doi.org/10.1098/rstb.2014.0111.
Tranter C, Lefevre L, Evison SEF, Hughes WOH. Threat detection: contextual recognition and response to parasites by ants. Behav Ecol. 2015;26:396–405.
Hoggard SJ, Wilson PD, Beattie AJ, Stow AJ. The effectiveness of antimicrobial defenses declines with increasing group size and genetic similarity. Ann Entomol Soc Am. 2013;106:53–8.
Donkersley P, Rhodes G, Pickup RW, Jones KC, Power EF, Wright GA, et al. Nutritional composition of honey bee food stores vary with floral composition. Oecologia. 2017;185:749–61. https://doi.org/10.1007/s00442-017-3968-3.
Walton A, Dolezal AG, Bakken MA, Toth AL. Hungry for the queen: honey bee nutritional environment affects worker pheromone response in a life-stage dependent manner. Funct Ecol. 2018. https://doi.org/10.1111/1365-2435.13222.
Genersch E, Ashiralieva A, Fries I. Strain- and genotype-specific differences in virulence of Paenibacillus larvae subsp. larvae, a bacterial pathogen causing American foulbrood disease in honeybees. Appl Environ Microbiol. 2005;71:7551–5. https://doi.org/10.1128/aem.71.11.7551-7555.2005.
Wedenig M, Riessberger-Gallé U, Crailsheim K. A substance in honey bee larvae inhibits the growth of Paenibacillus larvae larvae. Apidologie. 2003;34:43–51. https://doi.org/10.1051/apido:2002043.
Brødsgaard CJ, Hansen H, Ritter W. Progress of Paenibacillus larvae larvae infection in individually inoculated honey bee larvae reared singly in vitro, in micro colonies, or in full-size colonies. J Apic Res. 2000;39:19–27.
Dakin SC, Tibber MS, Greenwood JA, Kingdom FAA, Morgan MJ. A common visual metric for approximate number and density. Proc Natl Acad Sci. 2011;108:19552–7. https://doi.org/10.1073/pnas.1113195108.
Jiménez J. Effect of sample size, plot size, and counting time on estimates of avian diversity and abundance in a Chilean rainforest. J F Ornithol. 2000;71:66–87. https://doi.org/10.1648/0273-8570-71.1.66.
Mattila HR, Otis GW. Dwindling pollen resources trigger the transition to broodless populations of long-lived honeybees each autumn. Ecol Entomol. 2007;32:496–505. https://doi.org/10.1111/j.1365-2311.2007.00904.x.
Winston ML. The Biology of the honey bee. Cambridge: Harvard University Press; 1987.
vanEngelsdorp D, Hayes J, Underwood RM, Pettis JS. A survey of honey bee colony losses in the United States, fall 2008 to spring 2009. J Apic Res. 2010;49:7–14. https://doi.org/10.3896/IBRA.1.49.1.03.
Brodschneider R, Crailsheim K. Nutrition and health in honey bees. Apidologie. 2010;41:278–94. https://doi.org/10.1051/apido/2010012.
Dolezal AG, Toth AL. Feedbacks between nutrition and disease in honey bee health. Curr Opin Insect Sci. 2018;26:114–9. https://doi.org/10.1016/j.cois.2018.02.006.
Rosengaus RB, Maxmen AB, Coates LE, Traniello JFA. Disease resistance: a benefit of sociality in the dampwood termite Zootermopsis angusticollis (Isoptera: Termopsidae). Behav Ecol Sociobiol. 1998;44:125–34. https://doi.org/10.1007/s002650050523.
Döke MA, McGrady CM, Otieno M, Grozinger CM, Frazier M. Colony size, rather than geographic origin of stocks, predicts overwintering success in honey bees (Hymenoptera: Apidae) in the Northeastern United States. J Econ Entomol. 2019;112:525–33. https://doi.org/10.1093/jee/toy377.
Natsopoulou ME, McMahon DP, Paxton RJ. Parasites modulate within-colony activity and accelerate the temporal polyethism schedule of a social insect, the honey bee. Behav Ecol Sociobiol. 2016;70:1019–31. https://doi.org/10.1007/s00265-015-2019-5.
Stephan JG, Lamei S, Pettis JS, Riesbeck K, de Miranda JR, Forsgren E. Honeybee-specific lactic acid bacterium supplements have no effect on American foulbrood-infected honeybee colonies. Appl Environ Microbiol. 2019;85:1–12. https://doi.org/10.1128/AEM.00606-19.
Dynes TL, Berry JA, Delaplane KS, Brosi BJ, De Roode JC. Reduced density and visually complex apiaries reduce parasite load and promote honey production and overwintering survival in honey bees. PLoS ONE. 2019;14:1–16.
Delaplane KS, van der Steen J, Guzman-Novoa E. Standard methods for estimating strength parameters of Apis mellifera colonies. J Apic Res. 2013. https://doi.org/10.3896/IBRA/1.52.1.03.
Pettis JS, Rose R, Chaimanee V. Chemical and cultural control of Tropilaelaps mercedesae mites in honeybee (Apis mellifera) colonies in Northern Thailand. PLoS ONE. 2017;12:e0188063.
Pettis JS, Feldlaufer MF. Efficacy of lincomycin and tylosin in controlling American foulbrood in honey bee colonies. J Apic Res. 2005;44:106–8.
Kruschke JK, Aguinis H, Joo H. The time has come: Bayesian methods for data analysis in the organizational sciences. Organ Res Methods. 2012;15:722–52.
McElreath R. Statistical rethinking: a Bayesian course with examples in R and Stan. J Educ Behav Stat. 2015. https://doi.org/10.3102/1076998616659752.
Kruschke JK. Doing Bayesian data analysis: A tutorial with R, JAGS, and Stan. 2nd ed. Amsterdam: Elsevier; 2014. https://doi.org/10.1016/B978-0-12-405888-0.09999-2.
Stan Development Team. Stan Modeling Language Users Guide and Reference Manual, Version 2.17.0. 2017. http://mc-stan.org/.
R Core Team. R: A language and environment for statistical computing. 2017. https://www.r-project.org.
Stephan JG, de Miranda JR, Forsgren E. Data_Rcode for: American foulbrood in a honeybee colony: spore-symptom relationship and feedbacks between disease and colony development. BMC Ecology. 2020. https://doi.org/10.5281/zenodo.3672367.
We thank Sepideh Lamei and Jeffery S. Pettis for the help during the collection of the data.
Open access funding provided by Swedish University of Agricultural Sciences. This research was supported by Grant Dnr 222-2013-423 from the Swedish Research Council for Sustainable Development (FORMAS). The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Ethics approval and consent to participate
Ethics approval was not required for experimental work on honeybees when this work took place, as covered by the Animal Welfare Act (AWA) of the United States Code (Title 7 Agriculture; Chapter 54 Transportation, Sale, and Handling of Certain Animals; Sections 2131–2159), which restricts the requirement for ethical approval for animal research to certain warm-blooded mammals. All bees and pathogens were obtained locally. The experimental apiaries were located on USDA property outside the flight range of surrounding non-experimental apiaries. All material used in the experimental apiaries was kept separate from non-contaminated material, in accordance with local AFB apiary management recommendations. No personal data was collected. Consent to participate therefore does not apply.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The original version of this article was revised: The typo in article title occured during the production process has been corrected.
About this article
Cite this article
Stephan, J.G., de Miranda, J.R. & Forsgren, E. American foulbrood in a honeybee colony: spore-symptom relationship and feedbacks between disease and colony development. BMC Ecol 20, 15 (2020). https://doi.org/10.1186/s12898-020-00283-w
- Apis mellifera
- Paenibacillus larvae
- Group size
- Enzootic disease
- Host–pathogen dynamics
- Social immunity
- Host density